有冰海域溢油运动数值模型研究

提出了寒冷海域溢油运动短期数值模型,其中海冰运动预报模型针对不同海冰密集度的本构关系亦不同,在近岸及密集度较高的冰带采用基于连续介质理论的粘塑性本构关系;在冰边缘密集度较低的冰带采用基于离散介质理论的碰撞本构关系计算冰内力;溢油运动预报模型包括连续冰盖及高、中、低密集度冰场中溢油的漂移和扩展过程的预报.模拟计算结果表明模型的建立是合理的.该模型可用于事故溢油预报,为溢油的围收治理提供依据.     

基于Mohr-Coulomb准则的黏弹-塑性海冰动力学本构模型

针对中小尺度下海冰动力作用过程中的漂移、重叠和堆积特征,在连续介质力学基础上建立了一个黏弹塑性海冰动力学本构模型.该模型主要包括四部分,即海冰在屈服前的Kelvin Vogit黏弹模型、海冰塑性屈服的Mohr Coulomb准则、屈服后相关联的正交流动法则和影响海冰强度的静水压力.采用黏弹塑性本构模型对规则海域内的海冰堆积过程进行了数值试验,其结果与传统黏塑性模型、改进的黏塑性模型和经典冰坝理论的解析解相符,可合理地模拟海冰堆积的动力过程.为进一步检验该黏弹塑性本构模型的适用性,对渤海海冰的动力过程进行了48h的数值模拟,计算结果与海冰卫星遥感现场观测资料一致.通过与传统的黏塑性和改进的黏塑性模型结果的对比分析,进一步验证了黏弹塑性模型在中小尺度海冰动力学模拟中的可靠性.在以上海冰堆积的数值试验和渤海海冰动力学模拟中均采用了光滑质点流体动力学(SPH)方法.     

海冰动力学过程的数值模拟

讨论了海冰动力学性质并阐述决定海冰漂移的动量平衡,冰脊和水道形成及确定冰应力与形变、强度之间关系的海冰流变学.提出了模拟海冰动力学过程的数值模式,模式中冰厚分布由开阔水、平整冰和堆积冰3种要素表示.在这3要素的预报方程中引入形变函数,采用一种参数化方法模拟冰脊和水道.为了表示冰内相互作用,将海冰作为一种非线性粘性可压缩物质,采用粘-塑性本构关系.本文还概述和讨论了模式中所采用的数值方法,应用此模式模拟了渤海、波罗的海的波的尼亚湾和拉布拉多海的冰漂移.渤海冰漂移模拟结果明显地显示出潮周期变化,还模拟了渤海的冰脊和水道,进行了海冰流变学参数的敏感性试验.并将此冰模式与大气模式和边界层模式联接,给出渤海海冰预报结果.     

海冰模式和数值预报

根据渤海冰情,在海冰动力学和热力学研究基础上,提出一种海冰动力一热力模式。该模式包括平整冰、堆积冰和开阔水三要素,采用粘一塑性本构关系计算冰内应力。利用参数化方法处理变形函数和热力增长率。该冰模式与潮流模式耦合研究渤海冰－潮相互作用,并与ＥＣＯＭ－Ｓｉ海洋模式耦合进行预报试验。该冰模式与大气模式联接,或应用气象部门数值天气预报产品进行海冰业务数值预报。本文还给出了预报结果。     

第2讲 海冰动力学

简单回顾了海冰动力学的发展,介绍了海冰漂移、形变和堆积等动力学过程,对比分析了不同海冰本构关系在海冰流变学的应用。     

膨胀式自应力灌浆卡箍灌浆本构关系模型的研究

膨胀式自应力灌浆卡箍是一种典型的钢/灌浆组合结构,由于膨胀受限灌浆处于三向受压状态,在这种结构的强度校核中采用合适的灌浆本构模型是影响设计结果的重要因素。采用ABAQUS建模这种卡箍,构筑了三种灌浆本构关系模型,主要测试了莫尔库仑模型和塑性损伤模型的相关材料参数,然后通过一个节点灌浆卡箍的设计算例,比较和分析了三种灌浆本构模型的计算结果。仿真分析发现对于卡箍的灌浆强度校核,线弹性模型最保守,容易校核不合格,而塑性损伤模型最不保守。建议采用能反映灌浆的压碎和开裂行为的塑性损伤模型作为膨胀式灌浆卡箍强度校核分析的本构关系。研究结果可以为自应力灌浆卡箍的有限元强度校核方法提供有益的借鉴和指导。     

