自噴井井底壓力分布

1、重慶科技學院課程設(shè)計報告院（系）：石油與天然氣工程學院_ 專業(yè)班級：石油普 09-4_09-4_學生姓名：_張鵬_學號:_20094400462009440046_設(shè)計地點（單位）重慶科技學院_設(shè)計題目: :_ _自噴井井底壓力分布完成日期：20122012 年 7 7 月 5 5 日指導教師評語：_成績（五級記分制）：_指導教師（簽字）:_ _一2自噴井井底壓力分布張鵬（重慶科技學院，石油工程09 級，學號 2009440046）刖言：2012 年 6 月 25 日至 2012 年 7 月 5 日在徐春碧老師的指導下，我們進行了為期兩周的采氣工程課程設(shè)計， 主要設(shè)計內(nèi)容為氣水同產(chǎn)井井底壓力分

2、布情況，在老師的指導及整個團隊的合作下，我們完成了最初的尋找資料、尋求思路到編程、寫報告等一系列的設(shè)計。本報告利用 Beggs-Brill 方法和 Duns-Ros 方法求解井底壓力。主要目的是利用采油、采氣工 程和多項管流的知識來解決現(xiàn)實生產(chǎn)中某一油、氣井站具體問題，提高學生解決實際問題的 的能力。設(shè)計題目3已知水產(chǎn)量 為 38m /d,GLR=200,油壓為 1.8Mpa（表壓），油管尺寸為 62mm（21/ 2），1、計算方法1）、持液率H|（）的相關(guān)規(guī)律2）、混合速度Vm及n的相關(guān)計算油層中部深度為 1800m,w1-05,0.65,e=0.016.0.8mPa?s，天然氣視臨界溫度為

3、 225k，視臨界壓力為4.7Mpa,井口溫度為 20C（ 293k），井筒內(nèi)溫度變化為3.5 C/100m，表面張力為0.06N/M,求井底壓力分布。（用兩種方法做）二、設(shè)計方法Beggs-Brill法基本壓力公式：P106mg sinfmVm2n/(2D)mVmVsg/ P 106|H|(0)g(1 H)H|(?H|(0)1 C sin 1.81 .sin31.8C (1|)|n(de fN|vNFR9)VmVslVsgQL/A Qg/A3fnexp(s)In y0.0523 3.18ln y 0.8725(ln y)20.01853(1 n y)4gP1/ZRT12、根據(jù)各參數(shù)的取值范圍

4、，可以確定管內(nèi)液體的流態(tài)，具體的標準如表1)、確定管內(nèi)液體的流態(tài)分層L40.01 andNFR=0.01 andNFR=0.01 andL2NFRL1間隔0.01=i0.4 andL2NFR=0.4 andL3NFR=L4分散l=L|l=0.4 andNFRL42)、判斷流形的基本公式表1d、e、f、g的確定流型上/下坡defg分離流上坡0.011-3.768：3.539-1.614 I間歇流上坡2.960.305：-0.44730.0978分散流上坡0000各種流型下坡4.7-0.36920.12440.5053、算法框圖已知參數(shù)P猜測參數(shù)F2(0)64、設(shè)計程序源代碼（見附頁程序一）三、設(shè)

5、計方法二Duns-Ros法1、計算方法1）、壓縮因子的計算Z = 1 +(0.3051 -1.0467 / Tr -0.5783 / Tr3)r+(0.5353-0.6123 / Tr + 0.6816 / Tr3)Tr = T/TcPrP/Pcr0.27 pr/ ZTr2)、氣體密度的計算g3484.4gP(i)/(T z)Vsg Qsc (3.447 104Z T/P(i)/SVsl Qw / S3) 、液體粘度計算mexp1.0q0.01479.8( Z i 0.035T T)/2 32 1.9821051.8( z i 0.035T T)/21034) 、平均溫度的計算T二(z i 0

6、.035 + 2T) / 2 + 273.154)、四個基本準數(shù)的計算NVIvsl4zgNvg = Vsg4|/ gND= Vsg4I/gNI=i4l/ g2、流動形態(tài)的判斷由課本氣液兩項管流P36頁流動形態(tài)分布圖可知：I 區(qū)：液相是連續(xù)的，包括泡狀流、彈狀流和部分的沫狀流；II 區(qū)：液相和氣相是連續(xù)的，包括段塞流和沫狀流的其余部分;III 區(qū)：氣相是連續(xù)的，指霧狀流。其中71)第 I 區(qū)發(fā)生在NvgL1LQNVI的范圍內(nèi)。此時Ns2F1F2NVIF3答1 Nvi其中FI,F2, F3和 F4作為 N 的函數(shù)在課本 P37圖 3-5 圖給出。L!和L2為若斯因數(shù)作為ND的函數(shù)在圖 3-6 中給

7、出。2) 第 II 區(qū)發(fā)生在L1LzNviNvg50 36Nvi的范圍內(nèi)。此時其中Fs和 F6作為 N1的函數(shù)在課本 P38圖 3-7 圖給出。3) 2)第 III 區(qū)0 75發(fā)生在Nvg75 84Nvl的范圍內(nèi)。此時4、算法框圖已知參數(shù)P猜測參數(shù)F2(0)NS(1+F5)卄(1 F7NQ3、NSP的計算v Ns(vPf(i)Ph(i)=(P(i)(0.982vgF6N(1+F5)(T F7Nvl)2(9.81/Vsl Vsg) !FrPh(i)Vsl2+ (1x 0.25l)(v Vsl Vsg)2I(1 Vsg/Vsl)-)I)9.8pf(i)/(1 (1Vsl1/24 v Vsg) )/(2 v)z/dgVsg) Vsg)/P(i)/1000000)/10685、設(shè)計程序源代碼（見附頁程序二） 四、程序代碼1、Beggs-Brill法源代碼2、Duns-Ros法源代碼參考文獻：1、徐春碧等.氣液兩相管流理論與計算，重慶科技學院，2006.12,陳家朗石油