Внесены ВНИИГ им. Б.Е. Веденеева

Утверждены ВНИИГ им. Б.Е. Веденеева решением N 27 от 12.05.1986 г. и согласованы с Главтехуправлением

"Рекомендации по расчету длины полыньи в нижних бьефах ГЭС" содержат изложение методики расчета длины полыньи при отсутствии и наличии тепловых стоков в нижнем бьефе. Приводятся также список необходимых исходных данных, примеры расчета.

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

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

Настоящие Рекомендации посвящены вопросу расчета длины полыньи в нижних бьефах ГЭС в режимах наступления, отступления и стабилизации кромки льда; в Рекомендациях приводится также методика расчета длины полыньи с учетом влияния тепловых сбросов, расположенных в нижнем бьефе.

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

Рекомендации составлены сотрудниками ВНИИГ им. Б.Е. Веденеева: мл. науч. сотр., канд. техн. наук Е.Л. Разговоровой и мл. науч. сотр. Г.А. Трегуб. Научный руководитель работы - доктор техн. наук Л.И. Кудояров.

В работе над Рекомендациями учтены ценные замечания, высказанные Я.Л. Готлибом, Г.Н. Нисар-Мухамедовой (Гидропроект им. С.Я. Жука), С.Н. Назаренко (Ленгидропроект), Р.В. Донченко (ГГИ), А.Ж. Жулаевым (КазНИИВХ), М.М. Бейлинсоном (КазПИ им. Абая), В.В. Баланиным, Б.С. Бородкиным, М.И. Жидких (ЛИВТ), А.И. Бефани, Л.Е. Кресс (ОГМИ), И.И. Макаровым, В.А. Кякком, И.Н. Соколовым, И.Н. Шаталиной, А.Б. Векслером (ВНИИГ им. Б.Е. Веденеева).

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

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

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

1.4. В Рекомендациях приняты международная система единиц СИ и следующие обозначения:

a_л, a_с - температуропроводность льда и снега, соответственно, м²/с;

b - ширина нижнего бьефа по урезу воды, м;

b_з - ширина заберегов, м;

c, c_л - удельная теплоемкость воды и льда, соответственно, Дж/(кг·К);

C - коэффициент Шези, м^0,5/с;

h_л, h_с - толщина льда и снега, соответственно, м;

h_н - толщина льда на кромке (начальная толщина льда), м;

H - глубина нижнего бьефа, м;

q_ст - мощность тепловых стоков, Вт;

Q, Q_ш - расходы воды и шуги, м³/с;

R - гидравлический радиус;

S - интенсивность теплового потока, Вт/м²;

V - скорость течения, м/с;

x - координата по длине потока, м;

x₀, x_н - координаты створа нулевой изотермы и начала внутриводного ледообразования, соответственно, м;

x_к.о - длина полыньи в начале расчетного периода, м;

x_к - длина полыньи в конце расчетного периода, м;

z - координата по глубине потока; м;

- коэффициент теплоотдачи от воды к воздуху, Вт/(м²·К);

- коэффициент теплоотдачи от воды к нижней поверхности льда, Вт/(м²·К);

- коэффициент теплоотдачи от льда (снега) к воздуху, Вт/(м²·К);

t₀ - температура воды, поступающей из верхнего бьефа в нижний, °C;

- степень покрытия шугой водной поверхности;

- температура воздуха, °C;

, - эквивалентная температура воздуха над поверхностью воды и льда, соответственно, °C;

, - теплопроводность льда и снега, Вт/(м·К);

- коэффициент гидравлического трения;

n_ш, n_р - коэффициенты шероховатости нижней поверхности шуги и русла, с/м^0,33;

, , - плотность воды, льда, шуги, соответственно, кг/м³;

- удельная скрытая теплота кристаллизации, Дж/кг;

- время, с.

Безразмерные критерии и параметры:

Bi - критерий Био, ;

Mi - критерий Михеева, ;

- параметр расхода шуги, ;

X - параметр координаты, ;

- параметр времени, ;

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

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

2.2. Морфологические параметры включают в себя:

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