基于MODIS数据的渤海海冰厚度反演算法优化

海冰厚度是监测与研究渤海海冰的重要参数。为了获取更加可靠的渤海海冰厚度数据,本研究基于MODIS数据对海冰厚度反演中的冰水分离环节和冰厚计算方法都进行了改进。对于冰水分离环节,本文在Canny边缘检测算子提取海冰基础上,加入了二值化处理、阈值判别等步骤,实现了较高精度的渤海海冰范围自动化提取。通过试验确定了海冰厚度与反照率指数关系模型中的参数,包括海冰衰减系数和海水反照率参数,使其更加符合渤海海区的物理特征。将改进后算法的海冰厚度反演结果与渤海海上石油平台实测数据进行比较,并分析了误差来源。结果表明,经过对算法的改进,海冰厚度与反照率指数关系模型的反演结果与实测数据之间的平均绝对误差由7.05 cm缩小到2.74 cm,相关系数由0.434提高到0.485。     

海冰动力过程的改进离散元模型及在渤海的应用

海冰的断裂、重叠和堆积等离散分布特性广泛地存在于极区和副极区的不同海域,并对海冰的生消、运移过程有着重要影响。针对海冰在不同尺度下的离散分布特点,发展海冰动力过程的离散元方法有助于完善海冰数值模式,提高海冰数值模拟的计算精度。为此,本文针对海冰生消运移过程中的非连续分布和形变特性,发展了适用于海冰动力过程的改进离散元模型(MDEM)。不同于传统离散元方法,该模型将海冰离散为具有一定厚度、尺寸和密集度的圆盘单元。海冰单元设为诸多浮冰块的集合体,其在运移和相互接触碰撞过程中,依照质量守恒发生单元尺寸、密集度和厚度的相应变化。基于海冰离散性和流变性的特点,该模型采用黏弹性接触本构模型计算单元间的作用力,并依据Mohr-Coulomb准则计算海冰法向作用下的塑性变形及切向摩擦力。为验证该模型的可靠性,本文对海冰在规则水域内的运移和堆积过程进行了分析,离散元计算结果与解析值相一致;此外,对旋转风场下海冰漂移规律的模拟进一步验证了本文方法的精确性。在此基础上,对渤海辽东湾的海冰动力过程进行了48h数值分析,计算结果与卫星遥感资料和油气作业区的海冰现场监测数据吻合良好。在下一步工作中将考虑海冰离散元模拟中的热力因素影响,发展具有冻结、断裂效应的海冰离散元模型,更精确地模拟海冰动力-热力耦合作用下的生消和运移过程。     

滩海桶形基础平台整体三维弹塑性有限元地震动力响应分析

在经典的D-P本构模型基础上,通过改进其屈服准则。引入硬化曲线建立了改进的D-P模型,使该本构模型能很好地模拟大港滩海软土的应力应变关系。利用MARC程序作为平台,建立了滩海桶形基础平台整体有限元动力分析模型,并对我国第一座桶形基础平台进行了三维弹塑挂地震时程动力分析,得到了一些很有用的结论。     

辽东湾海冰单轴压缩长期强度的初步分析

冰的蠕变特性是其自身独特的一种性质,正是由于这一特性才导致了其本构关系的复杂性。室内确定海冰静压力强度特征的标准方法有:恒变形速率或加载速率下的压缩强度试验和恒荷载下的压缩蠕变试验。本文是根据近年来进行的海冰圆柱状试样单轴压缩强度和单轴压缩蠕变试验的结果,分析海冰的长期强度特征并给出统计结果。分析结果还表明海冰的弹性范围很小,即便是单轴压缩强度试验,也只有在加载时间很短的情况下,才表现出拟弹性阶段。     

The effects of resolving the Canadian Arctic Archipelago in a finite element sea ice model

