[優(yōu)秀畢業(yè)設(shè)計(jì)精品]采油方式優(yōu)化設(shè)計(jì)

1、 目 錄第一章 前言1第一節(jié) 設(shè)計(jì)概要1第二節(jié) 采油優(yōu)化目標(biāo)和基本原則2 第三節(jié) 工藝設(shè)計(jì)步驟3第二章 油井基本情況3第一節(jié) 油井基礎(chǔ)數(shù)據(jù)3第二節(jié) 油井假設(shè)條件4第三章 采油工藝選擇過程4 第一節(jié) 各單元流體物性4第二節(jié) 流出動(dòng)態(tài)曲線11第三節(jié) 流入動(dòng)態(tài)曲線19第四節(jié) 協(xié)調(diào)點(diǎn)的求取及采油方式的選擇21第四章 所選采油方法的方案設(shè)計(jì)22第五章 方案設(shè)計(jì)結(jié)果展示23第一章 前言第一節(jié) 設(shè)計(jì)概要一、采油工程方案設(shè)計(jì)編制的目的和意義油田開發(fā)是一項(xiàng)龐大而復(fù)雜的系統(tǒng)工程，因此，在油田投入正式開發(fā)之前，必須編制油田開發(fā)總體方案，作為油田開發(fā)工作的指導(dǎo)性文件。采油工程方案是總體方案的重要組成部份和總體方案實(shí)施

2、的核心。二、采油工程方案設(shè)計(jì)在油田開發(fā)總體方案中的地位采油工程設(shè)計(jì)是依據(jù)油藏地質(zhì)研究成果和油藏工程方案，制定應(yīng)采取的一系列工程措施方案。其中也包括根據(jù)油藏地質(zhì)條件和當(dāng)前工藝技術(shù)水平對實(shí)現(xiàn)油藏工程方案的開發(fā)指標(biāo)(主要是產(chǎn)量、注入量及壓力保持水平)的可行性進(jìn)行評價(jià)。由于油田產(chǎn)能建成之后，油田開發(fā)指標(biāo)和原油生產(chǎn)計(jì)劃，必須通過采油工程措施的實(shí)施來實(shí)現(xiàn)。所以，采油工程方案在總體方案的實(shí)施中處于十分重要的位置。三、采油工程方案設(shè)計(jì)編制的原則和要求采油工程方案編制的原則有：1、符合本油田開發(fā)的總體布署和技術(shù)政策；2、適合油藏地質(zhì)和環(huán)境特點(diǎn)；3、采用先進(jìn)而適用的工藝技術(shù)，具有良好的可操作性；4、體現(xiàn)少投入多產(chǎn)

3、出的經(jīng)濟(jì)原則。采油工程方案編制的原要求有：1、充分應(yīng)用油藏地質(zhì)研究和油藏工程提供的基本資料，并作為主要的設(shè)計(jì)依據(jù)；2、重點(diǎn)論證本油田(設(shè)計(jì)對象)開發(fā)的主要問題、基本工藝和關(guān)鍵技術(shù)；3、結(jié)合油藏特點(diǎn)開展必要的室內(nèi)和現(xiàn)場工藝試驗(yàn)，并充分借鑒同類油田的經(jīng)驗(yàn)；4、采用先進(jìn)的理論和設(shè)計(jì)方法，進(jìn)行科學(xué)論證和方案的優(yōu)選；5、具有科學(xué)性、完整性、適用性、可操作性和經(jīng)濟(jì)性。只有滿足采油方式優(yōu)選及工藝設(shè)計(jì)的要求，才能最大程度的減少油田在管理方面和 工藝設(shè)備方面的投入，從而最大程度的增加油田的收入。采油方式優(yōu)選是每個(gè)油田在開發(fā)階段必須完成的方案設(shè)計(jì)，要求設(shè)計(jì)人員在充分掌握油田相關(guān)數(shù)據(jù)的基礎(chǔ)之上，設(shè)計(jì)油田開發(fā)各個(gè)階段

