![FX_GRP_ID800006F7](/web/fileinfo/upload/magazine/12781/321271/image002.jpg)
DOI: 10.11817/j.issn.1672-7207.2020.05.018
基于广义Kelvin模型的非定常盐岩蠕变模型
韩伟民1,闫怡飞2,闫相祯1
(1. 中国石油大学(华东) 储运与建筑工程学院,山东 青岛,266580;
2. 中国石油大学(华东) 机电工程学院,山东 青岛,266580)
摘要:针对川东北油气田盐岩地层普遍出现的套损现象,对该地层段盐岩的蠕变特性进行研究。用非定常黏壶元件替换广义Kelvin模型中的线性黏壶元件,建立非定常广义Kelvin模型;将其与能够描述盐岩稳态蠕变特征的Heard模型进行串联,构建新的四元件NGKH盐岩非线性蠕变本构模型,可以描述盐岩的衰减蠕变阶段和稳态蠕变阶段特征。基于室内三轴蠕变试验结果,利用Levenberg-Marquardt算法对NGKH模型蠕变参数进行识别。基于FLAC3D的UDM接口程序,采用C++编程语言对其进行二次开发;并通过三轴蠕变试验的数值模拟对开发NGKH模型进行验证。研究结果表明:NGKH模型的拟合结果与室内试验结果吻合良好,该蠕变模型能够准确描述盐岩的非线性蠕变特性;通过黏弹性与非线性计算结果,验证了开发NGKH模型程序的正确性和合理性。
关键词:蠕变本构模型;非定常广义Kelvin模型;Heard稳态蠕变模型;二次开发;有限差分法
中图分类号:TU45 文献标志码:A 开放科学(资源服务)标识码(OSID)
文章编号:1672-7207(2020)05-1337-13
Non-stationary creep model for salt rock based on generalized Kelvin model
HAN Weimin1, YAN Yifei2, YAN Xiangzhen1
(1. College of Pipeline and Civil Engineering, China University of Petroleum, Qingdao 266580, China;
2. College of Mechanical and Electronic Engineering, China University of Petroleum, Qingdao 266580, China)
Abstract: The creep characteristics of salt rock were studied in northeastern Sichuan oil and gas field in view of the casing damage commonly occurring in salt rock formation. A non-stationary generalized Kelvin model was established by replacing the linear viscous component with a nonlinear viscous component. By connecting it with the Heard model, which could describe the steady creep stage of salt rock, a new non-linear four-element NGKH creep constitutive model for salt rock was constructed, which could describe the characteristics of decay creep stage and steady creep stage of salt rock. Based on the results of triaxial creep tests in laboratory, the creep parameters of NGKH model were identified by Levenberg-Marquardt algorithm. The NGKH model was redeveloped using UDM interface program of FLAC3D with C++ programming language, and verified by numerical simulation of triaxial compression tests. The results show that the fitting results of NGKH model are in good agreement with the laboratory test results, and the creep model can accurately describe the non-linear creep characteristics of salt rock. The correctness and rationality of developing NGKH model can be verified by the viscoelastic and non-linear calculation results.
Key words: creep constitutive model; non-stationary generalized Kelvin model; Heard steady creep model; secondary development; finite difference method
四川盆地作为油气资源的重要赋存地,分布着大量厚度不一的盐岩地层。而位于川东北三叠系嘉陵江组盐岩地层的多口油气田生产井出现了套管变形损坏情况[1],因此,研究该地层段盐岩的流变特性就显得尤为重要。当应力水平低于岩石长期强度时,其蠕变曲线仅存在衰减蠕变阶段和稳态蠕变阶段;当应力水平高于岩石长期强度时,其蠕变曲线呈现完整的3个阶段,即衰减蠕变阶段、稳态蠕变阶段和加速蠕变阶段。室内试验结果表明盐岩的蠕变曲线也存在上述3个阶段[2],但对于工程问题而言,处于高应力状态下的深层盐岩体蠕变过程中并不会进入加速蠕变阶段而发生损伤破坏,因此,研究衰减蠕变阶段和稳态蠕变阶段盐岩的非线性蠕变特征对工程实际更有意义。其中,衰减蠕变阶段的存在时间相对短暂,在恒定载荷作用下稳态蠕变阶段则更为持久。国内外学者针对盐岩的蠕变特征开展了大量研究[3-6],但由于地质环境、矿物成分及含水率等影响因素的存在,不同地区的盐岩弹性参数及流变特性存在明显差异,因此,单一的蠕变本构模型形式并不具有普适性。目前针对盐岩的非线性流变问题,主要有以下几种常见研究方法[7]:1) 以组合的黏弹性或黏塑性元件模型为基础,在其后串联加入非线性元件,蠕变参数可根据室内试验结果来确定。刘开云等[8]将Maxwell模型与非线性黏塑性体串联组成新的非线性黏塑性模型以描述泥岩的蠕变过程。一般的工程问题可采用该方法进行近似处理,但对于复杂的岩石蠕变本构,组合的黏弹性或黏塑性元件模型并不能反映与蠕变时间和应力水平相关的非线性蠕变特征,此时,该方法的近似处理结果并不适用。2) 将元件模型中的线性元件视为时间的函数,并以非线性元件进行替代,改进后的非线性元件组合模型的蠕变参数可根据室内试验结果进行确定。王军保等[9-10]用非线性黏壶代替常规Burgers模型的线性黏壶,通过考虑损伤因子和蠕变对黏滞系数的影响等,建立了改进的非线性元件组合模型以描述盐岩的非线性蠕变特性。ZHOU等[11-12]用分数阶Abel黏壶元件替代西原模型的线性黏壶,建立了改进的西原模型以反映盐岩流变的3个阶段。该方法虽较为理想,但模型中蠕变参数辨识的计算处理相对繁琐和困难,且改进后的非线性元件模型并不能反映温度因素对岩石蠕变过程的影响。3) 基于室内试验结果,考虑损伤和断裂力学,建立经验本构关系式。高小平等[13]应用不同应力和温度条件下的盐岩变形机制,得到了盐岩的稳态蠕变应变率与偏应力、能量与温度之间的函数关系式。杨春和等[14-15]通过不同应力水平作用下盐岩蠕变变形与时间的关系,建立了盐岩瞬态与稳态蠕变的耦合本构方程。王安明等[16-17]采用Norton Power盐岩蠕变本构模型分别对盐穴储气库及盐岩井段井眼缩径进行了计算分析。该方法可充分利用试验数据进行回归分析得到较为可靠的本构关系式,但由于岩石取芯费用较高且室内蠕变时间较长,短期内较难获得大量可靠的室内试验数据。鉴于此,为了避免上述几种非线性蠕变模型构建方法的局限性,根据川东北油气田盐岩三轴蠕变试验结果,并参考其他学者在研究盐岩蠕变特性中采用的方法[7,10-12,18],本文作者将改进后的非线性元件组合模型与非线性元件进行串联,构建新的适用于川东北油气田嘉陵江组盐岩的非线性蠕变本构模型,以便能够更好地描述蠕变时间、温度和应力水平对盐岩衰减蠕变阶段和稳态蠕变阶段的非线性蠕变特征的影响。首先,在广义Kelvin模型的基础上用非定常黏壶替代线性黏壶构建非定常广义Kelvin模型,将其与能够描述盐岩稳态蠕变特征的Heard模型进行串联组合,构建新的NGKH非线性盐岩蠕变模型,以描述盐岩的瞬态蠕变和稳态蠕变过程;然后,基于套损现象集中发生的嘉陵江组盐岩地层的岩心室内蠕变试验结果,采用Levenberg-Marquardt算法分别对不同应力水平下的NGKH模型、非定常广义Kelvin模型和广义Kelvin模型参数进行辨识,通过比较验证新模型的适用性;最后,基于FLAC3D提供的接口程序对NGKH模型进行二次开发,并通过数值试验对其正确性进行验证。
1 盐岩蠕变本构模型
1.1 非定常广义Kelvin模型的构建
广义Kelvin模型是由胡克体与Kelvin模型串联而成,如图1所示。
![FX_GRP_ID80002BBA](/web/fileinfo/upload/magazine/12781/321271/image004.jpg)
图1 广义Kelvin模型
Fig. 1 Generalized Kelvin model
广义Kelvin模型的蠕变方程为[18]:
(1)
式中:
为串联胡克体的弹性模量;
为Kelvin体的黏弹性模量;
为Kelvin体的初始黏滞系数。
由于广义Kelvin模型只能描述岩石的线性蠕变特征,对于蠕变特性复杂的盐岩来说并不适用。将Kelvin体中牛顿黏壶的黏滞系数视为时间的函数,建立非定常黏壶[19]:
(2)
式中:
为非定常黏壶的黏滞系数;
为待定系数。
对式(2)求导可得
(3)
由式(3)可知:当
>0时,
≤0,黏滞系数随时间衰减;当
=0时,
,非定常黏壶退化为线性黏壶;当
<0时,
≥0,黏滞系数随时间增大。因此,可通过变化
来实现非定常Kelvin体的非线性蠕变特征。用上述非定常黏壶替换图1中广义Kelvin模型的线性黏壶,构建而成的非定常广义Kelvin蠕变模型如图2所示。
![FX_GRP_ID800032C0](/web/fileinfo/upload/magazine/12781/321271/image026.jpg)
图2 非定常广义Kelvin模型
Fig. 2 Non-stationary generalized Kelvin model
非定常Kelvin体的状态方程为
(4)
式中:
为蠕变速率。
分离变量求定积分,并由
可得其蠕变柔量
为
(5)
由式(1)和式(5)可得非定常广义Kelvin模型的蠕变柔量为
(6)
则非定常广义Kelvin模型的一维本构方程为
(7)
在三维应力状态下,岩石的应力张量可分解为偏应力张量
和球应力张量
。通常考虑应力偏量仅引起蠕变变形,球应力引起弹性变形,同时假定岩石为各向同性体,由式(1)可得广义Kelvin模型的三维蠕变方程为[19]
(8)
式中:
为应变张量;
和
分别为串联胡克体的体积模量和剪切模量;
为Kelvin体的剪切模量。
在等围压三轴压缩试验中
,代入式(8)可得常规三轴压缩试验中广义Kelvin模型的轴向应变方程:
![](/web/fileinfo/upload/magazine/12781/321271/image058.gif)
(9)
同理,由式(7)和式(8)可得非定常广义Kelvin模型的三维蠕变方程和轴向应变方程分别为[19]:
![](/web/fileinfo/upload/magazine/12781/321271/image062.gif)
(10)
![](/web/fileinfo/upload/magazine/12781/321271/image058.gif)
(11)
1.2 Heard模型的构建
研究盐岩的稳态蠕变本构方程即确定其稳态蠕变速率与温度、偏应力之间的变化关系,盐岩的蠕变机制和蠕变本构方程随其所处地质环境的不同而变化[20]。对于川东北油气田的深部盐岩地层来说,其蠕变机制属于位错滑移机制[6],符合Heard本构模型,本构方程如式(12)所示。曾德智等[21-22]利用Heard蠕变模型对盐岩地层套管蠕变外载、盐岩层井壁稳定性及井眼缩径量进行了计算分析;LI等[6]通过对盐岩和盐膏岩进行的室内三轴蠕变试验,分析了不同矿物成分和围压对盐岩蠕变速率的影响,并采用Heard模型对试验结果进行了拟合。
(12)
式中:
为稳态蠕变速率;
为激活能;
为摩尔气体常量,
=8.314 J/(mol·K);
为偏应力;
为热力学温度,K;
和
为流变常数。
按照GB/T 50266—2013“工程岩体试验方法标准”,将取自川东北某气田嘉陵江组地层埋深3 838~3 858 m处的岩石岩芯,加工成高为100 mm、直径为50 mm的标准圆柱体岩石试样共2组。该地层段盐岩层厚度为174 m,地层岩性以灰白色硬石膏岩、白色盐岩、灰白色膏盐岩为主,夹灰色灰质白云岩、深灰色、灰色白云岩、膏质灰岩、泥灰岩。本试验采用MTS材料试验机进行,该装置由轴向加压、侧向加压、温控设备和微机系统组成。对第一组岩石试样,在室温条件下将围压加载到15 MPa并保持不变,轴向采用单级加载方式分别施加30,35,40和45 MPa的均布载荷并保持不变,加载时间持续40 h左右。由前人研究成果可知:温度因素对盐岩的蠕变特性影响非常显著,因此,对于第二组岩石试样,首先将三轴室的温度分别加到50 ℃和100 ℃并保持3 h左右,再将围压加载至15 MPa并保持不变,轴向仍然采用单级加载方式施加30MPa均布载荷并保持不变,加载时间持续30 h左右。图3(a)和图3(b)所示分别为围压15 MPa时,不同偏应力及不同温度下的盐岩应变-时间曲线。基于Heard蠕变模型,采用Origin分析软件对图3中不同条件下的稳态蠕变率进行非线性拟合,得到的参数
,
和
见表1。
![FX_GRP_ID80003B4A](/web/fileinfo/upload/magazine/12781/321271/image084.jpg)
图 3 盐岩室内试验蠕变曲线
Fig. 3 Creep curves of salt rock
表1 拟合得到的Heard模型蠕变参数
Table 1 Creep parameters of fitted Heard model
![](/web/fileinfo/upload/magazine/12781/321271/image085.jpg)
1.3 NGKH模型的构建
王军保等[18]基于芒硝的室内三轴蠕变压缩试验,将分数阶广义Kelvin模型与未考虑温度因素的Heard模型相结合,构建了能够反映芒硝非线性蠕变特性的蠕变本构模型。而在不同地层深度和地层温度下,盐岩的蠕变过程受温度影响较为显著,因此,本文采用的Heard模型中仍然考虑温度影响因素。
上述建立的非定常广义Kelvin模型能够描述盐岩的瞬时变形及衰减蠕变状态,而Heard蠕变模型能够反映盐岩的偏应力、温度与稳态蠕变速率的关系,因此,参考岩石非线性蠕变模型的经典构建思路,将非定常广义Kelvin模型与Heard蠕变模型进行串联,构成新的盐岩非线性蠕变模型,简称NGKH蠕变模型(见图4),下式即为NGKH蠕变模型的轴向应变方程:
![](/web/fileinfo/upload/magazine/12781/321271/image058.gif)
![](/web/fileinfo/upload/magazine/12781/321271/image087.gif)
(13)
式(13)中第1项和第2项反映初始蠕变阶段的瞬时应变,第3项反映衰减蠕变阶段的黏弹性应变,第4项反映稳态蠕变阶段的黏性流动应变。
1.4 蠕变参数敏感性分析
在应力状态、温度相同的条件下,对式(13)中的待定系数
和
进行敏感性分析,图5(a)和图5(b)所示分别为系数
和
对NGKH模型蠕变曲线的影响。对式(3)进行分析可知
的变化能改变Kelvin体中非线性黏壶的蠕变性质。由图5(a)可以看出
的变化对蠕变曲线的形状尤其是衰减蠕变阶段影响较为显著:当
>0时,随着
增大,衰减蠕变阶段的蠕变速率衰减速度加快,进入稳态蠕变阶段的时间提前,此后各蠕变曲线几乎重合;当
<0时,随着
减小,衰减蠕变阶段的蠕变速率衰减速度同样加快,进入稳态蠕变阶段的时间同样提前,但此后各蠕变曲线近似平行,轴向应变呈近似线性减小。从图5(b)可知:
的变化主要影响稳态蠕变阶段的曲线形态,对衰减蠕变阶段则影响很小。随着
增大,稳态蠕变率明显增大,轴向应变也明显增大。蠕变参数
对蠕变曲线的影响规律与
的影响规律相似。
![FX_GRP_ID8000593A](/web/fileinfo/upload/magazine/12781/321271/image091.jpg)
图4 NGKH蠕变模型
Fig. 4 NGKH creep model
![FX_GRP_ID80000BA6](/web/fileinfo/upload/magazine/12781/321271/image093.jpg)
图5 蠕变参数
和
对NGKH模型蠕变曲线的影响
Fig. 5 Influence of parameters
and
on creep curves of NGKH model
将围压固定为15 MPa,在其他蠕变参数不变的情况下,分别观察应力
的变化对NGKH模型及广义Kelvin模型蠕变曲线的影响,如图6所示。
![FX_GRP_ID800046BB](/web/fileinfo/upload/magazine/12781/321271/image101.jpg)
图6 轴向载荷对NGKH模型和广义Kelvin模型蠕变曲线的影响
Fig. 6 Effect of axial load on creep curves of NGKH model and generalized Kelvin model
从图6可以看出:随着轴向载荷增大,NGKH模型及广义Kelvin模型的轴向应变均随之增大,但广义Kelvin模型的蠕变曲线呈线性变化,而NGKH模型的稳态蠕变速率和蠕变曲线均呈非线性变化。因此,在不同应力状态下,盐岩的蠕变曲线形态可通过变化NGKH模型中
,
和
来调整,使NGKH模型更加适用于岩石的流变分析。
2 蠕变模型参数辨识
针对本文提出的NGKH蠕变模型,式(13)中需要给出的参数分别是
,
,
,
,
,
,
和
,其中
,
和
已经通过室内蠕变试验曲线拟合分析确定。令
,再利用
,
与弹性模量
和泊松比
之间的关系求得
,进而反求出
和
,此时,需要拟合确定的参数变为
,
,
和
。根据图3中不同应力状态下盐岩室内三轴蠕变试验结果,运用Origin分析软件的Levenberg-Marquardt算法[23],通过自定义函数模型,对
,
,
和
进行辨识。表2、表3和表4所示分别为NGKH模型、非定常广义Kelvin模型和广义Kelvin模型的蠕变参数辨识结果。将拟合得到的各模型蠕变曲线与室内试验结果进行对比,如图7所示。其中图7(a)~7(d)所示分别为不同偏应力下的蠕变曲线对比,图7(a)、图7(e)和图7(f)所示分别为不同温度下的蠕变曲线对比。从表2~4可见:在相同应力状态下,不同蠕变模型的
、
和
较为接近,而
和
却存在明显差异。从图7(a)和图7(b)可知:
和
组成的非定常黏壶对衰减蠕变阶段起着控制作用,而Heard体的加入直接影响了稳态蠕变阶段的蠕变曲线走向,同时反映出NGKH模型的灵活性和适用性。
![FX_GRP_ID8000669E](/web/fileinfo/upload/magazine/12781/321271/image109.jpg)
图7 NGKH模型、非定常广义Kelvin模型和广义Kelvin模型拟合曲线对比
Fig. 7 Comparison of fitting curves of NGKH model, non-stationary generalized Kelvin model and generalized Kelvin model
表2 NGKH模型蠕变参数辨识
Table 2 Creep parameters identification for NGKH model
![](/web/fileinfo/upload/magazine/12781/321271/image111.jpg)
表3 非定常广义Kelvin模型参数辨识
Table 3 Creep parameters identification for non-stationary generalized Kelvin model
![](/web/fileinfo/upload/magazine/12781/321271/image113.jpg)
表4 广义Kelvin模型参数辨识
Table 4 Creep parameters identification for generalized Kelvin model
![](/web/fileinfo/upload/magazine/12781/321271/image115.jpg)
从图7可知:与非定常广义Kelvin模型及常规广义Kelvin模型相比,由于引入了非定常黏壶及Heard蠕变体,本文提出的NGKH模型的拟合效果相对较好,理论计算结果与室内试验结果吻合度明显更高,能够较好地对不同应力状态下盐岩的瞬时应变、瞬态蠕变阶段和稳态蠕变阶段进行描述,更能反映盐岩的非线性蠕变特征,由此验证了该模型的合理性和可靠性。
3 蠕变模型二次开发
FLAC3D软件内置的蠕变模型均有特定的使用范围,在一定程度上影响了FLAC3D的广泛应用[24]。因此,部分学者利用其提供的接口程序对新的蠕变模型进行了二次开发,例如熊良宵等[25-27]就基于FLAC3D分别对改进Burgers模型和河海模型进行了二次开发。为了能够将本文提出的NGKH蠕变模型程序化,首先将其本构方程改写成三维差分形式。
3.1 NGKH模型的有限差分形式
3.1.1 模型各部分的偏应力与偏应变关系式
由图4和式(13)可知:NGKH模型各部分之间为串联关系,因此应力相等,应变相加,则有[24,28]
(14)
(15)
对于串联弹簧胡克体,偏应力与偏应变速率的关系为
(16)
对于非定常Kelvin体,偏应力
和偏应变
有如下关系[28]:
(17)
式中:
为非定常Kelvin体的黏滞系数,且
。
对于串联Heard体,有
(18)
引入应力强度
,且
[28],在FLAC3D中可通过相应指针读取应力张量的各分量,然后求出
。由经典弹塑性力学可知
,因此,对于Heard体,偏应力
与偏应变速率
之间的关系为
(19)
3.1.2 模型各部分的增量形式
为了利用FLAC3D软件进行二次开发,将NGKH本构方程的上述内容改写成有限差分形式。将式(15)写成增量形式:
(20)
式(16)的增量形式可写为
(21)
将式(17)写成增量形式:
(22)
式中:
;
;
;
和
分别为非定常Kelvin体平均偏应力和平均偏应变;
和
分别为第i步和第i
1步的偏应力张量;
和
分别为第i步和第i
1步的偏应变张量。
将式(22)整理后,得到非定常Kelvin体第i步偏应变更新公式为
(23)
将式(23)进一步整理后可得:
(24)
对于非线性Heard体,偏应变张量
可表示为
(25)
式中:
为蠕变乘子,随应力和蠕变应变的变化而变化。
三轴等围压应力状态下,
,因此,
(26)
由此式(19)的增量形式可写为
(27)
将式(21)、式(24)和式(27)代入式(20),有
![](/web/fileinfo/upload/magazine/12781/321271/image185.gif)
![](/web/fileinfo/upload/magazine/12781/321271/image187.gif)
(28)
整理可得NGKH模型第i步应力更新公式为
![](/web/fileinfo/upload/magazine/12781/321271/image191.gif)
![](/web/fileinfo/upload/magazine/12781/321271/image193.gif)
(29)
令
,
,![](/web/fileinfo/upload/magazine/12781/321271/image201.gif)
,则式(29)可简化为
(30)
3.2 二次开发程序的验证
基于C++语言,利用Visual Studio2010平台将上述NGKH蠕变模型本构方程的差分形式进行编译写入动态链接库.dll文件中,利用FLAC3D的UDM接口程序进行二次开发。图8所示为NGKH模型的二次开发流程。
![FX_GRP_ID80000D04](/web/fileinfo/upload/magazine/12781/321271/image207.jpg)
图8 自定义NGKH本构模型二次开发流程
Fig. 8 Secondary development process of self-defined NGKH constitutive model
为了验证开发NGKH模型程序的正确性,建立与室内试验规格相同的直径为50 mm、高度为100 mm的盐岩圆柱体数值模型,共划分1 440个单元和1 595个节点,如图9所示。模型底面施加法向约束,顶面分别施加30,35,40和45 MPa的均布载荷,侧面施加15 MPa围压,分别进行三轴蠕变压缩试验。
![FX_GRP_ID8000062C](/web/fileinfo/upload/magazine/12781/321271/image209.jpg)
图9 盐岩的有限差分数值试样
Fig. 9 Finite difference numerical sample for salt rock
3.2.1 黏弹性特征验证
由式(11)可知:当
=0时,非定常广义Kelvin模型退化为广义Kelvin模型,因此,在NGKH模型中将
和
同时设置为0,此时NGKH模型即退化为广义Kelvin模型。在FLAC3D程序内置的常规Burgers模型中,对Maxwell体串联牛顿黏壶的黏性系数不进行赋值,即退化为广义Kelvin模型。退化后的NGKH模型和退化后的Burgers模型均属于三元件模型,按照表4所示的广义Kelvin参数分别对这2个退化模型进行赋值,计算不同轴向载荷下的盐岩试样蠕变曲线。图10(a)~10(d)所示分别为偏应力25 MPa时由退化NGKH模型计算得到的蠕变10,20,30和40 h后的盐岩数值试样位移分布云图,图11所示为不同偏应力下退化NGKH模型与退化Burgers模型的蠕变曲线对比图。
![FX_GRP_ID80000274](/web/fileinfo/upload/magazine/12781/321271/image211.jpg)
图10 不同蠕变时刻退化Burgers模型的轴向位移分布
Fig. 10 Axial deformation distributions of degenerated Burgers model at different creep time
![FX_GRP_ID80002579](/web/fileinfo/upload/magazine/12781/321271/image213.jpg)
图11 退化NGKH模型和退化Burgers模型的黏弹性蠕变曲线对比
Fig. 11 Comparison of viscoelastic creep curves between degraded NGKH model and degraded Burgers model
由图10可以看出:盐岩试样的顶部轴向位移最大,往底部方向逐渐减小;随着蠕变时间的增加,盐岩试样的最大轴向位移随之增大,但在20 h后轴向应变增速开始放缓。从图11可知:在非定常参数
及Heard体不起作用的情况下,退化NGKH模型与退化Burgers模型的瞬时应变一致,曲线吻合良好,表明开发的NGKH模型程序进行的黏弹性计算结果是可靠的。
3.2.2 非线性特征验证
按照表2中蠕变参数对NGKH蠕变模型赋值并进行三轴蠕变数值计算,图12所示为开发的NGKH模型程序模拟结果与室内试验结果的对比。其中图12(a)所示为不同偏应力下数值模拟结果与室内试验结果对比,图12(b)所示为不同温度下数值模拟结果与室内试验结果对比。从图12可以看出开发NGKH模型的程序模拟结果与室内试验结果吻合良好,证明该模型的非线性蠕变计算结果是可靠的。
![FX_GRP_ID800064C8](/web/fileinfo/upload/magazine/12781/321271/image215.jpg)
图12 NGKH模型的非线性计算结果与试验结果的对比
Fig. 12 Comparison between nonlinear computation results from NGKH model and creep test results
上述分析验证了本文提出的将包含非定常黏壶的非定常广义Kelvin模型与稳态蠕变Heard体进行串联,最终构建NGKH蠕变模型以描述盐岩的非线性蠕变特征思路的正确性以及该蠕变模型在FLAC3D中二次开发的正确性。开发得到的FLAC3D蠕变本构模型为盐岩地层井壁围岩稳定性分析及套损问题研究提供了参考。
4 结论
1) 将牛顿黏壶的黏滞系数视为时间的函数,采用非定常黏壶替换常规广义Kelvin模型中的线性牛顿黏壶,建立了非定常广义Kelvin模型,将其与非线性Heard体进行串联,组成能够描述盐岩瞬时蠕变阶段和稳态蠕变阶段的四元件非线性NGKH蠕变模型。
2) 基于盐岩的室内三轴蠕变试验结果,对NGKH模型的蠕变参数进行辨识。拟合结果与室内试验曲线吻合良好,表明与广义Kelvin模型及非定常广义Kelvin模型相比,本文提出的NGKH蠕变模型更加适合描述川东北油气田嘉陵江组盐岩的非线性蠕变特征。
3) 基于FLAC3D二次开发NGKH模型程序的黏弹性计算结果与非线性计算结果都是可靠的,从而证实了开发NGKH模型的正确性与合理性。
参考文献:
[1] 夏竹君, 杨锡兵, 曾敏梅, 等. 普光气田套管变形规律研究[J]. 石油天然气学报, 2010, 32(3): 271-273.
XIA Zhujun, YANG Xibing, ZENG Minmei, et al. Study on the rule of casing deformation in Puguang gas field[J]. Journal of Oil and Gas Technology, 2010, 32(3): 271-273.
[2] 王贵君, 张磊, 张昱文, 等. 盐岩流变-损伤-破坏特性的实验研究[J]. 工程力学, 2013, 30(4): 288-293.
WANG Guijun, ZHANG Lei, ZHANG Yuwen, et al. Experimental investigation on the creep-damage-rupture characteristics of salt rock[J]. Engineering Mechanics, 2013, 30(4): 288-293.
[3] 刘江, 杨春和, 吴文, 等. 盐岩蠕变特性和本构关系研究[J]. 岩土力学, 2006, 27(8): 1267-1271.
LIU Jiang, YANG Chunhe, WU Wen, et al. Study on creep characteristics and constitutive relation of rock salt[J]. Rock and Soil Mechanics, 2006, 27(8): 1267-1271.
[4] 陈锋, 李银平, 杨春和, 等. 云应盐矿盐岩蠕变特性试验研究[J]. 岩石力学与工程学报, 2006, 25(增1): 3022-3027.
CHEN Feng, LI Yinping, YANG Chunhe, et al. Experimental study on creep behaviors of rock salt in Yunying salt mine[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(S1): 3022-3027.
[5] CORNET J S, DABROWSKI M, SCHMID D W. Long term creep closure of salt cavities[J]. International Journal of Rock Mechanics and Mining Sciences, 2018, 103: 96-106.
[6] LI Ping, DENG Jingen, ZHAO Wenlong, et al. An experimental study on creep characteristics of salt rock and gypsum-salt rock in Puguang gas field[J]. Petroleum Science and Technology, 2012, 30(16): 1715-1724.
[7] 孙钧. 岩石流变力学及其工程应用研究的若干进展[J]. 岩石力学与工程学报, 2007, 26(6): 1081-1106.
SUN Jun. Rock rheological mechanics and its advance in engineering applications[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(6): 1081-1106.
[8] 刘开云, 薛永涛, 周辉. 参数非定常的软岩非线性黏弹塑性蠕变模型[J]. 中国矿业大学学报, 2018, 47(4): 921-928.
LIU Kaiyun, XUE Yongtao, ZHOU Hui. A nonlinear viscoelastic-plastic creep model of soft rock with unsteady parameters[J]. Journal of China University of Mining and Technology, 2018, 47(4): 921-928.
[9] 王军保, 刘新荣, 郭建强, 等. 盐岩蠕变特性及其非线性本构模型[J]. 煤炭学报, 2014, 39(3): 445-451.
WANG Junbao, LIU Xinrong, GUO Jianqiang, et al. Creep properties of salt rock and its nonlinear constitutive model[J]. Journal of China Coal Society, 2014, 39(3): 445-451.
[10] 易其康, 马林建, 刘新宇, 等. 考虑频率影响的盐岩变参数蠕变损伤模型[J]. 煤炭学报, 2015, 40(S1): 93-99.
YI Qikang, MA Linjian, LIU Xinyu, et al. A variable parameters-based creep damage model considering the effect of frequency for rock salt[J]. Journal of China Coal Society, 2015, 40(S1): 93-99.
[11] ZHOU Hongwei, WANG Chunping, HAN Bingbing, et al. A creep constitutive model for salt rock based on fractional derivatives[J]. International Journal of Rock Mechanics and Mining Sciences, 2011, 48(1): 116-121.
[12] 丁靖洋, 周宏伟, 陈琼, 等. 盐岩流变损伤特性及本构模型研究[J]. 岩土力学, 2015, 36(3): 769-776.
DING Jingyang, ZHOU Hongwei, CHEN Qiong, et al. Characters of rheological damage and constitutive model of salt rock[J]. Rock and Soil Mechanics, 2015, 36(3): 769-776.
[13] 高小平, 杨春和, 吴文, 等. 盐岩蠕变特性温度效应的实验研究[J]. 岩石力学与工程学报, 2005, 24(12): 2054-2059.
GAO Xiaoping, YANG Chunhe, WU Wen, et al. Experimental studies on temperature dependent properties of creep of rock salt[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(12): 2054-2059.
[14] 杨春和, 陈锋, 曾义金. 盐岩蠕变损伤关系研究[J]. 岩石力学与工程学报, 2002, 21(11): 1602-1604.
YANG Chunhe, CHEN Feng, ZENG Yijin. Investigation on creep damage constitutive theory of salt rock[J]. Chinese Journal of Rock Mechanics and Engineering, 2002, 21(11): 1602-1604.
[15] 杜超, 杨春和, 马洪岭, 等. 深部盐岩蠕变特性研究[J]. 岩土力学, 2012, 33(8): 2451-2456.
DU Chao, YANG Chunhe, MA Hongling, et al. Study of creep characteristics of deep rock salt[J]. Rock and Soil Mechanics, 2012, 33(8): 2451-2456.
[16] 王安明, 李小根, 杨春和, 等. 层状盐岩蠕变变形相互作用研究[J]. 岩土力学, 2010, 31(12): 3964-3970.
WANG Anming, LI Xiaogen, YANG Chunhe, et al. Study of interaction between creep deformation of bedded salt rock[J]. Rock and Soil Mechanics, 2010, 31(12): 3964-3970.
[17] 黄小兰, 杨春和, 李银平. 蠕变作用下层状盐岩界面剪切应力变化规律研究[J]. 地下空间与工程学报, 2014, 10(3): 547-551.
HUANG Xiaolan, YANG Chunhe, LI Yinping. Interface shear stress distribution of bedded salt rock considering rheology effect[J]. Chinese Journal of Underground Space and Engineering, 2014, 10(3): 547-551.
[18] 王军保, 刘新荣, 张倩倩, 等. 芒硝蠕变特性及本构模型研究[J]. 四川大学学报(工程科学版), 2015, 47(5): 78-85.
WANG Junbao, LIU Xinrong, ZHANG Qianqian, et al. Study on creep properties and constitutive model of thenardite[J]. Journal of Sichuan University(Engineering Science Edition), 2015, 47(5): 78-85.
[19] 韩阳, 谭跃虎, 李二兵, 等. 岩石非定常Burgers蠕变模型及其参数识别[J]. 工程力学, 2018, 35(3): 210-217.
HAN Yang, TAN Yuehu, LI Erbing, et al. Non-stationary Burgers creep model of rock and its parameter identification[J]. Engineering Mechanics, 2018, 35(3): 210-217.
[20] DUAN Kang, ZHANG Qiangyong, CAI Bing, et al. Experimental study on creep property of similar material for Jintan salt rock gas storage medium[J]. Applied Mechanics and Materials, 2012, 170/171/172/173: 772-776.
[21] 曾德智, 林元华, 卢亚锋, 等. 复合岩盐层井眼蠕变缩径的数值模拟[J]. 地质力学学报, 2012, 18(2): 158-164.
ZENG Dezhi, LIN Yuanhua, LU Yafeng, et al. Numerical simulation of hole creep shrinkage for composite rock salt layers[J]. Journal of Geomechanics, 2012, 18(2): 158-164.
[22] 葛伟凤, 陈勉, 金衍, 等. 深部盐膏岩地层套管磨损后等效应力分析[J]. 中国石油大学学报(自然科学版), 2013, 37(1): 75-79.
GE Weifeng, CHEN Mian, JIN Yan, et al. Analysis of equivalent stress on casings after casing wear in deep salt-gypsum formation[J]. Journal of China University of Petroleum(Edition of Natural Science), 2013, 37(1): 75-79.
[23] 何志磊, 朱珍德, 朱明礼, 等. 基于分数阶导数的非定常蠕变本构模型研究[J]. 岩土力学, 2016, 37(3): 737-744.
HE Zhilei, ZHU Zhende, ZHU Mingli, et al. An unsteady creep constitutive model based on fractional order derivatives[J]. Rock and Soil Mechanics, 2016, 37(3): 737-744.
[24] 杨文东, 张强勇, 张建国, 等. 基于FLAC3D的改进Burgers蠕变损伤模型的二次开发研究[J]. 岩土力学, 2010, 31(6): 1956-1964.
YANG Wendong, ZHANG Qiangyong, ZHANG Jianguo, et al. Second development of improved Burgers creep damage constitutive model of rock based on FLAC3D[J]. Rock and Soil Mechanics, 2010, 31(6): 1956-1964.
[25] 熊良宵, 杨林德, 张尧. 硬岩的复合黏弹塑性流变模型[J]. 中南大学学报(自然科学版), 2010, 41(4): 1540-1548.
XIONG Liangxiao, YANG Linde, ZHANG Yao. Composite viscoelasto-plastic rheological model for hard rock[J]. Journal of Central South University(Science and Technology), 2010, 41(4): 1540-1548.
[26] 张涛, 何利军. 含分数阶导数元件非线性蠕变模型的二次开发[J]. 华东交通大学学报, 2017, 34(5): 21-28.
ZHANG Tao, HE Lijun. Secondary development of non-linear creep model with fractional order derivative elements[J]. Journal of East China Jiaotong University, 2017, 34(5): 21-28.
[27] 徐卫亚, 杨圣奇, 褚卫江. 岩石非线性黏弹塑性流变模型(河海模型)及其应用[J]. 岩石力学与工程学报, 2006, 25(3): 433-447.
XU Weiya, YANG Shengqi, CHU Weijiang. Nonlinear viscoelasto-plastic rheological model (Hohai model) of rock and its engineering application[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(3): 433-447.
[28] 王军保. 不同加载路径下盐岩蠕变力学特性与盐岩储气库长期稳定性研究[D]. 重庆: 重庆大学土木工程学院, 2012: 100-102.
WANG Junbao. Study on the creep mechanical properties of salt rock under different loading paths and long-term stability of salt rock gas storage[D]. Chongqing: Chongqing University. College of Civil Engineering, 2012: 100-102.
(编辑 杨幼平)
收稿日期: 2019 -07 -14; 修回日期: 2019 -09 -18
基金项目(Foundation item):国家自然科学基金资助项目(51804330);国家科技重大专项(2016ZX05017-003);中央高校基本科研业务费专项资金资助项目(18CX02154A) (Project(51804330) supported by the National Natural Science Foundation of China; Project(2016ZX05017-003) supported by the National Science and Technology Major Program; Project(18CX02154A) supported by Fundamental Research Funds for the Central Universities)
通信作者:闫怡飞,博士,讲师,从事油气安全工程研究;E-mail:yanyf163@163.com