Особенности термического режима рек

Особенности термического режима рек

3. Физические закономерности изменения температуры воды по глубине, ширине и длине рек

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

3.1 Уравнение теплопроводности

В общем случае к элементарному объему воды V поступают тепловые потоки Qi через все грани и вдоль каждой координатной оси (рис. 3.1). Потоки тепла через грани Q1, Q2, Q3, Q4, Q5, Q6 соответственно равны адвекции тепла с вышележащих участков реки и его выносу ниже по течению, поступлению тепла в результате дисперсии, его удаление в объеме поперечного переноса речной воды, удаление и поступление тепла в результате процесса конвективного теплопереноса. Перераспределение тепла по всем направлениям в пределах объема dV также связано с процессами турбулентного перемешивания и физической теплопроводности. В общем случае Q1Q2, Q3Q4, Q5Q6. Это приводит к изменению теплосодержания в этом объеме воды и ее температуры.

Рис. 3.1 Схема к поступлению и удалению тепла на гранях элементарного объема воды

Изменение потоков тепла вследствие физической (молекулярной) теплопроводности учитывается уравнением (Караушев, 1969)

                   (3.1)

где - тепловой поток по i-му координатному направлению, обусловленный физической теплопроводностью, V – объем воды,  – интервал времени. В соответствии с законом Фурье поток теплоты (Вт/м2), обусловленный этим механизмом теплопередачи, пропорционален градиенту температуры по направлению i и коэффициенту физической теплопроводности (Вт/м2 ×0С):

.                                             (3.2)

Замена в уравнении (3.1) соотношением (3.2) приводит к выражению:

.                  (3.3)

Считая, что температурное поле изотропно (т.е. ) получаем:

.                                  (3.4)

Поскольку сумма изменений частных потоков тождественно равна изменению теплосодержания (в соответствии с уравнение (2.1)), то

,                                     (3.5)

где С – удельная теплоемкость, ρ – плотность воды. Раскрытие полной производной dq/dt преобразует уравнение (3.5) к уравнению теплопроводности (уравнению Фурье-Кирхгофа)

,                           (3.6)

где v, u, w – продольная, поперечная и вертикальная компоненты скорости соответственно. Члены, связанные с ними учитывают вклад процессов адвекции, дисперсии и конвекции в изменение температуры воды. Отношение  называется коэффициентом температуропроводности (м2/с).

В водных потоках изменение теплового состояния в основном зависит от турбулентного теплопереноса. Суммарный эффект влияния физической теплопроводности и турбулентного теплопереноса с учетом осреднения всех членов уравнения (3.6) дает (Алексеевский, 2006):

,                (3.7)

где  – температура; , , , – осредненные, а ,  – пульсационные продольная, поперечная и вертикальная компоненты скорости течения. Условия переноса тепла в турбулентных потоках характеризует коэффициент турбулентной температуропроводности Он интегрально учитывает роль конвективной, адвективной, дисперсионной, а также турбулентной теплопередачи в суммарном изменении температуры объема воды. Роль физической теплопроводности несущественна по сравнению с турбулентным теплопереносом, поэтому уравнение (3.7) трансформируется к виду (Алексеевский, 2006):

,            (3.8)

в котором использовано условие изотропности температурного поля (). Более точным является выражение:

                     (3.9)

3.2 Эпюры вертикального распределения температуры воды

Закономерности вертикального изменения температуры воды q в реках изучены недостаточно. Первый способ теоретического описания распределения температуры воды по глубине реки был предложен В.А. Бергом (1962). Теоретические эпюры температуры хорошо согласуются с реальным изменением температуры воды по вертикали. Однако их получение трудоемко и ограничено условиями постановки решаемой задачи. В общем случае формулу для расчета теоретической эпюры температуры можно получить из уравнения (3.9). Для случая установившегося равномерного движения потока (), отсутствия поперечных составляющих осредненной скорости, неизменных по длине потока x и его ширине z температур:

                                       (3.10)

Производные в этом уравнении полные, поскольку учитывается изменение q лишь по одному координатному направлению. Согласно А.В. Караушеву (Караушев, 1977), коэффициент турбулентной диффузии

                                 ,        (3.11)

