CN114442481B - 针对模型不确定的海洋柔性立管边界振动控制方法 - Google Patents

针对模型不确定的海洋柔性立管边界振动控制方法 Download PDF

Info

Publication number
CN114442481B
CN114442481B CN202210371917.7A CN202210371917A CN114442481B CN 114442481 B CN114442481 B CN 114442481B CN 202210371917 A CN202210371917 A CN 202210371917A CN 114442481 B CN114442481 B CN 114442481B
Authority
CN
China
Prior art keywords
flexible riser
control
boundary
riser
vibration
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.)
Expired - Fee Related
Application number
CN202210371917.7A
Other languages
English (en)
Other versions
CN114442481A (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.)
Southwest Petroleum University
Original Assignee
Southwest Petroleum University
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 Southwest Petroleum University filed Critical Southwest Petroleum University
Priority to CN202210371917.7A priority Critical patent/CN114442481B/zh
Publication of CN114442481A publication Critical patent/CN114442481A/zh
Application granted granted Critical
Publication of CN114442481B publication Critical patent/CN114442481B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Earth Drilling (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了针对模型不确定的海洋柔性立管边界振动控制方法，对柔性立管建模、得到柔性立管系统的控制方程，确定柔性立管系统的边界条件和初始条件；基于自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律；获取柔性立管的横向振动实际偏移量；将所述实际偏移量代入所述自适应模糊反演控制律，得到实际控制输出；基于所述实际控制输出，在顶端对柔性立管进行控制。本发明用以解决现有技术中的控制方法依赖海洋柔性立管的精确数学模型所导致的控制精度较低的缺陷，实现不依赖精确模型，能够很好的补偿柔性立管系统的不确定性和外部荷载的不确定性，同时还能有效地跟踪和调节振动状况、主动抑制立管振动的目的。

Description

针对模型不确定的海洋柔性立管边界振动控制方法
技术领域
本发明涉及海洋油气开采领域，具体涉及针对模型不确定的海洋柔性立管边界振动控制方法。
背景技术
随着海洋油气资源勘探力度的逐步加强与钻井深度的不断加深，海洋油气开发环节中所需解决的技术难题也越来越多。对于浮式平台/船舶或半潜式平台/船舶而言，柔性立管是连接海底油井和海面平台的重要通道，也是海上石油输送的薄弱环节。由于其特殊的海洋工作环境，柔性立管的振动现象难以避免，而振动将使立管疲劳，并缩短其工作寿命、提高其生产成本并可能带来致命的安全生产风险。
对于海洋油气开采而言，其柔性立管系统因较大的长径比通常被简化为典型的欧拉伯努利梁模型，其动力学特性由无限维偏微分方程和有限维常微分方程共同表示。其柔性结构由无穷多个模态表示，其结构参数也可能随立管振动发生变化，因此对其控制设计具有较大难度。
为解决本质上是分布式参数系统的柔性立管系统的振动抑制问题，现有技术逐渐开始从对柔性立管系统主动控制的角度出发进行研究。其中边界控制方案既能避免模型截断引起的溢出效应，又具有经济和可执行性。
然而现有的控制算法依赖海洋柔性立管的精确数学模型，由于该系统实际存在的复杂性、非线性、时变性、不确定性等，在实际工程中难以对柔性立管系统状态分布式荷载和边界干扰进行有效建模，无法获得精确的数学模型，导致难以克服系统的各种不确定性干扰，控制精度始终较低。
发明内容
本发明提供针对模型不确定的海洋柔性立管边界振动控制方法，以解决现有技术中的控制方法依赖海洋柔性立管的精确数学模型所导致的控制精度较低的缺陷，实现不依赖精确模型，能够很好的补偿柔性立管系统的不确定性和外部荷载的不确定性，同时还能有效地跟踪和调节振动状况、主动抑制立管振动的目的。
本发明通过下述技术方案实现：
针对模型不确定的海洋柔性立管边界振动控制方法，包括：
对柔性立管建模、得到柔性立管系统的控制方程，确定柔性立管系统的边界条件和初始条件；
基于自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律；
获取柔性立管的横向振动实际偏移量；
将所述实际偏移量代入所述自适应模糊反演控制律，得到实际控制输出；
基于所述实际控制输出，在顶端对柔性立管进行控制。
针对现有技术中控制方法依赖海洋柔性立管的精确数学模型，导致控制精度较低的问题，本发明提出一种针对模型不确定的海洋柔性立管边界振动控制方法，本方法首先建立柔性立管模型，得到整个系统的控制方程，并确定其边界条件和初始条件。然后引入自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律。本案发明人在研究过程中发现，海洋柔性立管系统可看作是包含各种不确定性的非线性系统，因此本申请采用自适应模糊反演算法，该算法能够将反演方法与模糊控制方法相结合，可以有效解决海洋柔性立管在特殊工况下的非线性的跟踪控制难题，并且基于该算法所建立的自适应模糊反演控制律，能够估计与系统性能有关的不确定因素，以此克服或补偿系统的不确定性和各种干扰。本方法在得到自适应模糊反演控制律后，可直接获取柔性立管的横向振动实际偏移量，并将其代入，得到实际所需的控制输出、定义为实际控制输出，以该实际控制输出来对海洋钻采平台的柔性立管从顶端进行控制即可，其中柔性立管的顶端即是与钻采平台或钻采船所连接的一端，即是立管远离海底水下井口的一端，对该端部进行控制，可方便对控制器的操作安装甚至是后期维护。
综上，本申请基于反演方法并引入自适应模糊系统，以此来逼近对海洋柔性立管控制过程中的未知非线性项，补偿系统的不确定性和海洋环境干扰，并且本方法的自适应参数有且只有一个，减小了计算量，提高了计算速度，能有效的抑制海洋柔性立管的振动。并且，本申请不依赖海洋柔性立管系统的精确模型，具有很好的鲁棒性和自适应性，从而能很好的补偿立管系统的不确定性和外部荷载的不确定性，同时还能很好地跟踪和调节振动状况。此外，本申请采用自边界主动控制的控制方式，被控驱动装置只需在海洋柔性立管顶端实施控制作用，即可实现抑制立管振动的效果，与现有的被动控制方式相比，安装维护方便，性价比显著提高。
进一步的，得到所述控制方程的方法包括：
基于动力学分析得到柔性立管系统的动能E k、势能E P、以及非保守力对海洋柔性立管所做的总虚功W 0；
基于哈密顿原理、通过如下公式得到柔性立管系统的控制方程：
Figure 339715DEST_PATH_IMAGE001
；
式中，δ为变分操作符，t 1、t 2分别为两个端点时刻，t为时间变量。
本方案所提出的控制方程依据动力学分析结果得出，充分考虑了柔性立管在海洋环境中的特殊工况，具有较强的代表性。
进一步的，所述柔性立管系统的控制方程为：
Figure 187805DEST_PATH_IMAGE002
所述边界条件为：
Figure 802457DEST_PATH_IMAGE003
所述初始条件为：
Figure 646916DEST_PATH_IMAGE004
；
式中：t为时间变量、x为空间变量，L为立管长度，w(x,t)为横向振动实际偏移量，f(x,t)为作用在立管上的分布式海流载荷，d(t)为边界扰动，M s为船舶质量，ρ为单位立管质量，EI为立管弯曲刚度，T为立管张力，c为结构阻尼系数，d s为船舶阻尼系数，δ为变分运算符；u (t)为作用在立管顶端的边界控制量；
Figure 598692DEST_PATH_IMAGE005
为w(x,t)对t的一阶偏导，
Figure 477786DEST_PATH_IMAGE006
为w(x,t)对t的二阶偏导；
Figure 856815DEST_PATH_IMAGE007
为w(x,t)对x的二阶偏导，
Figure 454149DEST_PATH_IMAGE008
为w(x,t)对x的四阶偏导；
Figure 209616DEST_PATH_IMAGE009
为w(x,t)中x=L时对x的一阶偏导，
Figure 208796DEST_PATH_IMAGE010
为w(x,t)中x=L时对x的二阶偏导，
Figure 899671DEST_PATH_IMAGE011
为w(x,t)中x=L时对x的三阶偏导，
Figure 843356DEST_PATH_IMAGE012
为w(x,t)中x=L时对t的一阶偏导，
Figure 12301DEST_PATH_IMAGE013
为w(x,t)中x=L时对t的二阶偏导，
Figure 990621DEST_PATH_IMAGE014
为w(x,t)中x=0时对x的一阶偏导，w(x,0)为w(x,t)中t=0时的值，
Figure 852398DEST_PATH_IMAGE015
为w(x,t) 中t=0时对t的一阶偏导。
进一步的，基于自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律的方法包括：
建立自适应模糊逻辑系统φ 2：
Figure 283379DEST_PATH_IMAGE016
；
式中，
Figure 250155DEST_PATH_IMAGE017
为模糊基向量的转置，θ2为逼近向量；
确定边界控制律：
Figure 817403DEST_PATH_IMAGE018
式中，u(t)为作用在立管顶端的边界控制量，λ 2为第一控制参数，e 1为海洋柔性立管边界实际位置与边界理想位置的误差，e 2为海洋柔性立管边界横向振动偏移量变化率与虚拟控制信号之间的误差，x 2为海洋柔性立管边界横向振动偏移量变化率，y为海洋柔性立管边界实际位置，y d为立管理想位置，α 11为虚拟控制信号；
根据所述边界控制律，确定自适应律：
Figure 115660DEST_PATH_IMAGE019
式中，
Figure 909303DEST_PATH_IMAGE020
为自适应参数，r 2为第二控制参数，k 2为第三控制参数，
Figure 544684DEST_PATH_IMAGE021
为模糊基向量，θ2为逼近向量。
本方法通过模糊系统φ 2来逼近未知非线性函数，并基于确定的边界控制律得到最终所需的自适应律。
进一步的，所述自适应模糊逻辑系统在其模糊隶属函数的论域两端采用Z形隶属函数，在论域中间区域采用三角形隶属函数。
本方案在设置模糊隶属函数时，需满足以下规则：
①隶属函数应是对称的；
②隶属函数应遵守日常的语义顺序；
③论域中的每个点至少属于一个并且不能超过两个隶属函数的区域;
④当两个隶属函数互相重叠时，重叠的区域中不能同时包含两个隶属函数的最大隶属度；
⑤当两个隶属函数互相重叠时，在重叠区域中的任意元素所对应的两个隶属函数值的和必须小于或等于1。
本方案在模糊隶属函数的论域两端采用Z形隶属函数，中间区域采用三角形隶属函数，不仅可满足上述规则，还充分考虑了模糊化、模糊推理与解模糊过程的精确性与简便性，显著提高了运算效率。
进一步的，在建立所述自适应模糊反演控制律后，还通过李雅普诺夫稳定性理论，验证系统是否稳定有界：若是，进入下一步；若否，重新调整第一控制参数和/或第二控制参数和/或第三控制参数。
本方案提出了在柔性立管应用领域内对系统稳定进行验证的有效方法，保证了在采用自适应模糊反演算法时，整个柔性立管系统的稳定性，为后续的控制精确度提供了充分保障。
进一步的，验证系统是否稳定有界的方法包括：
对于整个柔性立管系统，设计李雅普诺夫函数；
确定李雅普诺夫函数的上下界；
根据李雅普诺夫函数对时间的导数得到如下不等式：若满足该不等式，认为系统稳定有界；若不满足该不等式，认为系统不是稳定有界：
Figure 107384DEST_PATH_IMAGE022
式中，t为时间变量、x为空间变量，w(x,t)为立管横向振动偏移量，L为立管长度，T为立管张力，υ、υ 1均为大于0的常量，ε 0为有界干扰。
进一步的，所述柔性立管的横向振动实际偏移量通过激光位移传感器实时获取。此种获取方式在作业现场具有成本低廉、传感器数量少、工程量低、实用性高等优点。
进一步的，在顶端对柔性立管进行控制的方法为：根据实际控制输出，驱动位于柔性立管顶端的执行器，由所述执行器向柔性立管施加控制作用力。其中执行器的具体输出方式在此不做限定，本领域技术人员能够实现的执行机构均可适用于本方案中。
进一步的，作用在立管上的分布式海流载荷f(x,t)通过如下公式表达：
Figure 435597DEST_PATH_IMAGE023
式中，ρ s为海水密度，C D为阻力系数，D为立管外径，A D 为阻力的振荡部分幅值，f v 为涡旋脱落频率，U(x,t)为海洋洋流速度与深度关系，θ为分布式海流的相位角。
本发明与现有技术相比，具有如下的优点和有益效果：
1、本发明针对模型不确定的海洋柔性立管边界振动控制方法，可以有效解决海洋柔性立管在特殊工况下的非线性的跟踪控制难题，能够估计与系统性能有关的不确定因素，以此克服或补偿系统的不确定性和各种干扰；并且本方法的自适应参数少，从而减小了计算量、提高了计算速度、能快速抑制海洋柔性立管的振动。
2、本发明针对模型不确定的海洋柔性立管边界振动控制方法，不依赖海洋柔性立管系统的精确模型，具有很好的鲁棒性和自适应性，从而能很好的补偿立管系统的不确定性和外部荷载的不确定性，同时还能很好地跟踪和调节振动状况。
3、本发明针对模型不确定的海洋柔性立管边界振动控制方法，采用自边界主动控制的控制方式，被控驱动装置只需在海洋柔性立管顶端实施控制作用，即可实现抑制立管振动的效果，与现有的被动控制方式相比，安装维护方便，性价比显著提高。
4、本发明针对模型不确定的海洋柔性立管边界振动控制方法，在模糊隶属函数的论域两端采用Z形隶属函数，中间区域采用三角形隶属函数，不仅可满足模糊隶属函数的规则，还充分考虑了模糊化、模糊推理与解模糊过程的精确性与简便性，显著提高了运算效率。
5、本发明针对模型不确定的海洋柔性立管边界振动控制方法，提出了在海洋柔性立管领域内对系统稳定进行验证的有效方法，保证了在采用自适应模糊反演算法时，整个柔性立管系统的稳定性，为后续的控制精确度提供了充分保障。
附图说明
此处所说明的附图用来提供对本发明实施例的进一步理解，构成本申请的一部分，并不构成对本发明实施例的限定。在附图中：
图1为本发明具体实施例中海洋柔性立管的示意图；
图2为本发明具体实施例的流程示意图；
图3为本发明具体实施例中未施加控制作用下海洋柔性立管的三维振动偏移量示意图；
图4为本发明具体实施例中施加控制后海洋柔性立管的三维振动偏移量示意图；
图5为本发明具体实施例中海洋柔性立管500m处振动偏移量幅值对比示意图；
图6为本发明具体实施例中海洋柔性立管1000m处振动偏移量幅值对比示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下面结合实施例和附图，对本发明作进一步的详细说明，本发明的示意性实施方式及其说明仅用于解释本发明，并不作为对本发明的限定。在本申请的描述中，需要理解的是，术语“前”、“后”、“左”、“右”、“上”、“下”、“竖直”、“水平”、“高”、“低”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本申请保护范围的限制。
实施例1：
针对模型不确定的海洋柔性立管边界振动控制方法，如图2所示，包括如下步骤：
对柔性立管建模、得到柔性立管系统的控制方程，确定柔性立管系统的边界条件和初始条件；
基于自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律；再通过李雅普诺夫稳定性理论，验证系统是否稳定有界：若是，进入下一步；若否，重新调整第一控制参数和/或第二控制参数和/或第三控制参数。
获取柔性立管的横向振动实际偏移量；
将所述实际偏移量代入所述自适应模糊反演控制律，得到实际控制输出；
基于所述实际控制输出，在顶端对柔性立管进行控制。
实施例2：
本实施例中海洋柔性立管系统模型的示意图如图1所示，图1中虚线箭头表示分布式海流，坐标系原点位于立管底部，控制器位于立管顶端（边界）做控制输入，以产生横向力来抑制立管振动，且控制方向向右。其中：立管长度为L，立管横向振动偏移量为w(x,t)，作用在立管上的分布式海流载荷为f(x,t)，边界扰动为d(t)，空间变量为x；时间变量为t，在本实施例中表示操作间隔时间段。
本实施例的控制过程如下：
一、建模，对柔性立管系统模型进行动力学分析：
将图1所示柔性立管系统的动能E k表达为：
Figure 716537DEST_PATH_IMAGE024
式中，M s为船舶质量，ρ为单位立管质量，
Figure 421187DEST_PATH_IMAGE025
为w(x,t)中x=L时对t的一阶偏导，
Figure 838393DEST_PATH_IMAGE026
为w(x,t)对t的一阶偏导。
将图1所示柔性立管系统的势能E P表达为：
Figure 337508DEST_PATH_IMAGE027
式中，EI为立管弯曲刚度，T为立管张力，为
Figure 371323DEST_PATH_IMAGE028
w(x,t)对x的二阶偏导，w’(x,t)为w(x,t)对x的一阶偏导。
再将系统所受总虚功W 0表达为：
Figure 223872DEST_PATH_IMAGE029
式中，c为结构阻尼系数，d s为船舶阻尼系数，δ为变分运算符，w(L,t) 为w(x,t)中x=L时的值，
Figure 885798DEST_PATH_IMAGE030
为w(x,t)对t的一阶偏导，
Figure 696759DEST_PATH_IMAGE031
为w(x,t)中x=L时对t的一阶偏导，u(t)为作用在立管顶端的边界控制。
基于哈密顿原理
Figure 686712DEST_PATH_IMAGE032
，得到系统控制方程：
Figure 467586DEST_PATH_IMAGE002
通过
Figure 124963DEST_PATH_IMAGE033
得到系统边界条件为：
Figure 965881DEST_PATH_IMAGE003
通过
Figure 708709DEST_PATH_IMAGE034
得到系统初始条件为：
Figure 168640DEST_PATH_IMAGE035
。
式中，
Figure 539578DEST_PATH_IMAGE005
为w(x,t)对t的一阶偏导，
Figure 426763DEST_PATH_IMAGE006
为w(x,t)对t的二阶偏导；
Figure 515942DEST_PATH_IMAGE007
为w(x,t)对x的二阶偏导，
Figure 717247DEST_PATH_IMAGE008
为w(x,t)对x的四阶偏导；
Figure 942692DEST_PATH_IMAGE009
为w(x,t)中x=L时对x的一阶偏导，
Figure 283935DEST_PATH_IMAGE010
为w(x,t)中x=L时对x的二阶偏导，
Figure 735776DEST_PATH_IMAGE011
为w(x,t)中x=L时对x的三阶偏导，
Figure 927723DEST_PATH_IMAGE012
为w(x,t)中x=L时对t的一阶偏导，
Figure 820724DEST_PATH_IMAGE013
为w(x,t)中x=L时对t的二阶偏导，
Figure 439924DEST_PATH_IMAGE014
为w(x,t)中x=0时对x的一阶偏导，w(x,0)为w(x,t)中t=0时的值，
Figure 113482DEST_PATH_IMAGE015
为w(x,t) 中t=0时对t的一阶偏导。
本实施例还对其给出证明假设条件进行了验证，验证结果表明f(x,t)和d(t)的能量是有限的，因此是有界的，假设合理，验证了本实施例所建立的控制方程的准确性。
二、设计海洋柔性立管系统的自适应模糊反演控制律：
使立管稳定在原位置的小邻域内，设计边界控制律为：
Figure 577961DEST_PATH_IMAGE036
式中，u(t)为作用在立管顶端的边界控制量，λ 2为第一控制参数，e 1为海洋柔性立管边界实际位置与边界理想位置的误差，e 2为海洋柔性立管边界横向振动偏移量变化率与虚拟控制信号之间的误差，x 2为海洋柔性立管边界横向振动偏移量变化率，y为海洋柔性立管边界实际位置，y d为立管理想位置，α 11为虚拟控制信号；
其中φ 2为用于逼近未知非线性函数的自适应模糊系统，表达为：
Figure 918944DEST_PATH_IMAGE037
式中，
Figure 584412DEST_PATH_IMAGE038
为模糊基向量的转置，θ2为逼近向量；
然后设计自适应律：
引入模糊系统来逼近非线性函数，采用单值模糊器、乘积推理机和重心平均反模糊器，确定模糊规则，则能够得到模糊系统的输出；
在逼近过程中，考虑最优逼近向量，对于给定的任意小的常量ε＞0，通过计算可求得本实施例所需的自适应律为：
Figure 604320DEST_PATH_IMAGE019
式中，
Figure 279015DEST_PATH_IMAGE039
为自适应参数，r 2为第二控制参数，k 2为第三控制参数，
Figure 333559DEST_PATH_IMAGE040
为模糊基向量，θ2为逼近向量。
需要说明的是，本实施例中第一控制参数λ 2、第二控制参数r 2、第三控制参数k 2均为大于0的常数。
三、系统稳定性验证：
对于整个系统，首先设计出李雅普诺夫（Lyapunov）函数，然后确定李雅普诺夫函数的上下界；再根据李雅普诺夫函数对时间的导数，得到如下不等式：
Figure 169928DEST_PATH_IMAGE041
该不等式表明本实施例中柔性立管横向振动的范围始终不会超过右侧的值，因此在满足该不等式的条件下，可认为系统稳定且有界。
当系统不满足如上不等式时，认为系统不稳定，此时需要根据实际情况重新调整第一控制参数和/或第二控制参数和/或第三控制参数，直至系统满足稳定有界为止。具体调整方法可采用：
首先调整第一控制参数λ 2，重新验证系统稳定性，若系统无法满足稳定有界，调整第二控制参数r 2并再次验证系统稳定性，若系统依然无法满足稳定有界，再调整第三控制参数k 2。按此调整顺序能够快速高效的得到稳定有界的海洋柔性立管系统，显著降低运算量。
四、获取海洋柔性立管系统的实际振动偏移量：
在t时刻利用激光位移传感器测量信号w(L,t)，再通过有限差分法得到
Figure 818078DEST_PATH_IMAGE042
。
五、将得到的w(L,t)、
Figure 889939DEST_PATH_IMAGE043
代入自适应模糊反演控制律，得到实际控制输出，根据该实际控制输出，驱动执行器向海洋柔性立管施加控制作用力。
其中执行器位于柔性立管顶端。
实施例3：
在实施例2的基础上，本实施例针对模型不确定的海洋柔性立管边界振动控制方法进行实验验证：
系统参数选择：λ 2=7.37×108，r 2=1，k 2=1，初始值θ 2(0)为0。模糊隶属函数μ Fi j (x i )在论域两端采用Z形隶属函数，中间区域采用三角形隶属函数。二维模糊控制器两输入量各设置为7个语言值等级，因此可得49条控制规则。柔性立管系统其余主要参数如表1所示：
表1 柔性立管系统其余主要参数
Figure 674356DEST_PATH_IMAGE044
作用在立管上的分布式海流载荷f(x,t)通过如下公式表达：
Figure 806260DEST_PATH_IMAGE045
式中，ρ s为海水密度，C D为阻力系数，D为立管外径，f v 为涡旋脱落频率，U(x,t)为海洋洋流速度与深度关系，θ为分布式海流的相位角；A D 为阻力的振荡部分幅值、其值取f(x,t)中第一项的20%；
其中：
涡旋脱落频率
Figure 941706DEST_PATH_IMAGE046
，St为斯特劳哈尔数；
海洋洋流速度与深度关系
Figure 817258DEST_PATH_IMAGE047
，U(t)为洋面海流速度，x为深度；
洋面海流速度：
Figure 721760DEST_PATH_IMAGE048
; i=1,2,3,4;
w i =(w 1, w 2, w 3, w 4)=（0.867，1.827，2.946，4.282）。
Figure 899932DEST_PATH_IMAGE049
为内流平均流速，U’=0.2 为内流速度波动振幅。
边界扰动d(t)为：
Figure 647308DEST_PATH_IMAGE050
图3所示为在本实例的系统基础上，未施加控制时海洋柔性立管三维振动偏移量示意图；图4所示为采用本申请控制方法施加控制时，海洋柔性立管振动三维图。
图5所示为在有控制作用下、无控制作用下海洋柔性立管在500m深度处的振动偏移量幅值对比。
图6所示为在有控制作用下、无控制作用下海洋柔性立管在1000m深度（边界）处的振动偏移量幅值对比。
需要说明的是，附图3-6中：
标记w(x,t)/m的坐标轴即为本申请中的横向振动实际偏移值w(x,t)，单位为米；
标记x/m的坐标轴即为本申请中的空间变量x，单位为米；
标记Time/s的坐标轴即为本申请中的时间变量t，单位为秒。
对比图3-6可知，本发明所提出的控制方法，能明显抑制模型不确定情况下海洋柔性立管的振动。与现有控制算法相比，本发明提出的边界主动控制算法结构简单，只有一个自适应参数，计算量小，显著提高了控制速度。
以上所述的具体实施方式，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施方式而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
需要说明的是，在本文中，术语“包括”、“包含”或者其任何其它变体，意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。此外，在本文中使用的术语“连接”在不进行特别说明的情况下，可以是直接相连，也可以是经由其他部件间接相连。

Claims (8)

1.针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，包括：
对柔性立管建模、得到柔性立管系统的控制方程，确定柔性立管系统的边界条件和初始条件；
基于自适应模糊反演算法，建立柔性立管系统的自适应模糊反演控制律；
获取柔性立管的横向振动实际偏移量；
将所述实际偏移量代入所述自适应模糊反演控制律，得到实际控制输出；
基于所述实际控制输出，在顶端对柔性立管进行控制；
其中，所述柔性立管系统的控制方程为：
Figure DEST_PATH_IMAGE001
所述边界条件为：
Figure 546276DEST_PATH_IMAGE002
所述初始条件为：
Figure DEST_PATH_IMAGE003
；
式中：t为时间变量、x为空间变量，L为立管长度，w(x,t)为横向振动实际偏移量，f(x,t)为作用在立管上的分布式海流载荷，d(t)为边界扰动，M s为船舶质量，ρ为单位立管质量，EI为立管弯曲刚度，T为立管张力，c为结构阻尼系数，d s为船舶阻尼系数；u(t)为作用在立管顶端的边界控制量；
Figure 650367DEST_PATH_IMAGE004
为w(x,t)对t的一阶偏导，
Figure 898946DEST_PATH_IMAGE005
为w(x,t)对t的二阶偏导；
Figure 28576DEST_PATH_IMAGE006
为w(x,t)对x的二阶偏导，
Figure 311790DEST_PATH_IMAGE007
为w(x,t)对x的四阶偏导；
Figure 852362DEST_PATH_IMAGE008
为w(x,t)中x=L时对x的一阶偏导，
Figure 588236DEST_PATH_IMAGE009
为w(x,t)中x=L时对x的二阶偏导，
Figure 521557DEST_PATH_IMAGE010
为w(x,t)中x=L时对x的三阶偏导，
Figure 659278DEST_PATH_IMAGE011
为w(x,t)中x=L时对t的一阶偏导，
Figure 121483DEST_PATH_IMAGE012
为w(x,t)中x=L时对t的二阶偏导，
Figure 593921DEST_PATH_IMAGE013
为w(x,t)中x=0时对x的一阶偏导，w(x,0)为w(x,t)中t=0时的值，
Figure 799775DEST_PATH_IMAGE014
为w(x,t) 中t=0时对t的一阶偏导；
其中，建立柔性立管系统的自适应模糊反演控制律的方法包括：
建立自适应模糊逻辑系统φ 2：
Figure 57581DEST_PATH_IMAGE015
；
式中，
Figure 956267DEST_PATH_IMAGE016
为模糊基向量的转置，θ2为逼近向量；
确定边界控制律：
Figure 916001DEST_PATH_IMAGE017
式中，u(t)为作用在立管顶端的边界控制量，λ 2为第一控制参数，e 1为海洋柔性立管边界实际位置与边界理想位置的误差，e 2为海洋柔性立管边界横向振动偏移量变化率与虚拟控制信号之间的误差，x 2为海洋柔性立管边界横向振动偏移量变化率，y为海洋柔性立管边界实际位置，y d为立管理想位置，α 11为虚拟控制信号；
根据所述边界控制律，确定自适应律：
Figure 659966DEST_PATH_IMAGE018
式中，
Figure 37858DEST_PATH_IMAGE019
为自适应参数，r 2为第二控制参数，k 2为第三控制参数，
Figure 841866DEST_PATH_IMAGE020
为模糊基向量，θ2为逼近向量。
2.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，得到所述控制方程的方法包括：
基于动力学分析得到柔性立管系统的动能E k、势能E P、以及非保守力对海洋柔性立管所做的总虚功W 0；
基于哈密顿原理、通过如下公式得到柔性立管系统的控制方程：
Figure 39629DEST_PATH_IMAGE021
；
式中，δ为变分操作符，t 1、t 2分别为两个端点时刻，t为时间变量。
3.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，所述自适应模糊逻辑系统在其模糊隶属函数的论域两端采用Z形隶属函数，在论域中间区域采用三角形隶属函数。
4.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，在建立所述自适应模糊反演控制律后，还通过李雅普诺夫稳定性理论，验证系统是否稳定有界：若是，进入下一步；若否，重新调整第一控制参数和/或第二控制参数和/或第三控制参数。
5.根据权利要求4所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，验证系统是否稳定有界的方法包括：
对于整个柔性立管系统，设计李雅普诺夫函数；
确定李雅普诺夫函数的上下界；
根据李雅普诺夫函数对时间的导数得到如下不等式：若满足该不等式，认为系统稳定有界；若不满足该不等式，认为系统不是稳定有界：
Figure 102132DEST_PATH_IMAGE022
；
式中，t为时间变量、x为空间变量，w(x,t)为立管横向振动偏移量，L为立管长度，T为立管张力，υ、υ 1均为大于0的常量，ε 0为有界干扰。
6.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，所述柔性立管的横向振动实际偏移量通过激光位移传感器实时获取。
7.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，在顶端对柔性立管进行控制的方法为：根据实际控制输出，驱动位于柔性立管顶端的执行器，由所述执行器向柔性立管施加控制作用力。
8.根据权利要求1所述的针对模型不确定的海洋柔性立管边界振动控制方法，其特征在于，作用在立管上的分布式海流载荷f(x,t)通过如下公式表达：
Figure 68951DEST_PATH_IMAGE023
；
式中，ρ s为海水密度，C D为阻力系数，D为立管外径，A D 为阻力的振荡部分幅值，f v 为涡旋脱落频率，U(x,t)为海洋洋流速度与深度关系，θ为分布式海流的相位角。
CN202210371917.7A 2022-04-11 2022-04-11 针对模型不确定的海洋柔性立管边界振动控制方法 Expired - Fee Related CN114442481B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210371917.7A CN114442481B (zh) 2022-04-11 2022-04-11 针对模型不确定的海洋柔性立管边界振动控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210371917.7A CN114442481B (zh) 2022-04-11 2022-04-11 针对模型不确定的海洋柔性立管边界振动控制方法

Publications (2)

Publication Number Publication Date
CN114442481A CN114442481A (zh) 2022-05-06
CN114442481B true CN114442481B (zh) 2022-06-07

Family

ID=81360555

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210371917.7A Expired - Fee Related CN114442481B (zh) 2022-04-11 2022-04-11 针对模型不确定的海洋柔性立管边界振动控制方法

Country Status (1)

Country Link
CN (1) CN114442481B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114578697B (zh) * 2022-05-09 2022-07-12 西南石油大学 一种电机驱动机械手的多约束自适应控制方法
CN114741991B (zh) * 2022-06-13 2022-08-23 西南石油大学 用于海洋立管的仿生鲨鱼皮减阻控制方法
CN114895570B (zh) * 2022-07-13 2022-09-20 西南石油大学 一种海洋柔性立管的多约束自适应控制方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537144A (zh) * 2014-11-27 2015-04-22 电子科技大学 带死区输入的海洋柔性立管边界控制方法
CN108729862A (zh) * 2018-06-26 2018-11-02 广州大学 针对海洋柔性立管系统的鲁棒自适应三维振动抑制方法
CN109828466A (zh) * 2019-02-28 2019-05-31 广州大学 一种针对海洋柔性立管的边界鲁棒自适应控制方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB0207382D0 (en) * 2002-03-28 2002-05-08 Holland Heating Uk Ltd Computer cabinet
CN110647104B (zh) * 2019-09-02 2021-05-11 华南理工大学 一种基于边界扰动观测器的柔性立管反步边界控制方法
WO2021162726A1 (en) * 2020-02-12 2021-08-19 Halliburton Energy Services, Inc. Identifying anomalies in well-environment flexible pipes
CN112966418B (zh) * 2021-03-11 2022-08-26 西南石油大学 一种海洋天然气水合物开采立管动力特性分析方法
CN114036857A (zh) * 2021-11-26 2022-02-11 天津大学 一种基于量化分析和概率推理的柔性立管风险评价方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104537144A (zh) * 2014-11-27 2015-04-22 电子科技大学 带死区输入的海洋柔性立管边界控制方法
CN108729862A (zh) * 2018-06-26 2018-11-02 广州大学 针对海洋柔性立管系统的鲁棒自适应三维振动抑制方法
CN109828466A (zh) * 2019-02-28 2019-05-31 广州大学 一种针对海洋柔性立管的边界鲁棒自适应控制方法及装置

Also Published As

Publication number Publication date
CN114442481A (zh) 2022-05-06

Similar Documents

Publication Publication Date Title
CN114442481B (zh) 针对模型不确定的海洋柔性立管边界振动控制方法
Gao et al. Velocity free leader-follower formation control for autonomous underwater vehicles with line-of-sight range and angle constraints
Fortuna et al. A roll stabilization system for a monohull ship: modeling, identification, and adaptive control
How et al. Control of coupled vessel, crane, cable, and payload dynamics for subsea installation operations
Hu et al. Robust adaptive NN control of dynamically positioned vessels under input constraints
CN111552182B (zh) 一种基于干扰观测器的船舶动力定位系统的抗干扰饱和控制方法
CN113238567B (zh) 一种基于扩展状态观测器的底栖式auv弱抖振积分滑模点镇定控制方法
Xia et al. Adaptive energy-efficient tracking control of a X rudder AUV with actuator dynamics and rolling restriction
Mahapatra et al. Nonlinear matrix inequality approach based heading control for an autonomous underwater vehicle with experimental realization
Quan et al. A geometrically exact formulation for three-dimensional numerical simulation of the umbilical cable in a deep-sea ROV system
Jensen Offshore pipelaying dynamics
Mu et al. State-unknown single parameter learning adaptive output feedback control for ship dynamic positioning
CN115686034B (zh) 考虑速度传感器失效的无人潜航器轨迹跟踪控制方法
Liang et al. Fully coupled time-domain simulation of dynamic positioning semi-submersible platform using dynamic surface control
Dworak et al. Improving quality of regulation of a nonlinear MIMO dynamic plant
CN114895570B (zh) 一种海洋柔性立管的多约束自适应控制方法
Wang et al. Nonsingular fast terminal sliding mode-based robust adaptive structural reliable position-mooring control with uncertainty estimation
Huang et al. Neural network adaptive backstepping control of multi-link underwater flexible manipulators
GB2351373A (en) Dimensioning or adapting an elastic structure subjected to a fluid in motion
Candeloro et al. Analysis of a multi-objective observer for uuvs
CN112904719B (zh) 一种适用于水下机器人位置环形区域跟踪控制方法
Chen et al. Vibration Control of Marine Flexible Riser with Uncertain Model
Liu et al. Dynamic positioning of AUVs in shallow water environment: Observer and controller design
CN113406884B (zh) 一种基于滑模自适应的多点系泊系统定位控制方法
Chen et al. Adaptive neural control for coordinated rolling motion of vessels with zero speed via composite learning

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20220607

CF01 Termination of patent right due to non-payment of annual fee