含蜡原油触变性测试方法

贾邦龙, 张劲军

贾邦龙, 张劲军. 含蜡原油触变性测试方法[J]. 油气储运, 2012, 31(4): 254-259. DOI: 10.6047/j.issn.1000-8241.2012.04.004
引用本文: 贾邦龙, 张劲军. 含蜡原油触变性测试方法[J]. 油气储运, 2012, 31(4): 254-259. DOI: 10.6047/j.issn.1000-8241.2012.04.004
Jia Banglong, Zhang Jingjun. Thixotropy test method of waxy crude oil[J]. Oil & Gas Storage and Transportation, 2012, 31(4): 254-259. DOI: 10.6047/j.issn.1000-8241.2012.04.004
Citation: Jia Banglong, Zhang Jingjun. Thixotropy test method of waxy crude oil[J]. Oil & Gas Storage and Transportation, 2012, 31(4): 254-259. DOI: 10.6047/j.issn.1000-8241.2012.04.004

含蜡原油触变性测试方法

基金项目: 

北京市优秀博士学位论文指导教师科技项目 YB20081141401

国家自然科学基金重点项目“深水环境下易凝高粘原油-天然气输送系统流动保障基础问题研究” 51134006

详细信息
    作者简介:

    贾邦龙, 在读博士生, 1984年生, 2006年毕业于中国石油大学(北京)油气储运工程专业, 现主要从事油气储运工程相关技术的研究工作。电话: 13810560726;Email: jiabanglong_112@126.com

  • 中图分类号: TE89

Thixotropy test method of waxy crude oil