4、的采油優(yōu)選方案，供基層工作人員在采油過程中參考和應(yīng)用。四、采油工程方案設(shè)計(jì)編制的前期準(zhǔn)備采油工程方案設(shè)計(jì)編制的前期準(zhǔn)備有：1、早期介入，熟悉和掌握油藏特點(diǎn)，參與油田開發(fā)的總體布署工作；2、收集和分析油藏地質(zhì)和試油、試采及鉆井完井資料；3、結(jié)合油藏特點(diǎn)和總體布署的要求進(jìn)行技術(shù)調(diào)研；4、擬訂采油工程方案編制大綱，確定重點(diǎn)，提出試采工藝要求和方案編制過程中需要開展的研究課題；5、針對油田的某些特殊問題開展專題研究或現(xiàn)場工藝試驗(yàn)，為方案編制進(jìn)行必要的技術(shù)準(zhǔn)備。五、采油工程方案設(shè)計(jì)的基本構(gòu)成采油工程方案設(shè)計(jì)的基本構(gòu)成有：1、完井工程與投產(chǎn)2、注水工藝3、采油方式選擇4、油層改造5、配套工藝6、動(dòng)態(tài)監(jiān)測五

5、、采油工程方案設(shè)計(jì)的內(nèi)容課程設(shè)計(jì)內(nèi)容選擇的是采油方式優(yōu)選及工藝設(shè)計(jì)。采油方式選擇的基本內(nèi)容：確定舉升方式選擇的原則、依據(jù)及要求、油井產(chǎn)能預(yù)測及分析。包括預(yù)測不同含水階段的油井產(chǎn)液、產(chǎn)油指數(shù)；油井產(chǎn)能分布研究：按油井產(chǎn)能進(jìn)行油井分類，并確定不同含水階段各類井井?dāng)?shù)的變化；分析各類油井在油藏上的分布情況。采油方式是貫穿油田開發(fā)全過程的基本技術(shù)。因此，確定采油方式便成為采油工程方案的主要組成部分，而油井供液能力是采油方式選擇的主要依據(jù)。優(yōu)選實(shí)質(zhì)上是根據(jù)油藏工程及油田生產(chǎn)的地面條件，選擇能夠發(fā)揮油井潛在產(chǎn)能的經(jīng)濟(jì)有效的舉升手段。由于在油田開發(fā)過程中，油田生產(chǎn)動(dòng)態(tài)和油井生產(chǎn)條件的變化，以及單井生產(chǎn)能力的差

6、異及分布的變化，因此構(gòu)成了采油方式選擇具有系統(tǒng)性、預(yù)測性、綜合性的特點(diǎn)，既要考慮不同開發(fā)階段油井的供液能力變化，又要考慮單井生產(chǎn)能力的分布狀況及其變化。 目前，在現(xiàn)場應(yīng)用最廣泛的采油方式有：自噴、氣舉、電潛泵、有桿泵和水力泵。在各種采油方式總，自噴是管理最簡單，也是最為經(jīng)濟(jì)高效的舉升方式。凡是能夠自噴完成配產(chǎn)的井，均應(yīng)采用自噴采油。不能自噴或自噴產(chǎn)量達(dá)不到陪產(chǎn)要求的井才考慮機(jī)械采油。目前，應(yīng)用最為廣泛的機(jī)械采油方式是有桿泵，其次是氣舉、電潛泵和水力泵。由于油田投入開發(fā)以后，機(jī)采設(shè)備和管理維護(hù)費(fèi)用是主要投入，不同采油方式的投入產(chǎn)出比差別較大。因此必須在技術(shù)經(jīng)濟(jì)綜合評價(jià)的基礎(chǔ)上選擇采油方式。第二節(jié)

