一种水下连接器密封性能退化随机模型建模方法与流程

本发明属于机械工程技术领域，尤其涉及一种水下连接器密封结构性能退化模型的建模方法。

背景技术：

水下连接器是海洋油气水下生产系统的关键设备，主要用于水下生产管道之间的连接，密封能力是其最重要的性能指标。由于连接器服役期间长期受到生产流体与海水的腐蚀，且还受到洋流等随机载荷的作用，连接器性能退化现象客观存在，需要建立连接器性能退化模型，获得服役期间密封性能指标的变化规律，便于后续开展连接器的可靠性分析工作。

目前，通常以水下连接器的密封面接触应力值为水下连接器性能指标，当接触应力值大于水下连接器内工质压力时即认为密封可靠。因此，在建立机械密封结构的退化模型时，通常以密封接触应力值为表征量，采集密封结构服役期间密封面接触应力值后，应用数据统计方法建立密封结构性能退化模型(王少萍,洪葳,梁贤赓等.一种基于接触应力的橡胶密封件性能退化监测方法[p].中国,cn201410074590.2,2016-05-04.)。而实际上，密封性能不但受密封面接触应力大小的影响，而且还受接触应力在整个有效接触长度上的分布情况影响(brownmw.sealsandsealinghandbook[m].elseviersciencepublishers,2007.)，仅仅以接触应力值的大小作为密封性能表征量，一定程度上影响密封性能分析的准确性，导致建立的密封结构性能退化模型不准确。

技术实现要素：

针对上述现有技术中存在的问题，本发明的目的是通过考虑水下连接器密封面有效接触长度对密封性能的影响，提供一种水下连接器密封结构性能退化模型的建模方法，以提高建立的水下连接器密封结构性能退化随机模型的准确性。

为了达到上述目的，本发明采用的技术方案是，一种水下连接器密封结构性能退化模型的建模方法，包括以下步骤：

1)根据水下连接器的原始结构尺寸和水下连接器的实际工况条件，通过有限元软件建立水下连接器密封性能热-固有限元分析参数化模型；

2)为考虑水下连接器的密封面有效接触长度对其密封性能的影响，以密封接触强度为密封性能表征量，使用步骤1)中建立的水下连接器密封性能热-固有限元分析参数化模型，对水下连接器密封性能进行分析，获得当前时刻水下连接器密封性能量化指标fsi，i为当前循环次数；首次循环获得初始t1时刻水下连接器密封性能量化指标fs1；

3)根据水下连接器的实际工况条件，采用腐蚀速率预测模型对水下连接器上毂座、密封件和下毂座的腐蚀速率进行预测，得到实际工况条件下的腐蚀速率，即可获得ti+1时刻对应的上毂座、密封件和下毂座腐蚀量，进而得到ti+1时刻水下连接器被腐蚀后上毂座、密封件和下毂座的结构尺寸；

4)使用步骤3)中ti+1时刻水下连接器被腐蚀后获得的上毂座、密封件和下毂座的结构尺寸，对步骤1)中建立的水下连接器密封性能热-固有限元分析参数化模型进行结构尺寸参数更新；

5)重复步骤2)至步骤4)3次；分别获得t2、t3和t4时刻水下连接器密封性能量化指标；fs2、fs3和fs4；

6)考虑水下连接器的密封性能退化过程具有随机性和不可逆性，选择gamma随机过程描述水下连接器的密封性能退化机理，使用步骤2)和步骤5)中获得的不同时刻水下连接器密封性能量化指标，即t1时刻对应的fs1、t2时刻对应的fs2、t3时刻对应的fs3、t4时刻对应的fs4，对gamma随机过程特征参数进行估计，建立基于gamma随机过程的水下连接器密封性能退化模型。

进一步的，所述水下连接器的实际工况条件包括压力、温度、液相与气相流速、含水率和管道内径。

进一步的，所述步骤2)中，所述密封接触强度的计算公式如下：

其中，fs为密封接触强度；σs(l)为密封面接触应力关于各单元节点对应有效长度l的函数；lss为密封面有效接触长度；n为相关性指数，当存在密封剂时，n＝1.2；当不存在密封剂时，n＝1.4；对于水下连接器，其工作介质温度较高，密封结构中密封剂对密封性能的影响较小，因此取相关性指数n＝1.4。

进一步的，所述步骤3)中，所述腐蚀速率预测模型包括经验型腐蚀速率模型、半经验型腐蚀速率模型和机制型腐蚀速率模型，可根据实际情况进行选择。

进一步的，所述步骤6)中，所述gamma随机过程的定义如下：

可称随机过程{d(t),t≥0}服从形状参数为v(t)、尺度参数为u的gamma过程，若其满足如下条件：

③d(t)＝0；

④d(t)在非相交区域中具有独立增量d(t)-d(s)(0≤s＜t)；

③当0≤s＜t时，独立增量d(t)-d(s)服从gamma分布ga(v(t)-v(s),u)；

gamma分布的概率密度函数为：

