|
楼主 |
发表于 2007-9-4 13:07
|
显示全部楼层
来自: 中国广西梧州
(转载) 渤海浅水油田FPSO运动性能及触底分析
李欣 杨建民 肖龙飞
(海洋工程国家重点实验室,上海交通大学)
摘 要 浮式生产储油系统(FPSO)在浅水中的运动响应特性和触底问题,是直接关系到系统安全的重大技术问题,备受设计部门和海洋石油公司的关注.本文在线性三维势流理论的基础上,采用时域计算方法对BZ25-1油田16万吨级FPSO不同吃水条件下的波频运动和触底情况进行了数值模拟,同时在海洋工程国家重点实验室水池进行了模型试验.计算和试验结果都表明,在浅水中,随着吃水的增加,即船底与海底间隙的减小,FPSO的垂荡,横摇和纵摇运动趋于缓和;满载FPSO在遭遇百年一遇风浪时会发生触底现象,须减少载油量.这对我国渤海等浅水海域油田应用FPSO系统进行开发具有一定的意义.
关键词:FPSO,浅水,时域计算,触底
1 前言
近年来,FPSO(Floating Production, Storage and Offloading System)在解决海上油田的开采和储存问题上越来越受到生产企业和研究机构的重视,发展迅速.我国自行研制的第一艘FPSO'渤海友谊'号于1989年4月交付使用至今,在渤海油田已有'渤海常青'号,'渤海明珠'号,'渤海世纪'号投产.我国的渤海地区有丰富的油田资源,到目前为止,探明地质储量6亿吨,其大部分为浅水油田.在水深较浅的渤海海域建造FPSO时,浅水深和船型尺度以及环境条件对超大型FPSO运动性能的影响非常重要.随着西非等深海海域陆续发现大型油田,国际上海洋石油开发和研究热潮正朝着深水方向发展,对于各种系泊形式FPSO的研究大多数为深水研究,国内外有关浅水FPSO的研究极少,特别是浅水软刚臂系泊超大型FPSO的数值模拟,由于船底与海底间的小间隙带来特殊复杂的水动力特性,国内外罕有相关研究报道.研究FPSO在浅水中的特性,对浅水油田的开发有着重要而深远的意义
BZ25-1油田位于渤海湾中南部,是近年来在渤海海域发现的5个亿吨级油田之一.根据该油田工程开发方案,将采用16万吨级FPSO系统,并直接引用满载吃水14.5m的QHD32-6 FPSO的设计结果.该油田所在地点水位极浅,极端水位为16.7m,而所采用的FPSO设计吃水为14.5m,间隙仅为2.2m,非常容易出现碰底的现象.如果按照一般最小安全作业间隙要求,则BZ25-1 FPSO的吃水只能达到12.34m,这将比原设计的原油装载量少近3万吨.因此,通过在上海交通大学风浪流水池进行理论计算和模型试验,确定BZ25-1 FPSO在百年一遇环境条件生存工况下,船底与海底之间的最小安全作业间隙,可为工程设计与施工提供可靠依据.如果试验结果使得BZ25-1FPSO的安全作业间隙得以减小,将使其原油舱容得到最大程度的利用,提高原油装载量和经济效益.
2 研究对象
本文对BZ25-1南油田FPSO在生存工况下的运动和碰底情况进行研究.试验和计算的内容包括在百年一遇海洋环境条件下的生存工况,考虑FPSO在不同吃水条件下其波频运动变化规律和触底情况.该FPSO系统由一艘16 kDWT FPSO,Yoke和塔构成.FPSO的主要参数如表1所示.百年一遇环境条件由表2列出.
表1 16万吨FPSO主要设计参数
总长:287.4 m
水线长:282.0 m
型宽:51.0 m
型深:20.6 m
设计吃水:14.5 m
载重量(DWT):160587.8 t
原油重:138006.0 t
满载排水量:201287.0 t
满载重心距基线高:13.002 m
表2 海洋环境条件
项目 单位 水深(m) 有效波高(m) 浪向( °deg ) 风速(m/s) 风向(°deg ) 流速(m/s) 流向(°deg)
1 16.7 5 180 33.9 180 1.73 135
2 16.7 5 180 33.9 180 1.73 90
波浪均采用PM谱,百年一遇工况波浪:有义波高5.0m,谱峰周期10.3s.所有工况均采用双参数PM谱,所选用的目标谱公式如下:
式中,
S(f) ———谱密度()
f ———波频(Hz)
g ———重力加速度()
Hs ———有义波高(m)
fp ———谱峰频率(Hz)=
Tp ———谱峰周期(s)
=1 ———形状参数
= 0.09 对于 f > fp
= 0.07 对于 f < fp
3 基础理论
3.1 三维势流理论
假设流体为理想流体,并且为无旋流.再Cartesian坐标系中,可以得到:
(1)
式中,,表示长峰波的频率,为自由表面起伏,为速度,为势函数, 下划线表示向量.线性化后的控制方程如下:
自由表面条件:
物面条件:
海底条件:
在线性化的条件下,可以将分解为绕射势,散射势和辐射势[4]:
(2)
式中,为绕射势,为散射势,为物体运动引起的势.
选择合适的格林函数,使它满足拉普拉斯方程,线性自由面条件,海底条件和辐射条件,应用FPSO船体的边界条件,可以计算出格林函数的数值解.积分方程可以简化为一系列的线性代数方程式.求解格林函数的过程在此就不作太多推导[4,5].
3.2 时域运动方程
FPSO的运动方程为:
(3)
.
这里,为波频运动,和分别为质量矩阵和静恢复力矩阵,为附加质量矩阵.为时延函数:
(4)
其中,是频域里的阻尼系数矩阵.为方向上的一阶波浪力,可根据Cummins(1962)提出的脉冲响应方法与频域计算中得到的波浪力联系起来,即设
(5)
其中为时刻的波形坐标,当时,是脉冲响应.而又与频率响应互为Fourier变换,即:
(6)
若已知整个频率范围内的,即可按上式求得,然后按波浪时历,据式(9)求得时域中的一阶波浪干扰力.
上述频域中的附加质量,阻尼系数以及波浪力均可通过的线性三维势流理论应用源汇分布法计算得到.因此求解时域运动方程(3),即可得到FPSO的波频运动时历.
4 数值模拟和模型试验
4.1 数值模拟
基于线性三维势流理论和时域计算方法,对该FPSO在百年一遇环境条件下的运动进行数值模拟.
建立FPSO的几何模型,对其外表面水下部分划分网格,划分时因为FPSO左右对称,所以只划分左舷部分即可,划分结果见图1,整条船共划分了个网格.
图1 FPSO的网格划分
用单根锚链组合模拟yoke单点系泊系统,使其刚度曲线与实验中的刚度曲线吻合,模拟结果刚度曲线见图2.
图2 静水回复力曲线
给定百年一遇风浪流和水深等环境条件,对FPSO单点系泊系统在风浪流下的运动进行时域模拟,得到运动时历曲线,对运动时历进行谱分析得到响应谱密度函数.选择与船重心处的roll和pitch,距船首柱14m距船底39.5m处点的heave这三个波频运动进行分析对比;选择四个参考点分别位于首柱后20m和尾柱前15m左右舷侧船底,即平行中体的起止点,计算四点距海底的距离,即船底与海底的间隙来判断FPSO是否碰底,如果间隙小于零,既说明发生了碰底.
4.2模型试验
模型试验在海洋工程国家重点实验室水池中进行,模型的缩尺比选为64.
(1)水深
实际水深16.7m,模拟水深0.261m.水深的模拟由调节可升降的假底来实现.
(2)流
百年一遇海况表层流速1.73m/s,中层流速1.63 m/s,底层流速1.29 m/s.流由造流系统模拟,以模拟表层流速为准,其相应模型试验流速为0.216m/s.
(3)风
风由造风系统模拟.造风系统由安装在拖车下方的轴流风扇组成.百年一遇海况定常风速为33.9m/s,其相应模型试验风速为4.24m/s.
(4)波
波浪的模拟由造波机来实现,水池中的不规则波浪根据给定的有义波高,周期及波谱进行模拟.模拟结果如下图所示.
图3 PM波浪测量谱与目标谱的比较
5 计算与实验结果
5.1 波频运动
图4为FPSO在百年一遇环境条件下,试验结果和计算结果垂荡,横摇和纵摇的时历曲线对比,从图中可以看出,垂荡和纵摇计算结果与实验的结果吻合良好.横摇运动符合程度稍差,主要原因是粘性横摇阻尼等因素的影响.
表3 列出了在百年一遇环境条件下不同吃水时,计算结果和试验结果垂荡,横摇和纵摇运动的统计值对比,可以看出,随着吃水的增加,船底与海底的间隙的减小,波频运动趋于缓和.图5表示吃水为14.5m,14.2m,14.0m,13.7m和12.34m时的波频运动计算结果进行谱分析后的对比,同样的趋势可以从图中看出,波频运动随着间隙的减小而减小.
5.2 触底分析
图5为百年一遇环境条件下船首右侧点船底与海底间隙的计算结果和试验结果时历曲线对比,可以看出时历曲线吻合较好.
表4列出了百年一遇环境条件下不同吃水时船底与海底间隙的最小值.在吃水为14.5m时,计算结果中出现了两个负值,分别为船首右侧和船尾左侧,说明在此吃水时这两处发生了碰底;而试验中其最小值虽然不小于零,但非常小,在测量的误差范围内可以认为发生了碰底现象.在14.2m计算结果出现一个绝对值很小的负值,即在此吃水时也发生了碰底,而试验结果没有发生碰底现象.在其他吃水时理论计算和试验均没有发生碰底.
6 结论
(1) 用数值计算来研究和预测在浅水中软钢臂系泊系统FPSO的运动是一种非常有效的方法.其中FPSO垂荡和纵摇运动的时域计算结果与相应的模型试验结果,包括时历曲线与统计值均基本吻合.横摇运动符合程度稍差,需要做进一步的研究.
(2) 在浅水中,随着吃水的增加,即船底与海底间隙的减小,FPSO的垂荡,横摇和纵摇这些波频振荡运动趋于缓和,从而不易导致碰底现象的发生.
(3) 百年一遇环境条件下,满载FPSO会发生触底现象,所以在遭遇百年一遇风浪时, FPSO须减少载油量,减少吃水.
|
|