2008/2009年冬季南海冷涌天气过程的海-气热通量交换及热量收支

冬季风期(11月—翌年3月)南海显著的气候特点是盛行东北季风并频繁地发生冷涌天气过程。使用2008年10月到2009年4月在西沙群岛永兴岛近海进行的海-气通量观测试验资料,分析了西沙海域冬季风期,尤其是冷涌时段的海-气通量交换和热量收支特征。结果表明:冬季风前期由于海-气温差增大,感热通量比西南季风期稍增加;潜热通量平均值与西南季风期接近;太阳总辐射明显降低,大气长波辐射减小,海洋热量净收入成为负值,使得秋季之后海面温度不断降低。冷涌期间海-气之间的感热通量高于冬季风期平均值,潜热通量大部分(1月份之前)也高于冬季风期平均值;由于潜热通量增大和太阳短波辐射减小,1月份之前的冷涌过程海洋热量净收支普遍出现较大负值,海洋失热量强于冬季风期,甚至强于2008年台风过程平均值。到了冬季后期太阳总辐射增强,海洋热量净收入转为正值,海水温度又逐渐升高。季节之间比较,观测区感热通量以冬季风期间最大,秋季次之,春季最小;而潜热通量夏季风期出现最大值,冬季次之,秋季最小。     

四维变分同化技术在风暴潮数值模拟中的应用

利用区域海洋模式ROMS(Regional Ocean Modelling System)及其四维变分同化模块,建立了具有资料同化能力的东中国海风暴潮数值模式,通过将海洋站水位观测资料同化到风暴潮模式中,提高了模式对风暴潮的模拟精度。四维变分同化技术能够在整个同化时间窗口保持动力协调,使模拟结果在该时间窗口内最大程度的靠近观测,同时,得到了最优预报初始场。利用该模式,对两次温带风暴潮过程进行了数值模拟,结果表明:在同化窗口内,同化对模拟精度有明显的提高;结束同化之后,得到的最优预报初始场对临近预报精度也有一定提高。     

2004年秋季冷空气活动对南海海表温度的影响

利用航次观测和网上的有关资料对南海2次强度不同的冷空气活动及其对南海SST的影响进行了分析.结果显示,9月22日弱冷空气过程南下速度慢,在陆地上变性较明显,未造成南海SST的明显变化;而10月2日前后的强冷空气过程南下速度快,陆地上变性比较弱,造成南海SST明显下降.通过对南海海表热收支分析,发现南海北部SST下降主要是冷空气造成净热通量急剧增加,海洋失去热量,而南部SST下降可能是南部海面气旋式风应力引起的下层冷水上涌.初步解释了2004年秋季冷空气活动对南海SST的影响.     

珠江口伶仃洋风暴潮传播特性及影响因素数值模拟研究——以台风“山竹”为例

基于FVCOM海洋模式对珠江口伶仃洋内及周边海域的风暴潮增水传播过程进行研究。首先,建立珠江口风暴潮模型,采用超强台风“山竹”作为典型案例进行风暴潮过程模拟,并对模型的计算结果进行验证,发现模拟结果和测站潮位结果比较吻合;然后,对伶仃洋在“山竹”登陆前后的风暴潮增水过程的时空分布及变化特征进行分析。从伶仃洋湾口到湾顶选取12个点进行时间序列分析,发现除了因距离外海远近不同导致的相位差异外,基本特征相似,符合国际上类似海湾内的风暴潮增水波动特征,可分为初振段、主振段和余振段。为了进一步研究台风参数差异对伶仃洋风暴潮增水的影响,本文基于“山竹”超强台风的特征参数,设计了一系列变化条件下的数值试验,结果发现：（1）台风的登陆时间会影响到风暴潮增水和天文潮之间的相位关系,进而影响到增水的大小。如果风暴潮增水极值正好在天文潮高潮位,风暴潮增水就会削弱,而风暴潮增水正好在天文潮低潮位,风暴潮增水就会增强。（2）台风中心压差决定了台风风力的大小,从而影响风暴潮增水。但是在同一海湾内的影响在空间上并不相同,在较浅区域影响大而较深区域影响小。（3）台风路径会对风暴潮增水产生较大影响。基于“山竹”的路径,...     

Numerical Simulation of Waves Generated by Seafloor Movements

Waves generated by vertical seafloor movements are simulated by use of a fully nonlinear two-dimensional numerical wave tank. In the souree region, the seafloor lifts to a designated height by a generation function. The numerical tests show that the linear theory is only valid for estimating the wave behaviors induced by the seafloor movements with a small amplitude, and the fully nonlinear numerical model should be adopted in the simulation of the wave generation by the large amplitude seafloor movements. Without the background surface waves, many numerical tests on the stable maximum elevations η0^max are carried out by beth the linear theory and the fully nonlinear model. The results of two models are compared and analyzed. For the fully nonlinear model, the influences of the amplitudes and the horizontal lengths on η^max are stronger than that of the characteristic duration times. Furthermore, results reveal that there are significant differences be- tween the linear theory and the fully nonlinear model. When the influences of the background surface waves are considered, the corresponding numerical analyses reveal that with the fully nonlinear model the η0^max near-linearly varies with the wave amplitudes of the surface waves, and the η0^max has significant dependences on the wave lengths and the wave phases of the surface waves. In addition, the differences between the linear theory and the fully nonlinear model are still obvious, and these differences are significantly affected by the wave parameters of the background surface waves, such as the wave amplitude, the wave length and the wave phase.     