б) поперечные разрезы русла по длине нижнего бьефа при тех же расходах.

2.3. Гидрологические параметры включают в себя:

а) расход воды: для строящихся и эксплуатируемых ГЭС при проведении поверочных расчетов расход воды в нижнем бьефе выбирается средним по пятидневкам (декадам), на которые разбивается расчетный период; для расчета длины полыньи в нижних бьефах проектируемых ГЭС выбираются расходы воды в реке для средневодного, маловодного и многоводного годов соответствующей обеспеченности по пятидневкам (декадам);

б) толщины льда и снега в нижнем бьефе по пятидневкам (декадам);

в) средние и средневзвешенные скорости течения по длине предполагаемого расчетного участка для заданных расходов воды;

г) даты ледостава и вскрытия;

д) предледоставный расход и плотность шуги; если пористость шуги неизвестна, то ее следует выбирать на основании приложения 1.

е) ширина заберегов по пятидневкам (декадам) за расчетный период; если неизвестна ширина заберегов, то ее рекомендуется рассчитывать по формуле [13]:

где m = 0,42V^0,2 + K; коэффициент K для русел с мелководными зонами и пологими берегами равен 0,525, при отсутствии мелководных зон и крутых берегах K = 0,465.

2.4. Метеорологические данные должны относиться к району расположения полыньи, охватывать расчетный период и содержать следующие сведения: температуру воздуха , скорость ветра W, облачность общую и нижнюю n_о, n_н, абсолютную влажность воздуха e. Возможно осреднение метеорологических данных по пятидневкам или декадам.

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

а) Средняя скорость ветра на высоте 2 м над водной поверхностью определяется по формуле

где K₁ - коэффициент, учитывающий степень защищенности метеорологической станции на суше (табл. 1); K₂ - коэффициент, учитывающий характер рельефа в пункте наблюдений (табл. 2); K₃ - коэффициент, учитывающий среднюю длину разгона воздушного потока над водной поверхностью при различной защищенности станции (табл. 3); W_ф - скорость ветра на высоте расположения флюгера, м/с.

б) Средняя влажность воздуха на высоте 2 м над водной поверхностью должна рассчитываться по формуле

где e' - средняя влажность воздуха, измеренная на континентальной метеостанции, мб; e₀ - максимальная упругость водяного пара, определенная по температуре поверхности воды (табл. 4), мб; M₁ - коэффициент трансформации, учитывающий изменение влажности и температуры воздуха над водной поверхностью, находится по табл. 5 в зависимости от среднего расстояния от подветренного берега и разности температур воды и воздуха.

в) Средняя температура воздуха над водной поверхностью на высоте 2 м должна рассчитываться по формуле

где - температура воздуха по данным метеорологической станции, расположенной на суше, °C; t_п - температура поверхности воды для расчета длины полыньи, t_п ~= 0,5t₀, °C.

г) Для участков рек шириной не более 500 м значения , e, W следует принимать без введения поправок.

2.6. При проведении расчета должны быть известны:

а) температура воды в начале нижнего бьефа t₀ или температура воды в створе нижнего бьефа, ближайшем к ГЭС;

б) при наличии тепловых сбросов в нижнем бьефе - их мощность или расходы и температуры воды тепловых сбросов.

2.7. Если по рассматриваемому объекту необходимые исходные данные отсутствуют, то их подбор должен производиться по методу аналогов.

2.8. Необходимые для расчета длины полыньи физические постоянные даны в приложении 1.

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

3.2. Расчет коэффициента теплоотдачи от воды к воздуху и эквивалентной температуры воздуха над поверхностью воды следует проводить на основании п. 2 "Рекомендаций по термическому расчету водохранилищ" [7].

а) Эквивалентная температура воздуха и коэффициент теплоотдачи от воды к воздуху рассчитываются по формулам:

где S_R - интенсивность радиационного теплообмена, Вт/м²; S_и - интенсивность теплообмена при испарении, Вт/м²; - функция, учитывающая увеличение интенсивности испарения за счет разности температур воды и воздуха, определяется по табл. 6.

б) Интенсивность радиационного теплообмена должна определяться зависимостью

