【正文】 圖46 單向應(yīng)力安全系數(shù)（1）以井深4500米處作為校核點(diǎn)進(jìn)行校核。 （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ù)第四強(qiáng)度理論，相當(dāng)應(yīng)力：（432）井深m剪切應(yīng)力Mpa彎曲正應(yīng)力Mpa軸向應(yīng)力Mpa徑向應(yīng)力Mpa周向應(yīng)力Mpa相當(dāng)應(yīng)力Mpa0005000010000015000020000025000030000035004000450030表49 不同井深處管柱相當(dāng)應(yīng)力圖47 第四強(qiáng)度理論校核安全系數(shù)，滿足安全要求。第五章 結(jié)論（1）整理了計(jì)算采油過(guò)程中流體參數(shù)的計(jì)算方法，已知井底，采用壓力按壓力增量迭代的方法得到井筒內(nèi)溫度、壓力隨井深的非線性分布曲線。將計(jì)算過(guò)程編程。（2）分析了溫度，壓力變化引起的四種基本效應(yīng)對(duì)深井帶封隔器的油管柱變形和載荷的影響。給出了相應(yīng)的計(jì)算公式。綜合考慮管柱內(nèi)外所對(duì)應(yīng)的“虛構(gòu)力”，用等效軸力來(lái)對(duì)管柱進(jìn)行分析計(jì)算。建立了深井測(cè)試管柱力學(xué)模型。（3）給出了深井油井管柱的載荷、變形、應(yīng)力計(jì)算的簡(jiǎn)化公式，分析了深井油管柱的應(yīng)力和變形。分別從單向應(yīng)力和第四強(qiáng)度了理論的角度對(duì)深井油管柱進(jìn)行了強(qiáng)度校核。通過(guò)計(jì)算管柱的安全系數(shù)來(lái)優(yōu)化管柱和施工參數(shù)。提高深井生產(chǎn)作業(yè)的安全性和成功率（4）深井管柱的受力分析和計(jì)算是設(shè)計(jì)和校核管柱強(qiáng)度、選擇作業(yè)工具、確定施工參數(shù)的前提和基礎(chǔ)，是保證深井生產(chǎn)作業(yè)的安全和成功的必做的工作。（5）深井管柱力學(xué)的分析和計(jì)算要綜合考慮井身結(jié)構(gòu)、管柱組合、管柱載荷、溫度壓力變化以及井況等因素。計(jì)算帶封隔器的管柱的載荷和變形要綜合考慮管柱屈曲和封隔器的約束。（6）對(duì)實(shí)際帶封隔器的自噴井進(jìn)行載荷、變形計(jì)算和強(qiáng)度校核。證明了模型的可靠和實(shí)用性。參考文獻(xiàn)參考文獻(xiàn)[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].石油機(jī)械,2002(12):3033+60.[6]楊新克,余杰,[J].新疆石油天然氣,2010(4):894+90.[7]陳明亮,王越支,[J].鉆采工藝,1998(6):4547+50+6.[8]蘇華,張學(xué)鴻,：鉆柱靜力學(xué)[J].大慶石油學(xué)院學(xué)報(bào),1994(1):4352.[9]李子豐,蔡雨田,李冬梅,[J].石油學(xué)報(bào),2011(4):163170.[10]張?jiān)x,蔣海,[J].海洋石油,2006(3):9698.[11][J].油氣井測(cè)試,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].石油機(jī)械,2010(5):9193+98.[17]劉琦,蔣建勛,石慶,[J].斷塊油氣田,2007(1):6466+96.[18]杜現(xiàn)飛,王海文,王帥,[J].石油礦場(chǎng)機(jī)械,2008(8):3338.[19]居迎軍,劉曉光,高永亮,[J].中國(guó)石油和化工,2008(10):5558.[20]高寶奎,[J].石油大學(xué)學(xué)報(bào)(自然科學(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] 于建梅. 深井測(cè)試管柱力學(xué)研究[D]. 中國(guó)石油大學(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。計(jì)算壓力梯度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