基于油田规模和油气的分布特征，半潜式生产平台往往拥有数量不等的立管，通常包含生产立管，注水立管，气举立管以及油气输出立管。为了维持立管的正常作业，需要采用锚泊系统对生产平台的偏移进行约束控制，以避免偏移过大对立管构成威胁。

立管系统的数值分析方法与锚泊系统类似，主要分为静态方法、准静态方法、准动态方法以及动态方法。准静态方法中，波浪引起的运动作为静态力来移动平台，锚索的质量、阻尼和流体加速度有关的垂向运动和动力影响予以忽略，该方法在预报缆绳张力时，精度往往不能满足要求。准动态分析是一种时域模拟方法，在计算船体的运动时，将其分离成低频部分和波频部分；低频响应直接由时域运动方程求出，波频响应通过运动RAO线性叠加得到；锚索的特征通过悬链线方程求解获得。动态耦合方法能准确计及流载荷，运动引起锚泊系统和立管系统上的阻尼、惯性力、海底摩擦的影响以及立管与船体接触的影响，在整体运动响应分析中，可考虑系泊/立管系统对浮式结构物所产生的耦合效应，但耦合计算耗费的时间长。

在悬链线理论中，假定缆绳的重力远大于水流力，从而忽略流体力的影响，因此无论是静态方法、准静态方法还是准动态方法，均不能考虑水流对缆绳及立管系统的影响。

急 性 毒 性LD50大 于5 g/kg（大 鼠 经 口）。ADI 0.15～ 2 mg/kg（FAO/WHO，1994）。根据国家标准《食品添加剂使用标准（GB 2760）》的规定，可作为营养强化剂添加在食品之中[12]。

对于浅水系泊系统而言，系泊缆绳往往较短，采用忽略流力影响的悬链线理论模型进行求解即可满足精度要求。当水深逐渐增加时，缆绳的长度急剧增大，此时缆绳受到的惯性力、海流拖曳力和水动力阻尼等的作用大大增加，对平台的水平偏移和缆绳的低频张力将产生较大的影响；此外，水深较大时，浮式平台与锚泊系统、立管系统之间的相互作用明显，会导致大幅的动态张力。Lee和Devlin等[1]、Kim等[2]研究表明，在深水条件下，准静态、准动态或者低频和波频分离的计算方法将产生不可靠的结果，此时需要采用动态计算方法对半潜式生产平台主体/锚泊系统/立管系统进行数值模拟。

随着水深的不断增加，平台主体/系泊系统/立管系统之间的耦合作用问题逐渐突出，一些学者对该问题开展了深入的研究。Garrett[3]采用频域计算和时域耦合分析方法分别研究了浮式生产系统和系泊缆、立管的相互干扰，该系统包括16根系泊缆和20根立管；Kim[4]等采用时域全耦合动态分析方法对FPSO整体运动响应进行模拟；zhang等[5]采用DEEPC软件模拟了Cell-truss Spar平台主体/系泊系统/立管系统的耦合效应。杨敏冬[6]采用非线性耦合动态分析方法对平台和系泊缆绳、立管之间的相互作用进行了计算，数值模型考虑各系统的相互作用以及缆绳质量、阻尼等影响。

AMVT早期临床诊断困难，无明显的体征表现，常因首发腹痛就诊，但因阑尾术后腹痛症状常被忽视，造成误诊或拖延诊断，有文献报到此病误诊率达到90%～95%，彩色多普勒对于肠系膜上静脉血栓可实现早期诊断，具有较大临床价值。有报到表明超声对肠系膜上静脉血栓整体的灵敏度、特异度，分别达到89.3%，90.0%。彩色多普勒对肠系膜上静脉的内径变化，管内低回声血栓的范围，大小，血流动力学改变以及血管再通情况都能很好的显示，随访观察，还可监测受累肠壁因缺血水肿增厚的情况、腹腔积液情况，细菌通过门静脉入肝后肝内炎性改变情况，是否形成脓肿，超声可多次动态观察后期用药疗效，无辐射，安全，经济。