7、 采油優(yōu)化目標(biāo)和基本原則一、采油方式選擇的目標(biāo)采油方式擇的目標(biāo)是：少井、高產(chǎn)、經(jīng)濟(jì)、實(shí)用。因此，必須在充分研究油田開發(fā)階段動(dòng)態(tài)變化，分析最大潛能的基礎(chǔ)上，進(jìn)行綜合經(jīng)濟(jì)評價(jià)優(yōu)選采油方式。二、采油方式優(yōu)選的基本原則：采油方式優(yōu)選的基本原則：1、能夠充分發(fā)揮油井的生產(chǎn)能力，滿足開發(fā)方案規(guī)定的配產(chǎn)任務(wù)。2、所選擇的舉升設(shè)備工作效率較高。3、所選擇的采油方式可靠，對油井的生產(chǎn)狀況具有較強(qiáng)的適應(yīng)性。4、立足于減少井下作業(yè)工作量，選擇的采油方式油井維修與管理方便。5、適合油田野外工作環(huán)境和動(dòng)力供應(yīng)條件。6、所選擇的采油方式投入少、效益高。第三節(jié) 工藝設(shè)計(jì)步驟一、采油方式的選擇及工藝設(shè)計(jì)的步驟：1、按單井產(chǎn)量

8、高低將油井分類，選出代表井，根據(jù)各種機(jī)采方式的特點(diǎn)初選幾種機(jī)才方式作為預(yù)備案，準(zhǔn)備所需要的基本數(shù)據(jù)。2、用計(jì)算機(jī)模擬各種備選采油方式在各代表井的生產(chǎn)可能性，即求得各種采油方式在預(yù)測開發(fā)期間內(nèi)各生產(chǎn)階段的最大產(chǎn)油量。3、對各種可行的采油方式進(jìn)行經(jīng)濟(jì)分析對比，求出各種采油方式的相對經(jīng)濟(jì)效果。4、與其它影響采油方式適應(yīng)性的非量化因素一起進(jìn)行綜合分析、選出適用于本油田的最佳采油方式。第二章 油井基本情況第一節(jié) 油井基礎(chǔ)數(shù)據(jù)所選井為x井，井別是開發(fā)井，井型是垂直井。通過對該井鉆完井資料、試井資料、生產(chǎn)測試資料和實(shí)驗(yàn)室分析資料，得到了設(shè)計(jì)所需要的井身數(shù)據(jù)、井動(dòng)態(tài)數(shù)據(jù)、流體物性數(shù)據(jù)見表 2-1。一、 井身數(shù)

9、據(jù)油管尺寸為27/8（內(nèi)徑62.0mm），套管尺寸為85/8（內(nèi)徑198.8mm），油層中部深度為2000m。 二、井動(dòng)態(tài)數(shù)據(jù)產(chǎn)液量為60m3/d,對應(yīng)井口壓力為3.3mpa，地層壓力為22.0mpa，生產(chǎn)氣液比為100m3/m3,油井含水率為20%，井底溫度為70，井口溫度為14。 三、流體物性數(shù)據(jù)原油相對密度為0.85，地層水相對密度為1.02，天然氣相對密度為0.65，原油泡點(diǎn)壓力為12.0mpa，水的粘度為1mpas，溶解氣油比為100m3/m3。表 2-1 設(shè)計(jì)所需基本數(shù)據(jù)表參數(shù)開采方向 設(shè)計(jì)油壓3.3mpa地層靜壓22.0mpa原油泡點(diǎn)壓力12.0mpa含水率20%天然氣原油的界面

10、張力=0.045n/m油管內(nèi)徑62.0mm套管內(nèi)徑198.8mm產(chǎn)層中部2000m井口溫度15地層溫度70測試產(chǎn)液量60.0m3/d水的粘度1mpas氣油比100m3/m3原油相對密度0.85地層水相對密度1.02氣體相對密度0.65第二節(jié) 油井假設(shè)條件在基礎(chǔ)數(shù)據(jù)恒定的情況下，生產(chǎn)井中不同井?dāng)嗟臏囟群蛪毫s各不相同，從而使各段流體的物性參數(shù)不同。計(jì)中我們采用orkiszewski方法來計(jì)算井筒內(nèi)的油壓分布曲線。為了適當(dāng)簡化設(shè)計(jì)，使所測得的數(shù)據(jù)能夠用于理論計(jì)算，需要做如下假設(shè)：1、油管下到油層中部；2、井筒流體溫度分布呈線性分布；3、水中不溶解天然氣。第三章 采油工藝選擇過程采油工藝選擇過程是指