где h – глубина потока, v – скорость течения в данной точке потока, а

M = 0,7Cш+6                                   (3.12)

при и M = 48 = const при . Этот параметр, как и коэффициент Шези, имеет размерность м0,5×с-1. Коэффициент Шези

.                                          (3.13)

Подстановка в уравнение (3.10) значение коэффициента турбулентной диффузии (3.11) и соответствующие преобразования дают

                                        (3.14)

Решение этого уравнения имеет вид:

                                         (3.15)

где q1 и q2 – постоянные интегрирования. Замена в этом уравнении глубины потока относительной глубиной потока , а также введение константы a1 = С/g = 427 м/0К приводит к уравнению

                                                (3.16)

В качестве константы интегрирования q1 примем придонную температуру потока, а q2 – разность температуры воды в поверхностном слое qn реки и температурой у дна q1, т.е. q2 = qn – q1. Относительную глубину будем учитывать со знаком «–» для получения прямой температурной стратификации в период весеннего и летнего нагревания водной массы. Такая необходимость связана с выбором начала координат. Относительная глубина у поверхности, а необходимая для этого коррекция соответствует вынесению знака «–» в показатель степени при экспоненте в уравнении (3.16). В этом случае эпюра температуры воды описывается уравнением:

.                                   (3.17)

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

3.3 Поперечное и продольное распределение температуры воды

Оценим поперечное распределение температуры воды для условий, когда изменение температуры воды по длине потока стационарно и неизменно, течение установившееся и равномерное, поперечные и вертикальные составляющие осредненной скорости равны нулю. Указанные условия означают, что процессы адвекции тепла на локальном участке реки отсутствуют и, следовательно, уравнение теплопроводности (3.8) имеет вид

(3.18)

Аналитическое решение этого уравнения в общем случае отсутствует. Оно появляется при использовании полученного выше теоретического распределения температуры воды по глубине потока. Такой подход (по аналогии с методом плоских сечений при построении поля скоростей на участке реки) можно назвать «1,5D», так как решение производится «одномерными» методами (Великанов, 1954).

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

(3.19)

Решение этого уравнения дает распределение поверхностной температуры воды по ширине потока. Для решения воспользуемся схемой обозначений для прямоугольного сечения русла (рис. 3.2), где В-ширина реки b=B/2 – половина ширины реки, z – расстояние от берега, y – отметка горизонта воды от дна, h – глубина потока. Использование прямоугольной схематизации русла позволяет предположить, что распределение температуры воды в поперечном сечении такой формы при прочих равных условиях симметрично, тепловое влияние обоих берегов – одинаково, влияние поверхностей раздела «вода – воздух» и «вода – ложе» также одинаково по всей ширине потока. В этом случае можно рассматривать распределение температуры воды только для одной, например, правой половины русла (считая распределение температуры в левой половине русла симметричным). В центре потока значения температуры максимально отличаются от прибрежной температуры воды.

В естественных условиях русло чаще бывает несимметричным. Поэтому заменим b на bп – расстояние от берега до «середины» потока (точки, в которой температура воды максимально отличается от прибрежной), а координату z в уравнении (3.19) – на относительное удаление от берега z/bп = (Гончаров, 1962).

Рис. 3.2. Схема принятых обозначений для прямоугольной формы поперечного сечения русла

В этом случае решение уравнения (3.19) (с учетом коэффициента турбулентной диффузии по уравнению (3.11)) имеет вид:

                                                (3.20)

При замене a1 = С/g = 427 м/0К

                                          (3.21)

где константа интегрирования q1n равна поверхностной температуре воды на середине потока, а q2n – разность поверхностной температуры воды у берега qбn и в центральной части русла q1n т.е. q2n = qбn – q1n. Показатель степени в уравнении (3.21) должен включать знак «минус» для воспроизводства экспоненциальной функцией реального распределения температуры воды по ширине потока

.                                        (3.22)

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

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

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

, (3.23)

где B/h – относительная глубина, v* – динамическая скорость. Ее величина

   (3.24)