本文以某深海半潜式生产平台为研究对象，平台最大作业水深为2 300 m，配置63根型式多样、功能不同的立管和脐带缆组成的立管系统。采用HYDROD程序(基于三维频域势流理论)，计算了一阶波浪力、二阶低频慢漂载荷、船体运动RAO、静水恢复力、附加质量和势流阻尼等船体的水动力参数；建立主船体、锚泊系统和立管系统的数值模型；采用DEEPC软件(基于动态耦合方法)分析了各系统之间的相互作用特征，研究了立管系统对锚泊系统定位能力的影响，重点考察立管系统的刚度、惯性力、拖曳力、阻尼等参数对生产平台偏移和锚索张力的影响。

1 分析方法和理论

1.1 平台的频域运动响应

半潜式生产平台/锚链系统/立管系统组成的整个系统主要由“水动力频域计算”和“系泊定位计算”两部分组成。HYDROD程序主要计算不同浪向和频率作用下，运动RAO、附加质量、辐射阻尼、一阶波浪力等船体的水动力性能。

各船体速度势Φ满足的边界条件如下所示[7]:

 

平台在微幅波作用下的六自由度运动[8]可表示为:

 

式中:m(ω)和mA(ω)分别为船体的质量、附加质量矩阵；B为船体阻尼矩阵；K为船体总刚度矩阵；和X分别为船体的加速度、速度以及位移矢量；F为系统所受的波浪激励力。

在HYDROD中，采用远场积分方法计算二阶低频慢漂载荷。其原理为:根据动量守恒定律，通过在远方控制面上积分，得到船体纵荡、横荡和艏摇三个方向上的平均二阶力，其计算公式为:

 

式中:Cg为波浪的群速度；C为波速；α为波浪入射角；A(θ)为绕射势在远场的幅值。

1.2 平台时域耦合运动

采用DNV船级社开发的DEEPC程序对平台系统的耦合运动开展计算。基于时域条件下，半潜式生产平台主体/锚泊系统/立管系统的运动方程式可写成[9]:

1.1 TOFD（衍射时差法）超声技术是利用缺陷部位的衍射波信号来检测和测定缺陷尺寸的一种超声波检测方法，通常使用纵波斜探头，采用一发一收模式。

 

式中:M为频率质量矩阵，包括结构物质量矩阵m和频率附加质量矩阵A；C为频率势流阻尼矩阵；D1和D2分别为线性阻尼和二次阻尼矩阵；f为速度矢量函数，可表示为K为静水刚度矩阵；x为位移矢量；q为激励力矢量。频域条件下，船体势流阻尼和附加质量通过适当转化可获得时延函数，从而考虑因自由面记忆效应所产生的影响。

而外部激励力q可以表示为如下形式:

 

表10显示了有无海流情况下，立管的附加质量对锚索张力的影响。结果表明，无论是否存在海流，立管的附加质量始终起到降低缆绳张力的作用。

1.3 系泊和立管系统的时域方程

立管的数学模型与系泊缆基本相同，通过系泊和立管系统的时域方程可求解出各单根缆绳的形状、张力随时间的变化情况，其中系统的张力为平台运动方程中的qext项。t时刻，在随机波浪和流速为U的水流作用下，缆绳将呈现一定形状。假设缆绳上的任一点与原点之间的长度为s，在该点锚索切线与水平面夹角为φ。取缆绳上的一微段ds，作用在此微段上的力包括重力、浮力、流体动力和张力，这些力可分解为沿轴线方向τ和与轴线的法线垂直的方向n。

微段ds在张力T的作用下将产生伸长应变ε，可以得到沿系泊缆的轴线方向τ、垂直于轴线的法线方向n的缆绳运动微分方程组。

  

图1 锚索微段上受力分析Fig.1 Force analysis at segment of mooring line

 