11、在自噴、氣舉、電潛泵、有桿泵、和水力泵，這五種油田最廣泛應(yīng)用的采油方式中選擇其中一種或幾種，以滿足油田在各個(gè)采油階段的采油需要。一般是以井底為解節(jié)點(diǎn)，通過求解此解節(jié)點(diǎn)的流入動(dòng)態(tài)ipr曲線和流出動(dòng)態(tài)曲線,把兩曲線繪制在同一坐標(biāo)圖中，在坐標(biāo)圖中求解協(xié)調(diào)點(diǎn)，根據(jù)協(xié)調(diào)點(diǎn)來選擇采油方式。在此次設(shè)計(jì)中，由于已知的是井口的油壓，因此需要從井口往井底進(jìn)行壓力梯度的計(jì)算，即以井底為解節(jié)點(diǎn)的流出動(dòng)態(tài)曲線計(jì)算。然后推出井底的壓力，以此為測試井底壓力，結(jié)合數(shù)據(jù)表中告訴的測試產(chǎn)量，計(jì)算出以井底為解節(jié)點(diǎn)的流入動(dòng)態(tài)ipr曲線。最后結(jié)合兩曲線求協(xié)調(diào)點(diǎn)，然后進(jìn)行采油方式的選擇。第一節(jié) 各單元流體物性一、原油粘度0=c表3-1

12、各段原油粘度深度 m平均溫度tcd原油粘度00-100289.525112.385-8.160623.5721322100-200292.275108.9891-8.120713.5231068200-300295.025105.7244-8.081163.4751782300-400297.775102.5849-8.041963.4283127400-500300.52599.56457-8.00313.3824778500-600303.27596.65805-7.964573.337642600-700306.02593.86009-7.926373.2937754700-800308

13、.77591.16577-7.888493.2508491800-900311.52588.57042-7.850933.2088351900-1000314.27586.06963-7.813683.16770691000-1100317.02583.65921-7.776743.12743851100-1200319.77581.33519-7.74013.08800511200-1300322.52579.09381-7.703763.0493831300-1400325.27576.9315-7.667713.01154891400-1500328.02574.84486-7.6319

14、52.97448071500-1600330.77572.83068-7.596472.9381571600-1700333.52570.8859-7.561282.90255721700-1800336.27569.00762-7.526352.86766131800-1900339.02567.19306-7.49172.833451900-2000341.77565.43961-7.457312.7999049二、原油溶解汽油比和體積系數(shù)rs=0.02124bo=式中:rs-m3/m3 p-mpa t-表3-2 各段原油溶解氣油比和體積系數(shù)深度h平均溫度t平均壓力prsbo0-10016

15、.3752.76514.042771.023251100-20019.1253.2116.600491.031078200-30021.8753.6819.328631.039423300-40024.6254.1822.255841.048353400-50027.3754.70525.347371.057802500-60030.1255.24528.534091.06762600-70032.8755.80531.843771.077869700-80035.6256.3835.239011.088463800-90038.3756.96538.681821.099314900-1000

16、41.1257.56542.20111.1104851000-110043.8758.17545.759241.1218831100-120046.6258.78549.284681.1333341200-130049.3759.452.807341.1449011300-140052.12510.02556.357791.1566521400-150054.87510.6659.933111.1685771500-160057.62511.29563.462821.1804971600-170060.37511.9366.944561.1923991700-180063.12512.5770

17、.410041.2043561800-190065.87513.2173.823661.2162721900-200068.62513.8577.18391.228138三、天然氣偏差因子式中：pr、tr擬對比壓力、溫度，p/ppc、t/tpc；ppc、tpc擬臨界壓力、溫度，mpa、k。表3-3 各段天然氣的偏差系數(shù)深度 hppctpc平均壓力p平均溫度tpprtpr偏差系數(shù)zg0-1004.603641206.27272.765289.5250.6006121.4036030.921827100-2004.603641206.27273.21292.2750.6972741.4169350