где I – уклон,‰. Преобразуем формулу (3.23) для получения выражения для расчета коэффициента турбулентной дисперсии в явном виде. Для этого запишем член -2,7, как lg0,002, а последний член – lg[(B/h)1,5], тогда

. (3.25)

Таким образом, коэффициент турбулентной дисперсии зависит от глубины и ширины потока, а также от величины его уклона. Подставляя значение DTy,z в уравнение (3.19), получаем:

 (3.26)

Объединив сомножители при втором члене уравнения (3.26), используя для этого формулы Шези и (3.24), в коэффициент a2 получаем выражение:

 (3.27)

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

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

                                       (3.28)

Решение этого уравнения имеет вид

                                            (3.29)

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

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

4. Особенности вертикального распределения температуры воды по глубине реки на локальных участках

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

Если температуру воды рассматривать в качестве консервативной примеси, то время ее выравнивания (после локального и мгновенного изменения на поверхности) по глубине потока  (Fischer, 1973), где  – коэффициент турбулентной диффузии, , - время Лагранжа, - время завершения процесса смешения по глубине потока. Источником изменения поля температур в локальной точке водного потока или в локальном створе может быть при впадении притоков, отведение подогретых вод, сброс воды из водохранилища, отличающейся по тепловым характеристикам. Неоднородность вертикальной структуры поля температуры существует на небольшом участке ниже места резкого изменения температуры. Его длина зависит от соотношения средней глубины, скорости потока и :  (Кондюрина, 2000). Условие  означает, что ниже источника примеси достигнута вертикальная однородность температуры воды.

Оценка масштабов неоднородности температуры воды по глубине реки затруднена несовершенством технических средств для измерения этой характеристики. Поэтому существование вертикальной неоднородности температуры воды подвергалось сомнению. В начале прошлого века обсуждался, например «закон равномерного распределения температуры в реках» во все сезоны года (Максимович, 1900; Альтберг, 1916). Однако измерения Ф.И. Быдина (1933), Л.А. Ячевского (1916), В.М. Сокольникова (1935), О.В. Ванеевой и М.Н. Панкратьевой (1941), Соколовой (1951), и др. показали, что основания для выделения такого закона нет, температура воды неодинакова по глубине рек.

4.1 Характеристика полевых данных

Тем не менее, данных о реальном распределении температуры воды по глубине потока очень мало, поскольку проведение таких исследований требует специального оборудования и значительных затрат времени. Наиболее детальны и убедительны данные о вертикальном распределении температуры воды, полученные Л.А. Ячевским (1916) на р. Неве с 4 июня до 2 июля 1915 г. (пример табл. 4.1, см. приложение №4) при изучении условий формирования донного льда. Измерения температуры воды проводились ртутными термометрами с точностью 0,010С, с плота, установленного на якорях, в 50 метрах от берега, на одной вертикали и на глубинах 0,05 м, 0,5, 1, 2 и 2,5 м. Измерения выполнялись на 1,8–1,9 км выше ответвления Большой Невы от р. Невы. Одновременно проводились наблюдения за облачностью и температурой воздуха. Значения температуры снимались в 1:00, 3:00, 5:00, 7:00, 13:00, 21:00 и 23:00, что связано с проверкой гипотезы Л.А. Ячевского о влиянии облачности и температуры воздуха на образование донного льда.

Таблица 4.1. Фрагмент базы данных по температуре р. Нева в 1915 г.

Достаточно детальные измерения вертикального распределения температуры проведены автором на р. Оке летом 2007 г. (табл. 4.2) в районе д. Трегубово (Озерский район Московской обл.). Измерения температуры воды на вертикалях производились вертушкой Нертока фирмы «Valeport», имеющей функции измерения скорости (м/с), направления течения, температуры воды (с точностью до 0,010С) и давления (мбар). Этот прибор позволяет выставлять время осреднения измерений от 1 с до 99 минут и устанавливать «ноль» давления, что обеспечивает более точно расположение вертушки на вертикали. Поскольку вертушка Нертока неинерционная, то осреднение, использованное во время измерений, составляет 1 мин. Вертикали делались на расстоянии 5–10 м от берега и далее по поперечному сечению реки через каждые 50 м. Точки измерений засекались с помощью GPS. На вертикали измерения температуры воды производились на глубинах 0,1 м, 0,5 м, 1 м, и далее через 1 м до дна. Измерения проводились в солнечную погоду, в различное время светлой части суток.    