式中:ρ为海水密度，g为重力加速度，A为锚索截面积，W为单位长度锚索重量，Cm为附加质量系数，D为缆绳截面直径，ε为伸长应变，CDτ和CDn为τ方向和n方向阻力系数，Vcτ和Vcn分别为锚索在τ方向和n方向的运动速度，Vrτ和Vrn为τ方向和n方向水质点和缆绳的相对运动速度。Vrτ和Vrn可表示为:

 

式中:Vwτ和Vwn分别为水质点在τ方向和n方向的速度。

第一，社区服刑人员的主观因素。一是存在文化程度差异，量表的表述不能完全理解，工作人员若协助解释则会形成暗示，影响测量效果；二是量表题量多且为客观题，社区服刑人员为节省时间随意填写，造成无效量表过多；三是多数社区服刑人员对心理常识知之甚少，都是“心盲”，主观配合程度差。

伸长应变ε可表示为:

 

式中:E为锚索弹性模数。

半潜式生产平台采用16点锚泊系统进行定位，该系统分为4组，每组4根，缆绳的布置沿中纵剖面和中横剖面对称，每一组内相邻两根系泊缆之间的夹角为5°，＃1与X轴的夹角为35°，定位绞车位于半潜式平台两舷侧，通过基线附近的导缆器将缆绳导入锚链舱内，调节锚链的长度来获取所需的预张力，如图3所示。

 

2 船体/锚泊系统/立管系统

2.1 船体水动力模型

研究作业水深为2 300 m时半潜平台的水动力性能，采用HYDROD软件计算平台的水动力系数、包括一阶波浪力、附加质量、势流阻尼、船体运动幅值响应等。其中，遭遇浪向为0°～180°，中间间隔15°，频率为0.02～2.0 rad/s，中间间隔0.02 rad/s。采用远场积分方法获得平台的二阶波浪漂移载荷。图2为用于水动力分析的船体湿表面网格。

  

图2 生产平台的水动力模型Fig.2 Hydrodynamic model of production platform

2.2 锚泊系统和立管系统

  

图3 系泊系统布置Fig.3 Arrangement of mooring system

作用在缆绳微段的流体拖曳力采用Morison公式求出:

每根锚索由三部分组成，上部与平台连接的为无档锚链(R4S无档锚链)，中部为聚酯缆(polyester rope)，下端为锚链(R4S无档锚链)，各海况下，抛出的锚索总长度为3 700 m，定位能力分析时，应考虑顶端和底端锚链的腐蚀，腐蚀量为0.5 mm/年，服役时间为30年，腐蚀总量为15 mm。分析时，假设锚索的破断负荷受腐蚀的影响，锚索的直径和单位长度重量保持不变。锚索参数见表1。

 

表1 锚索参数Tab.1 Parameters of mooring line

  

规格　直径/mm重量/(kg˙m－1)空气中　水中考虑腐蚀后的破断负荷/kN　长度m底端锚链　R4S无档锚链　165　545　474.2　21 888　300聚酯缆　279　54.6　14.0　22 563　3 200顶端锚链　R4S无档锚链　165　545　474.2　21 888　200

目标平台的立管和脐带缆数目较多(共63根)，布置方式如图2所示，其中，编号m1～m16为1号到16号系泊缆；编号r1～r63为1号到63号立管。

处理试验结果的目的在于确定各试验因素的优水平及试验范围内的优组合。正交表的选择一般是根据因素和水平的多少及试验工作量的大小而定，依据正交表的综合可比性，利用极差分析法可以非常直观简便地分析试验结果，确定因素的主次和最优组合[11-12]。

  

图4 系泊和立管系统布置Fig.4 Arrangement of mooring and riser system

立管的相关属性如表2所示。

海上风电场嵌岩区域的桩型选择，针对单桩基础，目前仍采用钢管桩。针对多桩基础，目前主要采用钢管桩加灌注桩的形式，即岩面以上部分采用钢管桩，岩面以下部分采用灌注桩。

 

表2 立管参数Tab.2 Parameters of riser

  