18、.91427200-3004.603641206.27273.68295.0250.7993671.4302670.907502300-4004.603641206.27274.18297.7750.9079771.4435990.901566400-5004.603641206.27274.705300.5251.0220171.456930.896708500-6004.603641206.27275.245303.2751.1393161.4702620.893197600-7004.603641206.27275.805306.0251.2609581.4835940.89104870

19、0-8004.603641206.27276.38308.7751.385861.4969260.890409800-9004.603641206.27276.965311.5251.5129331.5102580.891383900-10004.603641206.27277.565314.2751.6432651.523590.8939891000-11004.603641206.27278.175317.0251.7757681.5369220.8983071100-12004.603641206.27278.785319.7751.9082721.5502540.9043641200-

20、13004.603641206.27279.4322.5252.0418621.5635850.9121471300-14004.603641206.272710.025325.2752.1776241.5769170.9216911400-15004.603641206.272710.66328.0252.3155581.5902490.933041500-16004.603641206.272711.295330.7752.4534931.6035810.9461271600-17004.603641206.272711.93333.5252.5914271.6169130.9609241

21、700-18004.603641206.272712.57336.2752.7304471.6302450.9774841800-19004.603641206.272713.21339.0252.8694681.6435770.9957271900-20004.603641206.272713.85341.7753.0084881.6569081.015628四、天然氣粘度： 式中：天然氣的粘度，mpas ； 天然氣的密度，g/m3； mg天然氣的視分子量，取17kg/kmol； mair空氣的分子量，取29kg/kmol； t溫度，k； p壓力，mpa； 原油的相對密度 g天然氣的相對密度

22、，（空氣/air=10） 表3-4 各段天然氣粘度深度hkxy平均壓力p平均溫度t天然氣密度天然氣粘度0-100109.93195.5619881.2876022.765289.5250.0234620.011492100-200110.98075.5441861.2911633.21292.2750.0272050.0117200-300112.02715.5267171.2946573.68295.0250.0311280.011918300-400113.0715.5095691.2980864.18297.7750.0352610.012147400-500114.11255.4927

23、361.3014534.705300.5250.039540.012387500-600115.15165.4762081.3047585.245303.2750.043850.012632600-700116.18835.4599771.3080055.805306.0250.0482110.012885700-800117.22255.4440351.3111936.38308.7750.0525520.013143800-900118.25435.4283751.3143256.965311.5250.0568030.013402900-1000119.28375.4129891.317

24、4027.565314.2750.0609780.0136641000-1100120.31075.3978691.3204268.175317.0250.0650090.0139251100-1200121.33535.383011.3233988.785319.7750.0687950.014181200-1300122.35745.3684041.3263199.4322.5250.0723610.0144291300-1400123.37725.3540451.32919110.025325.2750.0757270.0146731400-1500124.39455.3399271.3

25、3201510.66328.0250.0788780.0149111500-1600125.40955.3260441.33479111.295330.7750.0817350.0151381600-1700126.4225.3123891.33752211.93333.5250.08430.0153541700-1800127.43225.2989581.34020812.57336.2750.0866030.0155591800-1900128.445.2857451.34285113.21339.0250.0886210.0157521900-2000129.44545.2727441.

26、34545113.85341.7750.0903610.015932五、地層水體積系數(shù)：表3-5 各段地層水體積系數(shù)深度h平均溫度t平均壓力p0-100289.5252.765-0.06924-0.000230.93055100-200292.2753.21-0.06889-0.000280.930853200-300295.0253.68-0.06851-0.000330.931178300-400297.7754.18-0.06811-0.00040.931523400-500300.5254.705-0.06768-0.000460.931888500-600303.2755.245-0