More Information
    Author Bio:

    Jia Banglong: Ph.D candidate, born in 1984, graduated from China University of Petroleum (Beijing), oil & gas storage and transportation engineering, in 2006, mainly engaged in the research of related technology of oil & gas storage and transportation engineering. Add: No.18, Fuxue Road, Changping District, China University of Petroleum (Beijing), Beijing, 102249, P. R. China. Tel: 13810560726; Email: jiabanglong_112@126.com

  • 摘要: 现有触变性测试方法(加载方式)包括恒剪切率、剪切率阶跃、控制剪切率滞回环、剪应力阶跃以及控制剪应力滞回环等。利用这些方法对同一条件下含蜡原油的触变性进行测试,并利用Houska触变模型和双结构参数触变模型进行拟合。利用一种方法测试数据拟合出的触变参数对其他测试方法加载条件下的应力进行预测,依据预测结果与实测值的吻合程度,对测试方法进行评价。结果表明:无论采用何种测试方法,利用剪切率高的加载条件得到的实验数据拟合的触变参数,均可成功地对较低剪切率加载条件下的应力进行预测,反之则会出现结构参数为负值进而导致预测失败的情况;利用同一测试方法的多组实验数据共同拟合触变参数,可以提高对其他加载方式下应力的预测值与实测值的吻合程度。
    Abstract: Existing thixotropy test methods (loading methods) include constant shear rate, shear rate stepping, control hysteresis loops of shear rates, shear stress stepping, and control hysteresis loops of shear stresses, etc.. These methods are used to test the thixotropy of waxy crude oil under the same conditions, and Houska thixotropy model and double-structural parameter thixotropy models are used for fitting. Thixotropy parameters from fitting of test data by one method are used to predict the stress of other test methods under loading conditions. Based on the coincidence degree between prediction results and measured values, the test method is evaluated. Results show that no matter which test method is used, the thixotropy parameters from fitting of test data obtained by the use of loading conditions with high shear rates can all successfully predict stress under loading conditions with lower shear rates, otherwise the structural parameter is negative so as to cause prediction failure; the thixotropy parameters from common fitting of multiple group of test data by the same testing method can improve the coincidence degree between predicted values and measured values of stress under other loading methods.
  • CNG地下储气井是一种将压缩天然气通过特殊工艺储存在地下50~200 m处的设施[1-2]。地下储气井具有占地面积小、运行费用低、操作维护简便以及事故影响范围小等特点,且具有恒温、抗静电等作用,现已成为国内加气站首选储气系统并大力推行。但储气井结构简单,井内未设置安全阀,地面上无检测系统,随着其使用年限的增长,由于构件的腐蚀和磨损,发生介质泄漏的可能性将不断上升。一旦储气井失效造成井下可燃气体泄漏,由于无法直接观测井底,很可能酿成重大安全事故,如四川某加气站储气井发生爆炸后套管冲出地面的恶性事故[3-4]

    净化后的天然气仍含有H2S、CO2、溶解氧、SRBC硫酸盐还原菌等介质,会造成管壁损伤[5-7]。储气井的损伤主要有套管腐蚀、气体腐蚀导致井筒均匀减薄、井筒局部腐蚀等,此外,裂纹也是影响储气井剩余强度的重要因素。由于韧性钢结构的应力集中不敏感性,相比裂纹形成过程,裂纹扩展阶段占总寿命的90%。考虑裂纹的存在,是分析储气井剩余强度更为合理、经济、合乎使用原则的方法[8-10]

    传统的输气管道剩余强度分析方法多从工作应力角度,分析服役储气井的最大Mises等效应力,或第一主应力[11-14]。祝效华等[15]就套管失效为顶事件的安全问题,通过故障树分析法(FTA),建立了套管失效模糊综合评判分析模型;刘颖等[16]对比了含腐蚀缺陷管道的剩余强度评价方法,指出采用有限元法进行数值模拟更加精确和有效;骆正山等[17]基于失效评定图(FAD)方法,使用标准给定的公式和失效评估曲线对含轴向裂纹管道进行了失效分析。考虑含裂纹的储气井,何家胜等[18]采用1/4节点法和二维奇异单元建立含轴向裂纹管道的有限元分析模型,计算了含不同深度裂纹的管道各节点应变和裂纹前沿最深处应力强度因子;姚安林等[19]针对高压输气管道副裂纹对主裂纹应力强度因子的影响,采用三维虚拟裂纹闭合技术(3D-VCCT),计算了裂纹表面点和最深点应力强度因子;刘卫国等[20]采用相同方法计算了在落石不同冲击速度下的裂纹前沿各点应力强度因子。

    在此,建立储气井井管-水泥环-地层组合模型,将三维虚拟裂纹闭合技术运用于CNG地下储气井裂纹分析中,分别建立了轴向和环向裂纹三维有限元模型。针对裂纹的所处位置和不同尺寸,通过CAE数值仿真模拟,研究裂纹对储气井井管剩余强度的影响,精确得出不同裂纹条件下最大工作应力和裂纹前沿各节点应力强度因子。

    天然气储气井结构简化为如下力学模型(图 1),储气井主要承受内压与地应力,地下100 m处地应力为4 MPa,从井口到地下,呈线性变化[21]。该井主要几何参数为:井深H=100 m,储气井内径D1=177.8 mm,井筒壁厚t1=10.36 mm,水泥环厚度t2=15 mm,地层厚度t3=200 mm,内压p=25 MPa。储气井结构各部分的材料参数如下(表 1)。

    图  1  储气井井管-水泥环-地层力学模型图
    表  1  储气井结构各部分的材料参数
    下载: 导出CSV 
    | 显示表格

    考虑到计算量和网格划分的问题,建立环向全尺寸、纵向截取高度h0=0.1 m研究段的局部储气井模型。并将h0研究段所受地应力简化为均布应力pj

    (1)

    式中:h为研究段纵向中点处距井口的高度,m,

    就载荷而言,由于CNG储气井地下深度一般仅100~200 m,相比地层压力,内压对储气井强度的影响更大,是储气井强度理论重点考虑的载荷因素。在计算裂纹尺寸对剩余强度的影响时,为了简化计算,将地应力忽略。

    管道内壁裂纹分别沿管道轴向分布和环向分布(图 2),裂纹长度为椭圆长半轴r1的2倍,裂纹深度为椭圆短半轴r2的2倍。由于环向裂纹需考虑轴向应力作用,利用薄壁容器轴向应力pz计算式(2),将内压p转换为开口容器横截面的轴向拉应力pz

    (2)
    图  2  储气井管道内壁裂纹模型示意图

    轴向和环向裂纹在长20 mm、深5 mm、内压25 MPa条件下,已知裂纹最大等效应力与距井口距离之间的关系(图 3)。对于轴向裂纹而言,随裂纹距井口距离的增加,最大等效应力逐渐增大。此时,由于地应力存在平行于裂纹面的分量,有弯曲作用效果,故地应力对固井是不利的,即处于井底的轴向裂纹更危险。然而对于环向裂纹,随裂纹距井口距离的增加,最大等效应力逐渐减小,与无裂纹存在的结构具有相同规律。这是因为地应力对井筒有径向的约束作用,对固井是有利的,即处于井口的环向裂纹更危险。

    图  3  含裂纹储气井最大等效应力随裂纹距井口距离的变化曲线

    分别考虑轴向和环向裂纹最危险的情况:①轴向裂纹处于井底位置,研究段内压25 MPa,地应力4 MPa;②环向裂纹处于井口位置,研究段内压25 MPa,地应力0。模拟得到裂纹长20 mm、深5 mm时,裂纹附近的Mises等效应力云图,最大等效应力位于裂纹长轴处,分别为490 MPa、368 MPa。同等尺寸的裂纹,位于井底的轴向裂纹比位于井口的环向裂纹更危险。CNG地下储气井井筒材料为N80石油套管,屈服强度552~758 MPa,安全系数取1.25,则许用应力为441~606 MPa。以许用应力441 MPa作为强度许可条件,长20 mm、深5 mm的轴向裂纹是不被允许的,环向裂纹相对安全。

    分析裂纹深度一定时最大等效应力随裂纹长度的变化。得到裂纹深度为5 mm,长度分别为10 mm、15 mm、20 mm、25 mm、30 mm时的含裂纹储气井最大等效应力(图 4)。以许用应力441 MPa作为强度许可条件,深5 mm、长30 mm以内的轴向裂纹不被允许,环向裂纹相对安全。需要注意的是,在以上工况下,等效应力最大值总出现在裂纹长轴处。轴向裂纹越长,最大等效应力越小;环向裂纹越长,最大等效应力越大。

    图  4  含裂纹储气井最大等效应力随裂纹长度的变化曲线

    分析裂纹长度一定时最大等效应力随裂纹深度变化的情况。得到裂纹长度为30 mm时,裂纹深度分别为1 mm、3 mm、5 mm、7 mm、8 mm、8.5 mm、8.75 mm、9 mm时含裂纹储气井的最大等效应力(图 5)。对于轴向裂纹,随着裂纹深度增加,最大等效应力先减少,当裂纹深度增加至一定程度后,最大等效应力将增大。对于环向裂纹,随着裂纹深度增加,最大等效应力增大。以许用应力441 MPa作为强度许可条件,长30 mm的裂纹,不允许深度大于7 mm的环向裂纹存在;不允许深度大于5 mm的轴向裂纹存在。

    图  5  含裂纹储气井最大等效应力随裂纹深度的变化曲线

    环向裂纹对裂纹长度与深度都很敏感,尺寸越大的环向裂纹越危险。深度5 mm的轴向裂纹,最大等效应力随长度增加呈下降趋势,即裂纹在长度方向的扩展,可以释放裂纹前沿附近的应力场;当裂纹长度为30 mm,深度发生变化时,很难以深度单参数变化描述最大等效应力的变化规律,表现出最大等效应力对裂纹形状(深长比)的依赖性。以下从断裂力学的角度,详细讨论裂纹尺寸、裂纹形状对储气井剩余强度的影响。

    将裂纹前沿轴向位移进行约束使其无法扩展,求得各点应力即内外压的反作用力之和;释放裂纹前沿轴向位移,裂纹因内外压作用而扩展,求得位移即扩展深度。虚拟加深后的裂纹在应力作用下闭合形成扩展前裂纹,可得内外压对裂纹自由面做功之和,由此求得产生单位面积自由表面的能量释放率[22-23]

    采用三维虚拟裂纹闭合法计算裂纹前沿应力强度因子,划分单元和节点(图 6,其中n1~n5为裂纹前沿,i为节点编号,bni为单元长度,Δ为单元宽度,δ为节点由于开裂产生的相对位移,可近似将扩展裂纹面等效为矩形)。

    图  6  基于三维虚拟裂纹闭合法的应力强度因子计算简图

    I型裂纹扩展δ时末端节点n1的单位面积自由表面能量释放率GI/n1为:

    (3)

    式中:Fyn为裂纹前沿上节点力,N。

    中间节点n2的单位面积自由表面能量释放率为:

    (4)

    其中:

    (5)

    对于平面各向同性线性弹性材料,平面应变能量释放率GI和应力强度因子KI关系如下:

    (6)

    由上式可得应力强度因子的表达式:

    (7)

    式中:E为弹性模量,MPa;ν为泊松比[20, 24-25]

    考虑轴向和环向裂纹最危险的情况:①环向裂纹处于井口位置,研究段内压25 MPa,地应力为0;②轴向裂纹处于井底位置,研究段内压25 MPa,地应力4 MPa。环向裂纹模型中的轴向拉应力pz由式(2)计算得107.26 MPa。当裂纹长度为20 mm,深度分别为1 mm、2 mm、4 mm、5 mm时,计算得到轴向和环向裂纹前沿应力强度因子(图 7)。

    图  7  裂纹前沿应力强度因子随各节点位置的变化曲线

    裂纹前沿应力强度因子在裂纹最深处(图 7a横坐标90°位置)出现峰值。当裂纹深长比较小时,该峰值即裂纹前沿应力强度因子的最大值。随裂纹深长比增大,该峰值应力强度因子先增后减,同时,出现峰值附近的应力强度因子变化趋势减缓。裂纹长轴处(图 7a横坐标为0°和180°位置)的应力强度因子不断增大。当深长比增大至一定程度后,裂纹长轴处应力强度因子突破裂纹最深处的峰值,成为裂纹前沿应力强度因子的最大值。由此认为,裂纹深长比较大时,在长度方向更易扩展;裂纹深长比较小时,在深度方向更易扩展,这与基于工作应力分析储气井剩余强度时,判定的危险点不同。利用断裂力学应力强度因子,可以更准确、合理地分析储气井的剩余强度。

    类似基于工作应力的强度分析方法,在获得储气井井管材料N80的断裂韧性KIC的条件下,可利用基于断裂力学的方法分析储气井的剩余强度,判定许可的裂纹尺寸及形状等。

    (1)裂纹的存在将削减储气井的剩余强度,轴向和环向裂纹的工作应力随距井口距离的增大呈线性变化。在同等裂纹尺寸条件下,位于井底的轴向裂纹,储气井工作应力更高,剩余强度更低;距离井口越近的环向裂纹,储气井工作应力越高,剩余强度越低。考虑N80许用应力441 MPa,长20 mm、深5 mm的轴向裂纹不允许出现,而环向裂纹相对安全。

    (2)深度为5 mm的裂纹,长度在10~30 mm变化时,轴向裂纹随裂纹长度增加,最大等效应力线性减小;环向裂纹随裂纹长度增加,最大等效应力线性增大。长度为30 mm的裂纹,深度在1~9 mm变化时,轴向裂纹随裂纹深度增加,最大等效应力先减后增;环向裂纹随裂纹深度增加,最大等效应力逐渐增加。考虑N80许用应力441 MPa,不允许深5 mm的轴向裂纹长度小于30 mm,不允许长30 mm的轴向裂纹深度大于5 mm;不允许深5 mm的环向裂纹长度大于30 mm,不允许长30 mm的环向裂纹深度大于7 mm。

    (3)轴向裂纹和环向裂纹在裂纹深长比较小的情况下,裂纹最深处应力强度因子最大,为危险点,即此处裂纹更容易发生深度方向的扩展,与刘卫国等的研究结论类似;当裂纹深长比增大至一定程度后,裂纹长轴处应力强度因子最大,为危险点,裂纹更容易在长度方向发生扩展。通过对轴向和环向裂纹在最危险情况下的深长比的讨论,得出更系统、完善的含裂纹储气井剩余强度变化规律。

  • 图  1   Houska模型和双结构参数触变模型对不同方法实验数据的拟合结果

    图  2   利用剪切率滞回环和剪应力阶跃方法所得触变参数分别对其他加载方式下的应力预测结果

    图  3   剪应力在20~50 Pa阶跃变化时的实验数据与采用Houska模型拟合的结构参数

    图  4   结构参数的预测值和拟合值

    图  5   同一测试方法单组和多组实验数据拟合的触变参数对18组触变测试加载方式下的应力预测偏差对比

    表  1   触变性测试方法

    下载: 导出CSV

    表  2   不同测试方法所得触变参数对所有试验加载方式下的应力预测结果

    下载: 导出CSV

    表  3   各组剪切率阶跃实验数据拟合的两种触变模型参数对所有触变测试加载方式下的应力预测结果

    下载: 导出CSV

    表  4   多组实验数据共同拟合的触变参数对全部触变测试方法加载方式下的应力预测结果

    下载: 导出CSV
  • [1]

    Zhang Fan, Zhang Jinjun, Yang Xiaojing. Comparison of thixotropic models of waxy crude oil: Proceedings of the XIV International Congress on Rheology, Seoul, Korea, 2004[C]. Seoul: Korea Society of Rheology, 2004.

    [2] 柳建军, 程华, 崔秀国, 等. 应用非线性隐式迭代法回归含蜡原油触变模型[J]. 油气储运, 2011, 30(1): 12-14. https://www.cnki.com.cn/Article/CJFDTOTAL-YQCY201101004.htm
    [3] 刘刚, 张国忠, 邓艳. 胶凝原油的粘弹—触变方程[J]. 中国石油大学学报(自然科学版), 2003, 27(6): 72-75. https://www.cnki.com.cn/Article/CJFDTOTAL-SYDX200306018.htm
    [4]

    Jinjun Zhang, Liping Guo, Houxing Teng. Evaluation of thixotropic models for waxy crude oils based on shear stress decay at constant shear rates[J]. Appl Rheol, 2010, 20(5): 53944.

    [5] 张劲军, 国丽萍. 基于滞回环的含蜡原油触变模型评价[J]. 石油学报, 2010, 31(3): 494-500. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB201003028.htm
    [6]

    Mewis J, Wagner N J. Thixotropy[J]. Advances in Colloid and Interface Science, 2009, 147-148: 214-227. doi: 10.1016/j.cis.2008.09.005

    [7]

    Barnes H A. Thixotropy-a review[J]. Journal of Non-Newtonian Fluid Mechanics, 1997, 70(1): 1-33.

    [8]

    Dullaert K, Mewis J. Thixotropy: build-up and breakdown curves during flow[J]. Journal of Rheology, 2005, 49(6): 1 213-1 230. doi: 10.1122/1.2039868

    [9]

    Dullaert K, Mewis J. A structural kinetics model for thixotropy[J]. Journal of Non-Newtonian Fluid Mechanics, 2006, 139(6): 21-30.

    [10] 马庆坤, 朱维耀, 高珉, 等. 可动凝胶体系渗流流变特性及其表征[J]. 石油学报, 2007, 28(5): 85-88. https://www.cnki.com.cn/Article/CJFDTOTAL-SYXB200705016.htm
    [11]

    Rønningsen H P. Rheological behaviour of gelled waxy North Sea crude oils[J]. J Petrol Sci Eng, 1992, 7(3-4): 177-213.

    [12]

    Houska M. Engineering aspects of the rheology of thixotrpic liquids[D]. Prague: Czech Technical University of Prague, 1981.

    [13] 赵晓东. 改性原油输送管道启动过程非稳态水力热力计算方法研究[D]. 北京: 中国石油大学, 1999.
    [14]

    Cheng D C. Yield stress: a time dependent property and how to measure it[J]. Rheological Acta, 1986, 25(5): 542-554.

    [15]

    Cawkwell M G, Charles M E. Characterization of Canadian Arctic thixotropic gelled crude oils utilizing an eight parameter model[J]. Journal of Pipelines, 1989, 7: 251-264.

图(5)  /  表(4)
计量
  • 文章访问数: 
  • HTML全文浏览量:  0
  • PDF下载量: 
  • 被引次数: 0
出版历程
  • 收稿日期:  2011-05-15
  • 网络出版日期:  2023-08-20
  • 刊出日期:  2012-04-24

目录

/

返回文章
返回