ЭКОЛОГИЯ
J
УДК 551.46:539.219.1:621.375.826
В. Г. Мерзликин
ОЦЕНКА ВЛИЯНИЯ НЕФТЯНЫХ ЗАГРЯЗНЕНИЙ НА ВОДНО-ТЕПЛОВОЙ РЕЖИМ АТЛАНТИЧЕСКОГО ТЕЧЕНИЯ ГОЛЬФСТРИМ
На основе концепции лучистого теплопереноса для деятельного 50-метрового слоя океана, верхний подповерхностный слой которого представляет собой водно-нефтяную эмульсию, рассматриваемую как светорассеивающую и поглощающую среду, исследован водно-тепловой режим поверхностного слоя, загрязненного нефтью. Установлено влияние увеличения потоков солнечного излучения отраженного и поглощенного слоем водно-нефтяной эмульсии, на максимальное значение температуры в приповерхностном слое и снижение прогрева глубинных слоев. Результаты проведенного анализа показывают, что нефтяное загрязнение, играющее роль теплоизоляции, ограничивает поверхностное испарение при штиле, но в штормовую погоду способствует зарождению и усилению ураганов, а процесс диссипации тепловой энергии перегретого поверхностного слоя океана влияет на формирование обезвоженных воздушных потоков, обуславливает дефицит осадков на европейском континенте и возможное похолодание, связанное с понижением средней температуры атлантического течения Гольфстрим. E-mail: [email protected]
Ключевые слова: нефтяные загрязнения, водно-нефтяная эмульсия, океан, климатические аномалии, излучение, отражение, рассеяние, поглощение, теплообмен, температурный максимум, скин-слой, испарение.
Анализ проблемы и постановка задачи. Климат европейского континента в большой степени характеризуется воздействием теплого течения Гольфстрим, определяющим мягкий климат Западной и Центральной Европы в сравнении с восточными регионами. Но проявившиеся с июля 2010 г. аномальные климатические изменения, никогда ранее не регистрировавшиеся за все время инструментальных наблюдений, вынуждают обратить особое внимание на возможные изменения водно-теплового режима знаменитого течения. К погодным катаклизмам можно отнести удивительно жаркое лето на Русской равнине (на фоне необычно холодного в Западной Сибири), рекордно высокие температуры поверхности океана (ТПО или SST по европейской терминологии) в северной Атлантике и последующую необычайно суровую и снежную зиму в Западной Европе 2010/11 г. Такие резкие температурные изменения за короткий период естественными климатическими колебаниями объяснить невозможно, и актуальным становится поиск причин антропогенного воздействия на природу. Так как
районом зарождения Гольфстрима, точнее усиления его энергии после прихода из Африки, является Мексиканский залив, то для многих наблюдателей стало логичным связать погодные аномалии с техногенной катастрофой в Мексиканском заливе — разливом нефти, несмотря на его относительно малую площадь. Во многих случаях антропогенный фактор играет роль малого параметра в неустойчивых термодинамических системах, в частности, природных катаклизмов [1-2]. Наблюдаемое увеличение техногенных катастроф, как правило, связано с интенсивным загрязнением окружающей среды, прежде всего с добычей и транспортом больших объемов углеводородного сырья. Это актуально для России с развитым нефтегазовым промышленным комплексом, нуждающимся в срочной модернизации и обновлении [3], особенно учитывая последствия катастрофической аварии в Мексиканском заливе [4-6]. Отдельные публикации указывают на то, что в настоящий момент до 20% Мирового океана уже было подвержено воздействию разливов нефти и нефтепродуктов. Суммарный выброс за год достигает по разным оценкам не менее 5 млн. тонн. Разлив нефти в Мексиканском заливе уже в августе 2010 г. способствовал появлению масштабных нефтяных загрязнений у восточного побережья США. При ширине Гольфстрима порядка 100-150 км и общей протяженности его маршрута до 5000 км перенос водных масс обусловил площадь загрязнения не менее миллиона квадратных километров, что проиллюстрировано на карте загрязнения поверхности Атлантического океана, полученной NASA (рис. 1, 4-я полоса обложки). Видно, что в центре Северной Атлантики образовалось громадное размываемое нефтяное пятно уже далеко вышедшее за пределы залива [5, 6].
Настоящая работа посвящена анализу крупномасштабных нефтяных загрязнений, их влияния на водно-тепловой режим системы нефтяное загрязнение - вода океана-атмосфера, представляемой многослойной полупрозрачной рассеивающей средой, и связанных с этим возможных климатических изменений.
Физическая модель и принятые допущения. Начиная с 70-х годов прошлого столетия был разработаны ряд комплексных физических моделей светорассеивающих и поглощающих сред, позволивших сформулировать новые математические модели природных процессов [7-9] и объяснить ряд наблюдаемых явлений для полярных и высокогорных ледников, а также решить некоторые прикладные задачи, связанные с созданием огнезащитных и теплоизолирующих материалов и покрытий различного назначения [10-17]. Отличительная особенность отмеченных исследований — изучение сложного теплообмена для оптических неоднородных сред путем совместного решения уравнений теплопроводности и переноса излучения на основе разработан-
Иллюстрации к статье В.Г. Мерзликина «Оценка влияния нефтяных загрязнений на водно-тепловой режим атлантического течения Гольфстрим»
Рис. 1. Карта нефтяного загрязнения Атлантики через четыре месяца после разлива нефти в Мексиканском заливе
Рис. 4. Тепловой баланс и распределение температуры у границы раздела атмосферы с чистым (левая часть) и загрязненным (правая часть) океаном
ных и апробированных оптических моделей материалов в условиях интенсивного конвективно-радиационного воздействия.
Однако в области океанологии и в России, и за рубежом десятилетиями не менялась методология оценки и анализа водно-теплового баланса, температуры поверхности океана и, в частности, интерпретации так называемого холодного скин-слоя [4, 18-23]. Несмотря на прецизионные измерения распределения температуры воды в океане по глубине с точностью до миллиметра, оценка тепловых режимов водных акваторий ограничивалась применением закона сохранения энергии, неполной краевой задачей лучистого теплообмена без детального анализа распространения проникающего солнечного излучения в толще воды. Нередко воздействующий на поверхность океана лучистый поток представлялся только аддитивной величиной потоков проникающего солнечного излучения и поглощаемого поверхностью океана излучения атмосферы. Исследование морской воды как среды, обладающей не только поглощением, но и рассеянием, ограничивалось так называемым параметрическим описанием формирующего теплового баланса [18, 19, 21, 22], что обусловило некорректную трактовку наблюдаемого подповерхностного максимума температуры. В качестве возможных причин называлось влияние конвективного и терморадиационного выхолаживания, турбулентных возмущений, но фактор лучистого объемного нагрева во всей подповерхностной области (включая скин-слой) не стал предметом исследования с позиций оптики рассеивающих сред [4, 18-23]. Здесь сказалась неосведомленность океанологов о подобных исследованиях в других областях науки и техники, в которых были достигнуты значительные результаты. Например, уже сейчас необходимо учитывать эффект формирования подповерхностного температурного максимума в океане при воздействии коллими-рованных потоков проникающего солнечного излучения в оптически неоднородной среде [12, 13].
Хотя оптика океана была предметом многолетних исследований в Институте океанологии им. П.П. Ширшова, в частности, известного советского ученого К.С. Шифрина, тем не менее задача сложного теплообмена на современном уровне в океанологии еще не решена. Это связано с различными причинами, одна из которых — необходимость дифференциации эффектов отражения и переизлучения водной средой — весьма существенна, на что, впрочем, указывал в своих работах сам К.С. Шифрин [24]. Но указанные разработки в России фактически завершились в 90-х гг. прошлого века [1, 23]. За рубежом проблема сложного теплообмена для океана также полностью не решена и не получила должного развития [4, 21, 22] из-за свойственных Западу, несмотря на тщательные инструментальные исследования, нередкого
застоя в теории и недостаточных междисциплинарных связей в науке. На это указывает, например, анализ работ немецких специалистов [4, 15], опубликованных без единой ссылки на работы российских исследователей, опубликованные как в российских, так и зарубежных изданиях.
Подобная ситуация складывалась и при изучении другой природной среды — ледников, в толще которых уже в 60-70-х гг. тоже регистрировался подповерхностный перегрев при сложном теплообмене для снежного слоя [25]. Теоретическое обоснование этого эффекта было получено автором совместно с М.С. Крассом в 1985 г. [8, 9], что было экспериментально подтверждено радиозондированием американскими учеными Ж.Кохом и Р.Йорданом в 1995 г. [26] (опять-таки без ссылок, например, на монографию [8]), а с позиций теории радиационно-кондуктивного теплообмена в многослойной рассеивающей среде проанализировано в работе [27].
Конечно, океан является более сложной средой, в которой необходимо учитывать дополнительно диффузию солей и турбулентное перемешивание, что существенно усложняет анализ процесс нагрева и его моделирование. Но в первом приближении этими явлениями можно пренебречь, чтобы выявить влияние основных процессов лучисто-кондуктивного теплообмена в рассеивающей и поглощающей морской среде. Присутствие какого-либо загрязнения только усиливает ослабление излучения, влияет на отражательную способность (альбедо) и снижает нагрев глубинных горизонтов океана.
В данной работе предлагаются новые физические модели водно-теплового баланса акваторий океана и атмосферных потоков, образующихся над водной поверхностью, загрязненной разливом нефтепродуктов. Загрязнения на поверхности воды рассматриваются как полупрозрачные водно-нефтяные эмульсии (ВНЭ), взвеси (пленки) инородных частиц с различными физическими и геометрическими параметрами, определяющие селективные свойства ослабления водной средой потоков излучения в видимом, ближнем и среднем ИК-диапазонах спектра солнечного излучения и атмосферы в условиях естественной конвекции. В проводимом анализе использованы данные об оптических свойствах морской воды из работ [23, 24].
Объект исследования представляется трехслойной физической моделью, содержащей (см. таблицу) нерассеивающий (показатель рассеяния а = 0) и слабопоглощающий (показатель поглощения к ~ 0,05 м-1) глубоководный слой акватории океана, изолированный полупрозрачным оптически неоднородным водно-нефтяным эмульсионным слоем с прилегающей надводной атмосферой, рассматриваемые как отражающие, поглощающие и излучающие среды.
Таблица
Оптические модели полупрозрачных ВНЭ (с объемной концентрацией частиц нефтепродуктов ^ 0,1 % с показателями рассеяния и поглощения <т « к = 0,05... 2,0 м нефтяной пленки (с коэффициентом пропускания т^ — 0,6 на поверхности океана) и чистой морской воды с затратами теплоты на испарение qv
Оптические модели рассеивающих М(к/<т-А%) и нерассеивающих сред Обозначение соответствующих оптической модели графиков на рис.2 и 3 ге,м 1 ст,м 1 Коэффициент объемного отражения (альбедо) полубесконечного слоя Ню А = г(Н со) Толщина полубесконечного плоского слоя с максимальным объемным коэффициентом отражения А (т < 1 %) Дх>> м Толщина слоя с коэффициентом пропускания г = 10% НТ=од, м
Поглощающая и рассеивающая ВНЭ с частичной экранировкой испарения — 12,5 Вт/м2
М1(1/1 - 27) кривая 1 (рис. 3) 1 1 0,268 1,83 1,28
М12(1/2 - 38) — 2 0,382 1,02 0,97
М2(0, 5/0,5 — 27) кривые 1, 2 и 2, 3 0,5 0,5 0,268 3,43 1,88
М21 (0,5/1 - 38) — 1 0,382 1,91 1,59
М3(0, 25/0,25-27) кривые 3 и 4 0,25 0,25 0,268 10,15 5,19
М4(0, 05/0,05-27) кривые 5 и б 0,05 0,05 0,268 49,52 26,15
Поглощающая и нерассеивающая сплошная нефтяная пленка, препятствующая испарению (qv)fw — 0
М/УК_| кривая 5 (рис. 3) | > 0,05 | 0 _0,268*_|_-_|_
Слабо поглощающая и нерассеивающая морская вода с естественными затратами теплоты на испарение (д„)сш = 25 Вт/м2 МС\¥ I кривые 7 и 4 I 0,05 I 0 I 0,07* I - I 44,53
— поверхностное отражение океана.
Этой модели соответствует задача об экранировании водной акватории теплоизолирующими (теплозащитными или терморегулирую-щими) покровами (покрытиями), состоящими из селективно рассеивающих, поглощающих и излучающих органических веществ. В связи с незначительным числом прямых экспериментальных измерений оптических и терморадиационных характеристик для морской воды и водно-нефтяной эмульсии принята их модельная оценка по данным об альбедо для чистого и загрязненного океанов. По имеющимся в литературе данным указанные величины не превышают соответственно 10% и 20... 30% [4, 19, 23, 24]. Тогда для чистого океана достаточно рассмотреть оптическую модель проникающего солнечного излучения, ослабляемого с глубиной по закону Бугера, с увеличенным граничным отражением за счет вспененной поверхности океана. Оптические модели полупрозрачных водно-нефтяных эмульсий с концентрацией частиц нефтепродуктов ~0,1 % определяются объемным отражением (при незначительном граничном отражении по закону Френеля), вызванным инородными рассеивающими и поглощающими частицами нефтепродуктов с определенной функцией их распределения по размерам, сравнимым с длиной волны проникающего излучения в видимом и ближнем ИК-диапазонах спектра, т.е. до Л ~ 1... 2 мкм. При этом необходимо учитывать и рассеянную компоненту, обусловленную взаимодействием солнечного излучения с флуктуациями плотности, определяемыми возмущениями водной среды. Тогда характеристики ослабления излучения за счет рассеяния и поглощения можно рассматривать как эффективные и использовать как феноменологические параметры для оценки функции поглощенной ВНЭ энергии излучения [8, 13-17, 27].
Принятые оптические характеристики морской воды приведены в таблице: для чистой воды показатели поглощения и рассеяния сравнимы по величине — а ~ к = 0,05 м-1 (0,05... 2,0 м-1 — с учетом загрязнения морской воды инородными частицами, например, ВНЭ) [23]. Сплошная нефтяная пленка характеризуется коэффициентами пропускания Tf = 0,6 и отражения г/ = 0,27. Для исследования теплового режима океана теплофизические характеристики чистой и загрязненной морской воды были приняты одинаковыми: р = 1100 кг/м3; ср = 3700 Дж/(кг-К); КТ = 0,5Вт/м-К.
Подповерхностный слой ВНЭ может частично препятствовать испарению вплоть до полной экранировки сплошной нефтяной пленкой невозмущенной водной поверхности с определенными (модельными) затратами теплоты на испарение. Последние были выбраны из оценок составляющих интегрального теплового баланса системы
атмосфера-земная поверхность. При плотности потока солнечного излучения 750Вт/м2 примерно 500Вт/м2 поглощается и отражается подстилающей поверхностью океана, а оставшаяся треть затрачивается на испарение [4, 18, 19, 22, 23]. Так как затраты на нагрев воды не менее чем на порядок превышают теплоту ее испарения, то в граничных условиях приняты энергетические затраты на испарение для чистой морской воды (¡,и = 25 Вт/м2, а по мере увеличения концентрации ВНЭ до состояния сплошной пленки значение (¡,и стремится к нулю.
Предлагаемые оптические модели классифицируются по выражению М(а/к — А) (см. таблицу). Вводимое значение объемных показателей рассеяния а и поглощения к имеет традиционную в оптике трактовку [8, 10-17] как доля рассеянной (поглощенной) энергии на единицу длины плоскопараллельного слоя ослабляющей излучение среды. В рамках микроскопической теории рассеивающих сред Г. Ми показатели рассеяния (поглощения) представляются интегральными величинами по нормированным сечениям (коэффициентам) рассеяния и поглощения при учете функции распределения полидисперсных частиц с относительным комплексным показателем преломления т = |ш2/ш1| для частиц наполнителя и его водного связующего. Моделированию оптических параметров в зависимости от структуры посвящены работы [8, 16].
Математическая модель и метод решения. Для расчета переноса излучения в поглощающей и рассеивающей среде, соответствующей принятым оптическим моделям чистой воды со слабым рассеянием и загрязненной светорассеивающей водно-нефтяной эмульсией, использовано двухпотоковое приближение для уравнения переноса излучения, согласно которому результирующий поток излучения представлен двумя потоками — в направлениях распространения падающего (дА1) и отраженного (дА2) излучения [8, 16, 23]:
а9А1(х) = — дА1(ж)(аА + кл) +--—— дл2(х); (1)
ах аА + кА
^А2(х) = ^(хХаА + ка)--аА— ?А1(х). (2)
ах аА + кА
В порядке приближения спектральной и температурной зависимостью оптических характеристик пренебрегается, а также не учитывается влияние турбулентного перемешивания верхнего слоя океана. Угловые характеристики рассеянного излучения характеризуются однонаправленной индикатрисой рассеяния с долей рассеяния в направлении, обратном распространяющемуся излучению, не менее 90-100%. Тогда под показателями поглощения и рассеяния понимаются нормированные сечения поглощения и рассеяния, проинтегрированные по
полному телесному углу, но рассматриваемые в одномерном приближении. В рамках двухпотокового приближения использованы следующие выражения для феноменологических параметров, характеризующих оптические свойства среды:
Ь =(к2 + 2ак)1/2 (3)
— для показателя ослабления светорассеивающей среды и
А = (Ь — к)/(Ь + к) (4)
— для альбедо (коэффициента отражения) полубесконечного слоя.
Решение уравнения переноса излучения в двухпотоковом приближении позволило рассчитать характеристики светового поля и коэффициенты отражения г(х), пропускания т(х) и поглощения а(х) плоского оптически неоднородного слоя чистой морской воды и водно-нефтяной эмульсии толщиной х [8, 9, 16]:
г (х) = [1 — ехр(—2Ьх)]А/[1 — А2 ехр(—2Ьх)]; (5)
т(х) = (1 — А2) ехр(—Ьх/[1 — А2 ехр(—2Ьх)]; (6)
а(х) = (1 — А)[1 — ехр(Ьх)]/[1 + А ехр(—Ьх)], (7)
удовлетворяющие условиям нормировки
а(х) + г(х) + т (х) = 1.
Модель для количественных расчетов и оценок температурного состояния рассматриваемой природной системы представляет замкнутую систему уравнений и граничных условий, содержащую:
— нелинейное одномерное дифференциальное уравнение теплопроводности с функцией теплового источника, коэффициенты которого включают в себя ряд оптических параметров морской воды согласно упомянутым выше физическим моделям;
— зависимость поглощения солнечного излучения в воде от глубины для так называемого в океанологии коротковолнового диапазона длин волн;
— сложное нелинейное граничное условие на поверхности воды, учитывающее отражение в коротковолновом (для непрозрачных объектов) и переизлучение в длинноволновом (до 8... 10 мкм) диапазонах спектра с учетом значений коэффициента пропускания солнечного излучения через сплошную нефтяную пленку, а также конвективную составляющую теплового потока.
-Kt— = ат(TA - -) + ee/co(-A - -4) + q0UR + q™ - pLxU, (10)
В общем виде лучистый теплоперенос с известной функцией ^(х, ¿, к, а, д0) поглощенной энергии и энергетическими затратами на фазовые переходы на поверхности облучаемого водно-эмульсионного слоя представляется нелинейным уравнением
дТ д ( дТ \
С^яТ = ТТ Кт-тр + ^ (8)
д£ дх \ дх у
где
^(х) = 1 — А^А-^^Я) {ехр (—Ьх) — ехр [Ь (х — 2Н)]} . (9)
Условие теплообмена на поверхности воды (х = 0) учитывает внешнюю конвекцию, воздействие длинноволновых составляющих потоков излучения океана и атмосферы, а также компонент солнечного излучения, для которых океан является непрозрачной средой:
дТ
дх
где ее/ = (1/е1 + 1/е2 — 1)-1 — эффективный коэффициент черноты для взаимно облучаемой воздушной среды и поверхности океана с коэффициентами черноты е1 и е2 в длинноволновом диапазоне; дО7^, дО7^ — компоненты суммарного потока излучения Солнца д0 соответственно в ближней УФ и ИК областях спектра; ат — коэффициент теплоотдачи для границе раздела с атмосферой; с0 — постоянная Стефана-Больцмана; КТ — коэффициент теплопроводности; р — плотность морской воды; хЦ, — скорость движения границы испарения; — теплота испарения воды.
На тыльной границе плоского слоя чистой морской воды при х = Н = 50 м полагаем
ат
—— = —0,04/Кт. (11)
ах
Температурное распределение в среде в начальный момент времени принимается линейным с градиентом —0,04 град/м и температурой на поверхности 22 °С.
Анализ результатов. Модельные расчеты проводились для нескольких сценариев сочетаний оптических характеристик пленки (эмульсии) на поверхности морской воды, ее оптических и тепло-физических параметров, а также условий теплообмена на границе раздела вода—воздух.
На рис. 2 приведены графики зависимости оптических характеристик — коэффициентов отражения и пропускания — для различных типов ВНЭ и чистого океана.
Для чистой морской воды имеет место ослабление проникающего солнечного излучения по закону Бугера (кривая 7); отражательная
О 20 40 60 х, м
Рис.2. Зависимость коэффициентов пропускания т(х) и отражения г(х) от толщины х подповерхностного слоя океана для оптических моделей (см. таблицу):
1,2 - М2(0,5/0,5 - 27); 3,4 - М3(0,25/0,25 - 27); 5,6 - М4(0,05/0,05 - 27); 7, 8 - МСШ
способность границы раздела двух сред определяется постоянной величиной 8... 10%, в которой не более 2-5 % (в зависимости от угла падения) соответствует френелевскому отражению, а остаток — отражению от вспененной поверхности океана. Это поверхностное отражение не зависит от толщины слоя вещества (кривая 8). Если в такой воде имеется рассеивающая примесь с показателем рассеяния, равным показателю поглощения, К4 = а4 « кош = 0,05 м-1, то отражение уже является функцией толщины, возрастает до ^0,27 на глубине ~0,2 м (так называемая толщина полубесконечного слоя) и с ростом толщины больше не изменяется. Для других вариантов рассеивающих слоев ВНЭ за счет подбора показателей поглощения и рассеяния удается изменить поглощение в слое морской воды, но альбедо толстого полубесконечного слоя в данных расчетах одинаковое, характерное для загрязненного океана.
На рис. 3 приведены расчетные распределения температуры по глубине, формируемые в толще морской воды на третьи сутки = 24 ч) в астрономические полдень % = 2,25 • %&ау и полночь % = 2,75 • %&ау при синфазных колебаниях температуры атмосферы от 17 до 27 0С и потока солнечного излучения от 500Вт/м2 до 0. Проведены расчетные оценки температурных распределений в толще океана, загрязненного
т °с
25
20
30
25
20
1 1 1 1
f\l
1 3 \
fr 2 vk
f'.5 1 ••.......
/ ^ i i 1 1
0 1 2 3 4 х, а
т,°с
1 1 1 1 1
! ?
// <К .....
- /iiii
0
0.2
0.4
0.6
м
0.8 х, м
Рис.3. Расчетные нестационарные распределения температуры загрязненной и чистой морской воды в штиль при радиационно-конвективном воздействии для 50-метровой толщи в мелком (а) и укрупненном (б) масштабах на третьи сутки в астрономические временные периоды:
полночь (рп = 2,75 • ¿¿ау) — для оптических моделей загрязненного ВНЭ М 1(1/1 — 27) (кривая 1) и М2(0,5/0,5 — 27) (кривая 3); моделей чистого океана МСЖ (кривая 4) и покрытого нефтяной пленкой М/Ж (кривая 5); полдень (¿а = 2,25 • ¿¿ау) — для модели М2(0, 5/0, 5 — 27) (кривая 2) с заданным начальным распределением температуры (кривая 6)
ВНЭ (кривые 1-3); чистой морской воды (кривая 4) и покрытой тонкой нефтяной пленкой (кривая 5).
На графиках представлены температурные распределения для двух оптических моделей ВНЭ. Для сильно поглощающей ВНЭ с показате-
лем поглощения к = 1 м-1 оптической модели М 1(1/1 — 27) рассчитан температурный профиль Т(х, 2,75 • Ьё,ау) (кривая 1). Для слабопогло-щающей ВНЭ (к = 0,5 м-1) оптической модели М2(0,5/0,5 — 27) рассчитаны температурные профили Т(х, 2,25 • Ьё,ау) и Т(х, 2,75 • Ыау) (кривые 2 и 3).
Для оптической модели М2(0,5/0,5 — 27) наибольшая температура «24 0С поверхности (х = 0) в полдень характеризуется температурным распределением Т(х, 2,25%^ау) (кривая 2) с присутствием так называемого холодного скин-слоя (толщиной до 0,2 м) относительно подповерхностного максимума температуры «26 0С. Для той же модели температура поверхности океана ТПО (ББТ) снижается до «21 0С на температурном профиле Т(х, 2,75 • Ьйау) (кривая 3) к полуночи % = 2,75 сут. Толщина скин-слоя практически не меняется, но сам скин-слой более холодный (на несколько градусов), а подповерхностный температурный максимум даже возрастает примерно на один градус до «27 0С за счет стока к поверхности океана кондуктивного потока теплоты, вызванного ночным выхолаживанием за счет излучения в длинноволновой области спектра и конвекции, что полностью соответствует экспериментальным данным [4, 18, 19, 23].
Рост поглощения для ВНЭ модели М 1(1/1 — 27) приводит к более существенному росту температуры (кривая 1) и перегреву в подповерхностной области в сравнении с ТПО модели М2(0,5/0,5 — 27) (кривая 2). Но проникающий в толщу более загрязненной морской воды поток излучения уменьшается сильнее в сравнении с ВНЭ модели М2(0,5/0,5 — 27) с меньшим ослаблением излучения, чем в модели М 1(1/1 — 27), что приводит к снижению температуры уже на глубине ~1 м. Таким образом, для более чистого океана температура на поверхности будет ниже, а прогрев глубинных слоев — сильнее.
В отличие от традиционной океанологии, где не рассматривалось влияние на температурный режим одновременно и поглощения и рассеяния, предлагаемый подход позволяет по-новому рассмотреть физический механизм нагрева океана, что имеет принципиальное значение. Проведенные расчеты показывают, что усиление рассеяния вызывает смещение максимума температурного профиля к поверхности при слабом изменении самой температуры, а рост поглощения способствует ее увеличению.
Суточные колебания температуры атмосферы и солнечного потока обуславливают осцилляции температуры на поверхности океана, которые меняют не только значение градиента температуры, но и его знак. К ночи имеет место сток кондуктивный сток теплоты к поверхности из перегретой поверхностной области, а к полудню наблюдается естественная инверсия температурных профилей Т(х, 2,25 • Ьдлуу) и
Т(х, 2,75 • Ыау) (модель М2(0,5/0,5 — 27), кривые 2 и 3) в подповерхностной области.
Температурный профиль для толщи морской воды, теплоизолированной сплошной нефтяной пленкой, оптической модели М/№ качественно повторяет структуру температурного поля для ВНЭ, но при слабом рассеянии нижележащих горизонтов чистой морской воды обуславливает меньшие температуры разогрева подповерхностной области и ТПО (БЭТ).
В случае сплошной полупрозрачной нефтяной пленки (кривая 5) имеет место рост поглощенного (до 30%) и отраженного (до 60%) потоков солнечного излучения; при этом поток солнечного излучения, воздействующий на глубинные слои морской воды, уменьшается до «200 Вт/м2 (что существенно меньше поглощенного потока для чистой воды « 470 Вт/м2). Для этого случая характерен интенсивный прогрев океана вблизи поверхности. При загрязнении океана ВНЭ (кривые 1, 2) температура морской воды может достигать 30 °С после трехсуточного прогрева в штиль, причем пропорционально росту не только поглощения, но и рассеяния. Для водно-нефтяной эмульсии и нефтяной пленки, частично экранирующих проникающий поток солнечного излучения, на глубинах х > 3 ... 5 м температура океана (см. рис.3,а, кривые 1-3, 5 меньше, чем в случае прогрева чистой воды (Т(х) < Т0(х)). Таким образом, имеет место снижение прогрева глубинных слоев океана на фоне перегрева подповерхностной области и роста ТПО (БЭТ).
Выводы. При длительном нагреве океана солнечным излучением температура воздуха также возрастает за счет увеличенного альбедо загрязненной воды и конвективного стока теплоты в атмосферу. Как следует из полученных результатов, если предлагаемая оптическая модель верна, то зона аномального прогрева загрязненной морской воды (и зарождения антициклонов) не только активно расширялась, но и продвигалась на восток вместе с Гольфстримом (см. рис. 1, 4-я полоса обложки). Это могло привести к формированию над обширной поверхностью океана обезвоженных конвективных потоков — воздушных атмосферных течений.
Гольфстрим — это источник огромного количества влаги (за счет испарения) и теплоты для атмосферы. В зоне его действия можно выделить так называемые энергоактивные зоны океана, в которых суммарное выделение теплоты в атмосферу на порядок превышает среднеширотные оценки. Еще в советские годы академик Г.ИМарчук сформулировал концепцию, согласно которой именно эти зоны играют ключевую роль в проблеме долгосрочного прогноза погоды и корот-копериодных колебаниях климата [1, 2, 5].
На рис.4 (см. 4-ю полосу обложки) показана схема глобального перераспределения тепловых потоков в атмосфере и океане при масштабном техногенном загрязнении нефтью водной поверхности. Основной вклад в нагрев океана дает поток солнечного излучения (желтые стрелки), до 95 % которого поглощается в толще чистой морской воды до глубины 60... 80 м [4, 18, 23]. Процессы конвективного (синие стрелки) и терморадиационного теплообмена в длинноволновом диапазоне (красные стрелки) в течение суток или подогревают океан, или его выхолаживают, что вызывает изменение температуры поверхности и появление скин-слоя, изменяющегося по размерам и перепадам температуры. Вклад терморадиации и конвекции в зависимости от времени суток и температуры океана и атмосферы может изменяться в пределах до ±30 %.
Для загрязненной морской воды (см. правую часть рис. 4) имеет место существенный рост рассеяния и поглощения солнечного излучения в приповерхностном слое. Это обуславливает как увеличение отраженного (поверхностью и объемом водной толщи) солнечного потока (до 30%), так поглощенного (до 35 % в загрязненной подповерхностной области). Но вглубь океана распространяется не более 60 % солнечного потока, вследствие чего прогрев морской воды на глубине значительно уменьшается. Таким образом, появление нефтяного загрязнения на поверхности океана приводит к кардинальной перестройке теплового режима системы океан-атмосфера, а именно:
— происходит дневное экранирование солнечного прогрева морской воды по глубине и значительный прогрев ее поверхностного слоя при интенсивном терморадиационном и конвективном выхолаживании (особенно в ночные часы), усиливающемся в штормовую погоду;
— рост температуры по глубине морской воды с загрязненной поверхностью уменьшается, что приводит к формированию более холодной водной толщи в океане;
— интенсивность испарения воды с загрязненной поверхности существенно снижается, резко возрастая при сильном волнении моря и штормовой погоде;
— имеет место дополнительный нагрев атмосферы за счет увеличенного отражения солнечной радиации от загрязненной и перегретой воды, а также роста терморадиационных и конвективных потоков.
За последние 10-20 лет оценки водно-теплового баланса отдельных участков мирового океана и суши показывают, что катастрофические потери или существенный избыток среднемесячных осадков не согласуются с традиционными моделями колебаний климата и погодных условий как результатов действия известных циклонов или антициклонов [1, 2, 22, 23]. На эти несоответствия еще указывал известный
русский ученый А.С. Монин: "... Расшифровка данных об изменениях климата все в возрастающей мере требует использования физических методов" [2].
Могут реализоваться два сценария циркуляции конвективных воздушных потоков: с дефицитом или избытком влаги, что обусловит либо продолжительную засуху на пространстве распространения обезвоженных воздушных масс, либо в случае возникновения циклона — усиление его энергии и подпитку водной массой.
В первом случае перегретые пограничные слои океана (приповерхностный воздушный и поверхностный водный) могут быть причиной формирования в атмосфере специфических аномально перегретых конвективных потоков, что в свою очередь способствует поддержанию очага тепловыделения в средней и верхней тропосфере (термодинамические условия зарождения антициклонов). При этом отрицательный водный баланс на морских акваториях может вызываться заторможенным механизмом испарения воды океана или внутренних морей.
Во втором случае перегретый слой воды в результате теплообмена с атмосферой постепенно охлаждается и при зарождении урагана отдает ему свою массу, что существенно усиливает мощность урагана. Происходит интенсивная диссипация (рассеяние, сброс) энергии перегретой воды с поверхности океана посредством ураганов, набирающих силу над акваторией океана, которые увеличивают свою кинетическую энергию и уносят громадную массу воды из атмосферы, в частности из Атлантики, на американский континент. Усиление ураганов по частоте и энергии сопровождается сбросами из атмосферы значительного объема водной массы, которая в обычных условиях достается другим регионам суши в виде циклонов.
Если предлагаемая модель получит свое подтверждение (при будущем проведении экспериментальных измерений оптических характеристик и температурного режима), то, значит, приход антициклонов на равнинные и удаленные от моря и гор территории суши Европы действительно способствовал установлению там аномально жаркой погоды — антициклональной устойчивой во времени и пространстве зоны с высоким атмосферным давлением, препятствующим продвижению циклонов из северной части Атлантики и Ледовитого океана. Именно здесь наблюдались (в том числе и в России летом 2010 г.) жаркая погода и летняя засуха.
В свою очередь, "запертые" циклоны из северной Атлантики и Ледовитого океана один за другим за короткое время обрушили все запасы своей влаги на ограниченную площадь от Западной Европы до Среднего Востока, вызвав там разрушительные наводнения.
Перегрев верхнего слоя морской воды под нефтяной пленкой ведет к интенсивному оттоку теплоты из воды в атмосферу, вследствие чего происходит охлаждение нижележащей водной толщи океана. Охлаждение океана приведет уже в ближайшее время к значительному уменьшению испарения воды в Атлантике.
В обычных условиях это уменьшение составляет около 120 тыс. км3/год; теперь, судя по приближенным оценкам, оно может существенно понизиться. Более охлажденные воздушные потоки не смогут достичь восточных регионов Европы, что будет вызывать в Западной Европе обильные дожди и снежные циклоны, а в Восточной Европе дефицит влаги летом и малоснежные зимы, в период которых вся надежда будет на аккумулирующие осадки с северо-восточных регионов. С высокой степенью вероятности это повлечет за собой дефицит осадков и недобор осенней влаги в почве в Восточной Европе. Уже сейчас наблюдается катастрофическое обмеление рек в России (и даже Волги). При уже образовавшемся летнем дефиците влаги в почве это неминуемо приведет к повышению среднелетних температур воздуха и засухам в будущие летние периоды, а также к малоснежным и суровым зимам на обширных территориях Европы.
Дефицит теплоты в глубине океана вдоль по Гольфстриму может привести к понижению его температуры. В обычных условиях Гольфстрим, как и океан, прогревается почти равномерно до глубины 80... 100 м (так называемый квазиоднородный слой), и приносимый им огромный объем нагретой воды в северную Атлантику и западный сектор Ледовитого океана формирует там относительно мягкий климат. Охлаждение Гольфстрима чревато уже в ближайшем будущем понижением среднезимних температур воздуха в странах Западной Европы и Скандинавии. Может быть, даже в течение нескольких лет зимы будут не холоднее обычного, но с обильными осадками в Западной Европе и будут малоснежными в ее континентальных районах. Это может способствовать увеличению ледового покрова на севере РФ в обычно незамерзающем Баренцевом море и изменить структуру эколого-экономических характеристик обширных регионов.
Автор признателен профессору М.С. Крассу и доценту О.В. Сидорову за полезное обсуждение различных аспектов настоящей работы.
СПИСОК ЛИТЕРАТУРЫ
1. НигматулинР. И. Океан: климат, ресурсы, природные катастрофы // Вестник РАН. - 2010. - Т. 80. - № 8. - С. 675-693.
2. М о н и н А. С. Климат как проблема физики // Успехи физических наук. -2000. - Т. 70. - № 4. - С. 419-445.
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
Стратегический прогноз изменений климата Российской Федерации на период до 2010-2015 гг. и их влияния на отрасли экономики России. - М.: Федеральная служба по гидрометеорологии и контролю окружающей среды (Росгидромет). 2006.
KeenlysideN.TsuangB.-J. Final report: Scale interactions in a coupled climate model / Advancing understanding of the upper ocean diurnal cycle and its relevance to climate. Joint DFG-NSC project: 446 TAI 113/33/0-1. http: // www.ifm-geomar.de/fileadmin/personal /fb1/me/nkeenlyside/Reports/DFG_finalreport.pdf/ Малинин В. Н. Разлив нефти в Мексиканском заливе. РГГМУ. - Санкт-Петербург. 2010. http://online812.ru/2010/09/28/012/pdf/
N O A A (National Oceanic @ Atmospheric Administration) / ESRL (Earth System Research Laboratory) Physical Science Division. http://blog.agu.org/geospace/2010/08/25/blocked-up- weather/ К р а с с М. С. Математические модели и численное моделирование в гляциологии. - М.: МГУ, 1981. - 140 с.
К р а с с М. С.,МерзликинВ. Г. Радиационная теплофизика снега и льда.
- Л.: Гидрометеоиздат, 1990. - 261 с.
К р а с с М. С., Геворкян С. Г., М е р з л и к и н В. Г., Товстоног В. А. Методика расчета радиационных и температурных полей снежных и ледяных массивов // Материалы гляциологических исследований. - 2000. - Вып. 90. - С. 142-147.
Howe J. T., Green M. J., Weston K. C. Сублимирующие объемно отражающие материалы как тепловая защита при совместном воздействии конвективного и интенсивного лучистого потоков тепла // AAIA Journal. - 1973. -Vol. 11, no. 7. NASA Ames research center. Moffet Field, Calif. USA. Dabby F. W. High intensity laser induced vaporization and explosion of solid material // IEEE Journal of Quantum Electronics. - 1972. - Vol. QE-8, no. 8. USA. Петров В. П., Рябцев Е. И., Сутугин В. Г., Мерзликин В. Г. Оптические свойства высокоотражающей керамики // Тематический сб. Вопросы авиационной науки и техники. Сер. Aвиационные материалы. - 1989. -С. 39-43.
Товстоног В. А. Анализ теплообмена в светорассеивающих материалах, нагреваемых излучением // Физика и химия обработки материалов. - 1985. -№ 3. -C. 35-40.
Товстоног В. А. Теплофизика рассеивающих материалов: прикладные проблемы и решения // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. - 2000. - № 3. - С. 67-85.
Manara J., Arduini-Schuster M., R a t z e r - S c h e i b e H.-J., S c h u l z U. Infrared-optical properties and heat transfer coefficients of semitransparent thermal barrier coatings // Surface and Coatings Technology. - 2009.
- Vol. 203. - Issue 8, 15 January. - P. 1059-1068.
Merzlikin V., Gutierrez Ojeda M., Sidorov O., TimoninV. New selectively absorbing and scattering ^at-insulating coatings of the combustion chamber for the low-heat-rejection diesel. SAE Techn. Pap. Ser., No. 07M-171, 2007. S i e g e l R. Internal Radiation Effects in Zirconia Thermal Barrier Coatings // AIAA J. Thermophysics Heat Trans. - 1996. - Vol. 10, no. 4. - P. 707-709. Федоров К. Н., Гинзбург А. И. Приповерхностный слой океана. - Л.: Гидрометеоиздат, 1988. - 304 с.
Ч а в р о А. И. Физические основы и методы определения температуры поверхности океана. - М., 1990. - 242 с.
ВершинскийН. В., Волков Ю. А., Соловьев А. В. O вертикальной микроструктуре тонкого поверхностного слоя океана // ДАН СССР. - 1981. -Т. 256. -№ 3. - С. 694-698.
W e 11 s A. J., C e n e d e s e C., F a r r a r J. T., Z a p p a C. J. Variations in ocean surface temperature due to near-surface flow: straining the cool skin layer // J. of Phys. Oceanography. - 2009. Vol 39. - P. 2685-2710.
22. G r a s s l H. The dependence of the measured cool skin of ocean on wind stress and total heat flux // Boundary-Layer Meteorology. - 1976. - Vol. 10. No. 4. - P. 465-474.
23. ДоронинЮ. П. Физика океана. СПб., 2000. - 287 с.
24. Шифрин К. С., И о н и н а С. Н. Тепловое излучение и отражение волнующейся поверхности моря в микроволновой области // Труды ГГО. - 1968. -Вып. 222. - С. 22-48.
25. А р т е м ь е в А. Н. Взаимодействие атмосферы и подстилающей поверхности на антарктическом плато // Тр. САЭ, 1976. - Т. 66. - 71 с.
26. K o h G., Jordan R. Sub-surface melting in seasonal snow cover // J. of Glaciology. - 1995. - Vol. 41. No. 139. - P. 474-482.
27. ТовстоногВ. А.,Мосалов Ф. Ф., М е р з л и к и н В. Г. Постановка и решение задач радиационно-кондуктивного теплообмена в многослойных рассеивающих средах // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. -2008. -№ 1. -C. 12-29.
Статья поступила в редакцию 17.01.2011
Владимир Гаврилович Мерзликин родился в 1952 г., окончил Московский инженерно-физический институт в 1975 г. Канд. техн. наук, доцент МГТУ "МАМИ". Автор более 70 научных работ в области теплообмена излучением, тепловой защиты, процессов теплообмена в природной среде.
V.G. Merzlikin (b. 1952) graduated from the Moscow Institute for Engineering and Physics in 1975. Ph. D. (Eng.), assoc. professor of the Moscow State Technical University "MAMI". Author of more than 70 publications in the field of study of radiation heat exchange thermal protection, heat exchange processes in nature.