27、.06722-0.000540.932274600-700306.0255.805-0.06674-0.000620.93268700-800308.7756.38-0.06623-0.000710.933106800-900311.5256.965-0.06569-0.000810.933553900-1000314.2757.565-0.06513-0.000910.9340191000-1100317.0258.175-0.06453-0.001030.9345051100-1200319.7758.785-0.06392-0.001140.9350131200-1300322.5259

28、.4-0.06327-0.001270.9355411300-1400325.27510.025-0.0626-0.00140.9360881400-1500328.02510.66-0.0619-0.001540.9366541500-1600330.77511.295-0.06117-0.001690.9372411600-1700333.52511.93-0.06042-0.001840.9378481700-1800336.27512.57-0.05964-0.002010.9384751800-1900339.02513.21-0.05883-0.002170.9391231900-

29、2000341.77513.85-0.058-0.002350.93979第二節(jié) 流出動(dòng)態(tài)曲線一、井筒溫度分布曲線由溫度假設(shè)條件知，井筒流體溫度分布呈線性分布。因此用線性關(guān)系求解溫度分布曲線。已知井口、井底溫度如下表 3-6。表 3-6 井口、井底溫度表井口溫度，15地層溫度，70根據(jù)表 3-1計(jì)算溫度梯度，。設(shè)每一100m為一個(gè)單元，根據(jù)溫度梯度，求出中每100m處的溫度，然后將其繪制在一坐標(biāo)圖上，及得出井筒壓力溫度分布曲線圖見圖 3-1.圖 3-1 井筒溫度分布曲線圖二、井筒壓力分布曲線用orkiszewski方法計(jì)算井筒內(nèi)的油管壓力分布，其vb編程如下：private sub comma

30、nd1_click()pwh = 3.3 * 10 6: psc = 0.101 * 10 6: twh = 15for i = 1 to 20dp = 0.1 * 10 6: dth = (70 - 15) / 2000: h = 100rg = 0.65: ro = 0.85: mg = 17: tsc = 298: oo = 0.045: fw = 0.2: q = 60 / 86400qw = fw * q: rp = 100: qosc = q * (1 - fw): wor = qw / qosc: di = 0.062: mdw = 1020: g = 9.81: vs = 0.