Таблица 4.2. Фрагмент базы данных автора по измерениям температуры воды в р. Ока

Измерения проводились в 4-х поперечных сечениях реки. Один из них располагался на участке одиночного разветвления (рис. 4.1, 1), в правом рукаве с 13:11 по 13:52, в левом рукаве – с 18:50 по 19:41. В левом рукаве (В=95 м, B/h=63), вертикали назначались в среднем каждые 16 м. Измерения проведены на 6 вертикалях. Русло левого рукава заросшее, больше половины ширины рукава занято водной растительностью. В правом рукаве (В=231 м, B/h=195), вертикали располагались каждые 21 м. Измерения здесь выполнены на 9 вертикалях. Весь день (17.06.2007) была солнечная погода, θвозд=24–300С, вода в реке за день нагрелась на 10С (по данным наблюдений на водомерном посту).

На втором, третьем и четвертом поперечных профилях измерения выполнены в один день (18.06.07) в фазу дневного нагревания (с 16:39 до 20:11) при солнечной погоде и температуре воздуха θвозд=24–300С.

Рис. 4.1. Схема участка р. Ока с расположением створов наблюдений

Второй профиль располагался на 100 м выше подводного руслового карьера (В=270 м, B/h=144). Водное пространство по левому берегу покрыто кувшинками, занято тростником. По ширине профиля вертикали назначались в среднем через каждые 27 м. Наблюдения проведены на 7 вертикалях.

Третий профиль располагался на верхней границе карьера: на одной из вертикалей глубиной до 10 м (В=260, B/h=67) измерялись температуры воды. Измерения на вертикалях проводились в среднем через 33 м. Всего изучались эпюры 6 вертикалей.

Четвертый профиль пересекал карьер в створе, ширина которого составляла В=250 м (B/h=55,3). Измерения на вертикалях проводились в среднем через 28 м. Всего анализировалось 6 вертикалей.

Дополнительные материалы автором получены летом 2008 на р. Протва в районе д. Сатино (Калужская обл.) с помощью комплексного зонда YSI (рис. 4.2). Наблюдения проводились в различных морфологических и динамических условиях. Измерения температуры выполнены на нескольких вертикалях в пределах каждого створ реки. Измерения на всех вертикалях осуществлялись на разных глубинах (0,1, 0,5, 1 м и далее через 1 м до дна). Один створ расположился в плесе (табл. 4.3), другой – на перекате. Измерения проводились 02.07.2008 в дневные часы в солнечную погоду. Ширина плеса B составляла 20 м, средняя глубина hср= 1,05 м. Ширина переката B=45 м, средняя глубина hср= 0,35 м.

Рис. 4.2. Комплексный зонд YSI

Несколько распределений температуры по глубине были изучены 03.07.2008. Они располагались в зоне слияния р. Исьма и р. Протва в соответствии со схемой на рис. 4.3. Работы проводились в солнечную погоду при температуре воздуха 24–280С с 10:30 до 18:00. Результаты наблюдений приведены в Приложении №4. Температура воды в Протве в период измерений в узле ее слияния с р. Исьма (рис. 4.4) составляла θ = 20,3–20,9 0С, в Исьме – θ = 16,6 – 17,10С. Поэтому особенности температурных эпюр изначально соответствовали условиям формирования теплового состояния воды Исьмы, Протвы и зоны смешения водных масс этих рек.

Измерения температуры воды на вертикалях в узле слияния Исьмы и Протвы проводились в основном на стандартных горизонтах: 0,1 м, 0,5 м, 1 м и далее до дна через 1 м. На профиле 6 и в точках 4,5,6 профиля 5 (рис. 4.2) измерения выполнялись также на глубине 0,03 см с целью оценить изменение температуры в приповерхностном слое воды.

Рис. 4.3. Район проведения исследования поля температуры в области слияния рр. Протва и Исьма.

Рис. 4.4. Узел слияния рр. Протвы и Исьмы в июле 2008 г.