При отсутствии данных о распределении облачности по ярусам для расчета следует использовать формулу

Здесь (Q + q)₀ - суммарная солнечная радиация на уровне моря при альбедо, равном нулю, Вт/м², находится по табл. 2-1 приложения 2; k_e - коэффициент, учитывающий отклонение влажности воздуха от ее среднесуточного значения, определяется по формуле k_e = 1 + n(e_пр - e), где e - влажность воздуха, наблюдаемая на высоте 2 м над подстилающей поверхностью, мб; величины n и e_пр даны в табл. 2-2 приложения 2; k_z - коэффициент, учитывающий влияние высоты местности над уровнем моря, определяется по табл. 2-3 приложения 2; n_о, n_н - облачность общая и нижняя, доли единицы; k_н и k_в+с - коэффициенты, учитывающие задержание суммарной радиации облаками нижнего и совместно верхнего и среднего ярусов, определяются по табл. 2-4, 2-5 приложения 2; A - среднесуточное альбедо поверхности воды, определяется по табл. 2-6 приложения 2; - доля рассеянной по направлению к поверхности воды радиации, которую следует определять по формуле

b', b" - величины, зависящие от влажности воздуха и облачности, определяются по табл. 2-7, 2-8 приложения 2; - постоянная Стефана - Больцмана, равная 5,777·10^-8 Вт/(м²·К⁴); - абсолютная температура воздуха на высоте 2 м, равная , К; - температура воздуха на высоте 2 м, °C; значения определяются по табл. 2-9 приложения 2; k - коэффициент, определяемый по табл. 2-10 приложения 2; n - наблюденная облачность, доли единицы; c₁ - коэффициент, определяемый по табл. 2-11 приложения 2; I - эффективное излучение при безоблачном небе, Вт/м², определяется по табл. 2-12 приложения 2; T_п - абсолютная температура поверхности воды T_п = (273,16 + t_п), К.

в) Интенсивность теплообмена при испарении следует рассчитывать по формуле

3.3. Расчет коэффициента теплоотдачи от поверхности льда (снега) к воздуху рекомендуется производить по формуле [11]

где W - скорость ветра, м/с; коэффициент B определяется по табл. 7.

3.4. Эквивалентную температуру воздуха над ледяным покровом следует рассчитывать по формуле

где S_R.л, S_и.л - соответственно, интенсивность радиационного теплообмена и теплообмена при испарении с поверхности льда.

Интенсивность радиационного теплообмена следует рассчитывать, используя зависимости (7) и (8), но вместо альбедо поверхности воды подставлять альбедо льда или снега (табл. 8).

Температура поверхности льда должна определяться по формуле

при наличии снежного покрова температура поверхности снега:

где h_{л. 0} - толщина льда в начале расчетного периода; h_с - толщина слоя снега.

Интенсивность теплового потока при испарении с поверхности льда следует находить, используя зависимость [20]

где K₄ определяется в зависимости от разности температур поверхности льда и воздуха (табл. 9).

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

4.2. Условия существования того или иного режима движения кромки льда представлены в табл. 10.

4.3. Расчет положения створа нулевой изотермы должен проводиться по формуле

4.4. Если имеют место забереги по длине нижнего бьефа, то вместо ширины русла b в формулу (15) подставляется величина b₁:

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

4.6. Длина участка между створом нулевой изотермы и створом начала внутриводного ледообразования рассчитывается в соответствии с приложением 3 [16].

4.7. Толщина льда на кромке при ее наступлении определяется из условий предельно-напряженного состояния шугового ковра и должна рассчитываться по формуле [4, 6, 21]

а) Коэффициент Шези следует определять по формуле

где R - гидравлический радиус, м; R = F₀/P₀; F₀ - площадь поперечного сечения русла при заданном расходе, м²; P₀ - смоченный периметр, м.

Если русло прямоугольное, то

где

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

где n_р, n_ш - коэффициенты шероховатости русла при открытой водной поверхности и нижней поверхности шуги, соответственно.

в) Определение коэффициентов шероховатости русла и нижней поверхности шуги необходимо проводить по следующим формулам [8]:

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