立管类别　生产管　气举管　注水管　外输气管　外输油管　脐带缆规定外径D/m　0.218　0.142　0.218　0.325　0.457　0.137规定壁厚t/m　0.023　0.020　0.013　0.033　0.028换算内径d/m　0.175　0.104　0.193　0.262　0.401抗拉伸刚度EA/kN　2 817 610　1 499 852　1 687 969　5 972 814　7 888 010　422 000抗弯刚度EI/(kN˙m2)　13 821　2 869　9 020　64 454　181 952　5抗扭转刚度/(kN˙m2)　10 632　2 207　6 938　49 581　139 965　243管内流体密度/(kg˙m－3)　970　202　1 026　202　970

3 环境条件和数值模型

该平台服役于南海海域，工作水深约为2 300 m。在进行锚泊定位分析计算时，波浪频谱形式为JONSWAP谱，风频谱形式为NPD谱，风、浪、流方向相同，服役海况环境条件参数如表3所示。

基于时域耦合方法，采用DeepC软件分析目标平台主体、立管和锚泊系统间的相互作用对平台的运动响应和锚索动力响应的影响。

计算了200年一遇生存工况，4种不同浪向下的平台偏移、导缆孔处锚索张力。在时域耦合分析方法中，风浪流的时历曲线均通过随机种子数生成，由计算机程序生成的随机序列来模拟随机的环境条件，取随

机种子数为1，100，1 000的三组随机数列进行模拟，对三组数据进行统计平均，作为最终的数值计算结果。

式(1)中，MI为某时段相对湿润度；P为某时段的降水量，单位为毫米(mm);PET为某时段的潜在蒸散量，用FAO Penman-Monteith方法计算，单位为毫米(mm)，算法参见文献[19]。

 

表3 目标平台的计算环境条件Tab.3 Environmental condition of platform

  

设计工况　回归期　Hs/m　Hmax/m　Tmax/s　Tp/s　Tz/s　Gama Vw/(m˙s－1)1min.average Vc(surface)/(m˙s－1)生存工况　200年一遇台风工况　14.3　24.5　14.7　15.6　11.6　2.42　56.3　2.20

采用表4所示的流剖面计算系泊和立管系统所受的流载荷，表面流计算船体受到的流载荷。

 

表4 海流剖面分布形式Tab.4 Distribution of current profile

  

深度200 年一遇/(cm˙s－1)表层(0 m)219.5中层(200 m)183.6底(2 300 m)　106.7

为考虑立管系统对系泊定位能力的影响，建立无立管的船体系统模型，采用时域耦合方法对目标平台的锚泊系统进行分析，并与包含立管系统的完整模型展开对比。

与锚泊系统相似，立管系统在平台产生漂移的情况下也会为平台提供一定的回复力，进而从一定程度上增加系统的刚度。表5为有/无立管两种数值模型下的系统刚度对比。

770 Surface modification of orthopaedic titanium implant for anti-infection and promoting fracture healing: an update

  

图5 船体和系泊系统耦合分析模型(M1)Fig.5 Coupling analysis model of hull and mooring system

  

图6 船体、系泊系统和立管系统耦合分析模型(M2)Fig.6 Coupling analysis model of hull，mooring and riser systems

4 计算结果和分析

4.1 立管对系统刚度的影响

数值模型共分为两种，M1:由船体和16根系泊缆组成；M2:由船体、系泊缆以及全部的立管和脐带缆组成，如图5和图6所示。

通过对比发现，由于立管系统包含的立管和脐带缆的数目较多(共63根)，其对锚泊系统整体刚度的影响不可忽略，在纵荡和横荡方向上的刚度相差分别为6%和7%，在艏摇方向上更是达到了11%。因此开展系统刚度的计算时，需要计及立管系统的作用。

 

表5 有/无立管线性刚度对比Tab.5 Linear stiffness with/without risers

  

自由度(DOF)刚度有立管　无立管有立管与无立管相差纵荡/(kN˙m－1)　349.75　329.4　6.18%横荡/(kN˙m－1)　299.55　279.5　7.17%艏摇/(kN˙m˙rad－1)　3.22E＋06　2.89E＋06　11.42%