Таблица 4.3. Результаты измерений температуры в плесе р. Протвы 2008 г.

4.2 Анализ фактических данных

Распределение температуры воды на вертикали зависит от сезона года, синоптических условий, глубины рек, морфологии русла, расстоянии от берега до точки измерений, скорости течения и времени суток. Характер распределения отражает изменение температуры воды по глубине рек (рис. 4.5). Обычно оно равно разнице температур поверхностного и придонного слоя воды (Δθэ) (Соколова, 1951). Измерения автора и других специалистов (Ячевский, 1916) показывают, что максимальная величина разницы температур может возникать не только в этих слоях водной массы, но и в ее промежуточных слоях.

Рис. 4.5. Распределение температуры воды на вертикали в устье р. Невы

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

Для графика на рис. 4.6: θ1 = θ0,5 – 0,04 (коэффициент корреляции r = 0,998), а для графика на рис. 4.7 θ2 = 1,02θ0,5 -0,45 (коэффициент корреляции r = 0,998).

Коэффициент корреляции между рядами температур на смежных горизонтах водной массы r свидетельствует о сильном перемешивании водного потока, отсутствии выраженных градиентов температуры воды на вертикали. Тем не менее, она неодинакова на разных глубинах. Разность температуры воды изменяется в разные дни от 0 до 0,720С. Отсутствие детальных данных об условиях полевых экспериментов не позволяет обосновано определить причины возникновения больших или меньших отличий температуры воды по глубине потока. Можно лишь отметить, что изменение вертикального распределения температур зависит от общего тренда изменения температуры во времени, имеет некоторые внутрисуточные особенности. Это следует из анализа эпюр температуры воды для двух дней (рис. 4.8–4.9).

Рис. 4.6. Связь температуры воды в р. Нева на глубине 0,5 м и 1 м в 13:00 за весь период наблюдений

Рис. 4.7. Связь температуры воды на р. Нева на глубине 0,5 м и 2 м в 13:00 за весь период наблюдений

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

Рис. 4.8 Трансформация эпюры температур воды р. Нева во времени (1 – 1:00, 2 – 3:00, 3 – 5:00, 4 – 7:00, 5 – 13:00, 6 – 21:00 и 7 – 23:00) в условиях синоптического нагревания водной массы

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

поверхности обусловлено процессами теплообмена, поскольку воздух теплее воды (θв=180С, θ=17,930С) и поток тепла направлен в сторону водной массы. До 13:00 форма приповерхностной части эпюры сохраняется, хотя температура водывозрастает почти на 10. Это связано с активным поглощением солнечной радиацией и теплообменом между атмосферой и речной водой (θв=24,30С, θ=18,760С). В 21:00 температура воды повышается по сравнению с условиями в 13:00. Она увеличивается от поверхности до глубины 0,5 м и с 1 м до дна. В слое с 0,5 до 1 м температура воды понижается. Такое распределение может быть связано с сочетанием влияния двух факторов: излучения тепла водной поверхностью и теплообменом с атмосферой. Температура воздуха в это время суток выше температуры воды (θв=20,80С, θ=18,730С), поэтому приповерхностный слой теплее, чем вода на глубине 0,5 м. Можно предположить, что пониженные (по сравнению с более глубокими горизонтами) температуры воды в нижерасположенном слое связаны с ролью этого слоя в эффективном излучении воды. За счет эффективного излучения реки, средняя температура этой толщи воды понизилась. Одновременно турбулентный теплообмен с атмосферой способствует увеличению температуры воды по направлению к поверхности в горизонте 0 – 0,5 м. К 23:00 температура воды в целом снижалась. При этом в верхнем метровом слое температура воды понизилась на 0,10, а с увеличением глубины температура стала выше в сравнении с поверхностной. Изменение температуры воды в поверхностном слое связано с усилившимся эффективным излучением воды и небольшой разницей температуры воды и воздуха (θв=18,70С, θ=18,540С).

