CN107463754A - 用于flng装置液化工艺的海上边界条件仿真模拟方法 - Google Patents

用于flng装置液化工艺的海上边界条件仿真模拟方法 Download PDF

Info

Publication number
CN107463754A
CN107463754A CN201710701931.8A CN201710701931A CN107463754A CN 107463754 A CN107463754 A CN 107463754A CN 201710701931 A CN201710701931 A CN 201710701931A CN 107463754 A CN107463754 A CN 107463754A
Authority
CN
China
Prior art keywords
msub
mrow
equation
hull
mover
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.)
Granted
Application number
CN201710701931.8A
Other languages
English (en)
Other versions
CN107463754B (zh
Inventor
李焱
喻西崇
王春升
谢彬
李玉星
王清
朱小松
程兵
韩旭亮
王世圣
赵晶瑞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201710701931.8A priority Critical patent/CN107463754B/zh
Publication of CN107463754A publication Critical patent/CN107463754A/zh
Application granted granted Critical
Publication of CN107463754B publication Critical patent/CN107463754B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于包括以下步骤:(1)建立海况条件下工质在FLNG的液化工艺系统内部流动的三维控制体单元模型,分析得到工质流动方向的动量方程;(2)基于实际作业海域的海况条件基础数据,得到船体或平台的三维运动方程;(3)计算得到工质重心的平移运动方程和角运动方程;(4)建立一维流体流动方程,用于表征实际海况条件下工质在FLNG的液化工艺系统内的流动;(5)根据一维流体流动方程计算海况边界条件参数,并加载入FLNG装置液化工艺动态模型中,实现实际海况环境在FLNG装置液化工艺中的应用。本发明可以广泛应用于提高FLNG装置液化工艺仿真模拟的真实性和预测应急工况的准确性。

Description

用于FLNG装置液化工艺的海上边界条件仿真模拟方法
技术领域
本发明属于海洋石油工程领域,特别是关于一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法。
背景技术
在深水气田开发中,由于距离海岸远,传统管道输送模式投资建造存在难度大、流动安全风险较高等问题,海洋工程界提出一种新型的深水气田开发模式FLNG装置,该装置将天然气的预处理、液化、储存、装卸等功能融为一体,成为解决深远海气田开发以及边际气田开发的首选设备。
然而,在设计FLNG装置上部液化工艺模块过程中,由于海上的特殊环境,在进行液化工艺仿真设计过程中无法准确描述海况条件对工艺和设备性能造成的影响,外界晃动在工艺动态仿真过程中也没有有效的方法对其进行边界定义。现有方法设计出的FLNG工艺流程为了满足海上安全性的要求,其工艺设计裕量大,空间占用大,从而导致FLNG上部液化工艺模块的投资占比高,总体费用不够经济有效。而且,单纯依靠试验进行数据采集和分析,需要耗费大量时间,经济成本过高。
发明内容
针对上述问题,本发明的目的是提供一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法,可以有效实现液化工艺的动态仿真研究和海况适应性研究,进而提高FLNG装置液化工艺仿真模拟的真实性和预测应急工况的准确性。
为实现上述目的,本发明采取以下技术方案:一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于包括以下步骤:(1)建立海况条件下工质在FLNG的液化工艺系统内部流动的三维控制体单元模型,分析海况条件在工质流动过程中的作用方式,得到海况环境及船体或平台运动条件下工质流动方向的动量方程;(2)基于实际作业海域的海况条件基础数据,得到船体或平台的三维运动方程;(3)对步骤(2)中得到的船体或平台的三维运动方程进行六自由度求解,得到工质重心的平移运动方程和角运动方程;(4)将得到的工质重心的平移运动方程和角运动方程代入步骤(1)中得到的海况环境下工质流动方向的动量方程进行惯性项求解,并考虑工质流动方向与船体或平台运动方向的夹角因素,建立一维流体流动方程,用于表征实际海况条件下工质在FLNG的液化工艺系统内的流动;(5)根据得到的一维流体流动方程,计算海况边界条件参数,并加载入FLNG装置液化工艺动态模型中,实现实际海况环境在FLNG装置液化工艺中的应用。
所述步骤(1)中,海况环境下工质流动方向的动量方程的建立过程,包括以下步骤:
(1.1)建立海况条件下工质在液化工艺系统内部流动的三维控制体单元模型;
(1.2)根据建立的三维控制体单元模型,分析工质所受重力、摩擦力和压力的作用力的影响,并计算得到海况条件及船体或平台无运动条件下工质流动方向的动量方程;
船体或平台无运动条件下工质流动方向的动量方程为:
式中,v为无运动条件下的工质流动方向x的速度,t为时间,P为压强,ρ为工质密度,g为重力加速度,f为达西摩擦系数,α为工质流动方向x与水平方向夹角;
(1.3)根据步骤(1.2)得到的船体或平台无运动条件下工质流动方向的动量方程,以及船体或平台运动对海况条件下工质流动的作用影响,得到海况条件及船体或平台运动条件下工质流动方向的动量方程:
式中,vm为工质因船体或平台运动的流动速度,v为工质流动方向x的速度,△η为船体或平台运动产生的工质流动方向位移量,△t为单位时间间隔,其中η=f(ω,t),ω为船体或平台运动速度,t为时间,γ为船体或平台运动方向与工质流动方向的夹角。
所述步骤(2)中,所述船体或平台的三维运动方程为:
式中,i,j=1,2,...,6,mij为船体或平台等系泊浮体质量矩阵,μij(ω)为附加质量系数矩阵,λij(ω)为阻尼系数矩阵,(Cij)为静水恢复力系数矩阵,FW为波浪力,ηij为船体或平台等系泊浮体六自由度运动响应,分别表示纵荡、横荡、垂荡、横摇、纵摇和艏摇运动。
所述步骤(3)中,工质重心的平移运动速度方程为:
式中,为工质重心的平移运动速度,m为工质的质量,是重力的受力矢量。
所述步骤(3)中,工质重心的角运动速度方程为:
式中,为工质重心的角运动速度,L为惯性张量,为角速度矢量,为工质力矩矢量,且
式中,是重心力矩矢量,R是一个转化矩阵,展开为
式中,Cx=cosx,Sx=sinx;φ、θ、ψ分别代表六自由度中x、y、z方向的旋转。
所述步骤(4)中,建立的一维流体流动方程为:
式中,vm为工质因船体或平台运动的流动速度,v为工质流动方向x的速度,η为系泊浮体六自由度运动响应,γ为工质流动方向与船体或平台运动方向夹角。
本发明由于采取以上技术方案,其具有以下优点:1、本发明有效搭建了三维船体或平台运动与一维液化工艺仿真之间的桥梁,能够将复杂的海况环境进行大幅度简化用于快速计算运动对FLNG装置天然气液化工艺的动态影响,通过实际作业海域的海况环境基本参数,能够快速计算液化工艺在该海域不同环境条件下的工艺适用性,缩短了计算时间并提高了计算可靠性。2、本发明为FLNG等浮式液化工艺的设计和关键设备选型提供了有利的环境条件参考,能够为工艺设计的环境影响参数提供可靠的计算依据和结果预测。本发明可以使FLNG装置的工艺动态仿真与海上真实环境相关联,可以快速高效的对海况作用条件对工艺性能的影响进行评估,能够较好的解决工程中遇到的现实问题,为FLNG装置工艺类型比选提供有效的海上适应性依据,非常适用于前期研究。
附图说明
图1是本发明的流程结构示意图;
图2是本发明海上环境工质流体流动受力图。
具体实施方式
下面结合附图和实施例对本发明的进行详细的描述。
如图1所示,本发明提供的一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法,包括以下步骤:
(1)建立海况条件下工质在FLNG的液化工艺系统内部流动的三维控制体单元模型,分析海况条件在工质流动过程中的作用方式,得到海况环境及船体或平台运动条件下工质流动方向的动量方程;
(2)基于实际作业海域的复杂海况条件基础数据(如风、浪、流等因素),获得船体或平台的三维运动情况,得到船体或平台的三维运动方程;
(3)对步骤(2)中得到的船体或平台的三维运动方程进行六自由度求解,得到工质重心的平移运动方程和角运动方程;
(4)将步骤(3)中得到的工质重心的平移运动方程和角运动方程代入步骤(1)中得到的海况环境下工质流动方向的动量方程进行求解,并考虑工质流动方向与船体或平台运动方向的夹角因素,建立一维流体流动方程,用于表征实际海况条件下工质在FLNG的液化工艺系统内的流动;
(5)根据步骤(4)中得到的一维流体流动方程,计算海况边界条件参数,并加载入FLNG装置液化工艺动态模型中,实现实际海况环境在FLNG装置液化工艺中的应用。
上述步骤(1)中,海况环境下工质流动方向的动量方程的建立过程,包括以下步骤:
(1.1)建立海况条件下工质在液化工艺系统内部流动的三维控制体单元模型。
(1.2)根据建立的三维控制体单元模型,分析工质所受重力、摩擦力和压力的作用力的影响,并计算得到海况条件及船体或平台无运动条件下工质流动方向的动量方程。
如图2所示,为海况条件下工质的流体流动受力图。图中工质所受的重力作用G为
式中,ρ为工质密度,g为重力加速度,x为工质流动方向,A为工质的截面积。
工质截面上所受的压力F为:
式中,P为压强。
工质外表面上所受的摩擦力Ff为:
Ff=-τ0πDδx (3)
式中,D为工质流动当量直径,τ0为剪应力,且
式中,f为达西摩擦系数,V为工质的流动速度。
根据式(1)~(4)即可得到船体或平台无运动条件下工质流动方向的动量方程为:
式中,v为无运动条件下的工质流动方向x的速度,t为时间,α为工质流动方向x与水平方向夹角。
(1.3)根据步骤(1.2)得到的船体或平台无运动条件下工质流动方向的动量方程,以及船体或平台运动对海况条件下工质流动的作用影响,得到海况条件及船体或平台运动条件下工质流动方向的动量方程。
由图2可知,工质受船体或平台运动的作用力为Fη,因而得到海况环境及船体或平台运动条件下工质流动方向的动量方程为:
式中,vm为工质因船体或平台运动的流动速度,△η为船体或平台运动产生的工质流动方向位移量,△t为单位时间间隔,其中η=f(ω,t),ω为船体或平台运动速度,t为时间,γ是船体或平台运动方向与工质流动方向x的夹角。
上述步骤(2)中,基于复杂海况条件(风、浪、流等因素)基础数据所获得的船体或平台的三维运动方程为:
式中,i,j=1,2,...,6,mij为船体或平台等系泊浮体质量矩阵,μij(ω)为附加质量系数矩阵,λij(ω)为阻尼系数矩阵,(Cij)为静水恢复力系数矩阵,FW为波浪力,η为船体或平台等系泊浮体六自由度运动响应,分别为η的一阶形式和二阶形式,分别表示纵荡、横荡、垂荡、横摇、纵摇和艏摇运动。
上述步骤(3)中,对船体或平台六自由度运动进行简单的方程求解,获得以工质重心的平移运动与角运动方程,其中平移运动时平移运动速度方程为:
式中,为工质重心的平移运动速度,m为工质的质量,是重力的受力矢量。
其中角运动时角运动速度方程为:
式中,为工质重心的角运动速度,L为惯性张量,为角速度矢量,为工质力矩矢量,且
式中,是重心力矩矢量,R是一个转化矩阵,展开为
式中,Cx=cosx,Sx=sinx。φ、θ、ψ分别代表六自由度中x、y、z方向的旋转。
在上述步骤(4)中,考虑工质流动方向方向与船体或平台运动方向的夹角因素,将带入方程(5)和方程(6)中,得到一维流体流动方程为:
式中,vm为工质因船体或平台运动的流动速度,v为工质流动方向x的速度,η为系泊浮体六自由度运动响应,γ为工质流动方向x与船体或平台运动方向的夹角,根据实际运动情况结合方程(8)和(9)进行求解。
上述实施例中,步骤(2)和步骤(3)中,采用计算机语言FORTRAN编写程序进行计算。步骤(4)中,采用计算机语言C#编写程序进行计算。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (6)

1.一种用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于包括以下步骤:
(1)建立海况条件下工质在FLNG的液化工艺系统内部流动的三维控制体单元模型,分析海况条件在工质流动过程中的作用方式,得到海况环境及船体或平台运动条件下工质流动方向的动量方程;
(2)基于实际作业海域的海况条件基础数据,得到船体或平台的三维运动方程;
(3)对步骤(2)中得到的船体或平台的三维运动方程进行六自由度求解,得到工质重心的平移运动方程和角运动方程;
(4)将得到的工质重心的平移运动方程和角运动方程代入步骤(1)中得到的海况环境下工质流动方向的动量方程进行惯性项求解,并考虑工质流动方向与船体或平台运动方向的夹角因素,建立一维流体流动方程,用于表征实际海况条件下工质在FLNG的液化工艺系统内的流动;
(5)根据得到的一维流体流动方程,计算海况边界条件参数,并加载入FLNG装置液化工艺动态模型中,实现实际海况环境在FLNG装置液化工艺中的应用。
2.如权利要求1所述的用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于:所述步骤(1)中,海况环境下工质流动方向的动量方程的建立过程,包括以下步骤:
(1.1)建立海况条件下工质在液化工艺系统内部流动的三维控制体单元模型;
(1.2)根据建立的三维控制体单元模型,分析工质所受重力、摩擦力和压力的作用力的影响,并计算得到海况条件及船体或平台无运动条件下工质流动方向的动量方程;
船体或平台无运动条件下工质流动方向的动量方程为:
<mrow> <mfrac> <mrow> <mi>d</mi> <mi>v</mi> </mrow> <mrow> <mi>d</mi> <mi>t</mi> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <mi>&amp;rho;</mi> </mfrac> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>+</mo> <mi>g</mi> <mi> </mi> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mfrac> <mi>f</mi> <mn>2</mn> </mfrac> <mi>v</mi> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>=</mo> <mn>0</mn> <mo>;</mo> </mrow>
式中,v为无运动条件下的工质流动方向x的速度,t为时间,P为压强,ρ为工质密度,g为重力加速度,f为达西摩擦系数,α为工质流动方向x与水平方向夹角;
(1.3)根据步骤(1.2)得到的船体或平台无运动条件下工质流动方向的动量方程,以及船体或平台运动对海况条件下工质流动的作用影响,得到海况条件及船体或平台运动条件下工质流动方向的动量方程:
<mrow> <msub> <mi>v</mi> <mi>m</mi> </msub> <mo>=</mo> <mi>v</mi> <mo>+</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>&amp;eta;</mi> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>=</mo> <mi>v</mi> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;gamma;</mi> <mo>;</mo> </mrow>
式中,vm为工质因船体或平台运动的流动速度,v为工质流动方向x的速度,△η为船体或平台运动产生的工质流动方向位移量,△t为单位时间间隔,其中η=f(ω,t),ω为船体或平台运动速度,t为时间,γ为船体或平台运动方向与工质流动方向的夹角。
3.如权利要求1所述的用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于:所述步骤(2)中,所述船体或平台的三维运动方程为:
<mrow> <mo>(</mo> <msub> <mi>m</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;mu;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>(</mo> <mi>&amp;omega;</mi> <mo>)</mo> <mo>)</mo> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <msub> <mi>&amp;lambda;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>(</mo> <mi>&amp;omega;</mi> <mo>)</mo> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mo>(</mo> <msub> <mi>C</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> <mo>)</mo> <mi>&amp;eta;</mi> <mo>=</mo> <msub> <mi>F</mi> <mi>W</mi> </msub> <mo>;</mo> </mrow>
式中,i,j=1,2,...,6,mij为船体或平台等系泊浮体质量矩阵,μij(ω)为附加质量系数矩阵,λij(ω)为阻尼系数矩阵,(Cij)为静水恢复力系数矩阵,FW为波浪力,ηij为船体或平台等系泊浮体六自由度运动响应,分别表示纵荡、横荡、垂荡、横摇、纵摇和艏摇运动。
4.如权利要求1所述的用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于:所述步骤(3)中,工质重心的平移运动速度方程为:
<mrow> <mover> <mover> <msub> <mi>v</mi> <mi>G</mi> </msub> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mfrac> <mn>1</mn> <mi>m</mi> </mfrac> <mi>&amp;Sigma;</mi> <mover> <msub> <mi>f</mi> <mi>G</mi> </msub> <mo>&amp;RightArrow;</mo> </mover> <mo>;</mo> </mrow>
式中,为工质重心的平移运动速度,m为工质的质量,是重力的受力矢量。
5.如权利要求1所述的用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于:所述步骤(3)中,工质重心的角运动速度方程为:
<mrow> <msub> <mover> <mover> <mi>&amp;omega;</mi> <mo>&amp;RightArrow;</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>B</mi> </msub> <mo>=</mo> <msup> <mi>L</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mrow> <mo>(</mo> <mi>&amp;Sigma;</mi> <msub> <mover> <mi>M</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>B</mi> </msub> <mo>-</mo> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>B</mi> </msub> <mo>&amp;times;</mo> <mi>L</mi> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;RightArrow;</mo> </mover> <mi>B</mi> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式中,为工质重心的角运动速度,L为惯性张量,为角速度矢量,为工质力矩矢量,且
<mrow> <mover> <msub> <mi>M</mi> <mi>B</mi> </msub> <mo>&amp;RightArrow;</mo> </mover> <mo>=</mo> <mi>R</mi> <mover> <msub> <mi>M</mi> <mi>G</mi> </msub> <mo>&amp;RightArrow;</mo> </mover> <mo>;</mo> </mrow>
式中,是重心力矩矢量,R是一个转化矩阵,展开为
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>C</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <msub> <mi>S</mi> <mi>&amp;theta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> <mo>-</mo> <msub> <mi>C</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> <mo>+</mo> <msub> <mi>C</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;theta;</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> <mo>+</mo> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;theta;</mi> </msub> <msub> <mi>S</mi> <mi>&amp;psi;</mi> </msub> <mo>-</mo> <msub> <mi>S</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;psi;</mi> </msub> </mrow> </mtd> <mtd> <mrow> <msub> <mi>C</mi> <mi>&amp;phi;</mi> </msub> <msub> <mi>C</mi> <mi>&amp;theta;</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>;</mo> </mrow>
式中,Cx=cosx,Sx=sinx;φ、θ、ψ分别代表六自由度中x、y、z方向的旋转。
6.如权利要求1所述的用于FLNG装置液化工艺的海上边界条件仿真模拟方法,其特征在于:所述步骤(4)中,建立的一维流体流动方程为:
<mrow> <msub> <mi>v</mi> <mi>m</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&amp;rho;</mi> </mfrac> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>P</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>-</mo> <mi>g</mi> <mi> </mi> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mfrac> <mi>f</mi> <mn>2</mn> </mfrac> <mi>v</mi> <mo>|</mo> <mi>v</mi> <mo>|</mo> <mo>+</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>w</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&amp;part;</mo> <mi>t</mi> </mrow> </mfrac> <mi>cos</mi> <mi>&amp;gamma;</mi> <mo>;</mo> </mrow>
式中,vm为工质因船体或平台运动的流动速度,v为工质流动方向x的速度,η为系泊浮体六自由度运动响应,γ为工质流动方向与船体或平台运动方向夹角。
CN201710701931.8A 2017-08-16 2017-08-16 用于flng装置液化工艺的海上边界条件仿真模拟方法 Active CN107463754B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710701931.8A CN107463754B (zh) 2017-08-16 2017-08-16 用于flng装置液化工艺的海上边界条件仿真模拟方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710701931.8A CN107463754B (zh) 2017-08-16 2017-08-16 用于flng装置液化工艺的海上边界条件仿真模拟方法

Publications (2)

Publication Number Publication Date
CN107463754A true CN107463754A (zh) 2017-12-12
CN107463754B CN107463754B (zh) 2020-06-02

Family

ID=60549833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710701931.8A Active CN107463754B (zh) 2017-08-16 2017-08-16 用于flng装置液化工艺的海上边界条件仿真模拟方法

Country Status (1)

Country Link
CN (1) CN107463754B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117272692A (zh) * 2023-11-21 2023-12-22 中国石油大学(华东) 一种海上天然气处理工艺晃荡适应性评价方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120296616A1 (en) * 2011-05-16 2012-11-22 Masaya Tsunoda Three-dimensional fluid simulation method
CN103745063A (zh) * 2014-01-16 2014-04-23 中国海洋石油总公司 一种基于液货船舱内剧烈晃荡荷载的液舱优化设计方法
CN103902751A (zh) * 2012-12-27 2014-07-02 中国地质大学(北京) 油藏数值模拟的边界处理方法及装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120296616A1 (en) * 2011-05-16 2012-11-22 Masaya Tsunoda Three-dimensional fluid simulation method
CN103902751A (zh) * 2012-12-27 2014-07-02 中国地质大学(北京) 油藏数值模拟的边界处理方法及装置
CN103745063A (zh) * 2014-01-16 2014-04-23 中国海洋石油总公司 一种基于液货船舱内剧烈晃荡荷载的液舱优化设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
韩旭亮等: "大型FLNG 液舱晃荡和船体耦合运动研究现状及展望", 《中国海上油气》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117272692A (zh) * 2023-11-21 2023-12-22 中国石油大学(华东) 一种海上天然气处理工艺晃荡适应性评价方法及系统
CN117272692B (zh) * 2023-11-21 2024-02-13 中国石油大学(华东) 一种海上天然气处理工艺晃荡适应性评价方法及系统

Also Published As

Publication number Publication date
CN107463754B (zh) 2020-06-02

Similar Documents

Publication Publication Date Title
Belibassakis et al. Hydrodynamic performance of flapping wings for augmenting ship propulsion in waves
Shen et al. RANS computations of added resistance and motions of a ship in head waves
CN103387038B (zh) 减小船舶横摇运动的分析方法
Faltinsen et al. Water entry of a wedge with finite deadrise angle
Zhuang et al. Numerical study on ship motion fully coupled with LNG tank sloshing in CFD method
Masuyama et al. Tacking simulation of sailing yachts with new model of aerodynamic force variation during tacking maneuver
CN111709151A (zh) 一种船舶快速数学模型的创建方法
Liang et al. Dynamic Modeling and Computer Simulation for Autonomous Underwater Vehicles with Fins.
CN104155043A (zh) 一种动力定位系统外界环境力测量方法
Qiu et al. 3D motion model for the freefall lifeboat during its launching from a moving ship
CN101615352A (zh) 航海模拟器用KaMeWa型喷水推进水翼船运动数学模型
CN107463754A (zh) 用于flng装置液化工艺的海上边界条件仿真模拟方法
CN108460206B (zh) 一种波浪滑翔器运动预测方法
Ullah et al. Effect of hydrostatic nonlinearity on the large amplitude response of a spar type floating wind turbine
Hu et al. Computing the hydrodynamic coefficients of underwater vehicles based on added momentum sources
Algarín et al. Numerical computation of the hydrodynamic coefficients on planing hulls in the six degrees of freedom
Hess et al. Naval maneuvering research and the need for shear stress measurements
Dalane et al. Nonlinear coupled hydrostatics of arctic conical platforms
Zhuang et al. The fully coupled effects of FPSO with different filling ratio tanks in CFD method
Zhang The study of 6-DOF motion simulation platform for passenger ship
Qiu et al. Numerical investigation and its application on the falling motion of freefall lifeboat
Ku et al. Development of a multibody dynamics kernel for motion analysis of a floating wind turbine
Raman-Nair et al. A model for deployment of a freefall lifeboat from a moving ramp into waves
Sen et al. Developing mathematical model for calculating forces affecting to Ship motions
Armudi et al. Numerical analysis of ship motion under linear and nonlinear waves

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information
CB02 Change of applicant information

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Applicant after: China Offshore Oil Group Co., Ltd.

Applicant after: CNOOC research institute limited liability company

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Applicant before: China National Offshore Oil Corporation

Applicant before: CNOOC Research Institute

GR01 Patent grant
GR01 Patent grant