- высота выступов нижней поверхности шуги (льда), м:

h_{н. 0} - начальная толщина льдин, из которых формируется кромка; при h_{н. 0} = 0,025 м; при h_{н. 0} = 0,05 м; при h_{н. 0} = 0,2 м [13].

г) Для нижних бьефов действующих ГЭС, когда имеются натурные наблюдения над толщиной льда, величину h_н следует рассчитывать по формуле

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

4.8. Для действующих ГЭС начальное положение кромки льда определяется по данным натурных наблюдений.

4.9. Для проектируемых ГЭС или в случае отсутствия данных о начальном положении кромки при проведении поверочных расчетов для действующих ГЭС начальное положение кромки льда рассчитывается по формуле

где X_н - безразмерный параметр, определяющий соотношение между интенсивностью теплообмена воды с атмосферой и объемной скрытой теплотой кристаллизации, которая выделяется при образовании шуги, формирующей кромку, при толщине льда на кромке h_н, в единицу времени; параметр X_н находится по графику на рис. 1 в соответствии с п. 4.10 настоящих Рекомендаций.

4.10. Параметр X_н определяют либо по известному предледоставному расходу шуги Q_{ш. 0}, либо по ширине заберегов в предледоставный период b_{з. 0}, измеренной у створа, ближайшего к кромке льда.

а) Если задан предледоставный расход шуги Q_{ш. 0}, то следует рассчитать параметр :

и по известному значению параметра и по графику (рис. 1) найти X_н.

б) Если задана ширина заберегов b_{з. 0}, то следует определить степень покрытия шугой водной поверхности вблизи кромки льда в предледоставный период:

и по известному значению и графику на рис. 1 найти параметр X_н.

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

где x_{к. 0} - положение кромки льда в начале расчетного периода.

4.12. Параметр X_к должен находиться по графику (рис. 2) при значении параметра :

где определяется по графику (рис. 2) при X_к = X_{к. 0}. Параметр характеризует отношение количества тепла, отданное водой в атмосферу за расчетный период , к объемной скрытой теплоте кристаллизации, выделившейся при образовании слоя шуги толщиной h_н.

4.13. Параметр X_{к. 0} следует вычислять по формуле

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

где V_кр - скорость отступления кромки льда.

4.15. Скорость отступления кромки должна определяться на основании анализа теплового баланса на прикромочном участке таяния в зависимости от температур воздуха и воды, подходящей к кромке льда. Выбор расчетных формул для определения V_кр проводится в соответствии с табл. 11, исходя из условий теплообмена льда с водой и воздухом:

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

4.16. Количество тепла, приносимое к кромке водой в единицу времени, q₀, должно определяться зависимостью

где t_кр - температура воды на кромке льда:

Количество тепла, уносимое потоком в единицу времени от створа, в котором прекращается таяние на нижней поверхности льда, q_т, следует находить по формуле

Здесь t_т - температура воды в створе, где прекращается таяние льда с нижней поверхности:

где ; - коэффициент теплоотдачи от воды ко льду, рассчитываемый по формуле [11]

Приведенную толщину льда на участке таяния с учетом влияния слоя снега следует вычислять по формуле

h_л.н, h_с - толщины льда и снега в начале расчетного периода, м.

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

где - средняя температура воды на участке таяния:

Длину участка таяния x_т следует находить по формуле

Величина скрытой удельной теплоты ледообразования с учетом запаса холода во льду должна определяться зависимостью [9]

4.17. Для условий пп. 2, 3 в табл. 11 скорость отступления кромки должна определяться формулой

4.18. Для условий п. 4 в табл. 11 скорость отступления кромки следует рассчитывать, пользуясь формулами:

где

4.19. Расчет длины полыньи в нижнем бьефе должен проводиться в той последовательности, которая соответствует блок-схеме, представленной на рис. 3.

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

5.2. При впадении теплового стока мощностью q_ст выше створа нулевой изотермы координата створа нулевой изотермы с учетом влияния теплового стока на температуру воды должна определяться формулой

где x_с - координата створа впадения стока, м;

- суммарный расход основного потока и стока.