31、244: e = 0.01524 * 10 (-3)uw = 1: rp = 100start:p1 = pwh + pscdp1 = p1 + dppav = (p1 + dp1) / 2tav = (twh + twh + h * dth) / 2pav1 = pav / (10 6)tav1 = tav + 273.15rs = 0.02124 * rg * (pav1 * 10 (1.76875 / ro - 0.001638 * tav) 1.205bo = 0.972 + 1.47 * 10 (-4) * (5.6146 * rs * (rg / ro) 0.5 + 2.25 *

32、tav + 40) 1.175mdo = (1000 * ro + 1.2 * rg * rs) / boc = 205735 * (5.625 * 10 (-2) * tav1 + 1) (-3.444)d = 10.313 * (log(5.625 * 10 (-2) * tav1 + 1) / log(10) + 1.5051) - 36.447uo = c * (log(1.067 / ro - 1) / log(10) + 2.1189) dppc = 4.8677 - 0.3565 * rg - 0.07653 * rg 2tpc = 103.8889 + 183.3333 * r

33、g - 39.72222 * rg 2pr = pav1 / ppctr = tav1 / tpczg = 1 - 3.52 * pr / (10 (0.9813 * tr) + 0.274 * (pr 2) / (10 (0.9813 tr)mdg = 3.4841 * 10 (-3) * rg * pav / (zg * tav1)k = 2.6832 * 10 (-2) * (470 + mg) * tav1 1.5 / (116.1111 + 10.5556 * mg + tav1)x = 0.01 * (350 + 54777.78 / tav1 + mg)y = 0.2 * (12

34、 - x)ug = 10 (-4) * k * exp(x * (mdg / 1000) y)bg = psc * zg * tav1 / (pav * tsc)vwt = -5.7325 * 10 (-3) + 2.40104 * 10 (-4) * (tav1 - 273.15) + 1.78412 * 10 (-6) * (tav1 - 255.37) 2vwp = -5.0987 * 10 (-7) * pav1 * (tav1 - 255.37) - 6.54435 * 10 (-9) * pav1 2 * (tav1 - 255.37) - 5.20574 * 10 (-5) *

35、pav1 - 4.74029 * 10 (-6) * pav1 2bw = (1 + vwt) * (1 + vwp)qg = bg * qosc * (rp - rs)ql = bo * qosc + bw * qwqm = qg + qlwm = (1000 * ro + 1.2 * rg * rp + 1000 * rw * wor) * qosca = 3.1415926 * di 2 / 4vsg = qg / avsl = ql / avm = vsg + vslmdl = mdo * (1 - fw) + mdw * fwngv = vsg * (mdl / (g * oo) 0

36、.25lb = 1.071 - 0.7277 * vm 2 / dils = 50 + 36 * ngv * ql / qglm = 75 + 84 * (ngv * ql / qg) 0.75ul = (uw * fw + uo * (1 - fw) / 1000if (qg / qm) 2300 then f = (1.14 - 2 * log(e / di + 21.25 / (re 0.9) / log(10) (-2) else f = 64 / re end iftf = f * mdl * (vsl / (1 - hg) 2 / (2 * di)dphz = (mdm * g +

37、 tf) / (1 - wm * qg / (a 2 * pav)elseif (qg / qm) lb and ngv 3000 and reb 8000 then vbb = vbc goto finish1 end if re1 = vm * di * mdl / ul vbc = (0.546 + 8.74 * 10 (-6) * re1) * (g * di) 0.5 reb = vbc * di * mdl / ulif reb = 8000 then vbb = vbc goto finish1end iffinish1:if vm 3.048 then co = (0.0023

38、6 * log(10 3 * ul + 1) / log(10) / (di) 1.415 - 0.14 + 0.167 * log(vm) / log(10) + 0.113 * log(di) / log(10) if co -0.2132 * vm then co = -0.2132 * vmelse co = (0.00537 * log(10 3 * ul + 1) / log(10) / (di) 1.371 + 0.455 + 0.569 * log(di) / log(10) - (log(vm) / log(10) + 0.516) * (0.0016 * log(10 3

39、* ul + 1) / log(10) / (di) 1.571 + 0.722 + 0.63 * log(di) / log(10) if co 2300 then f = (1.14 - 2 * log(e / di + 21.25 / (re1 0.9) / log(10) (-2) else f = 64 / re1 end if tf = f * mdl * vm 2 * (ql + vbb * a) / (qm + vbb * a) + co) / (2 * di) dphz = (mdm * g + tf) / (1 - wm * qg / (a 2 * pav)end ifxx

40、 = abs(dphz * h - dp)if xx 1 thenlist1.additem 壓力梯度為 & dphz else dp = dphz * h goto start end if if (qg / qm) lb and ngv ls then list1.additem 此時(shí)流態(tài)為段塞流 else list1.additem 此時(shí)流態(tài)為純液流 end if p = p1 + dphz * h y = h * ilist1.additem 第 & y & m壓力為 & p pwh = pwh + dphz * h twh = twh + dth * hnext i end sub由

41、程序計(jì)算得井筒壓力分布見表 3-7。表 3-7 井筒壓力分布表深度 h/(*100m)壓力p/mpa03.31003.3842004.3433005.3224005.9625006.6256007.3127008.0198008.7439009.481100010.232110010.994120011.764130012.665140013.566150014.467160015.367170016.268180017.169190018.069200019.103由井筒壓力分布表做井筒壓力分布曲線見下圖 3-2。圖 3-2 井筒壓力分布曲線從上圖，我們可以看出，從0-2000m均為段塞流。三、流出動(dòng)態(tài)曲線由程序得出不同產(chǎn)量下的井底流壓的數(shù)據(jù)如下表3-8。 表3-8 流出動(dòng)態(tài)數(shù)據(jù)井底 19.3 18.618.920.721.822.722.9產(chǎn)液量 1