Impact of the Loading Method on Gas Seal Performance of Set Cement
-
摘要: 页岩气井水泥环的气密封性能至关重要,而水泥环在页岩气井生产过程中所受的载荷复杂,其气密封性必然会受到影响,进而会影响整个井筒的完整性。因此,为了了解水泥环在不同载荷条件下的密封失效形式和各因素的影响规律,在对水泥石进行三轴循环加载试验的基础上,对其渗透率和声发射特性进行监测,研究了水泥石气密封性在不同加载方式下的变化。结果表明:水泥石在加载过程中,有扩容损伤的特性,且渗透率会随着扩容点的出现而迅速升高;循环加卸载峰值高于扩容点应力时,循环加卸载对水泥石的损伤严重,渗透率升高较快,气密封性失效风险明显增大。研究表明,添加增韧材料可以提高水泥石的抗损伤能力,从而提高水泥环的气密封性。Abstract: Gas sealing performance of the cementsheath in natural gas wells is very important.Becausethe cement sheath is subjected to complicated loads during natural gas well production,gas sealing is inevitably affected by such loads and furtherimpacts wellbore integrity.In order to investigate the sealing failure modeof cement sheath under different loading conditions and the influence laws of various factors,the gas sealing performance of set cement under different loading modes was studied on the basis of triaxial cyclic loading tests on set cement while monitoring its permeability and acoustic emission characteristics.Experimental results show that dilatancydamageoccurs,and the permeability of set stone increases rapidly with the appearance of dilatancy point during the loading process.When the peak load is higher than the stress on dilatancy point,the cyclic loading-unloadingprocess will seriously damage theset cement,the permeability rises rapidly,and the risk of gas sealing failure increases significantly.The research demonstrates that the gas sealing performance of cement sheath can be improved by enhancing the set cement damage resistance capability.
-
Keywords:
- well cementing /
- cement sheath /
- damage /
- cyclic loading /
- gas sealing performance
-
在进行井眼轨道设计[1-7]、随钻测量数据处理时[8-15],需要计算北东坐标、垂深、水平投影长度等参数,如果井眼轨道上每一点的井斜角和方位角与井深呈线性关系,则这些参数可以通过解析计算公式计算;如果井斜角或者方位角与井深呈非线性关系,则北东坐标或者水平投影长度就需要使用数值积分法来计算[16]。刘修善[16]提出采用定积分表示这些参数,但是并没有详细说明如何计算这些定积分。鲁港等人[17-18]推导出了空间圆弧轨道的井斜角计算公式,并使用数值积分法来计算水平投影长度[19]。对于恒工具面模型,还没有查阅到讨论如何计算北东坐标数值积分的公开文献,而开发钻井软件过程中,也会遇到恒工具面模型北东坐标数值积分的问题,因此,笔者对积分参数进行了数学变换,使积分函数进一步化简,从而快速计算数值积分。
1. 井眼轨道模型
井眼轨道在数学上可以抽象为三维空间中的光滑曲线,这意味曲线上的每一点处都有切线,并且切线随切点的移动而连续转动。
以井口为坐标原点
O ,以正北方向为x轴,以正东方向为y轴,以指向地心方向为z轴,建立三维直角坐标系O-xyz,称为整体坐标系,也称为井眼坐标系。1.1 井眼轨道微分方程
井眼坐标系中井眼轨道参数可以写成下面的形式:
{N=f(l)E=g(l)(L1⩽ (1) 式中:
N 为北坐标,m;E 为东坐标,m;H 为垂深,m;l 为井深,m;{L_1} 为起始井深,m;{L_2} 为结束井深,m;f(l) ,g(l) 和h(l) 均为井深l 的函数。式(1)也可以写成更简洁的矢量形式:
{\boldsymbol{r}} = {\boldsymbol{r}}(l) = \left[ {\begin{array}{*{20}{c}} {f(l)} \\ {g(l)} \\ {h(l)} \end{array}} \right]{\text{ }}({L_1} \leqslant l \leqslant {L_2}) (2) 式中:
{\boldsymbol{r}} 为位置矢量;{\boldsymbol{r}}(l) 为井深l 的矢量函数。切线矢量
{\boldsymbol{t}} 定义为:{\boldsymbol{t}} = \frac{{{\boldsymbol{\dot r}}}}{{\left\| {{\boldsymbol{\dot r}}} \right\|}} = \frac{1}{{\sqrt {{{[f'(l)]}^2} + {{[g'(l)]}^2} + {{[h'(l)]}^2}} }}\left[ {\begin{array}{*{20}{c}} {f'(l)} \\ {g'(l)} \\ {h'(l)} \end{array}} \right] (3) 井眼轨道的光滑性要求
{[f'(l)]^2} + {[g'(l)]^2} + {[h'(l)]^2} \ne 0 。切线矢量是单位长度矢量,也称为井眼方向矢量,它与z轴的夹角为井斜角,记为
\alpha ;x轴顺时针转到切线矢量水平投影所转过的角度为井斜方位角,简称方位角,记为\phi 。因此切线矢量可写为[17-23]:{\boldsymbol{t}} = \left[ {\begin{array}{*{20}{c}} {\sin\alpha \cos\phi } \\ {\sin\alpha \sin\phi } \\ {\cos\alpha } \end{array}} \right] (4) 式中:
\alpha 为井斜角,rad;\phi 为井斜方位角,rad;{\boldsymbol{t}} 为切线矢量。将井眼轨道投影到水平面上,井眼轨道上任一点
P 在投影曲线上的对应点为P' ,则投影曲线上开始点到P' 的长度称为水平投影长度,记为S 。因此,得到井眼轨道参数的微分表达式:
\left\{ {\begin{array}{*{20}{l}} {\dfrac{{{\text{d}}N}}{{{\text{d}}l}} = \sin \alpha \cos \phi }&{} \\ {\dfrac{{{\text{d}}E}}{{{\text{d}}l}} = \sin \alpha \sin \phi }&{({L_1} \leqslant l \leqslant {L_2})} \\ {\dfrac{{{\text{d}}H}}{{{\text{d}}l}} = \cos \alpha }&{} \\ {{\text{d}}S/{\text{d}}l = \sin \alpha }&{} \end{array}} \right. (5) 式中:
S 为水平投影长度,m。如果已知井斜角和方位角随井深的变化规律,就可以通过式(5)计算出北东坐标、垂深、水平投影长度。
\left\{ {\begin{array}{*{20}{l}} {\Delta N = \displaystyle\int_{{L_1}}^L {\sin \alpha \cos \phi {\text{d}}l} } \\ {\Delta E = \displaystyle\int_{{L_1}}^L {\sin \alpha \sin \phi {\text{d}}l} } \\ {\Delta H = \displaystyle\int_{{L_1}}^L {\cos \alpha {\text{d}}l} } \\ {\Delta S = \displaystyle\int_{{L_1}}^L {\sin \alpha {\text{d}}l} } \end{array}} \right. (6) 式中:
\Delta N 为北坐标增量,m;\Delta E 为东坐标增量,m;\Delta H 为垂深增量,m;\Delta S 为水平投影长度增量,m。1.2 恒工具面模型
恒工具面模型的特征是井眼轨道上每一点的井眼曲率和工具面角保持为常数。
井斜角和方位角的变化规律为[16]:
\, \left\{ {\begin{array}{*{20}{l}} {\alpha = {\alpha _1} + (\kappa \cos \omega )(l - {L_1})} \\ {\phi = {\phi _1} + \tan \omega \left( {\ln \tan \dfrac{\alpha }{2} - \ln \tan \dfrac{{{\alpha _1}}}{2}} \right)} \end{array}} \right.{\text{ }}({L_1} \leqslant l \leqslant {L_2}) (7) 式中:
\kappa 为井眼曲率,rad/m;\omega 为工具面角,rad;{\alpha _1} 和{\phi _1} 分别为起始井深处的井斜角和方位角,rad。垂深、水平投影长度增量的积分形式有封闭形式的表达式,而北东坐标增量的积分一般没有封闭形式的表达式,需要采用数值积分法来计算。
将积分变量转化为
\alpha ,得到:\begin{split} \Delta N =& \frac{1}{{\kappa \cos \omega }}\int_{{\alpha _1}}^\alpha \cos \left[ {{\phi _1} + \tan \omega \left( {\ln \tan \frac{\alpha }{2} - \ln \tan \frac{{{\alpha _1}}}{2}} \right)} \right]\cdot\\ &\sin \alpha {\text{d}}\alpha \\[-10pt] \end{split} (8) 由于
\tan \dfrac{\alpha }{2} = \sqrt {\dfrac{{1 - \cos \alpha }}{{1 + \cos \alpha }}} 、{\text{d}}(\cos \alpha ) = - \sin \alpha {\text{d}}\alpha ,令t = \cos \alpha ,则式(8)变为:\Delta N = \frac{1}{{\kappa \cos \omega }}\int_{\cos {\alpha _1}}^{\cos \alpha } {\cos \left[ {{\phi _1} + \frac{{\tan \omega }}{2}\left( {\ln \frac{{1 - {t_1}}}{{1 + {t_1}}} - \ln \frac{{1 - t}}{{1 + t}}} \right)} \right]{\text{d}}t} (9) 类似地,有:
\Delta E = \frac{{ - 1}}{{\kappa \cos \omega }}\int_{\cos {\alpha _1}}^{\cos \alpha } {\sin \left[ {{\phi _1} + \frac{{\tan \omega }}{2}\left( {\ln \frac{{1 - {t_1}}}{{1 + {t_1}}} - \ln \frac{{1 - t}}{{1 + t}}} \right)} \right]{\text{d}}t} (10) 1.3 空间圆弧模型
空间圆弧模型是空间斜平面上的一段圆弧,井眼曲率为常数。它在水平面上的投影是一段椭圆曲线,所以水平投影长度(积分)没有封闭形式的表达式,需要使用数值积分法来计算。
{\boldsymbol{t}} = {{\boldsymbol{n}}_1}\sin \theta + {{\boldsymbol{t}}_1}\cos \theta (11) \,其中\qquad \qquad\qquad \theta = \kappa (l - {L_1})\quad\qquad (12) {{\boldsymbol{n}}_1} = {{\boldsymbol{h}}_1}\cos {\omega _1} + {{\boldsymbol{v}}_1}\sin {\omega _1} (13) {{\boldsymbol{h}}_1} = \left[ {\begin{array}{*{20}{c}} {\cos {\alpha _1}\cos {\phi _1}} \\ {\cos {\alpha _1}\sin {\phi _1}} \\ { - \sin {\alpha _1}} \end{array}} \right] (14) {{\boldsymbol{v}}_1} = \left[ {\begin{array}{*{20}{c}} { - \sin {\phi _1}} \\ {\cos {\phi _1}} \\ 0 \end{array}} \right] (15) 式中:
\theta 为圆心角,rad;{\omega _1} 为起始井深处的工具面角,rad;{{\boldsymbol{n}}_1} ,{{\boldsymbol{h}}_1} 和{{\boldsymbol{v}}_1} 分别为起始井深处的内法线矢量、井眼高边矢量和方位垂直矢量。如果令
{{\boldsymbol{n}}_1} = {[{a_1},{b_1},{c_1}]^{\text{T}}} ,{n_1} = \cos {\alpha _1} ,则由式(11)得到井斜角的变化规律:\cos \alpha = {c_1}\sin \theta + {n_1}\cos \theta (16) 2. 数值积分
数值积分法需要反复计算被积函数的值,因此在选定了数值积分公式后,要提高数值积分法的计算速度,主要就是减少被积函数所用的计算量。测试了C++语言中正弦函数、余弦函数、对数函数、乘法运算1亿次的用时,共进行了10次测试,结果见表1。测试软件环境为C++Builder 6.0,硬件为Inter i5-700U CPU 2.60 GHz。
表 1 主要数学函数的计算速度Table 1. Calculation speed of major mathematical functions测试序号 运算1亿次的用时/ms 余弦函数 正弦函数 对数函数 乘法 1 4 964 4 521 3 139 282 2 4 960 4 527 3 121 284 3 4 953 4 507 3 135 282 4 4 953 4 505 3 137 282 5 4 975 4 539 3 147 285 6 4 966 4 538 3 142 285 7 4 981 4 542 3 145 285 8 4 970 4 530 3 144 289 9 4 969 4 538 3 142 286 10 4 961 4 522 3 135 287 平均 4965 4527 3139 285 在不同的软硬件环境中,不同函数的计算速度略有不同,但存在这样的统计结论,乘法、对数函数、正弦函数和余弦函数的计算速度从快到慢的排序是乘法、对数函数、正弦函数和余弦函数。因此,进行数值积分计算时,要尽可能地用乘法代替对数和三角函数,能使用对数就尽量不使用三角函数,能使用正弦函数就不使用余弦函数。
2.1 变步长Simpson数值积分法
定积分为:
I = \int_a^b {F{\text{(}}x{\text{)d}}x} (17) 计算该定积分的步骤为[24]:
1)
n = 1 ,h = b - a ,用梯形公式计算。{T_n} = \frac{{F(a) + F(b)}}{2}h (18) 2)用变步长梯形法计算。
{x_k} = a + kh (19) {T_{2n}} = \frac{1}{2}{T_n} + \frac{h}{2}\sum\limits_{k = 0}^{n - 1} {F\left({x_k} + \frac{h}{2}\right)} (20) 3)用Simpson求积公式计算。
{I_{2n}} = \frac{{4{T_{2n}} - {T_n}}}{3} (21) 4)用
{I_n} 表示第n 步计算出来的积分值。如果\left| {{I_{2n}} - {I_n}} \right| < {\varepsilon _1} ,则结束,{I_{2n}} 即为所求积分的近似值。5)如果
\left| {{I_{2n}} - {I_n}} \right| \geqslant {\varepsilon _1} ,令n = 2n ,h = h/2 ,返回2)。{\varepsilon _1} 为给定的精度要求,例如{\varepsilon _1} = 1\times {10^{ - 6}} 。2.2 恒工具面法北东坐标增量的计算
\,令\qquad \qquad \quad{\phi _0} = {\phi _1} + \frac{{\tan \omega }}{2}\ln \frac{{1 - {t_1}}}{{1 + {t_1}}} \quad (22) \varphi = \frac{{\tan \omega }}{2}\ln \frac{{1 - t}}{{1 + t}} (23) 则式(9)中的被积函数可以写成:
\cos ({\phi _0} - \varphi ) = \cos {\phi _0}\cos \varphi + \sin {\phi _0}\sin \varphi (24) 式(10)中的被积函数可以写成:
\sin ({\phi _0} - \varphi ) = \sin {\phi _0}\cos \varphi - \cos {\phi _0}\sin \varphi (25) \,令\qquad \qquad s = \int_{\cos {\alpha _1}}^{\cos \alpha } {\cos \left( {\frac{{\tan \omega }}{2}\ln \frac{{1 - t}}{{1 + t}}} \right){\text{d}}t} \quad (26) c = \int_{\cos {\alpha _1}}^{\cos \alpha } {\sin \left( {\frac{{\tan \omega }}{2}\ln \frac{{1 - t}}{{1 + t}}} \right){\text{d}}t} (27) 则式(9)和式(10)可以写成:
\Delta N = \frac{{s\cos {\phi _0} + c\sin {\phi _0}}}{{\kappa \cos \omega }} (28) \Delta E = \frac{{c\cos {\phi _0} - s\sin {\phi _0}}}{{\kappa \cos \omega }} (29) 由于
\left| s \right| \leqslant 1 ,\left| c \right| \leqslant 1 ,必定存在2个参数r 和{\varphi _0} ,使s = r\sin {\varphi _0} ,c = r\cos {\varphi _0} ,于是式(28)和式(29)可以写成:\Delta N = r\frac{{\sin ({\varphi _0} + {\phi _0})}}{{\kappa \cos \omega }} (30) \Delta E = r\frac{{\cos ({\varphi _0} + {\phi _0})}}{{\kappa \cos \omega }} (31) 变步长Simpson积分算法并不只是用来计算
\Delta N 和\Delta E ,而是使计算\Delta N 和\Delta E 的算法做些稍许改变,使计算\Delta N 和\Delta E 时对数函数的计算次数减少。2.3 空间圆弧法水平投影长度增量的计算
\Delta S = \int_{{L_1}}^L {\sqrt {1 - {{\cos }^2}\alpha } {\text{d}}l} (32) 通过积分参数变换
\theta = \kappa (l - {L_1}) 、{\text{d}}\theta = \kappa {\text{d}}l ,得:\Delta S = \frac{1}{\kappa }\int_0^\varepsilon {\sqrt {1 - {{({c_1}\sin \theta + {n_1}\cos \theta )}^2}} {\text{d}}\theta } (33) 通过以下变换:
{b_2} = \frac{{c_1^2 + n_1^2}}{2} (34) {a_2} = 1 - {b_2} (35) {c_2} = {c_1}{n_1} (36) {r_2} = \sqrt {b_2^2 + c_2^2} (37) {e_2} = \frac{{{a_2}}}{{2{r_2}}} (38) {\gamma _2} = \arctan \frac{{{b_2}}}{{{c_2}}} (39) {\tau_1} = \tan \frac{{{\gamma _2} - 2\varepsilon }}{2} (40) {\tau_2} = \tan \frac{{{\gamma _2}}}{2} (41) 式(33)可以写成:
\Delta S = \frac{{\sqrt {2{r_2}} }}{\kappa }\int_{{\tau_1}}^{{\tau_2}} {\frac{1}{{1 + {\tau^2}}}\sqrt {{e_2} + \frac{\tau}{{1 + {\tau^2}}}} {\text{d}}\tau} (42) 式(33)需要进行2次三角函数、3次乘法、2次加法和1次开方运算;式(42)需要进行3次乘法、1次除法、2次加法和1次开方,不需要进行三角函数运算,计算量显然比式(33)要减少很多。
3. 算 例
算例1为文献[16]中的例5-4。假定井眼轨道符合恒工具面模型,开始点的井斜角为30°、方位角为135°,井眼曲率为12°/30m、工具面角为320°,计算距开始点100 m长井眼轨道的参数。
采用上文给出的数值积分算法进行计算,结果(见表2)与文献[16]中表5-9完全相同(小数点后保留2位数字),文献[16]在算例5-4中给出了
\Delta L = 60 m的计算细节,计算结果小数点后保留了4位数字,也与本文计算结果相同(见表2中标记为*60的行)。另外,也用Compass软件进行了验证,Compass软件没有计算\Delta S 的功能,在保留小数点后2位数字的情况下,该软件计算出的其他井眼轨道参数与表2中的相同。表 2 恒工具面模型井眼轨道参数Table 2. Borehole trajectory parameters in constant tool face model\Delta L /m \alpha /(°) \phi/(°) \Delta N /m \Delta E /m \Delta H /m \Delta S /m 0 30.00 135.00 0.00 0.00 0.00 0.00 10 33.06 130.08 −3.53 3.86 8.52 5.23 20 36.13 125.55 −7.00 8.34 16.75 10.91 30 39.19 121.34 −10.37 13.44 24.67 17.02 40 42.26 117.40 −13.56 19.13 32.25 23.54 50 45.32 113.68 −16.54 25.37 39.46 30.46 60 48.39 110.15 −19.26 32.14 46.30 37.75 *60 48.3851 110.1535 −19.2585 32.1413 46.3024 37.7532 70 51.45 106.79 −21.68 39.40 52.74 45.40 80 54.51 103.57 −23.77 47.10 58.76 53.39 90 57.58 100.47 −25.49 55.21 64.35 61.68 100 60.64 97.47 −26.83 63.69 69.48 70.26 算例2为文献[16]中的例5-2。假定空间圆弧井段开始点的井斜角为40°、方位角为60°,井眼曲率为10°/30 m、工具面角为315°,计算距开始点100 m长井眼轨道的参数。
采用空间圆弧的矢量描述法[22]和上文给出的数值积分算法进行计算,结果(见表3)与文献[16]中的表5-3完全相同(小数点后保留2位数字),文献[16]还在算例5-2中给出了
\Delta L = 60 m的计算细节,计算结果小数点后保留了4位数字,也与本文计算结果相同(见表3中标记为*60的行)。表 3 空间圆弧模型井眼轨道参数Table 3. Borehole trajectory parameters in space arc model\Delta L /m \alpha /(°) \phi/(°) \Delta N /m \Delta E /m \Delta H /m \Delta S /m 0 40.00 60.00 0.00 0.00 0.00 0.00 10 42.41 56.51 3.47 5.60 7.52 6.59 20 44.92 53.32 7.44 11.24 14.76 13.49 30 47.51 50.41 11.90 16.92 21.68 20.71 40 50.17 47.74 16.83 22.60 28.26 28.23 50 52.88 45.27 22.22 28.28 34.48 36.06 60 55.64 42.96 28.05 33.92 40.32 44.18 *60 55.6400 42.9644 28.0498 33.9239 40.3234 44.1769 70 58.44 40.81 34.30 39.52 45.76 52.57 80 61.28 38.78 40.94 45.06 50.78 61.21 90 64.14 36.87 47.96 50.50 55.37 70.10 100 67.03 35.04 55.33 55.85 59.50 79.20 4. 结 论
1)利用恒工具面法计算北东坐标,采用数值积分法,经过积分参数变换,被积函数从需要进行7次三角函数和1次对数函数运算减少为进行3次三角函数和1次对数函数运算,计算量减少了约一半。
2)利用空间圆弧法计算水平投影长度时,采用数值积分法,经过积分参数变换,将被积函数中的三角函数运算简化为普通四则运算和开方运算,计算量约为原来的十分之一。
3)定积分的数值计算采用变步长Simpson方法,计算精度可以控制到给定的精度,弥补了使用固定步长Simpson方法需要更多的计算量才能达到给定精度的不足。
4)算例表明,本文算法的计算结果在小数点后保留4位数字的情况下与公开数据完全相同,在小数点后保留2位数字的情况下与Compass软件完全相同。本文算法的可靠性得到了验证。
-
[1] 尹丛彬,叶登胜,段国彬,等.四川盆地页岩气水平井分段压裂技术系列国产化研究及应用[J].天然气工业,2014,34(4):67-71. YIN Congbin,YE Dengsheng,DUAN Guobin,et al.Research about and application of autonomous staged fracturing technique series for horizontal well stimulation of shale gas reservoirs in the Sichuan Basin[J].Natural Gas Industry,2014,34(4):67-71. [2] GOODWIN K J,CROOK R J.Cement sheath stress failure[J].SPE Drilling Engineering,1992,7(4):291-296.
[3] JACKSON P B,MURPHEY C E.Effect of casing pressure on gas flow through a sheath of set cement[R].SPE 25698,1993.
[4] WANG W,TALEGHANI A D.Three-dimensional analysis of cement sheath integrity around wellbores[J].Journal of Petroleum Science and Engineering,2014,121:38-51.
[5] XU Honglin,ZHANG Zhi,SHI Taihe,et al.Influence of the WHCP on cement sheath stress and integrity in HTHP gas well[J].Journal of Petroleum Science and Engineering,2015,126:174-180.
[6] GHOLAMI R,AADNOY B,FAKHARI N.A hermos-poroelastic analytical approach to evaluate cement sheath integrity in deep vertical wells[J].Journal of Petroleum Science and Engineering,2016,147:536-546.
[7] 于荣泽,卞亚南,张晓伟,等.页岩储层非稳态渗透率测试方法综述[J].科学技术与工程,2012,12(27):7019-7027,7035. YU Rongze,BIAN Yanan,ZHANG Xiaowei,et al.The review of non-steady permeability test in shale gas reservoir[J].Science Technology and Engineering,2012,12(27):7019-7027,7035. [8] BRIDGMAN P W.Volume changes in the plastic stages of simple compression[J].Journal of Applied Physics,1949,20(12):1241-1251.
[9] DZ/T 0276.19-2015岩石物理力学性质试验规程:第19部分:岩石单轴压缩变形试验[S]. DZ/T 0276.19-2015 Regulation for testing the physical and mechanical properties of rock:part 19:test for determing the deformability of rock in unizxial compression[S]. [10] 许东俊,李小春,蔡忠理,等.应力状态与岩石扩容特性[J].岩土力学,1992,13(2/3):37-44. XU Dongjun,LI Xiaochun,CAI Zhongli.Stress state and dilatation property of rock[J].Rock and Soil Mechanics,1992,13(2/3):37-44. [11] 张哨楠.致密天然气砂岩储层:成因和讨论[J].石油与天然气地质,2008,29(1):1-10. ZHANG Shaonan.Tight sandstone gas reservoirs:their origin and discussion[J].Oil Gas Geology,2008,29(1):1-10. [12] SHADRAVAN A,KIAS E,LEW R,et al.Utilizing the evolving cement mechanical properties under fatigue to predict cement sheath integrity[R].SPE 175231,2015.
[13] 李早元,郭小阳,罗发强,等.油井水泥环降脆增韧作用机理研究[J].石油学报,2008,29(3):438-441. LI Zaoyuan,GUO Xiaoyang,LUO Faqiang,et al.Research on mechanism of increasing flexibility and decreasing brittleness of cement sheath in oil well[J].Acta Petrolei Sinica,2008,29(3):438-441. [14] 张峰,姚晓.油井水泥增韧剂BF的室内研究[J].钻井液与完井液,2003,20(1):8-10. ZHANG Feng,YAO Xiao.Study on oilwell cement toughness reinforcing agent BF[J].Drilling Fluid Completion Fluid,2003,20(1):8-10. [15] 姚晓,樊松林,吴叶成,等.油井水泥纤维增韧材料的研究与应用[J].西安石油大学学报(自然科学版),2005,20(2):39-42. YAO Xiao,FAN Songlin,WU Yecheng,et al.Study and application of toughness-enhancing fiber material for the cement used for well-cementing[J].Journal of Xi’an Shiyou University(Natural Science Edition),2005,20(2):39-42. [16] 徐海民,赵林,湛峰.一种新型胶乳水泥浆性能研究[J].断块油气田,2011,18(6):803-804,808. XU Haimin,ZhAO Lin,ZHAN Feng.Performance of a new latex cement slurry[J].Fault-Block Oil Gas Field,2011,18(5):803-804,808. [17] 宋梅梅,张瑞萍,杨振杰,等.玻璃纤维对硫铝酸钙水泥石力学性能的影响[J].石油钻采工艺,2017,39(2):192-196. SONG Meimei,ZHANG Ruiping,YANG Zhenjie,et al.Effect of fiber glass on mechanical property of sulphoaluminate set cement[J].Oil Drilling Production Technology,2017,39(2):192-196. [18] 谭春勤,刘伟,丁士东,等.SFP弹韧性水泥浆体系在页岩气井中的应用[J].石油钻探技术,2011,39(3):53-56. TAN Chunqin,LIU Wei,DING Shidong,et al.Application of SFP elasto-toughness slurry in shale gas well[J].Petroleum Drilling Techniques,2011,39(3):53-56. [19] 何聪,金解放,周学进,等.静载荷与循环冲击组合作用下岩石损伤本构模型研究[J].有色金属科学与工程,2016,7(4):114-120. HE Cong,JIN Jiefang,ZHOU Xuejin,et al.Damage constitutive model of rock subjected to coupled static loadings and cyclic impacts[J].Nonferrous Metals Science and Engineering,2016,7(4):114-120. -
期刊类型引用(3)
1. 魏文绮,蒋敏,罗烈强,金志轩. 一种自研水泥浆在榆37区块储气库井中的应用. 当代化工研究. 2025(03): 143-145 . 百度学术
2. 完颜祺琪,毛涛,王云,黄盛,李康,李早元,苏东华. 储气库固井水泥石渗透率演化规律. 断块油气田. 2024(05): 930-935 . 百度学术
3. 关志刚,邓宽海,吴彦先,林志伟,林元华. 非常规油气井压裂过程中水泥环完整性实验评价. 西南石油大学学报(自然科学版). 2023(04): 121-132 . 百度学术
其他类型引用(2)
计量
- 文章访问数: 4277
- HTML全文浏览量: 78
- PDF下载量: 3983
- 被引次数: 5