5.3. Температуру воды всего потока в створе стока следует рассчитывать по формуле

5.4. После определения положения створа нулевой изотермы и температуры воды расчет длины полыньи должен производиться в соответствии с требованиями раздела 4 настоящих Рекомендаций, но за начальный створ принимается створ впадения стока, т.е. вместо x₀ и t₀ используются, соответственно, x_{0. с} и t_ст, за расчетный расход воды принимается расход , все гидравлические и морфометрические характеристики русла должны выбираться для расхода .

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

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

5.6. Формулы для расчета выбираются на основании табл. 12 в зависимости от значения параметра Ш, который представляет собой отношение объемной скрытой теплоты кристаллизации, необходимой для образования шуги расходом в основном потоке до впадения стока, к мощности теплового стока:

где - расход шуги основного потока выше створа стока:

параметр должен определяться по графику (рис. 1) при

5.7. При условии Ш > 1 (табл. 12) величина должна рассчитываться по формуле

где параметр X₁ находится по формуле (53), параметр X₂ находится по графику на рис. 1 при значении параметра :

5.8. При условии Ш < 1 (табл. 12) имеет место полное таяние шуги и нагревание воды в потоке выше нуля, при этом порядок расчета следующий.

а) Сначала необходимо определить температуру потока в створе стока:

б) Затем следует рассчитать длину полыньи в соответствии с требованиями раздела 4 настоящих Рекомендаций, но вместо t₀ используется значение температуры потока в створе впадения стока t_ст, рассчитанное по формуле (56).

5.9. Расчет длины полыньи при отступлении кромки льда должен проводиться в соответствии с требованиями раздела 4.

5.10. Температуру воды на кромке следует определять по формуле

где температура воды в створе стока рассчитывается по формуле (49).

Примечание. При n_о = 1 (сплошная облачность) b' = 0.

Начало процесса внутриводного ледообразования возможно при выполнении трех условий [10, 12] - переохлаждения воды, наличия устойчивых зародышей льда и отвода теплоты кристаллизации от границы фазового перехода. Участок от створа нулевой изотермы до створа, в котором начинается образование кристаллов внутриводного льда, может быть весьма протяженным и составлять несколько километров, поэтому длину этого участка необходимо оценивать при проведении расчета длины полыньи. Расчет длины участка между створами нулевой изотермы и начала внутриводного ледообразования следует проводить на основании зависимостей, представленных в [16]:

где

x_н - координата створа начала внутриводного ледообразования; x₀ - положение створа нулевой изотермы, вычисляемое по формуле (15) настоящих Рекомендаций. Номограмма, позволяющая рассчитывать величину x_н - x₀ и температуру воды в створе начала внутриводного ледообразования t_н, построенная по приведенной зависимости (1), представлена на рисунке.

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

Для Q = 500 м³/с, b = 280 м, V = 1 м/с, , параметр . Аппроксимируя график в области от 0 до 10^-6 зависимостью и используя эту зависимость, можно найти Mi_н = 2,8·10^-4, откуда длина участка между створами нулевой изотермы и начала внутриводного ледообразования x_н - x₀ = 175 м.

Пример 1. В нижнем бьефе ГЭС ширина заберегов в предледоставный период в районе кромки льда b_{з. 0} = 80 м, забереги наблюдаются только на шугообразующем участке. Ледостав наступает 10 ноября. Толщина льда, измеренная у берега в районе кромки 15 ноября, h_{л. 0} = 0,32 м. Расход воды в нижнем бьефе Q = 500 м³/с, скорость течения V = 1 м/с, ширина b = 280 м, температура воды в начале нижнего бьефа t₀ = 3 °C, коэффициенты теплоотдачи от воды к воздуху и льду и эквивалентная температура воздуха , , эквивалентная температура воздуха надо льдом , плотность шуги .

Найти начальное положение кромки льда.

Решение.

1. Положение створа нулевой изотермы определяется по формуле (15):

2. Толщина льда на кромке рассчитывается по формуле (25); продолжительность периода от даты ледостава до даты измерения толщины льда составляет 5 дней (с 10 по 15 ноября):