大亚湾风暴潮波浪传播数学模型

本文建立了大亚湾风暴潮涌浪传播数学模型,以任意多边形离散计算海区,每一多边形构成单元,以波浪运动学和动力学守恒方程模拟单元内能量传递,以风暴潮过程模拟边界入射波高过程,用风暴潮涌浪传播基本方程和波能缓坡方程结合模拟湾口巨浪向湾内的传播过程。通过分析大亚湾不同种类岸线反射系数的概率分布,并结合实测波高对模型进行率定,最终确定模型参数。将大亚湾特征点计算波高与统计推荐波高比较进行模型验证,结果显示SE向波高与H_(13%)推荐波高对应较好,可以用于大亚湾海区的波浪预报。计算当大亚湾口分别出现10年一遇、50年一遇及100年一遇的波高时,在E、ESE以及ES向入射波浪条件下大亚湾海域极值波高的分布。分别对风暴潮涌浪在不同类型岸线的爬高以及风暴潮涌浪传播对岸线的作用力进行计算。     

南黄海辐射沙洲近岸海域波浪特性研究

Errors will be caused in calculating the fatigue damages of details in liquid cargo tanks by using the traditional spectral analysis method which is based on linear system, for the nonlinear relationship between the dynamic stress and the ship acceleration. An improved spectral analysis method for the assessment of the fatigue damage in detail of a liquid cargo tank is proposed in this paper. Based on assumptions that the wave process can be simulated by summing the sinusoidal waves in different frequencies and the stress process can be simulated by summing the stress processes induced by these sinusoidal waves, the stress power spectral density (PSD) is calculated by expanding the stress processes induced by the sinusoidal waves into Fourier series and adding the amplitudes of each harmonic component with the same frequency. This analysis method can take the nonlinear relationship into consideration and the fatigue damage is then calculated based on the PSD of stress. Take an independent tank in an LNG carrier for example, the accuracy of the improved spectral analysis method is proved much better than that of the traditional spectral analysis method by comparing the calculated damage results with the results calculated by the time domain method. The proposed spectral analysis method is more accurate in calculating the fatigue damages in detail of ship liquid cargo tanks.     

Antiflotation design for water tank using pressure relief technique

A pressure relief technique has been proposed to reduce hydrostatic pressures by opening drainage holes around a water tank. This solution is evaluated by laboratory experiments and numerical parametric investigation. After seepage flow is allowed, the hydraulic head acting on the water tank is reduced due to seepage losses. The drainage holes should be opened on the base slab near the sidewalls, and at a certain height on the sidewalls. Given a fixed total opening area, the number of drainage holes with smaller diameter should be allowed near the edge of the water tank to increase the efficacy of the approach. The properties of the surrounding soil influence the results significantly, where a cushion layer with a higher hydraulic conductivity and a greater thickness is beneficial to the stability of the water tank, and a backfill layer with adequate thickness and hydraulic conductivity should be selected. An illustrative example is given in the end to demonstrate the advantage of the proposed antiflotation design strategy compared to the conventional enhanced self-weight method, and more economic design using less reinforcement and concrete can be achieved.     

东南沿海台风风暴潮增水过程中非线性机制和地形的作用研究:以1509号台风"灿鸿"为例

本文基于FVCOM（Finite Volume Community Ocean Model）构建了一个覆盖中国渤海、黄海和东海的数值模型,采用NCEP-CFSR风场数据对1509号台风“灿鸿”产生的风暴潮进行模拟,与实测水位数据的对比表明该模型可靠、模拟结果合理。基于此模型,本文对非线性作用和地形在风暴潮增水过程中的作用进行了研究。首先,重点分析了增水过程中潮汐与风暴潮的非线性作用,结果表明：高潮时非线性作用使增水值降低;低潮时非线性作用使增水值升高。另外,开边界处分别只添加M₂、S₂和K₁分潮,分析天文潮的潮高和周期对非线性作用的影响,结果表明：潮高越高,非线性作用越明显;半日潮的非线性作用较全日潮更明显;并且,增水极值附近出现的半日周期的波动也与非线性作用有关。其次,除了非线性作用,地形对风暴潮的增水也有一定影响,本文改变地形的实验结果表明：坡度越大,增水极值越小。琉球群岛的存在使得东南沿海出现风暴潮增水的面积减小,但使得风暴潮增水的高值区域扩大。     

基于GPU技术的风暴潮集合预报

台风预报的准确性在风暴潮预报中起着重要作用。台风强度和路径的不确定性意味着使用集合模式来预报风暴潮。本文利用中央气象台的最优路径台风参数驱动国家海洋环境预报中心业务化的水动力学模型,开展华南沿海的风暴潮模拟,模式模拟结果与实测吻合较好。为了改进计算效率,采用CUDA Fortran 语言对模型进行了改造,改造后的模型在计算结果与原模型基本一致的基础上,计算时间缩短了99%以上。通过融合欧洲中期天气预报中心（ECWMF）的50条路径与3种可能台风强度构造出了150个台风事件,并用150个台风事件驱动改进的风暴潮数值模型,计算结果可以提供集合预报产品和概率预报产品。通过“山竹”台风风暴潮过程可以发现集合平均预报结果和概率预报结果与实测吻合较好。改进的数值模型可以运行普通工作站上,非常适合风暴潮集合预报,并且可以提供更好的决策产品。