Though narrow straits may have a strong influence on the large-scale sea ice mass balance, they are often crudely represented in coarse resolution sea ice models. Unstructured meshes, with their natural ability to fit boundaries and locally increase the mesh resolution, propose an alternative framework to capture the complex oceanic areas formed by coasts and islands. In this paper, we develop a finite element sea ice model to investigate the sensitivity of the Arctic sea ice cover features to the resolution of the narrow straits constituting the Canadian Arctic Archipelago. The model is a two-level dynamic-thermodynamic sea ice model, including a viscous-plastic rheology. It is run over 1979–2005, forced by daily NCEP/NCAR reanalysis data. Confronting qualitatively numerical experiments with observations shows a good agreement with satellite and buoys measurements. Due to its simple representation of the oceanic interactions, the model overestimates the sea ice extent during winter in the southernmost parts of the Arctic, while the Baffin Bay and Kara Sea remain ice-covered during summer. In order to isolate the benefits from resolving the Canadian Arctic Archipelago, a numerical experiment is performed where we artificially close the archipelago. Focusing on the large-scale sea ice thickness pattern, no significant change is found in our model, except in the close surroundings of the archipelago. However, the local and short-term influences of the ice exchanges are nonnegligible. In particular, we show that the ice volume associated to the Canadian Arctic Archipelago represents 10% of the Northern Hemisphere sea ice volume and that the annual mean ice export towards Baffin Bay amounts to 125 km³ yr⁻¹, which may play an important role on the convective overturning in the Labrador Sea.     

Numerical simulation for dynamical processes of sea ice

ＮｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｆｏｒｄｙｎａｍｉｃａｌｐｒｏｃｅｓｓｅｓｏｆｓｅａｉｃｅＷｕＨｕｉｄｉｎｇ，ＢａｉＳｈａｎ，ＺｈａｎｇＺｈａｎｈａｉａｎｄＬｉＧｕｏｑｉｎｇ（ＲｅｃｅｉｖｅｄＭａｙ１６，１９９６；ａｃｃｅｐｔｅｄＪａｎｕａｒｙ...     

普里兹湾海冰季节性变化的高分辨率数值模拟

普里兹湾海冰以一年冰为主,海冰覆盖存在较大的季节性变化.海冰的分布及其季节性变化主要受当地大气环流及海流的影响.基于一个海洋-海冰耦和模式,模拟研究了该海区海冰的季节性变化特征.海洋模式基于MIT环流模式(MITgcm),海冰动力学模式参考Hibler类型的VP模型,热力学过程取自Winton三层模型.模式区域覆盖整个...     

质点-网格海冰模式在渤海的数值预报试验

在海冰动力学和热力学的研究基础上,开发应用于渤海的质点-网格海冰模式.该模式采用质点-网格法,有效地避免了传统海冰模式的数值扩散问题.该模式采用冰厚分布函数取代传统渤海业务海冰预报模式中平整冰、堆积冰和开阔水的三类海冰分型.应用质点-网格海冰模式对2003年～2004年冬季渤海冰情进行了业务化逐日数值预报试验,并对预报结果进行统计检验和分析比较.     

基于SMAP卫星雷达资料的海冰密集度反演技术研究

SMAP是美国于2015年初发射的一颗卫星，搭载了L波段的雷达。它采用圆锥扫描方式，具有固定的入射角、较大的幅宽和千米级的分辨率，在海冰监测方面具有独特的优势。本文利用SMAP卫星雷达资料分别与德国Bremen大学海冰密集度产品和美国国家冰雪数据中心（NSIDC）海冰密集度产品建立3.125 km和25 km匹配数据集，分析了L波段雷达后向散射系数、极化比和归一化极化差与海冰密集度之间相关性，建立基于人工神经网络的海冰密集度反演算法。为了验证SMAP卫星雷达资料反演海冰密集度的精度，本文选择德国Bremen大学和美国冰雪数据中心发布的海冰密集度产品分别与SMAP海冰密集度产品进行对比分析，SMAP海冰密集度与Bremen海冰密集度的偏差为0.07、均方根误差为0.14；与NSIDC海冰密集度的偏差为0.04、均方根误差为0.18，这表明SMAP海冰密集度产品与现有业务化海冰密集度产品具有很好的一致性。     

Ocean circulation and sea ice distribution in a finite element global sea ice–ocean model

