测试管柱的力学分析.doc

高温高压地层测试工程理论研究4 测试管柱的力学分析测试管柱在井筒中要受到各种外力的作用，如内外压力、重力、井壁的反力等的作用。这些作用力与温度共同作用在测试管柱上，造成管柱的变形，如拉伸变形和屈曲变形等，以及在测试管柱中产生内力，如轴向力、弯矩等。如果这些变形或内力过大，就可能对测试管柱产生损坏。在不同的操作中，这些外力是不同的。因而，各种工况所产生的内力也不尽相同。例如，下放测试管柱时，测试管柱受的外力为重力和完井液对管柱的浮力，上部则由钻机大钩吊着；在坐封时，大钩逐步加上钻压，即松弛力，使封隔器坐封；在开井时，测试管柱中有天然气流过，因而测试管柱内外压力会发生变化，此外，测试管柱的温度变化会使管柱伸长。因此，在分析时必须根据不同工况进行具体分析。管柱在受到外力作用时产生变形，根据不同的内力，变形有所不同。众所周知，当管柱的轴向力是受拉时，管柱只是伸长，而当管柱的轴向力是受压时，除了轴向缩短外，对于这种长细比很大的管柱，管柱还会产生屈曲变形。屈曲变形反过来又会影响内力。因此，对测试管柱在井筒中的力学分析有助于合理地设计测试管柱及其测试操作。在本章中，我们研究井眼中管柱的受力分析、受压部分的屈曲分析和测试管柱的强度分析。4.1 测试管柱各工况的受力分析在地层测试过程中，需要进行测试管柱的下放（简称为下钻）、用低比重流体替代测试管柱中的流体（简称为低替）、封隔器坐封（简称为坐封）、打开井口关井阀诱喷（简称为开井）、井下关井阀关井（简称为1关）、井口关井阀关井（简称为2关）、高比重泥浆循环压井（由井口油管将高比重泥浆压入，从环形空间流出；简称为循环）或高比重泥浆反循环压井（由井口环形空间将高比重泥浆压入，从油管流出；简称为反循环）和压裂与酸化（简称为高挤酸）等操作。在这些操作中，测试管柱受力是不一样的。下面我们根据不同工况分析测试管柱的受力情况。4.1.1 下钻完测试管柱在下放的过程中，井眼中存在有完井液。测试管柱此时受有重力、悬挂力和液体的作用力（浮力）。设井眼垂直，测试管柱全长为l，从井口沿井眼中心垂直向下建立坐标系，深度坐标为z，参见图41。重力均匀作用在测试管柱整个深度范围，微长dz的管柱所受的重力dQ为(41)式中：r为单位长度测试管柱的质量，kg/m；g为重力加速度，m/s2。悬挂力P直接作用在测试管柱的顶部。悬挂力的大小按下式计算(42)式中：g为完井液密度，kg/m3；L为测试管柱长度，m。V为测试管柱浸入井筒完井液的体积，m3。浮力与测试管柱的下放深度成正比。当下放到深度z = l时，测试管柱受的浮力Fb作用在测试管柱的底部，其大小为(43)式中：A为测试管柱底部的面积，m2。对于组合测试管柱，即由不同截面积组成的测试管柱，各段管柱分别计算该段管柱的浮力。假设组合管柱的某一分界面位于深度z = l1的地方，则应在此分界面上加上一个向上的集中力F1(44)式中：A1为组合管柱分界面以上段底部的面积，m2；A2为组合管柱分界面以下段顶部的面积，m2。4.1.2 低替进行低比重液体替换操作时，测试管柱仍然是由大钩悬挂，作用在测试管柱上的重力不变。悬挂力P由平衡方程计算出。低替时，从井口油管内将低比重液体压入井中，井中多余完井液从井口的环空处流出。低替载荷最大时是低替完成的时候。此时，在测试管柱内充满低比重液体，环空中则仍有完井液占据。测试管柱在井底的压力应等于环空在井底的压力。由于测试管柱中低替液比重小于完井液的比重，在测试管柱的井口必须加上压力p(45)式中：g 1为完井液的密度，kg/m3；g 2为低替液的密度，kg/m3；假设组合管柱的某一分界面位于深度z = l1的地方，则应在此分界面上加上一个向上的集中力F1(46)式中：，m2；和分别为测试管柱分界面上段的外半径和下段外半径，m。，m2；和分别为测试管柱分界面下段的外半径和上段外半径，m。4.1.3 坐封坐封是在低替之后进行的操作。坐封操作时，通过旋转测试管柱使封隔器动作，再将大钩往下放，封隔器坐封。大钩松弛力越大，封隔器坐封力也越大。作用在测试管柱上的载荷除了低替的载荷外，还有封隔器上作用有坐封力。在封隔器以下的测试管柱的重量已经通过封隔器传递到套管上，因此，对封隔器以上的测试管柱已经不再起作用。坐封后，封隔器成为固定支点，坐封以后各操作封隔器以上测试管柱的变形以此为基准计算。一般来说，封隔器坐封力大小与井口松弛力的大小并不相同，原因是大钩在下放的过程中，下部的管柱因受压而进入螺旋屈曲，测试管柱与井壁接触，一部分力被摩擦力抵消。在计算封隔器坐封力时，必须考虑测试管柱的螺旋屈曲。封隔器的坐封力与井口松弛高度的距离的关系在下节中导出。4.1.4 射孔坐封之后，通过油管加压或通过环空加压进行射孔，建立起油层与测试管柱之间的通道，诱导油气由测试管柱喷出。射孔操作与坐封操作不同的是，在射孔时，在井口测试管柱中或在环空中施加井口压力。该压力的大小必须能使射孔枪动作。4.1.5 开井射孔之后，建立起油气流的通道。储油层的油气流在地层压力的驱动下流入封隔器下套管口袋中，并经由测试管柱流出。如第三章所述，流体在流过测试管柱时，压力和温度要发生变化，变化的规律由一组非线性微分方程描述。开井后在测试管柱上，外部压力仍然维持坐封时的状态，测试管柱的内部压力和测试管柱的温度需要用数值方法求出（见第三章）。测试管柱内流体的压力和温度变化使得测试管柱的变形产生变化。4.1.6 井下关井阀关井开井以后第一次关井一般采用井下关井阀关井。通过在井口的环空加压使得井下关井阀动作，从而达到关井的目的。关井后，地层油气流通过测试管柱流出井口的途径被阻断，测试管柱中的流体密度仍然保持开井时的密度，井口测试管柱压力减为零。由于流量减为零，井底的流压逐渐增加，作用在关井阀上的压力将使封隔器受很大的向上的力。4.1.7井口关井阀关井当最后关井时，采用井口关井阀进行关井操作。关井后，井中测试管柱中的流体密度仍然保持开井时的密度。由于流量为零，井口的压力进一步升高，直到井底的流压恢复到地层压力。4.2 考虑摩擦井眼中杆柱屈曲方程的建立4.2.1 无摩擦井眼中管柱的屈曲方程假设测试管柱和井眼内壁摩擦力可以忽略，管柱是均匀的（无接头），井眼垂直，管柱无初始弯曲内力。设井眼中心线向下为x轴，y轴、z轴如图43所示，井眼的半径为，杆柱的半径为，杆柱受压发生螺旋屈曲后，其轴线上任一点必在以井眼中心为轴线，半径为的圆柱面上，圆柱面的矢径为(47)微元体作用有内力，内力矩及分布外力(48)其中：q为单位长度杆柱重量(N/m)；N为井壁对杆柱的法向支反力(N/m)。微元体的静力平衡方程为，(49)即，(410)(411)(412)(413)式(410)中：m为管柱与井筒内壁之间的滑动摩擦系数。井眼的约束使屈曲变形仍在弹性小变形范围内，梁假设仍起作用，因此(414)(415)上二式中：E为管柱材料的弹性模量，MPa；J为管柱的截面惯性模量，m4。根据方程(410)至(415)，用Maple软件，可导出杆柱在垂直井眼中的屈曲方程， Maple程序如下：mz:=-EJ*r*diff(sin(th(x),x$2);Fy:=-diff(mz,x)+Fx(x)*r*diff(sin(th(x),x);my:=-EJ*r*diff(cos(th(x),x$2);Fz:=diff(my,x)-Fx(x)*r*diff(cos(th(x),x);N:=diff(Fz,x)*cos(th(x)-diff(Fy,x)*sin(th(x);simplify(N);dFz:=diff(Fz,x);eq:=dFz-N*cos(th(x);simplify(eq);运行后，得到垂直井眼中的屈曲方程：(416)反力N的表达式为(417)以及弯矩的表达式(418)在方程(416)中，如果假设轴向力为常量，则方程左边的第二项为零，屈曲方程可简化为：(419)4.2.2 有摩擦屈曲方程的解根据方程(419)，Lubinski等得出了井眼中管柱的等重、与井壁无摩擦的螺旋屈曲解1。(420)式中：F0为管柱的轴向载荷，N；L为管柱的长度，m。摩擦力在井眼中杆柱的屈曲中起着重要作用，为了在计算中考虑摩擦力，我们进行以下的近似计算：设杆柱与井眼壁的接触力仍可用Lubinski解，即(421)设杆柱的中性点为坐标原点，向下为正，现考虑几种不同的加载工况。4.2.2.1 上部加载在由上部加载情况下，可以得到如下方程：(422)设，上述方程化为(423)分离变量，得(424)两边积分，得出(425)式中，C为积分常数，由边界条件确定。当时，即，解出y并代入，得(426)当时，即，解出y并代入，得(427)相应于时，轴向力引起的缩短为(428)相应于时，轴向力引起的缩短为(429)相应于时，屈曲引起的缩短为(430)管柱刚下到垂直井眼中时，整个管柱处于拉伸状态，管柱与井壁没有摩擦。当封隔器坐封后，如果给管柱加一定的钻压，管柱的一部分会因为受压而进入螺旋屈曲，此时，。管柱开始受压点为中性点LN，该点的计算公式为：(431)式中：LT为管柱的全长；WB为加的钻压。作为例子，假定套管内径Dci = 127mm，油管外径Dto = 88.9mm，油管内径Dti = 70.2mm，油管全长为4000m，井中的液体比重为1；在地面下放油管坐封并加上钻压，钻压为400KN，中性点以下的管柱的轴向压力如图44所示。随着摩擦系数的增加，能加到封隔器上的载荷越小。摩擦系数较大时，距离中性点较远的受压段油管的轴向力随深度变化不大；摩擦系数减小，加到封隔器上的载荷增大。当摩擦系数为零时，井口的松弛力全部由封隔器承担，这时，管柱的轴向力为一直线。图44 钻压状态下不同摩擦系数下的轴向载荷现考虑另一种情况，如果在井口管柱加大松弛力直到的悬挂力为负，即井口管柱不是悬挂而是在油管的顶部施加向下的载荷，此时，。假设油管顶部加压为200KN。油管轴向力如图45所示，可以看出，底部的轴向力与顶部载荷关系不大，只与摩擦系数有关；顶部的轴向力变化很快，顶部的压力很难传递到井中。摩擦系数越大，顶部压力传递的衰减越厉害。图45 油管顶部加压下不同摩擦系数下的轴向载荷4.2.2.2 底部加载在由底部加载情况下，可以得到如下方程：(432)设，上述方程化为(433)分离变量，得(434)两边积分，得出(435)轴向力引起的缩短为(436)屈曲引起的缩短为(437)如果从管柱底部向上加载200KN，油管轴向力如图46所示，可以看出，摩擦系数越大，中性点越深。图46 油管底部向上加载时不同摩擦系数下的轴向载荷如果从管柱底部向上加载分别为200KN、400KN、600KN、800KN；摩擦系数为0.2；油管轴向力如图47所示。图47 油管底部向上加不同载荷时的轴向载荷油管底部加200KN时，中性点深度为3500m，当加到800KN时，中性点深度仅为3400m。这表明，底部加压时，影响的是距离底部不远的地方。图 某井的套管程序从工程实际来看，底部加载类似于井下封隔器滑动，井底压力直接作用在管柱底部，巨大的压力可能使井底几百米的管柱出现永久变形。4.2.2.3 多段井筒和多段管柱轴向力计算完井与地层测试中，往往管柱是复合管柱，即由不同尺寸的油管组合的管串，油层套管也不是直接接到井口，而是用套管挂挂接在前一次下的套管上（右图），这样就形成了多段井筒和多段管串的情况。计算轴向力时，必须将井段分成很多段，每一段具有同一套管内径、同一油管尺寸。在每一段中就可以应用上面导出的公式计算轴向力。在两段连接处考虑连接条件。4.3 屈曲方程的精确解4.3.1 无因次化由于方程(419)是一个非线性的常微分方程，一般来说，可能没有解，也可能有一个解，也可能有许多解，2002年，Mitchell得到了方程(419)另外一些解析解2。其中一个解具有一定的物理意义，我们下面导出该解。为了求解这个方程，将方程作无因次处理，设无因次长度(413)方程(419)化为(414)反力(417)化为(415)设无因次反力为(416)类似的，设无因次弯矩为(417)4.3.2 精确解 2设，方程(414)变为(418)式中C为积分常数，方程(416)简化为(419)用雅可比椭圆函数(420)其中，和均为常数，代入方程(418)，得，(421)设，得解为(422)积分上式，得(423)式中的dn和cn均为雅可比椭圆函数。4.3.3 精确解的分析4.3.3.1 变形形状为了求出精确解的形状，假设，r = 0.7，管柱外半径r1 = 0.7；用Matlab计算并绘出杆柱在井眼中的立体图，程序如下：clear;M=0.99999;U1=0:0.5:30;r=0.7;r1=0.5;nd=40;U=U1/sqrt(1+M);SN,CN,DN = ellipj(U,M);Th=log(DN-sqrt(M)*CN);X,Y,Z = cylinder(r+r1,50);Z=Z*30/1.4;surf(X,Y,Z);alpha(.5);shading interp;axis off; m,n=size(U1);for i=1:n for j=1:nd x(i-1)*nd+j)=r*cos(Th(i)+r1*cos(2*j*pi/nd); y(i-1)*nd+j)=r*sin(Th(i)+r1*sin(2*j*pi/nd); z(i-1)*nd+j)=U(i); if(i=n) if(j = nd) f(i-1)*nd+j,:)=(i-1)*nd+j (i-1)*nd+j+1 i*nd+j+1 i*nd+j; else f(i-1)*nd+j,:)=(i-1)*nd+j (i-1)*nd+1 i*nd+1 i*nd+j; end end endendvt = x y z;hold on;patch(Faces,f,Vertices,vt,FaceVertexCData,0.5,FaceColor,flat);立体图中可以看出，与Lubinski解不同的是：Lubinski解在整个屈曲管柱长度都是螺旋形状，雅可比椭圆函数的解是周期性的，在井中左旋和右旋交替，在图中指出了由左旋改变成右旋的交界点，左右旋的长度与常数M的大小有关，M越大，螺旋圈数越多，但其周期随参数M的增加而增大，当M趋于1时，精确解趋于Lubinski解。因此，Lubinski解是精确解的特殊形式。4.3.3.2反力计算根据(422)式，结合(419)式，可用Maple导出反力，Maple程序如下：u:=k/sqrt(1+k2)*JacobiSN(z/sqrt(1+k2),k);N:=u2+4*u*diff(u,z$2)+3*(diff(u,z)2-u4;simplify(N);运算结果为：(424)当时，无因次反力如图49所示。反力在大部分区域为常值0.25，与Lubinski解一致。反力在一小部分区域为负值，这与实际假设是不符的，因此在这一部分区域，解是不正确。在另一小部分区域中，反力的值很大，约为0.75，是平均反力的三倍，这个区域是杆柱螺旋角变化为0的点的附近。2、杆柱缩短计算由于屈曲，可引起杆柱的缩短，在屈曲段的缩短量由下式确定：(425)与值有关，图3是不同值的缩短值与的关系，随值的增加，屈曲缩短值增加，Lubinski解的缩短值最大。这意味着当较小时，管柱比较硬，变形比较小，在计算变形时必须考虑这种影响。图49 无因次反力图图4-10 不同k值时管柱的屈曲引起的缩短参考文献1 Lubinski, A., Althouse, W.S.,