3. Степень покрытия шугой водной поверхности вблизи кромки перед установлением ледостава при полном покрытии поверхности заберегами и шугой рассчитывается по формуле (28):

4. Параметр координаты, соответствующий начальному положению кромки, находится по рис. 1, исходя из величины :

5. Начальное положение кромки рассчитывается по формуле (26):

Пример 2. На начало третьей декады ноября (20 ноября) длина полыньи в нижнем бьефе составляла x_{к. 0} = 100000 м. Температура воды, сбрасываемой с ГЭС в нижний бьеф, t₀ = 2,8 °C; расход воды Q = 500 м³/с; V = 1 м/с, b = 280 м, ширина заберегов по всему нижнему бьефу в среднем b_{з. 0} = 80 м, , , коэффициенты шероховатости русла и нижней поверхности шуги n_р = 0,05, n_ш = 0,015, плотность шуги .

Найти длину полыньи в конце третьей декады ноября (30 ноября) .

Решение.

1. Положение створа нулевой изотермы определяется по формуле (15) с учетом требований п. 4.3 настоящих Рекомендаций:

Так как x₀ < x_{к. 0}, то в соответствии с п. 1 табл. 11 имеет место наступление кромки.

2. Толщина льда на кромке определяется в соответствии с п. 4.7 настоящих Рекомендаций:

3. Параметр X_{к. 0} рассчитывается по формуле (31):

4. Параметр находится по графику на рис. 2 при X = X_{к. 0} = 0,129:

5. Параметр рассчитывается по формуле (30):

6. Параметр X_к находится по графику на рис. 2 при :

7. Положение кромки льда в конце третьей декады рассчитывается по формуле (29):

Пример 3. Исходные данные те же, что и в примере 2, но , , . Толщина слоя снега на льду h_с = 0,1 м, толщина льда в нижнем бьефе h_{л. 0} = 0,2 м.

Найти длину полыньи в конце третьей декады ноября (30 ноября).

Решение.

1. Положение створа нулевой изотермы определяется по формуле (15) с учетом требований п. 4.4 настоящих Рекомендаций:

Так как x₀ > x_{к. 0}, то в соответствии с п. 2 табл. 11 имеет место отступление кромки льда.

2. Температура воды на кромке рассчитывается по формуле (35) при x = x_{к. 0}:

3. Приведенная толщина льда находится по формуле (39):

4. Температура воды под ледяным покровом, при которой прекращается таяние льда на нижней поверхности, рассчитывается по формулам (37), (38):

5. По формуле (43) определяется величина :

6. Длина участка таяния рассчитывается по формуле (42):

7. Средняя температура воды подо льдом на участке таяния определяется по формуле (41):

8. Скорость отступления кромки находится по формулам (33), (34), (36), (40):

9. Длина полыньи в конце расчетного периода рассчитывается по формуле (32):

Пример 4. Исходные данные те же, что и в примере 2, но в створе x_с = 85000 м действует тепловой сброс мощностью q_ст = 10⁸ Вт. Расход воды ниже створа стока .

Найти длину полыньи в конце третьей декады ноября.

Решение.

1. Положение створа нулевой изотермы определяется по формуле (15) с учетом требований п. 4.4 настоящих Рекомендаций:

Так как x₀ < x_с, то сток расположен на шугообразующем участке.

2. Безразмерный параметр координаты при x = x_с определяется по формуле (53):

3. Безразмерный параметр расхода шуги в створе стока определяется по графику (рис. 1) при :

4. Объемный расход шуги в створе стока определяется по формуле (52):

5. Параметр Ш рассчитывается по формуле (51):

Так как Ш > 1, то согласно табл. 12 шуга в створе стока тает частично.

6. Параметр расхода шуги после ее таяния рассчитывается по формуле (55):

7. Безразмерный параметр координаты, соответствующий , находится по рис. 1 при :

8. Увеличение длины шугообразующего участка при таянии шуги за счет теплового стока определяется по формуле (54):

9. Длина полыньи в нижнем бьефе в данном примере определяется по формуле (50):