В период охлаждения водной массы (рис. 4.9) эпюры температур имеют одинаковый вид во все сроки на глубинах от 0,5 до 2,5 м: заметное понижение температуры при переходе от слоя 0,5 к слою 1 м, менее заметное увеличение температуры с глубиной на горизонтах 1–2 м (исключение – эпюры в 03:00 и 05:00, когда наблюдается понижение температуры с глубиной), заметное понижение температуры воды характерно для глубины 2–2,5 м. Заметное отличие эпюр есть в верхнем полуметровом слое. В 01:00, 03:00, 05:00 температура в этом слое понижается по направлению к поверхности воды. В 07:00, 13:00 и 21:00 температура воды в поверхностном слое выше, чем на глубине 0,5 м, а в 23:00 температура воды у поверхности вновь ниже, чем на глубине 0,5 м. Все эти трансформации эпюры температур сочетаются с колебанием средней на вертикали температуры воды в диапазоне 20,56 – 21,11 0С.

Таким образом, форма эпюры температур в течение суток меняется слабо, ее изменение в основном происходит в поверхностном слое воды, хотя средняя температура вертикали за день может изменяться на 10С и более. Изменение температуры в поверхностном слое воды связано с тепловым обменом между водной массой и атмосферой. Он обусловливает охлаждение воды ночью и ее нагревание днем. Общее же изменение температуры водной массы связано с суточным сочетанием составляющих теплового баланса. Ночью собственное излучение водной массы преобладает над притоком лучистой энергии, поэтому происходит охлаждение воды. При достаточной прозрачности воды, излучение и поглощение энергии характерно для всей водной массы, поэтому и охлаждение воды происходит по глубине потока приблизительно с одинаковой скоростью. В течение суток разность θ – θв меняется от 0,49 до 5,040С. Наличие относительно холодной воды в придонных слоях водного потока (рис. 4.9) может быть связано с охлаждением водной массы вследствие разгрузки грунтовых вод.

Перемешивание водной массы приводит к тому, что различия температуры на вертикали за все время измерений не превышали 0,720. Однако, эта величина достаточно большая, что позволяет говорить об умеренном характере влияния перемешивания воды на вертикальное распределение температуры воды. Можно предложить две гипотезы, объясняющие эти особенности распределения температуры воды с глубиной: 1) турбулентность водного потока не является определяющим фактором в перемешивании водной массы и 2) поскольку наблюдения производились в придельтовом участке реки, то влияние турбулентного перемешивания уменьшается вследствие снижения скорости течения. Водная масса в целом хорошо перемешана, поэтому отличия температуры воды на разном удалении от поверхности не столь велики по сравнению с водоемами. Скорость течения воды мала настолько, что турбулентное перемешивание не достаточно интенсивно для разрушения температурной стратификации.

Обработка данных по 145 эпюрам температуры на одной вертикали р. Нева показала (рис. 4.10), что изменение температуры воды по глубине реки (h = 2,5 м), в основном составляет 0,10С. За время измерений (22 дня) перепад температур воды по глубине потока колебался от 0,02 до 0,490С. Это довольно большая величина по сравнению с данными, полученными для рек с большими скоростями течения и большей интенсивностью турбулентного перемешивания.

Анализ данных по 34 вертикалям, находящимся на разном расстоянии от берега р. Оки, отличающимся по глубине и средней скорости течения, показывает, что изменение температуры воды Δθэ по глубине потока находится в пределах от 0 до 0,130С.

Рис. 4.10. Гистограмма значений Δθэ для р. Нева

Повторяемость величин максимальных перепадов температуры воды (рис. 4.11) существенно неодинакова. Наиболее часто максимальная разность температур воды по глубине потока составляет 0,010С. В 73,6% случаев эта разность и лишь в 26,4% она превышает 0,02. Максимальная величина .

Степень однородности значений θ по глубине мало зависит от расположения вертикали по отношению к берегам реки. В трех из 7 случаев полная однородность температуры воды наблюдалась возле берега, а в остальных случаях – близко к стрежню потока. Устойчивой связи, между величиной изменения температуры воды на вертикали и числом Рейнольдса не обнаружено (рис. 4.12).

