【正文】 圖46 單向應(yīng)力安全系數(shù)（1）以井深4500米處作為校核點進行校核。 （424） （425）（2）截面處接觸支反力： （426）（3）截面處的彎矩為： （427）（4）彎矩產(chǎn)生的彎曲正應(yīng)力： （428）（5）接觸支反力產(chǎn)生的剪切應(yīng)力： （429）（6）軸向力產(chǎn)生的軸向應(yīng)力：（429）（7）內(nèi)壓、外壓產(chǎn)生的徑向應(yīng)力： （430）（8）內(nèi)壓、外壓產(chǎn)生的軸向應(yīng)力：（431）（9）根據(jù)第四強度理論，相當應(yīng)力：（432）井深m剪切應(yīng)力Mpa彎曲正應(yīng)力Mpa軸向應(yīng)力Mpa徑向應(yīng)力Mpa周向應(yīng)力Mpa相當應(yīng)力Mpa0005000010000015000020000025000030000035004000450030表49 不同井深處管柱相當應(yīng)力圖47 第四強度理論校核安全系數(shù)，滿足安全要求。第五章 結(jié)論（1）整理了計算采油過程中流體參數(shù)的計算方法，已知井底，采用壓力按壓力增量迭代的方法得到井筒內(nèi)溫度、壓力隨井深的非線性分布曲線。將計算過程編程。（2）分析了溫度，壓力變化引起的四種基本效應(yīng)對深井帶封隔器的油管柱變形和載荷的影響。給出了相應(yīng)的計算公式。綜合考慮管柱內(nèi)外所對應(yīng)的“虛構(gòu)力”，用等效軸力來對管柱進行分析計算。建立了深井測試管柱力學(xué)模型。（3）給出了深井油井管柱的載荷、變形、應(yīng)力計算的簡化公式，分析了深井油管柱的應(yīng)力和變形。分別從單向應(yīng)力和第四強度了理論的角度對深井油管柱進行了強度校核。通過計算管柱的安全系數(shù)來優(yōu)化管柱和施工參數(shù)。提高深井生產(chǎn)作業(yè)的安全性和成功率（4）深井管柱的受力分析和計算是設(shè)計和校核管柱強度、選擇作業(yè)工具、確定施工參數(shù)的前提和基礎(chǔ)，是保證深井生產(chǎn)作業(yè)的安全和成功的必做的工作。（5）深井管柱力學(xué)的分析和計算要綜合考慮井身結(jié)構(gòu)、管柱組合、管柱載荷、溫度壓力變化以及井況等因素。計算帶封隔器的管柱的載荷和變形要綜合考慮管柱屈曲和封隔器的約束。（6）對實際帶封隔器的自噴井進行載荷、變形計算和強度校核。證明了模型的可靠和實用性。參考文獻參考文獻[1]劉延鑫,王旱祥,侯乃賀,[J].鉆采工藝,2012(4):1186+84.[2]竇益華,[J].鉆采工藝,2007(5):634+31+40.[3]曾志軍,胡衛(wèi)東,劉竟成,[J].天然氣工業(yè),2010(2):9294+152.[4]王安義,龍堯,[J].西部探礦工程,2010(8):5356.[5]李子豐,蔣恕,[J].石油機械,2002(12):3033+60.[6]楊新克,余杰,[J].新疆石油天然氣,2010(4):894+90.[7]陳明亮,王越支,[J].鉆采工藝,1998(6):4547+50+6.[8]蘇華,張學(xué)鴻,：鉆柱靜力學(xué)[J].大慶石油學(xué)院學(xué)報,1994(1):4352.[9]李子豐,蔡雨田,李冬梅,[J].石油學(xué)報,2011(4):163170.[10]張袁輝,蔣海,[J].海洋石油,2006(3):9698.[11][J].油氣井測試,2005(4):6768+82.[12]畢博,朱衛(wèi)新,[J].科學(xué)技術(shù)與工程,2010(10):1420.[13]李欽道,謝光平,——封隔器管柱受力分析系統(tǒng)討論之五[J].鉆采工藝,2002(2):6165+6.[14] Mitchell R F. Comprehensive analysis of buckling with friction[J]. SPE Drilling amp。 Completion, 1996, 11(3): 178184..[15]王祖文,林玉璽,[J].大慶石油地質(zhì)與開發(fā),2007(6):105109.[16]唐海雄,張俊斌,汪順文,[J].石油機械,2010(5):9193+98.[17]劉琦,蔣建勛,石慶,[J].斷塊油氣田,2007(1):6466+96.[18]杜現(xiàn)飛,王海文,王帥,[J].石油礦場機械,2008(8):3338.[19]居迎軍,劉曉光,高永亮,[J].中國石油和化工,2008(10):5558.[20]高寶奎,[J].石油大學(xué)學(xué)報(自然科學(xué)版),2002(6):5456+8.[21] Hammerlindl D J. Movement, forces, and stresses associated with bination tubing strings sealed in packers[J]. Journal of Petroleum Technology, 1977, 29(2): 195208.[22]高寶奎,[J].天然氣工業(yè),2002(6):6769+8.[23] Mitchell R F. Buckling behavior of well tubing: the packer effect[J]. Old SPE Journal, 1982, 22(5): 616624.[24]李欽道,謝光平,——封隔器管柱受力分析系統(tǒng)討論(之四)[J].鉆采工藝,2002(1):6973+7.[25] 于建梅. 深井測試管柱力學(xué)研究[D]. 中國石油大學(xué), 2007.附錄井筒多相流程序代碼：Private Sub Command1_Click()Dim G, qt, Kp, Fww, Bata, pw, po, fw, trr, too, tpj, ppj, pwf, ph, H, L, yapi, mo, xg, yng, Rs, F, Bo, pl, z, y, x, uod, uo, A, B, uw, ulDim sog, swg, sl, Tc, Tr, Pc, Pr, Prr, zz, pg, vo, vw, vg, Qvg, Qvl, qg, Ql, Qvt, qm, pll, vgg, vm, LB, LB1, Hgfw = rgg = roo = pwf = ph = 0dd = too = 15H = 4500trr = * H / 100qt = 25G = qt * 1000 / 24Kp = 1 / ( + * Exp(G / 1000))39。注意po和pw的取值Fww = 1 * fw / (1 * fw + roo * (1 fw))Bata = 2 * * Kp / G / (1 + Fww)L = 0tp = too + (trr too) / Bata / H * (Bata * L + 1 Exp(Bata * (H L))) = tpL = 500ttrr = too + (trr too) / Bata / H * (Bata * L + 1 Exp(Bata * (H L))) = ttrrFor iii = 1 To 5 Step 1tpj = (tp + ttrr) / 2 = tpjppj = (pwf + ph) / 2 = ppj39。Print t=。 tpj39。Print p=。 ppjppj = ppj * 10 ^ 6yapi = / roo 39。If yapi 15 Then39。AA = * yapi * ( * tpj + 32)39。Rs = * rgg * ( * ppj * 10 ^ AA) ^ 39。End If39。If yapi Then39。mo = 10 ^ ( * Log(1346 / (yapi )))39。Elsemo = ( yapi) / 39。End Ifxg = * ppj * rgg / 10 ^ 5 / (tpj + )39。Print xg=。 xg39。If xg Thenyng = * Log(xg / )39。ElseIf xg Then39。yng = * Log(xg / )39。Else 39。yng = * Log(xg / )39。End If39。Print yng=。 yng39。If yapi 15 ThenRs = 23650 * roo / mo * yng / (1 yng) Rs = (H L) / 500 * 1039。End If39。Print Rs=。 RsF = * Rs * (rgg / roo) ^ + * tpj + 40Bo = + * F ^ 39。Print Bo=。 Bopo = 1000 * (roo + * 10 ^ 3 * Rs * rgg) / Bo39。Print po=。 popl = po * (1 fw) + 1000 * fw = pl39。Print 混合物平均密度=。 plz = * yapiy = 10 ^ zx = y * (32 + * tpj) ^ uod = (10 ^ x 1) / 1000A = * ( * Rs + 100) ^ B = * ( * Rs + 150) ^ uo = A * (1000 * uod) ^ B / 100039。Print uo=。 uouw = Exp( * 10 ^ 2 * (32 + * tpj) + * 10 ^ 5 * (32 + * tpj) ^ 2) / 100039。Print uw。 uwul = uo * (1 fw) + uw * fw = ul39。Print ul。 ulsog = ( * ( * tpj + 32) * yapi) * Exp( * 10 ^ 7 * ppj) / 100039。Print sog=。 sogs2 = 76 * Exp( * 10 ^ 7 * ppj) / 1000s1 = ( * 10 ^ 7 * ppj) / 1000swg = ((248 * tpj) / 206 * (s2 s1) + s1)39。Print swg=。 swgsl = sog * (1 fw) + swg * fw39。Print sl=。 slTc = + * rggIf rgg = ThenPc = 10 ^ 6 * ( * rgg)ElsePc = 10 ^ 6 * ( * rgg)End IfTr = ( + tpj) / TcPr = ppj / PcPrr = * Pr / Trzz = 1 + ( / / ^ 3) * Prr + ( / ) * Prr ^ 2 + / * Prr ^ 239。Print zz=。 zzpg = * 10 ^ 3 * rgg * ppj / zz / (tpj + )39。Print pg=。 pgug = * 10 ^ 239。計算壓力梯度vo = Bovw = fw / (1 fw)39。生產(chǎn)汽油比55vg = zz * * 10 ^ 6 * ( + tpj) / ppj / ( + 20) * (100 Rs)39。是Fww還是fwQvg = vg * qt * (1 Fww) / 86400 / roo39。Print 氣相體積流量Qvg=。 QvgQvl = (vo + vw) * qt * (1 Fww) / 86400 / roo39。Print 液相體積流量Qvl=。 Qvlqg = Qvg * pg39。Print 氣相質(zhì)量流量qg=。 qgQl = (po + 1000 * vw + * rgg * Rs) * qt * (1 Fww) / 86400 / roo39。Print 液相質(zhì)量流量Ql=。 QlQvt = (vo + vw + vg) * qt * (1 Fww) / 86400 / roo39。Print 總體積流量Qvt=。 Qvtqm = qg + Ql39。Print 總質(zhì)量流量qm=。 qm39。液相平均密度pll = Ql / Qvl39。Print 液體密度（油水體積加權(quán)）pll=。 pllAp = * * dd ^ 2vgg = Qvg / Ap * (pll / / sl) ^ 39。Print vgg=。 vggvm = Qvt / Ap = vm39。Print 總的混合物流速vm=。 vmLB = * vm ^ 2 / dd39。Print LB=。 LBLB1 = Qvg / Qvt39。Print qg/qt。 LB1LS = 50 + 36 * vgg * Ql / qg39。Print Ls。 LSLM = 75 + 84 * (vgg * Ql / qg) ^ 39。Print Lm。 LMIf LB1 LB Thenaaa = 泡流ElseIf LB1 LB And vgg LS T