套管摩阻计算刚杆模型.doc

套管摩阻计算刚杆模型 （1）受力分析 从静力学分析入手，在井眼轴线坐标系上任取一弧长为ds的微元体AB，并对其进行受力分析。以A点为始点，其轴线坐标为s，B点为终点，其轴线坐标为s+ds，其单元体的受力如图2.2所示。 微元段套管柱受力分析 （2）各力的表达式 轴线坐标s处（A点）的集中力（即A点的内力）为： (1) 式中，，——套管轴向、主、次法线方向上的力，N； A点的内力矩为为： (2) 式中,,——套管轴向、主、次法线方向上的力矩，； B点的内力为： (3) B点的内力矩为： (4) 微元段ds上的均布接触力为： (5) 式中N——井壁对管柱的接触压力，； 、——主法线和副法线方向的均布接触压力，N； ——轴向摩擦系数，无因次量； “”代表上提、下放套管，上提套管时取“+”，下放套管时取“—”（下同）。 单位长度套管柱浮重为： (6) 式中 ——浮力系数，即 (7) ——钻井液密度； ——套管柱材料密度； ——套管单位长度自重，； ——套管重力的方向； （3）平衡方程 由微元段ds的受力平衡条件，即： (8) （4）推理过程 将(1)、（2）、（3）、（4）、（5）、（6）式代入（8）式得： （10） 略去微量的乘积可得： （11） 整理得： （12） 根据Frenetic-Serret公式： （13） 式中：K——微元段井眼曲线曲率； ——微元段井眼曲线挠率； 则有： （14） 将（14）代入（12）得： （15） 将力分别向切线方向、主法线和副法线投影可得到： （16） 另外，由自然曲线坐标的单元向量与直角坐标系之间的关系为： （17） 式中：——平均井斜角，； ——平均方位角，； 、和、分别为微元段两端的井斜角和方位角。 因此可以得到： （18） 将（18）代入（16）得： （19） 式??? ——井斜角，rad； ——方位角，rad； ——井斜变化率，rad/m； ——方位变化率，rad/m； 其中：； ； ； 由微元段ds 的力矩平衡方程（李子丰.油气井管柱动力学基本方程及应用【J】石油学报，1999.20（3）：87-90）整理得： （20） 其中：m——单位长度分布外力矩； m= （21） 由（1）、（2）、（4）和（20）整理得： （22） 将力矩向切线方向、主法线和副法线投影可得： （23） 其中：； 由于管柱的变形受到井壁的限制，故其变形仍属于弹性小变形范围，设管柱的抗弯刚度为EI，抗扭刚度为GJ，忽略剪切力的影响则本构方程为： （24） （李子丰.油气井管柱动力学基本方程及应用【J】石油学报，1999.20（3）：87-90） 推理： （25） 由（24）（25）可知： （26） 将（26）代入（23）可得： （27） （5）结果 将（27）代入（19）并整理可得全刚度管柱摩阻计算公式： （29） 其中： ；