Наибольшие вертикальные градиенты температуры воды чаще формируются на небольшом расстоянии от берега (15 м), при глубине потока h >2 м и скорости течения м/с. Гистограммы Δθэ при средних скоростях течения на вертикали м/с (рис. 4.13) обнаруживают выраженную неравномерность температуры воды для всех вертикалей. Тем не менее, пока данных недостаточно для однозначного утверждения об увеличении вертикальных аномалий температуры воды при переходе от поверхностных к придонным слоям воды на участках рек при уменьшении скоростей течения. Наибольший градиент 0,130С, например, соответствует средней скорости на вертикали 0,4 м/с.

Рис. 4.11. Изменчивость величины Δθэ на разных вертикалях р. Ока

Рис. 4.12. Соотношение максимальной разности температуры воды на вертикали и числом Рейнольдса (р. Ока, июнь 2007 г.)

Рис. 4.13. Гистограмма распределения значений Δθэ при средней скорости на вертикали v < 0,2 м/с

Характер распределения температуры воды по глубине потока учитывает и градиент температуры gradθ (в 0С/м). Изменчивость этой характеристики для р. Ока характеризует гистограмма на рис. 4.14. Из нее следует, что градиенты температуры воды меняются от 0 до 0,20С/м. Для большинства случаев gradθ равен < 0,030С/м. Чем больше градиент температуры воды, тем меньше вероятность его формирования.

Изменение градиентов температуры воды по глубине р. Оки, например, характеризует рис. 4.15. Из рисунка следует, что величина градиента больше 0,050С не встречается в слое воды y > 1 м, где y – расстояние от поверхности воды до точки измерений. При y > 2 м практически отсутствуют градиенты температуры, превышающие 0,020С. Единственная вертикаль, на которой нарушается это правило, находится в русловом карьере, где распределение температуры воды имеет более сложный характер (рис. 4.16). Здесь, вероятно, на глубинах 5–7 м формируется противотечение, поставляющее более теплую воду на значительные глубины. В переуглубленном русле реки возникают процессы, подобные процессам в водоемах, которые влияют на формирование типичных распределений гидрофизических характеристик по глубине водного объекта.

Рис. 4.14. Гистограмма распределения максимального вертикального градиента температур для условий р. Ока

Рис. 4.15. Изменчивость локальных значений градиента температуры по глубине р. Ока

Рис. 4.16. Распределение температуры воды по глубине на р. Ока в пределах руслового карьера

Данные, полученные на р. Протве, характеризуют разнообразие условий формирования эпюр температуры. Оценка величины изменения температуры воды на вертикалях (рис. 4.17) показывает, что наиболее часто встречаются эпюры температур с величиной Δθэ = 0,01 – 0,08. В большинстве случаев они формируются в пределах водной массы Протвы и Исьмы. Все вертикали, для которых изменение температуры воды Δθэ > 0,10С, находятся в зоне смешения. Следовательно, смешение вод сопровождается своеобразным «гашением» турбулентности, что и обусловливает большую вертикальную неоднородность водной массы в зоне смешения речных вод. Часть вертикалей, для которых величина Δθэ < 0,10С находится в зоне смешения водных масс рек (рис. 4.18). Средняя температура воды здесь выше, чем в р. Исьма (на 1–1,50С), и ниже по сравнению с Протвой. Это свидетельствует о формировании подзоны смешения с чертами процесса, соответствующими рр. Протва и Исьма.

Рис. 4.17. Изменчивость величины Δθэ на разных вертикалях р. Протва, июнь 2007 г.

Рис. 4.18. Распределение величины Δθэ на разных вертикалях р. Протва и р. Исьма

Рис. 4.19. Эпюры в водной массе р. Протва (1), р. Исьма(2) и в зоне смешения

Сравнение температурных эпюр в р. Протва, р. Исьма и в зоне их смешения (рис. 4.19) показывает, что форма эпюр в водных массах этих рек значительно слабее выражена, чем в зоне смешения. Величина максимального изменения температуры воды Δθэ в зоне смешения может достигать 1,380С, тогда как в Протве и Исьме эта величина обычно <0,10C.