其中，为gamma函数；ia(d)为示性函数，当d∈a时，ia(d)＝1；当时，ia(d)＝0。

进一步的，所述步骤1)中，所述的有限元软件为通用有限元分析软件，包括abaqus或ansys。

与现有技术相比，本发明至少具有以下有益效果，本发明提供的一种水下连接器密封结构性能退化模型的建模方法，使用密封接触强度作为密封性能表征量，可以同时考虑水下连接器的密封面接触应力和密封面有效接触长度对其性能的影响；使用gamma随机过程描述水下连接器的密封性能退化机理，可以描述水下连接器的密封性能退化时表现的随机性与不可逆性。由此可见，应用本发明方法建立的水下连接器密封结构性能退化随机模型，其准确性得到了提高。

附图说明

图1本发明水下连接器密封性能退化随机模型建模流程图；

图2本发明实施例中采油树井口连接器基本结构；

图3本发明实施例中连接器密封性能热-固有限元分析参数化模型；

图4本发明实施例中水下连接器密封性能退化随机模型。

图中：1-上毂座；2-密封件；3-下毂座。

具体实施方式

下面结合附图和具体实施方式，对本发明作进一步的说明：

如图2所示，为设计压力为34.5mpa(5000psi)的与水下井口20in套管头配套的水下采油树井口连接器基本结构，包括上毂座1、密封件2和下毂座3。

如图1所示，本发明的水下连接器的密封性能退化模型建模方法包括以下步骤：

s1，根据水下连接器的原始结构尺寸和实际工况条件，实际工况条件包括压力、温度、液相与气相流速、含水率和管道内径；使用abaqus或ansys有限元软件建立连接器密封性能热-固有限元分析参数化模型，如图3所示；

s2，为考虑水下连接器的密封面有效接触长度对密封性能的影响，以密封接触强度为密封性能表征量，使用步骤s1中建立的水下连接器密封性能热-固有限元分析参数化模型，对水下连接器密封性能进行分析，得到密封面各单元节点的接触应力σs1，σs2，…，σsi和各单元节点对应的密封面有效接触长度l1，l2，…，li，以此作为样本数据，使用插值方法建立密封面接触应力关于有效长度的函数σs(l)，代入密封接触强度公式进行计算，即获得连接器初始密封性能指标值fs1，如式(1)所示：

其中，fs为密封接触强度；σs(l)为密封面接触应力关于各单元节点对应有效长度l的函数；lss为密封面有效接触长度；n为相关性指数，当存在密封剂时，n＝1.2；当不存在密封剂时，n＝1.4；对于水下连接器，其工作介质温度较高，密封结构中密封剂对密封性能的影响较小，因此取相关性指数n＝1.4。

s3，根据水下连接器的实际工况条件，包括压力、温度、液相与气相流速、含水率和管道内径，采用如式(2)所示norsok腐蚀速率预测模型(经验型)对水下连接器上毂座1、密封件2和下毂座3的腐蚀速率进行预测，得到实际工况条件下密封结构的腐蚀速率为vcorrmm/a，即可获得任意时刻t1上毂座1、密封件2和下毂座3的腐蚀量vcorr×t1，从而得到腐蚀后各结构的具体尺寸；

其中，vcorr为腐蚀速率，mm/a；kt为与温度和腐蚀产物膜相关的常数；为co2逸度，bar；s为与流速相关的管道内壁切应力，pa；为管道内生产流体ph值对腐蚀速率的影响因子；

s4，使用步骤s3获得的腐蚀后上毂座1、密封件2和下毂座3的结构尺寸，对步骤s1中建立的水下连接器密封性能热-固有限元分析参数化模型进行结构尺寸参数更新；

s5，重复步骤s2至步骤s4共计3次，分别获得t2、t3和t4时刻水下连接器密封性能量化指标fs2、fs3和fs4；

s6，考虑连接器的密封性能退化过程具有随机性和不可逆性，选择gamma随机过程描述密封性能退化机理，所述gamma随机过程的定义如下：

可称随机过程{d(t),t≥0}服从形状参数为v(t)、尺度参数为u的gamma过程，若其满足如下条件：

⑤d(t)＝0；

⑥d(t)在非相交区域中具有独立增量d(t)-d(s)(0≤s＜t)；

③当0≤s＜t时，独立增量d(t)-d(s)服从gamma分布ga(v(t)-v(s),u)；

gamma分布的概率密度函数为：

其中，为gamma函数；ia(d)为示性函数，当d∈a时，ia(d)＝1；当时，ia(d)＝0。

然后，使用步骤s2和s5中获得的不同时刻连接器密封性能指标fs1，fs2，fs3，fs4对gamma随机过程特征参数(形状参数为v(t)、尺度参数为u)进行估计，得到基于gamma随机过程的连接器密封性能退化模型，如图4所示。由图4可知，随着时间增加，水下连接器的退化量呈现出逐步递增的趋势，而且在任意时刻密封性能退化量并非一个确定的值，而是具有明显的随机性，表明利用本发明方法构建的退化模型能够描述水下连接器密封性能退化过程的随机性和不可逆性。