CC BY
375
Ключевые слова:
морская и
атмосферная
циркуляция,
ветровое волнение,
морфо- и
литодинамические
процессы,
численное
моделирование.
Keywords:
sea and atmospheric circulation, wind waves, morpho- and lithodinamic processes,
numerical simulation.
УДК 551.5(571.121)
Н.А. Дианский, И.М. Кабатченко, В.В. Фомин, В.В. Архипов, А.С. Цвецинский
Моделирование гидрометеорологических характеристик в Карском и Печорском морях и расчет наносов у западного побережья полуострова Ямал
При проектировании и строительстве гидротехнических сооружений в береговой зоне необходима информация как о динамике поперечного профиля берега, так и о характеристике вдольберегового потока наносов на достаточно большом протяжении берега, поскольку вмешательство в природную среду может нарушить ее равновесное состояние и привести к катастрофическому изменению режима береговой зоны, а также к негативным последствиям при эксплуатации гидротехнических сооружений.
Основное воздействие на морфо- и литодинамические процессы и трансформацию рельефа береговой зоны оказывают ветровые волны и течения [1]. Природа волн и течений определяется метеорологическими факторами, в первую очередь штормовыми ситуациями. Картина штормовых деформаций дна и берегов может быть чрезвычайно многообразной вследствие сложности и нелинейности самой системы «волны - течения - транспорт наносов» [1]. Однако благодаря установлению важных количественных взаимосвязей внутри этой системы в настоящее время становится реальным количественный прогноз морфодинамических изменений конкретных участков побережий, основанный на численном моделировании штормовых ситуаций. При этом речь может идти о прогнозе не только в масштабе времени отдельного шторма, но и в более длительной перспективе, измеряемой годами и десятилетиями.
В свою очередь, моделирование штормовых ситуаций подразумевает комплексное воспроизведение атмосферной и морской циркуляции, а также ветрового волнения с помощью соответствующих численных моделей. Именно такому подходу посвящена настоящая статья. Полученные результаты воспроизведения течений и волнения используются для моделирования морфо- и литодинамических процессов в береговой зоне в районе Харасавэйского месторождения (далее - расчетная область).
Методика моделирования гидрометеорологических характеристик для расчета наносов у западного побережья полуострова Ямал
Общая методика расчета транспорта наносов подразумевает решение ряда задач, каждая из которых сама по себе является сложной гидродинамической проблемой. Расчет транспорта наносов базируется на методиках, для которых необходим расчет течений и полей ветрового волнения в расчетной области [1], для чего, в свою очередь, требуется воспроизведение атмосферного воздействия. В настоящей работе для расчета гидрометеорологических параметров используется вычислительный комплекс хорошо апробированных моделей. Расчет атмосферного воздействия проводится по модели WRF (англ. Weather Research and Forecasting model); расчет течений, уровня, температуры, солености моря и морского льда - по модели INMOM (англ. Institute Numerical Mathematics Ocean Model) [2, 3]; расчет параметров волнения - по Российской атмосферно-волновой модели (РАВМ) [4]. Описание этого комплекса моделей и его верификация применительно к Карскому и Печорскому морям представлены в работе [3].
Блок гидрометеорологического моделирования отличается от такового из работы [3] использованием ст-координатной версии INMOM [2]. Решение итоговой задачи расчета наносов в акватории расчетной области у западного побережья Ямала
№ 2 (22) / 2015
Современные подходы и перспективные технологии в проектах освоения нефтегазовых месторождений российского шельфа
99
требует исходных данных о течениях и характеристиках ветрового волнения с высоким пространственным разрешением. Для расчета течений с требуемой точностью применена версия INMOM для юго-западной части Карского моря и Печорского моря, реализованная на неравномерной сетке с детализацией в акватории Байдарацкой губы, получающейся в модельной сферической системе координат путем размещения ее полюса в точке с географическими координатами 71,151791° с.ш. и 67,198387° в.д. Таким образом, шаги по пространству составляют порядка 50 м в расчетной области и доходят до 5-10 км на периферии. По вертикали задаются 20 неравномерно распределенных по глубине уровней.
Общий размер сеточной области (рис. 1) составляет в плане 960*295 точек по модельным долготе и широте соответственно. Поскольку используется математически точная сферическая система координат, то никаких топологических искажений при расчете полей течений не вносится, в отличие от зачастую используемых декартовых систем координат, «натягиваемых» на сферичную поверхность Земли. Как видно, «жидкие границы» сеточной области находятся далеко от расчетной области. Описание способа задания условий на «жидких границах» приведено в [3]. В такой постановке течения через Карские Ворота, важные
для формирования циркуляции в Байдарацкой губе, будут автоматически рассчитываться в модели INMOM. Для расчета приливных течений на «жидких границах» ставится условие изменения уровня под действием крупномасштабной приливной волны.
Как уже было сказано, подробная верификация блока расчета гидрометеорологических характеристик приведена в работе [3], поэтому здесь качество расчетов демонстрироваться не будет. Отметим только, что качество воспроизведения гидрометеорологических характеристик отвечает уровню современных возможностей численного моделирования.
Таким образом, для расчета донных отложений и транспорта взвешенных наносов необходимо выполнить следующие виды численного моделирования: 1) расчет атмосферного воздействия по модели WRF на грубой атмосферной сетке; 2) расчет параметров волнения на глубокой воде на атмосферной сетке по РАВМ;
3) интерполяцию параметров глубоководного волнения на мелкую сетку основной модели течений INMOM и их трансформацию в параметры мелководного волнения согласно топографии дна; 4) расчет параметров напряжения трения, вызванного обрушением волн на малых глубинах; 5) расчет течений по INMOM с учетом приливов и предписанного атмосферного воздействия (см. п. 1); 6) расчет формирования
Рис. 1. Сеточная область версии INMOM для Печорского моря с детализацией в акватории Байдарацкой губы:
слева - область моделирования морской циркуляции целиком (показаны каждая 10-я точка по долготе и каждая 3-я точка по широте); справа - детализация сеточной области в расчетной области (показаны каждая 10-я точка по долготе и каждая 3-я точка по широте)
№ 2 (22 ) / 2015
100
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
пространственной структуры и среднемноголетних объемных величин потоков илистых наносов в толще морских вод; 7) расчет потоков донных песчаных наносов по методике [1] с использованием рассчитанных полей течений и ветрового волнения. Пункты 6 и 7 реализуются на сеточной области INMOM с использованием необходимых разностных операторов, уже реализованных в ней.
Метод расчета литодинамических процессов
В прибрежных акваториях полуострова Ямал в безледный период происходит перенос большого количества наносов. Эти наносы существенно различаются как по фракционному составу, так и по происхождению. Можно выделить 2 основных типа подобных наносов:
1) влекомые и сальтирующие наносы, транспортируемые в основном в придонном слое моря;
2) взвешенные наносы, которые транспортируются во всей толще воды.
Экспериментальные исследования, выполненные в 2011-2012 гг., показали, что в грунтах доминирует мелкозернистый песок со средним диаметром частиц ~ 0,16-0,17 мм. В то же время экспедиционные данные показали и наличие взвешенных илистых наносов, обязанных своим происхождением интенсивному штормовому разрушению береговых обрывов Карского моря. Средний линейный размер взвешенных частиц по результатам измерений гранулометрического состава взвешенных наносов на акватории Харасавэйского полигона во время экспедиции 2012 г. составил 0,0011-0,0012 мм.
В связи с тем что физические закономерности транспорта наносов в придонном слое моря и транспорта взвешенных наносов в толще воды отличны друг от друга, то их моделировали по двум различным методикам, которые приведены далее.
Методика расчета транспорта донных наносов. Для расчета потока придонных наносов использовалась энергетическая концепция Р. Бэгнольда [1]. В основу расчетов положено уравнение сохранение массы
dd=^х+^*1. (1)
dt дх ду
где qx - вдольбереговой поток наносов; qy - поперечный поток наносов; d - изменяющаяся во времени t толщина слоя наносов. В данном слу-
чае использовался алгоритм расчета вектора q вдольбереговых и поперечных потоков наносов [1]:
(2)
где е4 и es - отношение мощностей, поглощаемых механизмами донного трения, к мощностям, расходуемым потоком на влечение песка по дну и перемещение сальтирующих частиц, соответственно; р - плотность воды; fW - коэффициент трения для волн; um - амплитуда горизонтальной скорости частиц; Ф - угол естественного откоса песчаного грунта (в наших условиях tgФ ~ 0,6 [1]); WS - скорость осаждения частиц в воде (WS = 0,155ds - 0,0075[см/с] [1]), где ds - гидравлическая крупность песка, мм; В - добавочная мощность, связанная с процессом обрушения волн, благодаря которому турбулентность в водной толще и придонном слое усиливается и количество взвешенных частиц значительно возрастает; u - полная переносная скорость песчаных наносов у дна.
При введении обозначений
f=-3- fw p«m>
3п
G = 9п £4 F 8 tg<& um
+ -^- (4F + B),
Ws '
(3)
где F - потери мощности на донное трение; G - погруженный вес движущихся песчаных зерен на единицу площади, (2) примет вид:
q = Ga. (4)
Полная переносная скорость у дна u в общем случае представляет собой сумму скорости воды в придонном слое U5 и переносной скорости UB, обусловленной турбулентными механизмами за счет обрушения волн, действующими в донном пограничном слое:
u = U5+ U в. (5)
Разложив вектор расхода наносов в формуле (4) на поперечную qx = qcos© и продольную qy = qsin© составляющие, можно записать:
qx = G(Ub + UB),
qy = G(v8 + Vg). (6)
№ 2 (22) / 2015
Современные подходы и перспективные технологии в проектах освоения нефтегазовых месторождений российского шельфа
101
Соответствующие компоненты скоростей рассчитывались следующим образом:
(Us, V)
(Qx, Qy)
h\ ln—-1
(7)
D* - D u2
(UB, VB) = a——*--(cos 0, sin 0),
EC
D* =---£- cos 0,
X
(8)
здесь Qx и Qy - полные потоки из модели течений (9); h - глубина; za - шероховатость; D - скорость диссипации; D* - пороговая скорость диссипации; X - расстояние от берега; E - полная энергия волн; Cg - групповая скорость; для нахождения um, cos© и sin© используется ветроволновая модель РАВМ.
Методика расчета транспорта взвешенных наносов. Составляющие удельного расхода (транспорта) взвешенных наносов определяются из следующих выражений:
си, ц - молекулярная вязкость воды); D = Н + Z (H - глубина места в невозмущенном состоянии; Z - отклонение уровня моря от невозмущенного состояния); v - вертикальный коэффициент турбулентной диффузии, который берется таким же, как и для температуры и солености в INMOM; F - оператор боковой диффузии; qc(x, у, о, t) - источники взвешенных частиц, в общем случае меняющиеся в пространстве и времени. При численной реализации (10) используется монотонная схема переноса, не дающая ложных экстремумов [5].
Уравнение (10) описывает изменение концентрации взвешенных частиц в условиях, учитывающих перенос примеси течениями, горизонтальной и вертикальной турбулентной диффузией, а также физические характеристики частиц взвеси. Для решения уравнения (10) применительно к поверхности ставится условие отсутствия потока взвешенных частиц, применительно к дну - оседания/взвеши-вания частиц:
D D
Qx = | ucdz, Qy = I vcdz,
0 0
(9)
V dc D do
c\ o=1
(11)
где с - объемная концентрация взвешенных наносов; u, v - горизонтальные составляющие скорости течения, вычисляемые в трехмерной термодинамической модели; D - полная глубина.
Объемная концентрация взвешенных наносов вычисляется при помощи уравнения переноса-диффузии взвешенных частиц, которое в о-координатной системе INMOM имеет вид:
dDc 1 -----1--
dt rr
X L
d(rn+rnf )c do
— (r Duc) + — (r Dvc) dx y dy x
d V dc do D do
+ Fc + qc (x, y, o, t), (10)
где x, у, о, t - координатная система; с - концентрация взвешенных частиц; u, v - горизонтальные компоненты скорости течения; ю - вертикальная составляющая скорости в о-системе координат; rx ry - метрические коэффициенты сферической системы координат со сме-
щенными полюсами; =
1 gd 2(рс -Р) 18 ц
гра-
витационная вертикальная скорость оседания взвешенных частиц (скорость седиментации или гидравлическая крупность), рассчитываемая по формуле Стокса (d - средний диаметр частиц взвеси, рс - плотность частиц взве-
где fd - параметр оседания/взвешивания.
На боковых границах задается условие отсутствия потока взвешенных частиц, за исключением участков береговой границы, с которых осуществляется вынос взвеси. На этих участках ставится следующее условие (12) для потока концентрации взвешенных веществ:
sh2
Чс = a X&W, (12)
где а - безразмерный эмпирический коэффициент; g - ускорение свободного падения; h - высота волны; ros - гидравлическая крупность; Т - период волны; £, - суммарное изменение уровня (волновое и сгонно-нагонное); V - модуль вдольбереговой скорости. Функция
Хэвисайда %(£)
0, £ < 0 £ £> 0
в формуле (12) ис-
пользуется для того, чтобы учесть только нагонный суммарный уровень: как только уровень становится меньше 0, то размыв прекращается.
Верификация расчета транспорта взвешенных наносов необходима прежде всего для того, чтобы выбрать в формуле (12) значение безразмерного эмпирического коэффициента a. Последний обычно берется в пределах 10-11-10-13. Для верификации использовались
№ 2 (22 ) / 2015
102
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
Рис. 2. Положение створов отбора проб (а) и сравнение наблюденных (красная линия) и модельных (синяя линия) придонных профилей мутности, мг/л, на створах: б - 1; в - 5; г - 11
б
в
а
г
профили мутности, построенные по данным, полученным 13-14.09.2011 г. на 15 створах (рис. 2, см. а).
Для сравнения результатов расчета мутности с данными наблюдений (см. рис. 2б-г) данные расчетов показаны с суточным осреднением за 14 сентября 2011 г. Согласно рис. 2 модель транспорта взвешенных наносов достоверно воспроизводит измеренные профили мутности.
Расчет структуры донных наносов и величин потоков взвешенных наносов
Как уже было показано, в прибрежных акваториях полуострова Ямал в безледный период происходит перенос наносов разного гранулометрического состава, для моделирования перемещения которых применяются различные методики (см. ранее). В соответствии с ними проведены расчеты для безледных (август -октябрь) периодов 2009-2012 гг. (рис. 3, 4).
№ 2 (22) / 2015
Современные подходы и перспективные технологии в проектах освоения нефтегазовых месторождений российского шельфа
103
Следует отметить, что для расчета расходов для 2010 г. и 2012 г. в модели введен учет влекомых наносов, в то время как для 2009 г. и 2011 г. при расчетах учитывалось только перемещение наносов путем сальтации. Поэтому рассчитанные вдольбереговые расходы песчаных наносов для 2009 г. и 2011 г. меньше, чем для 2010 г. и 2012 г. Различие между результатами расчета расходов песчаных донных наносов для 2010 г. и 2012 г. объясняется различной продолжительностью безледного периода. В 2012 г. безледный период длился на 1 месяц дольше, чем в 2010 г., поэтому в 2012 г. и поток наносов оказался больше. По-видимому, в качестве характерных климатических величин вдольберегового переноса донных наносов следует опираться на расчеты для 2010 г., когда продолжительность безледного периода была близка к климатической. Значения, полученные для 2009 г. и 2011 г. (практически равные между собой), можно использовать в качестве характерных величин перемещения наносов путем сальтации, а разность между значениями 2010 г. и 2009 (2011) гг. - в качестве характерных величин перемещения наносов путем влечения.
Для безледных периодов 2010 г. и 2012 г. хорошо видно, что расходы через створы № 7 и № 8 (см. рис 3), лежащие к югу от подходного канала, существенно превышают расходы через створы № 1-6, лежащие к северу от канала. Так как для безледных периодов 2009 г. и
Рис. 3. Схема расположения створов, для которых проводились расчеты транспорта наносов за 2009-2011 гг. применительно к Харасавэйскому полигону в целом и к подходному каналу порта в устье реки Харасавэй в частности
2011 г. (для которых не учтен транспорт влекомых наносов) такой закономерности не наблюдается, можно заключить, что к югу от подходного канала резко возрастает транспорт влекомых наносов.
На рис. 5 цветом показаны расчетные мощности песчаных отложений в сантиметрах: положительные значения соответствуют накоплению наносов, отрицательные - их размыву. Для всех периодов отчетливо видны занесение подходного канала и размыв грунта подводного продолжения месторождения Харасавэй. Следует отметить, что процесс аккумуляции/размыва наносов как сильно нестационарный зависит от волнения, приливных и сгонно-нагонных течений. Поэтому схемы отложений (см. рис. 5) соответствуют итоговым состояниям на конец каждого безледного периода в течение четырех лет. Нестационарность заносимости канала проявляется в том, что в зависимости от синоптической ситуации он может и размываться. Но в целом литодинамическая ситуация способствует его занесению, что и отражено на рис. 5.
Высокая концентрация наносов на Хараса-вэйском полигоне наблюдается в основном вдоль берега (рис. 6, см. а), что хорошо согласуется с данными натурных наблюдений и объясняет максимумы расходов взвешенных наносов вблизи берега. Инструментальная погрешность определения концентрации взвешенных наносов составляет 0,0001 % , т.е. все концентрации ниже этого значения следует считать чистой водой, в которой не происходит транспорта взвешенных наносов. Транспорт взвешенных наносов происходит во вдольберего-вой полосе шириной не более 10 км.
120000
100000
80000
^ 60000 о
40000
20000
0
0
1
2 3 4 5 6 7
№ створа
Рис. 4. Результирующие через створы № 1-8 (см. рис. 3) вдольбереговые расходы донных песчаных наносов за 2009-2012 гг., м3
№ 2 (22 ) / 2015
104
Научно-технический сборник • ВЕСТИ ГАЗОВОЙ НАУКИ
-70 -15 -10 -5
а
б
-75 -70 -15 -10 -5 О
О 15
-50 -70 -10 -5
в
г
Рис. 5. Отложения донных наносов в районе Харасавэйского подходного канала за 4 безледных периода в 2009-2012 гг. (соответственно а, б, в, г): изолиниями показана топография дна, м
Необходимо отметить, что на рис. 5, 6 показаны мощности илистых отложений в «сухом» виде. Поэтому для определения реальной мощности илистых осаждений (они должны быть больше) нужно вводить поправку, которую, в свою очередь, необходимо выявить из данных наблюдений. В целом в канале отмечается отложение илистых наносов.
***
Таким образом, реализован современный комплекс моделирования, позволяющий рассчитывать характеристики атмосферного воздействия, ветрового волнения, гидротермодина-
мические характеристики морской циркуляции (скорости течений, температуру, соленость, уровень моря, толщину и сплоченность морского льда), транспорта и отложения наносов. Полученные характеристики хорошо согласуются с данными натурных наблюдений в ходе экспедиционных работ в Байдарацкой и Обской губах и в районе Харасавэйского месторождения. В частности, получены характеристики структуры и величины потока наносов для района Харасавэйского месторождения, а также выявлена изменчивость рельефа подводного берегового склона.
№ 2 (22) / 2015
Современные подходы и перспективные технологии в проектах освоения нефтегазовых месторождений российского шельфа
105
а
б в г
Рис. 6. Средняя за безледный период расчетная концентрация, %, взвешенных наносов на Харасавэйском полигоне (а) и расчетная мощность осажденных илистых наносов, мм, за безледные периоды 2009 г. (б), 2010 г. (в) и 2012 г. (г)
Список литературы
1. Леонтьев И.О. Прибрежная динамика: волны, течения, потоки наносов / И.О. Леонтьев. - М.: Геос, 2001.- 272 с.
2. Дианский Н.А. Моделирование циркуляции океана и исследование его реакции на короткопериодные и долгопериодные атмосферные воздействия / Н.А. Дианский. -М.: Физматлит, 2013. - 272 с.
3. Дианский Н.А. Воспроизведение циркуляции Карского и Печорского морей с помощью системы оперативного диагноза и прогноза морской динамики / Н.А. Дианский,
В.В. Фомин, И.М. Кабатченко, В.М. Грузинов // Арктика: экология и экономика. - 2014. -№ 1 (13). - С. 57-73.
4. Кабатченко И. М. Моделирование ветра и волн при вторичных термических циклонах на Черном море / И.М. Кабатченко,
Г В. Матушевский, М.В. Резников,
М.М. Заславский // Метеорология и гидрология. - 2001. - № 5. - С. 61-71.
5. Дианский Н.А. Особенности распространения загрязнений в северо-западной части Тихого океана / Н.А. Дианский, А.В. Гусев,
B. В. Фомин // Известия РАН. Физика атмосферы и океана. - 2012. - Т. 48. - № 2. -
C. 247-266.
№ 2 (22 ) / 2015