Вертикальные градиенты температуры чаще всего (в 60% случаев) меньше 0,10С/м. Градиенты более 20С/м встречаются в 7% случаев. Самые большие величины gradθ наблюдаются в приповерхностном слое воды при переходе от глубины 0,03 к 0,1 м, здесь градиенты температуры достигают 7,710С/м, однако встречаются и случаи gradθ = 0. Это свидетельствует о большом разнообразии теплового состояния поверхностного слоя, зависящим, в частности, от изменения глубины реки.

Анализ данных, полученных на плесе и перекате (рис. 4.20) в р. Протва свидетельствуют о подобии градиентов температуры воды и разности Δθэ на разных вертикалях. Для плеса gradθ для различных вертикалей меняется от 0 до 0,030С, а для переката эта величина меняется от 0 до 0,050С, т.е. различий между ними практически нет. В приповерхностном слое (0,03 – 0,1 м) величина градиента температур составляет 0,530С/м. Она отражает краткосрочную синоптическую обстановку теплового взаимодействия между водной массой реки и приводным слоем атмосферы.

Выявленные особенности изменения температуры воды по глубине рр. Ока, Протва и Исьма не являются характерными для всех рек. В частности, для многих рек свойственна большая вертикальная неоднородность температур. Экстремально большие различия температуры воды по глубине характерны, например, для истока Ангары (Верещагин, 1933). Максимальная разность температур по глубине реки здесь достигает 70С, что связано с озерным генезисом ее вод. Большие вертикальные температурные градиенты характерны и для некоторых устьев крупных северных рек. В районе Усть-Енисейского порта придонная и поверхностная температуры могут отличаться на 60С (Соколова, 1951). По ее же данным на р. Лена (с. Кюсюр) в период летнего нагревания водной массы разность температуры воды поверхностного и придонного слоя воды достигает 0,40С. В период ее охлаждения эта разность еще больше – 1,10С, при этом придонная температура воды оказывается больше поверхностной. На р. Лена у с. Солянка наибольшие положительные разности Δθэ наблюдаются в июне, июле и августе (период нагревания) и достигают 1,550С. В период последующего охлаждения температура воды в потоке выравнивается, а в конце сентября и начале октября возникает обратная температурная стратификация (Δθэ > -0,60С).

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

Рис. 4.21. Возможные типы (1,2,3,4,5,6,7,8,9) температурных эпюр на границе водной массы и русел рек, приземного слоя атмосферы

Форма нижней части эпюры зависит от направления потока тепла, формирующегося при поступлении грунтовых вод или тепловом взаимодействии речной водной массы с грунтами. Если поток тепла на границе «вода-ложе» отсутствует, то формируются эпюры типа 1,2,3. Если он направлен в сторону водной массы (ее теплосодержание возрастает), то преобладает формы эпюр типа 4,5,6. Такие типы эпюр могут возникать в зимний период при разгрузке относительно теплых грунтовых вод или в период весеннего нагревания, когда ложе потока нагревается быстрее, по сравнению с водной массой. Если поток тепла направлен от водной массы к грунтам, то происходит охлаждение водной массы. Это может происходить в летний период при поступлении в основной поток охлаждающих грунтовых вод и в период осеннего охлаждения, когда ложе потока охлаждается быстрее, чем водная масса. Типы эпюр 3,6,9 формируются при равенстве потоков тепла от грунта к водному объекту и в обратном направлении.

Форма верхней части эпюры зависит от многих факторов, влияющих на теплообмен на границе «река-атмосфера». Типы 1,4,7 соответствуют охлаждению водной массы за счет теплообмена с атмосферой. На реках такие формы эпюры формируются в период суточного, синоптического или сезонного охлаждения водной массы. Типы эпюр 2,5,8 в естественных условиях формируются в период суточного, синоптического или сезонного нагревания водной массы. Типы эпюр 3,6,9 возникают в том случае, если тепловое взаимодействие атмосферы и речной водной массы уравновешено. Они могут формироваться в утренние и вечерние часы, когда интенсивность солнечной радиации компенсируется эффективным излучением воды. Кроме того, нередко формирование такой формы эпюры у берегов при малых глубинах и малых скоростях течения. При очень больших скоростях течения происходит практически полное смешивание различных слоев водной массы, поэтому в любую погоду, независимо от сезона, эпюра температур выражена слабо.

Страницы: 1, 2, 3, 4