提取方法的选择。通过试验，溶剂提取法、超声波辅助提取法及微波辅助提取法的吸光值分别为0.231、0.412和0.378，故而得出，超声波辅助法提取枸杞叶黄素效果最佳。

1)由于45°条件下风浪流载荷最大，所以此时平台所受的偏移和张力最大。

2)当考虑立管系统时，系统刚度会相应增大，进而为平台提供了更大的回复力，可一定程度上降低平台的最大漂移和锚索的最大张力。

3)在0°浪向时，锚索最大张力会略大于无立管的情况。这是由于立管系统的设计取决于井口位置、立管的型式和数量以及水深等诸多因素，很难保证立管系统布置的对称性，当立管系统存在时，平台的初始平衡状态会存在稍小的偏移。

 

表6 平台最大漂移和导缆孔最大张力Tab.6 Maximum offset of platform and maximum force at fairleads

  

风浪流方向最大漂移/m有立管　无立管有立管与无立管相差最大张力/kN有立管　无立管有立管与无立管相差0°　64.47　65.34　1%　7 941.2　7 840.7　－1%15°　65.3　67.37　3%　8 246.7　8 512.9　3%30°　69.51　73.29　5%　8 803.1　9 114.8　3%45°　72.57　77.94　7%　9 038.3　9 462.3　4%

4.2 立管对系泊定位能力的影响

环境载荷对平台、锚索和立管产生的静态作用力主要包括平台上所受的风力、流力、波浪慢漂力以及锚索和立管上所受的拖曳力，如表7所示。对比可知，在该环境条件下，锚索和立管上的拖曳力占总静态力不可忽略，最大比例可达34%，其中，相比锚索，立管上的拖曳力所占的比例较为明显，最大比例可达总静态力的24%。

表8和图7显示的是各浪向条件下，锚索的张力和平台的偏移变化特征。综合考虑立管的刚度、附加质量、阻尼以及水流的拖曳载荷等因素的作用，通过对比可知:

1)考虑立管系统时，平台的最大偏移比无立管系统时增加了25%，最大张力增加了13%。

例2 (2016.烟台中考)美国心理学家桑代克通过对动物行为的大量观察和实验提出了“尝试与错误”学习理论。下列观点是在对动物行为研究基础上形成的，其中应用归纳法获得的结论是( )

忽略作用于立管系统上的流载荷，仅考察立管对系泊刚度的贡献，分析船体在风、浪、流载荷作用下，有/无立管系统对平台漂移和锚索张力的影响。通过表6的对比可知:

2)立管对系泊系统的影响主要集中在系统的刚度，附加质量和阻尼之上。立管增加了整个系统的刚度，增大了回复力，增大的回复力将降低平台的偏移和系泊缆的张力；在台风工况下，海流较大，相比系统刚度，作用于立管上的拖曳力占主导，引起平台朝来流方向产生较大运动。

 

表7 200年一遇环境条件下的风、浪、流载荷Tab.7 Wind，wave and current forces in environment condition of 200 year return period

  

环境方向/(°)平台风载荷/kN平台流载荷/kN波浪慢漂力/kN锚索和立管拖曳力/kN立管拖曳力/kN总静态力/kN锚索和立管拖曳力占百分比/(%)立管拖曳力占百分比/(%)纵荡　横荡　纵荡　横荡　纵荡　横荡　纵荡　横荡　纵荡　横荡　纵荡　横荡　纵荡　横荡　纵荡　横荡0　6 305　0　6 410　0　1 607　0　6 843　0　4 984　0　21 165　0　32.3　0.00　23.5　0.00 15　6 624 1 806 6 748 1 814 1 329　325　6 645 1 980 4 812 1 286 21 346 5 925　31.1　33.4　22.5　21.70 30　6 741 4 066 6 877 4 134　939　529　6 049 3 765 4 312 2 491 20 606 12 494 29.3　30.1　20.9　19.94 45　5 675 5 955 5 823 6 113　683　726　5 054 5 213 3 518 3 527 17 236 18 007 29.3　28.9　20.4　19.59

 

