CN104990660B - 一种应用于浮式平台的内波流载荷获取方法 - Google Patents
一种应用于浮式平台的内波流载荷获取方法 Download PDFInfo
- Publication number
- CN104990660B CN104990660B CN201510379297.1A CN201510379297A CN104990660B CN 104990660 B CN104990660 B CN 104990660B CN 201510379297 A CN201510379297 A CN 201510379297A CN 104990660 B CN104990660 B CN 104990660B
- Authority
- CN
- China
- Prior art keywords
- mrow
- wave current
- interior wave
- interior
- fluid
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
本发明涉及一种应用于浮式平台的内波流载荷获取方法，其包括以下步骤：1)建立内波流的描述方程；2)获取内波流的波面时程；3)获取内波流流速和加速度；4)求取内波流载荷。本发明由于在获取内波流载荷时将内波流流速和加速度随时间的变化考虑在内，从而能够考虑到作用在平台上的拖曳力和惯性力随时间的变化，使得最终得到的内波流载荷是一个随时间变化的时程，因而得到的内波流载荷更加准确。本发明可以广泛应用于浮式平台的内波流载荷的获取。
Description
技术领域
本发明涉及一种载荷获取方法，特别是关于一种应用于浮式平台的内波流载荷获取方法。
背景技术
随着我国深水海洋油气勘探开发的快速发展，近几年我国的大型深水半潜式钻井平台的数量大为增加。中国南海是我国深水油气开发的主战场，南海海洋环境条件十分恶劣，比较典型的特点是内波流频发，且强度大。系泊定位系统动力定位条件下，浮式平台受到内波流载荷的作用，会发生大范围的偏移，内波流载荷的获取是内波流作用下平台偏移计算的重要基础，当平台的偏移超出一定范围时，则会影响到深水钻井作业的进行，甚至威胁到平台系泊系统和钻井系统的安全。目前国内外海洋工程相关规范获取内波流载荷的方法是将内波流按照常规的定常流来处理，这种处理方法忽略了内波流流速的时变特征(即内波流的波动特征)，从而导致内波流载荷中的惯性力部分被忽略，这种传统计算方法会低估内波流的作用，从而影响平台的正常作业及安全。
发明内容
针对上述问题，本发明的目的是提供一种将内波流的波动性考虑在内的应用于浮式平台的内波流载荷获取方法。
为实现上述目的，本发明采取以下技术方案：一种应用于浮式平台的内波流载荷获取方法，其包括以下步骤：
1)建立内波流的描述方程为：
其中，η(x,t)为内波流的波面时程，t为时间，x为目标浮式平台的坐标，C0为内波流的线性速度，α为非线性参数，β为弥散参数；
2)获取内波流的波面时程：
根据式(1)可以推导得出内波流的波面时程表示为：
其中，Cp为内波流的非线性速度，l为内波流的特征波长，η0为内波流的波幅；
3)获取内波流流速和加速度：
假设内波流有两个分层，上层流体流速u1与下层流体流速u2相反，上层流体厚度为h1,下层流体厚度为h2，上层流体密度为ρ1，下层流体密度为ρ2，对式(2)进行求导得到内波流的上层流体流速u1和下层流体流速u2：
然后对式(3)进行求导即可得到内波流的上层流体加速度a1和下层流体加速度a2：
4)求取内波流载荷：
根据内波流的描述方程得到内波流流速、加速度的分布，将浮式平台的将内波流载荷F分成拖曳力FD和惯性力FI两部分，即：
F＝FD+FI (5)
其中，拖曳力FD可表示为：
式中，A为受流面积；ρ为海水密度，u为内波流流速，当浮体处在上层流体中时：ρ＝ρ1，u＝u1,当浮体处在下层流体中时，ρ＝ρ2，u＝u2；Cd为拖曳力系数，根据模型试验方法或者根据规范的方法来计算获得；
惯性力FI可表示为：
FI＝CmρVa (7)
其中，V为目标浮式平台的排水量；a为内波流加速度，当浮体处在上层流体中时：a＝a1,当浮体处在下层流体中时，a＝a2；Cm为附加质量系数，根据势流理论的水动力软件计算获得。
本发明由于采取以上技术方案，其具有以下优点：1、本发明由于在获取内波流载荷时将内波流的波动性考虑在内，因而能够更加真实的评估内波流的作用，从而能够提高平台设计的可靠性和安全性。2、本发明由于在获取内波流载荷时将内波流流速和加速度随时间的变化考虑在内，从而能够考虑到作用在平台上的拖曳力和惯性力随时间的变化，使得最终得到的内波流载荷是一个随时间变化的时程，因而得到的内波流载荷更加准确。本发明可以广泛应用于浮式平台的内波流载荷的获取。
附图说明
图1是内波流关键参数示意图；
图2是深水半潜式钻井平台的结构示意图；
图3是内波流流速时程图；
图4是内波流加速度时程图；
图5是平台纵向内波流载荷时程图；
图6是平台横向内波流载荷时程图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提供的一种应用于浮式平台的内波流载荷获取方法，其包括以下步骤：
1)建立内波流的描述方程为：
其中，η(x,t)为内波流的波面时程，t为时间，x为目标浮式平台的坐标，C0为内波流的线性速度，α为非线性参数，β为弥散参数。
2)获取内波流的波面时程：
根据式(1)可以推导得出内波流的波面时程表示为：
其中，Cp为内波流的非线性速度，l为内波流的特征波长，η0为内波流的波幅。
3)获取内波流流速和加速度：
假设内波流有两个分层(如图1所示)，上层流体流速u1与下层流体流速u2相反，上层流体厚度为h1,下层流体厚度为h2，上层流体密度为ρ1，下层流体密度为ρ2，对式(2)进行求导得到内波流的上层流体流速u1和下层流体流速u2：
然后对式(3)进行求导即可得到内波流的上层流体加速度a1和下层流体加速度a2：
4)求取内波流载荷：
根据内波流的描述方程得到内波流流速、加速度的分布，可以将浮式平台的将内波流载荷F分成拖曳力FD和惯性力FI两部分，即：
F＝FD+FI (5)
其中，拖曳力FD可表示为：
式中，A为受流面积；ρ为海水密度，u为内波流流速，当浮体处在上层流体中时：ρ＝ρ1，u＝u1,当浮体处在下层流体中时，ρ＝ρ2，u＝u2；Cd为拖曳力系数，可以根据模型试验方法或者根据规范的方法来计算获得。
惯性力FI可表示为：
FI＝CmρVa (7)
其中，V为目标浮式平台的排水量；a为内波流加速度，当浮体处在上层流体中时：a＝a1,当浮体处在下层流体中时，a＝a2；Cm为附加质量系数，可以根据势流理论的水动力软件计算获得。
本发明考虑了内波流流速u以及加速度a随时间的变化，从而也考虑到作用在目标浮式平台上的拖曳力FD和惯性力FI随时间的变化，因而，最终得到的内波流载荷F＝FD+FI也是一个随时间变化的时程(时间序列)，该时程可用于计算目标浮式平台在内波流作用下的偏移。
本发明由于同时考虑了内波流载荷F的拖曳力FD和惯性力FI，两者分别与内波流流速u和加速度a有关系，而又由于内波流流速u和加速度a的相位差为90°，通过拖曳力FD和惯性力FI时间序列的叠加，因此考虑了两者的不同相位。而传统方法要么不考虑惯性力FI(该方法会低估内波流载荷F)，即使考虑了惯性力FI，也只是简单地将惯性力FI的最大值和拖曳力FD的最大值叠加(该方法会高估内波流载荷F)，因此本发明提出的方法可以更为准确的获取内波流载荷F。
下面结合一个实施例来说明本发明的方法。
如图2所示，本实施例涉及的目标浮式平台为一座深水半潜式钻井平台，表1为该深水半潜式钻井平台的主要参数，表2为内波流实测数据。
表1深水半潜式钻井平台主要参数
| 名称 | 尺度 |
| 下浮体(m) | 114.07×20.12×8.54R2.13 |
| 立柱(m) | (15.86～17.39)×(15.86～17.39)×21.50/R3.96 |
| 横撑(m) | 2.438×1.83(4根) |
| 吃水(m) | 19 |
| 下浮体中心线距(m) | 58.56 |
表2内波流实测数据
已知内波流的上层流体最大流速为2.0m/s,内波流周期为1200s，波长为2000m，上层深度h1为180m，下层深度h2为420m，则根据式(3)分别得到内波流的上、下层流体流速时程：图3(a)为上层流体流速时程，图3(b)为下层流体流速时程；内波流的上、下层流体加速度时程：图4(a)为上层流体加速度时程，图4(b)为下层流体加速度时程。
在周期1200s的内波流载荷下，根据势流理论得到平台的附加质量系数纵荡附加质量系数Cm是1.302，横荡附加质量系数Cm是1.754。
当内波流来向为船艏向，即平台遭遇纵浪向内波流时，根据式(5)～(7)获取内波流载荷时程、拖曳力时程和惯性力时程(如图5(a)-(c)所示)，内波流载荷最大值为：1583.7kN；拖曳力成分最大值为：1435kN；惯性力最大值为：450kN。
考虑到拖曳力和惯性力最大值不会同时出现，在时域叠加后，本发明方法的计算结果(1583.7KN)比传统方法的计算结果(即拖曳力成分最大值1435KN)大10.36％。
当内波流来向为横浪向，即平台遭遇横浪向内波流时，根据式(5)～(7)获取内波流载荷时程、拖曳力时程和惯性力时程(如图6(a)-(c)所示)，内波流载荷最大值为3209.8kN；拖曳力成分最大值为：3075kN；惯性力最大值为：607.22kN。
考虑到拖曳力和惯性力最大值不会同时出现，在时域叠加后，本发明方法的计算结果(3209.8KN)比传统方法的计算结果(即拖曳力成分最大值3075KN)大4.38％。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
Claims (1)
1.一种应用于浮式平台的内波流载荷获取方法，其包括以下步骤：
1)建立内波流的描述方程为：
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mn>0</mn>
</msub>
<mo>+</mo>
<mi>&alpha;</mi>
<mi>&eta;</mi>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mi>&beta;</mi>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>3</mn>
</msup>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>x</mi>
<mn>3</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mn>0</mn>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，η(x,t)为内波流的波面时程，t为时间，x为目标浮式平台的坐标，C0为内波流的线性速度，α为非线性参数，β为弥散参数；
2)获取内波流的波面时程：
根据式(1)可以推导得出内波流的波面时程表示为：
<mrow>
<mi>&eta;</mi>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>&PlusMinus;</mo>
<msub>
<mi>&eta;</mi>
<mn>0</mn>
</msub>
<msup>
<mi>sech</mi>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mi>x</mi>
<mo>-</mo>
<msub>
<mi>C</mi>
<mi>p</mi>
</msub>
<mi>t</mi>
</mrow>
<mi>l</mi>
</mfrac>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，Cp为内波流的非线性速度，l为内波流的特征波长，η0为内波流的波幅；
3)获取内波流流速和加速度：
假设内波流有两个分层，上层流体流速u1与下层流体流速u2相反，上层流体厚度为h1,下层流体厚度为h2，上层流体密度为ρ1，下层流体密度为ρ2，对式(2)进行求导得到内波流的上层流体流速u1和下层流体流速u2：
然后对式(3)进行求导即可得到内波流的上层流体加速度a1和下层流体加速度a2：
4)求取内波流载荷：
根据内波流的描述方程得到内波流流速、加速度的分布，将浮式平台的将内波流载荷F分成拖曳力FD和惯性力FI两部分，即：
F＝FD+FI (5)
其中，拖曳力FD可表示为：
<mrow>
<msub>
<mi>F</mi>
<mi>D</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>C</mi>
<mi>d</mi>
</msub>
<mn>2</mn>
</mfrac>
<msup>
<mi>&rho;Au</mi>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
式中，A为受流面积；ρ为海水密度，u为内波流流速，当浮体处在上层流体中时：ρ＝ρ1，u＝u1,当浮体处在下层流体中时，ρ＝ρ2，u＝u2；Cd为拖曳力系数，根据模型试验方法或者根据规范的方法来计算获得；
惯性力FI可表示为：
FI＝CmρVa (7)
其中，V为目标浮式平台的排水量；a为内波流加速度，当浮体处在上层流体中时：a＝a1,当浮体处在下层流体中时，a＝a2；Cm为附加质量系数，根据势流理论的水动力软件计算获得。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379297.1A CN104990660B (zh) | 2015-07-01 | 2015-07-01 | 一种应用于浮式平台的内波流载荷获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510379297.1A CN104990660B (zh) | 2015-07-01 | 2015-07-01 | 一种应用于浮式平台的内波流载荷获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104990660A CN104990660A (zh) | 2015-10-21 |
| CN104990660B true CN104990660B (zh) | 2017-08-25 |
Family
ID=54302502
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510379297.1A Active CN104990660B (zh) | 2015-07-01 | 2015-07-01 | 一种应用于浮式平台的内波流载荷获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104990660B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109572930B (zh) * | 2019-01-31 | 2020-12-04 | 烟台中集来福士海洋工程有限公司 | 一种半潜式平台流载荷的估算方法 |
| CN112417726B (zh) * | 2020-11-20 | 2022-12-20 | 中国船舶工业集团公司第七0八研究所 | 一种远程智能浮式平台锚泊定位安全评估系统 |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101368382B (zh) * | 2008-10-09 | 2011-04-06 | 中国海洋石油总公司 | 一种深水半潜式钻井平台甲板可变载荷的计算与验证方法 |
| CN101475048B (zh) * | 2009-01-13 | 2011-08-10 | 中国海洋大学 | 一种新型海洋深水浮筒平台 |
| JP5440772B2 (ja) * | 2009-09-29 | 2014-03-12 | 国立大学法人三重大学 | アンカー荷重計着脱システム、アンカー荷重計共設置セット、アンカー荷重計着脱ユニット、およびアンカーへの荷重計着脱方法 |
| CN101949404B (zh) * | 2010-08-27 | 2012-11-21 | 中国计量学院 | 一种圆柱体波流载荷的消减方法 |
| CN102680270B (zh) * | 2012-06-07 | 2014-08-06 | 哈尔滨工程大学 | 不平衡激励力及不平衡弯矩联合作用时设备对船体激励载荷的间接测量方法 |
-
2015
- 2015-07-01 CN CN201510379297.1A patent/CN104990660B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104990660A (zh) | 2015-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103295078A (zh) | 深水悬链线系泊缆的结构设计优化方法 | |
| CN104343445B (zh) | 识别气层与水层、高气油比储层及含气饱和度的方法 | |
| CN104990660B (zh) | 一种应用于浮式平台的内波流载荷获取方法 | |
| Simpson et al. | Internal seiche modes and bottom boundary‐layer dissipation in a temperate lake from acoustic measurements | |
| CN102393873A (zh) | 一种内河船舶航行安全距离的计算方法 | |
| CN105005695A (zh) | 一种用于时域疲劳分析的波浪散布图组块等效方法 | |
| Aristizábal et al. | Mechanisms driving stratification in Delaware Bay estuary | |
| Xing et al. | The influence of mixing length formulation and stratification upon tidal currents in shallow seas | |
| Yen et al. | Investigation of cylinder planing on a flat free surface | |
| CN102998083B (zh) | 通过自升式钻井平台风洞试验获得平台实际风载荷的方法 | |
| Yang et al. | Viscous drift force and motion analysis of semi-submersible in storm sea states compared with model tests | |
| CN112946737A (zh) | 一种利用纵横波速度增量交汇图识别天然气水合物的方法 | |
| Song et al. | Study on vessel–pipe coupling dynamic behavior under regular waves in deep sea mining process | |
| Fujiwara et al. | VIM model test and assessment on a semi-submersible type floater with different column intervals | |
| CN106802978A (zh) | 一种自升式钻井平台环境图谱绘制方法 | |
| Sun et al. | CFD analysis on hydrodynamic characteristics for optimizing Torpedo anchors | |
| CN104359700B (zh) | 一种fpso原油外输系统耦合特性试验方法 | |
| CN106991228B (zh) | 重力锚水平承载力的设置方法 | |
| CN110096744B (zh) | 一种基于船长波长比的波浪富裕深度计算方法 | |
| CN109726421A (zh) | 基于相长相消的圆柱阵列波浪力幅值包络线的获取方法 | |
| CN110096745B (zh) | 一种基于波浪遭遇周期的波浪富裕深度计算方法 | |
| Anundsen | Operability comparison of three ultra-deepwater and harsh environment drilling vessels | |
| Zhai et al. | Numerical simulation of the dynamic behavior of deep-water semi-submersible platform under wind and waves | |
| CN102424096A (zh) | 预测深水半潜式钻井平台慢漂运动的方法 | |
| Zeng et al. | Shallow water effects on ship-generate d waves |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| CB03 | Change of inventor or designer information |
Inventor after: Xie Bin Inventor after: Li Yang Inventor after: Wang Junrong Inventor after: Xie Wenhui Inventor after: Su Jing Inventor after: Zeng Hengyi Inventor after: Wang Shisheng Inventor after: Du Qinggui Inventor after: Zhu Xiaosong Inventor after: Zhang Wei Inventor before: Xie Bin Inventor before: Wang Junrong Inventor before: Xie Wenhui Inventor before: Zeng Hengyi Inventor before: Wang Shisheng Inventor before: Du Qinggui Inventor before: Zhu Xiaosong Inventor before: Zhang Wei Inventor before: Li Yang |
|
| CB03 | Change of inventor or designer information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee after: CNOOC research institute limited liability company Patentee after: China Offshore Oil Group Co., Ltd. Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No. Co-patentee before: CNOOC Research Institute Patentee before: China National Offshore Oil Corporation |
|
| CP01 | Change in the name or title of a patent holder |