1. Пехович А.И., Трегуб Г.А. Расчет шугообразования и движения кромки ледяного покрова в нижних бьефах ГЭС//Известия ВНИИГ им. Б.Е. Веденеева: Сборник научных трудов. - 1980. - Т. 143. - С. 87 - 91.

2. Трегуб Г.А. Метод расчета длины полыньи в нижних бьефах ГЭС//Матер. конф. и совещ. по гидротехнике: Борьба с ледовыми затруднениями на реках и водохранилищах при строительстве и эксплуатации гидротехнических сооружений/ВНИИГ им. Б.Е. Веденеева. - 1984. - С. 18 - 23.

3. Чижов А.Н. Образование внутриводного льда и формирование шугохода на горных реках//Труды ГГИ. - 1962. - Вып. 93. - С. 3 - 23.

4. Донченко Р.В. Модель процесса замерзания рек//Труды ГГИ. - 1980. - Вып. 270. - С. 3 - 11.

5. Донченко Р.В. Физические свойства внутриводного льда (шуги)//Труды ГГИ. - 1956. - Вып. 55 (109). - С. 5 - 40.

6. Берденников В.П. Физические характеристики льда заторов и зажоров//Труды ГГИ. - 1965. - Вып. 129. - С. 19 - 43.

7. Рекомендации по термическому расчету водохранилищ: П 78-79/ВНИИГ. - Л.: Гидрометеоиздат, 1977. - 30 с.

8. Рекомендации по расчету зажорных явлений в нижних бьефах ГЭС. - Л.: Гидрометеоиздат, 1977. - 30 с.

9. Пехович А.И., Жидких В.М. Расчеты теплового режима твердых тел. - Л.: Энергия, 1976. - 352 с.

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

11. Шаталина И.Н. О расчете коэффициента теплоотдачи при обтекании поверхностей достаточно большой длины//Известия ВНИИГ им. Б.Е. Веденеева. - 1975. - Т. 109 - С. 140 - 148.

12. Лед в водохранилищах и нижних бьефах ГЭС//Я.Л. Готлиб, Р.В. Донченко, А.И. Пехович, И.Н. Соколов. - Л.: Гидрометеоиздат, 1983. - 197 с.

13. Нежиховский Р.А. Расчеты и прогнозы стока шуги и льда в период замерзания рек//Труды ГГИ. - 1963. - Вып. 103. - С. 3 - 40.

14. Чижов А.Н. О расчетах толщины ледяного покрова на реках и водохранилищах//Труды ГГИ. - 1980. - Вып. 270. - С. 40 - 55.

15. Ляпин В.Е., Трегуб Г.А., Разговорова Е.Л. Ледотермический режим нижних бьефов ГЭС и влияние на него тепловых стоков//Гидрофизические процессы в реках и водохранилищах. - М. 1985. - С. 263 - 269.

16. Ляпин В.Е., Трегуб Г.А., Разговорова Е.Л. Методы прогноза и регулирования ледотермических явлений в нижних бьефах высоконапорных ГЭС//Матер. конф. и совещ. по гидротехнике: Инженерное мерзлотоведение в гидротехническом строительстве/ВНИИГ им. Б.Е. Веденеева. - 1984. С. 158 - 163.

17. Крицкий С.Н., Менкель М.Ф., Россинский К.И. Зимний термический режим водохранилищ, рек и каналов. - М.-Л.: Госэнергоиздат, 1947. - 154 с.

18. Назаренко С.Н., Кожевникова Т.Е., Сулимова Л.И. Опыт прогнозирования элементов ледотермического режима нижних бьефов ГЭС//Гидротехническое строительство. - 1980. - N 9. - С. 40 - 43.

19. Руководство по гидрологическим расчетам при проектировании водохранилищ. - Л.: Гидрометеоиздат, 1983. - 283 с.

20. Донченко Р.В. Экспериментальные исследования испарения с поверхности льда//Труды ГГИ. - 1968. - Вып. 159. - С. 56 - 60.

21. Берденников В.П. Динамические условия образования заторов льда на реках//Труды ГГИ. - 1964. - Вып. 110. - С. 3 - 11.