表8 平台最大漂移和锚索最大张力对比Tab.8 Maximum offset of platform and maximum force at fairleads

  

风浪流角度/(°)平台最大漂移/m　锚索最大张力/kN有立管　无立管　有立管与无立管相差/(%)　有立管　无立管　有立管与无立管相差/(%)0 105.29　84.47　24.6　10 340.75　9 194.88　12.5 15　108.18　89.96　20.3　11 088.85　10 350.62　7.1 30　116.17　100.54　15.5　12 166.85　11 398.42　6.7 45　121.32　107.26　13.1　12 654.58　12 002.70　5.4

  

图7 立管对锚索张力的影响Fig.7 Effect of risers on mooring tension

表9显示了有无海流情况下，立管粘性力对锚索张力的影响。通过分析可知，不考虑海流影响时，立管的粘性效应起到正阻尼的作用，降低了缆绳上受到的张力；当计及海流的影响且流速较大时，立管粘性效应反而导致缆绳上的张力增大。

造成上述现象的原因在于:立管上的粘性力作用效果取决于水流与缆绳运动的相对速度，当相对速度与平台偏移方向相同时，立管上的粘性力效应将使缆绳张力更大，反之，当相对速度与平台偏移方向相反时，阻碍平台的运动并降低缆绳的张力。

 

表9 立管粘性力系数对锚索张力的影响Tab.9 Effect of riser drag coefficient on mooring tension

  

来流角度/(°)不考虑海流的影响　考虑海流的影响考虑粘性力的缆绳最大张力/kN不考虑粘性力的缆绳最大张力/kN差距/(%)考虑粘性力的缆绳最大张力/kN不考虑粘性力的缆绳最大张力/kN差距/(%)0 7 808.8　8 148.8　－4.35　10 245.2　9 029.5　11.9 15　8 057.8　8 315.4　－3.20　11 004.8　9 523.4　13.5 30　8 959.6　9 185.3　－2.52　12 210.8　10 453.3　14.4 45　9 301.8　9 552.6　－2.70　12 764.0　10 788.9　15.4

式中:qWI为风力；q(1)WA为一阶波浪力；q(2)WA为二阶慢漂波浪力；qCU为海流作用力；qext为多种外力的综合作用力，包括波浪慢漂阻尼、系泊回复力以及特定力等其它外力。

 

表10 立管附加质量系数对锚索张力的影响Tab.10 Effect of riser added mass coefficient on mooring tension

  

来流角度/(°)不考虑海流的影响　考虑海流的影响考虑粘性力的缆绳最大张力/kN不考虑粘性力的缆绳最大张力/kN差距/(%)考虑粘性力的缆绳最大张力/kN不考虑粘性力的缆绳最大张力/kN差距/(%)0 10 245.2　10 251.6　－0.06　7 808.8　7 814.2　－0.07 15　11 004.8　11 019.7　－0.13　8 057.8　8 074.8　－0.21 30　12 210.8　12 229.9　－0.16　8 959.6　8 983.1　－0.26 45　12 764.0　12 783.1　－0.15　9 301.8　9 327.4　－0.28

图8和图9显示的是带立管系统状态下，平台的最大偏移和缆绳的最大张力随时间变化情况。从对比可知:无论是偏移还是缆绳张力，其时历曲线均表现为较明显的低频大幅度振荡，其频谱在频率ω为0 rad/s附近处的能量最大，低频载荷的影响较为明显；此外，在大的振荡曲线上，还出现了众多周期较短、幅值较小的振荡运动，同时，频谱曲线上，频率ω为0.4 rad/s附近再次出现了能量较小的聚集区域，这是由平台在波频范围内的运动响应引起的。因此，尽管二阶波浪载荷相比一阶波浪载荷是一个小量，但其对半潜式锚泊系统的定位性能影响起到主导作用，是平台发生较大偏移、系泊缆受较大张力的主因。

  