A newly developed global Finite Element Sea Ice–Ocean Model (FESOM) is presented. The ocean component is based on the Finite Element model of the North Atlantic (FENA) but has been substantially updated and extended. In addition to a faster realization of the numerical code, state-of-the-art parameterizations of subgrid-scale processes have been implemented. A Redi/GM scheme is employed to parameterize the effects of mesoscale eddies on lateral tracer distribution. Vertical mixing and convection are parameterized as a function of the Richardson number and the Monin–Obukhov length. A finite element dynamic-thermodynamic sea ice–model has been developed and coupled to the ocean component. Sea ice thermodynamics have been derived from the standard AWI sea ice model featuring a prognostic snow layer but neglecting internal heat storage. The dynamic part offers the viscous-plastic and elastic-viscous-plastic rheologies. All model components are discretized on a triangular/tetrahedral grid with a continuous, conforming representation of model variables. The coupled model is run in a global configuration and forced with NCEP daily atmospheric reanalysis data for 1948–2007. Results are analysed with a slight focus on the Southern Hemisphere. Many aspects of sea ice distribution and hydrography are found to be in good agreement with observations. As in most coarse-scale models, Gulf Stream transport is underestimated, but transports of the Kuroshio and the Antarctic Circumpolar Current appear realistic. The seasonal cycles of Arctic and Antarctic sea ice extents and Antarctic sea ice thickness are well captured; long- and short-term variability of ice coverage is found to be reproduced realistically in both hemispheres. The coupled model is now ready to be used in a wide range of applications.     

海冰中侧向传播光衰减系数的解析验证

The attenuation of lateral propagating light(LPL) in sea ice was measured using an artificial light source in the Canadian Arctic during the 2007/2008 winter. The apparent attenuation coefficient μ(λ) for lateral propagating light was obtained from the measured logarithmic relative variation rate. In this study an analytical solution based on the strict optical theories is developed to validate the measured result. There is a good consistency between theoretical solution and measured data, by which a quite simple but very rigorous relationship among the light source, measurement geometry, and measured irradiance is established. The attenuation coefficients acquired by measurement and theory are the diffusion attenuation as an apparent optical property of ice, independent of the light source and shining condition. The attenuation ability of sea ice should be caused by the microstructure of sea ice, such as crystal size, ice density, brine volume, air inclusion, etc. It also includes the leak from both interfaces by directional scattering. It is verified that the measuring approach is operational and accurate to measure the attenuation of the LPL. The solution from this study did not tell the connection among the extinction and the inclusions of sea ice theoretically because of insufficient understanding.     

海冰对海浪影响研究综述

受全球气候变化的影响,极区海浪尤其是北极海浪在过去几十年发生了显著的变化,使得海冰边缘区海冰与海浪的相互作用愈发显著。本文从物理海洋学的角度出发,较系统地总结了海冰对海浪作用研究的国内外现状,从理论和实测的角度分别探讨了海冰对海浪能量的耗散及其引起的波动频散关系的变化,同时分析了当前海冰覆盖海域海浪的数值模拟与现场观测研究,指出了未来开展有冰海域海浪数值模拟与预报所面临的主要问题,并对该方向今后的研究做出展望。总体来看,尽管海冰对海浪作用的机理复杂且与海冰类型高度相关,但是海冰对海浪能量的衰减与传播距离基本呈指数关系,并且海冰会一定程度上影响海浪的传播速度。未来依然需要更多不同海冰类型下海浪的观测数据以开展进一步的机理分析、模型检验和参数校准,进而实现高精度的业务化预报。     

Polynyas in a high‐resolution dynamic–thermodynamic sea ice model and their parameterization using flux models

This paper presents an analysis of the solutions for a steady state latent heat polynya generated by an applied wind stress acting over a semi‐enclosed channel using: (a) a dynamic–thermodynamic sea ice model, and (b) a steady state flux model. We examine what processes in the sea ice model are responsible for the maintenance of the polynya and how sensitive the results are to the choice of rheological parameters. We find that when the ice is driven onshore by an applied wind stress, a consolidated ice pack forms downwind of a zone of strong convergence in the ice velocities. The build‐up of internal stresses within the consolidated ice pack becomes a crucial factor in the formation of this zone and results in a distinct polynya edge. Furthermore, within the ice pack the across‐channel ice velocity varies with the across‐channel distance. It is demonstrated that provided this velocity is well represented, the steady state polynya flux model solutions are in close agreement with those of the sea ice model. Experiments with the sea ice model also show that the polynya shape and area are insensitive to (a) the sea ice rheology; (b) the imposition of either free‐ slip or no‐slip boundary conditions. These findings are used in the development of a simplified model of the consolidated ice pack dynamics, the output of which is then compared with the sea ice model results. Finally, we discuss the relevance of this study for the modelling of the North Water Polynya in northern Baffin Bay.