图8 风浪流为0°时平台的最大偏移时历和频谱曲线(有立管)Fig.8 The maximum offset and corresponding spectrum at 0°of incident direction(with risers)

  

图9 风浪流为0°时缆绳的最大张力时历和频谱曲线(有立管)Fig.9 The maximum tension and corresponding spectrum at 0°of incident direction(with risers)

5 结 语

基于三维频域势流理论，计算船体的水动力参数，采用动态耦合方法分析了船体/锚泊系统/立管系统之间的相互作用特征，重点考察立管系统的刚度、惯性力、拖曳力等参数对生产平台偏移和锚索张力的影响，可获得以下结论:

是夜，虚弱的、奄奄一息的我看见西边林场边的树丛里，飘起缕缕的亮光，那亮光像北极光，它们飘散，星星点点地朝我飘来，在我旁边，那团白光聚集成范峥峥的人形，她幽幽地对我说，带我回家……

在风、浪、流载荷同时存在的情况下，立管对锚索张力的影响很大，立管系统将增大整个系统的刚度，降低平台的运动。

立管上粘性力的作用效果取决于水流与缆绳运动的相对速度，当相对速度与平台偏移方向相同时，立管上的粘性力效应将使缆绳张力更大，反之，阻碍平台的运动并降低缆绳的张力。而立管的附加质量始终起到降低缆绳张力的作用。

平台服役于南海海域，该海域的水流速度较大，立管的拖曳载荷随之增加，尤其对“配备众多立管的半潜式生产平台”而言，海流拖曳力的影响较为明显，立管的存在导致平台最大偏移量增大了25%，系泊缆的最大张力增大了近13%。

因此，对于流速较大的海域，带有多立管的平台系泊系统设计不能忽略立管的影响。

参考文献:

[1] LEE M Y，DEVLIN P.Development of API RP2SM for synthetic fiber rope moorings[C]//Proceeding of the Offshore Technology Conference.2000:OTC12178.

[2] KIM M H，TAHAR A，KIM Y B.Variability of spar motion analysis against various design methodologies/parameters[C]//Proceedings of the Twentieth Offshore Mechanics and Artic Engineering Coference.OMAE-1063.

[3] GARRETT D L.Coupled analysis of floating production systems[J].Ocean Engineering， 2005， 32:802-816.

[4] KIM M H，KOO B J，MERCIER R M，et al.Vessel/mooring/riser coupled dynamic analysis of a turret-moored FPSO compared with OTRC experiment[J].Ocean Engineering， 2008， 35:1 676-1 685.

[5] ZHANG F，YANG J M，LI R P，et al.Coupling effects for cell-truss spar platform:comparison of frequency-and time-domain analyses with model tests[J].Journal of Hydrodynamics， 2008， 20(4):424-432.

叶霭玲看我不说话，又多心了。她用肩膀蹭了我一下，说，嗨，嗨，瞧你魂不守舍的样子，小心叫人把小魂儿勾走了吧？

[6] 杨敏冬.深水浮式结构与系泊/立管系统的全时域非线性耦合动态分析[D].大连:大连理工大学，2012.(YANG M D.Fully time-domain nonlinear coupled dynamic of floating structure/mooring/riser system in deep water[D].Dalian:Dalian University of Technology，2012.(in Chinese))

[7] 戴遗山.舰船在波浪中运动的频域与时域势流理论[M].北京:国防工业出版社，1998.(DAI Y S.Potential flow theory of ship motions in water[M].Beijing:National Defernce Press，1998.(in Chinese))

[8] Bureau Veritas.Hydrostar for experts user manual[R].2012.

[9] Marintek Report.SIMO-theory manual version 3.6[R].

由表4可看出，基于13个形态指标划分的缺齿蓑藓8个聚类群枝叶均为长条形，先端渐尖至锐尖，少数钝尖，枝叶中上部疣发达，基部光滑无疣。但方差分析结果显示，各聚类群在叶片长度、叶形、上中部细胞疣的发达程度、基部细胞是否微疣、中肋是否突出叶尖等方面存在明显分化(P < 0.01)。