井筒压力分布计算设计与实现

题目井筒压力分布计算目录第1章概述111设计的目的和意义112设计的主要内容1第2章基础数据2第3章能量方程理论331能量方程的推导332多相垂直管流压力分布计算步骤6第4章气液多相垂直管流压力梯度的摩擦损失系数法841基本压力方程842平均密度平均流速的确定方法843摩擦损失系数的确定1144油气水高压物性参数的计算方法1245井温分布的的计算方法1646实例计算17第5章设计框图及结果2151设计框图2152设计结果22结束语29参考文献30附录31第1章概述11设计的目的和意义目的确定井筒内沿程压力损失的流动规律，完成自喷井系统从井口到井底的所有相关参数的计算，运用深度迭代方法计算多相垂直管流的压力分布。意义利用所学的专业知识，结合已有的基础数据，最终计算井筒内的压力分布。对于油气井的优化设计、稳产高产及测试技术的预测性与精确性具有重要的现实意义。12设计的主要内容根据已有的基础数据，利用所学的专业知识，完成自喷井系统从井口到井底的所有相关参数的计算，最终计算井筒内的压力分布。计算出油井温度分布；确定平均温度压力条件下的参数；确定出摩擦阻力系数；确定井筒内的压力分布；详见第四章。第2章基础数据数据表见下表（表21）表21基础数据表地面脱气原油密度（KG/M3）841地层水比热（J/KG）4400天然气密度（KG/M3）0929天然气比热（J/KG）2227水密度（KG/M3）1000天然气分类（贫气或富气）富气水油比（M3/M3）01井号B1112P56井口温度（）15井深（M）1082地温梯度（/100M）315油管内径（MM）62传热系数（W/M）275油压（MPA）064饱和压力（MPA）628日产油量（T/D）405原油比热（J/KG）2200日产气量（M3/D）24447第3章能量方程理论31能量方程的推导流体流动系统都可根据能量守恒定律写出两个流动断面间的能量平衡关系进入断面1的流体能量在断面1和2之间对流体额外所做的功在断面1和2之间耗失的能量从断面2流出的流体的能量根据流体力学及热力学,对质量为M的任何流动的流体，在某一状态参数下（P、T）和某一位置上所具有的能量包括内能U；位能MGH；动能；压缩或膨胀2V能。V据此，就可以写出多相管流通过断面1和断面2的流体的能量平衡关系。为了得到各种管流能量平衡的普遍关系，选用倾斜管流。（31）22112MVMVUZPVQUZPVGSINGSIN式中流体质量，公斤；流体体积，；V3米压力，帕；P重力加速度，；G2米秒管子中心线与参考水平面之间的夹角，度；液流断面沿管子中心线到参考水平Z面的距离，米；HSIN图11流体流动示意图流体的内能，包括分子运动所具有的内部动能及分子间引力引起的内部位U能以及化学能、电能等，焦尔；流体通过断面的平均流速，米/秒。V（31）式中，除了内能外，其他参数可用测量的办法求得。内能虽然不能直接测量和计算其绝对值，但可求得两种状态下的相对变化。根据热力学第一定律，对于可逆过程或DQUPVDQPV式中DQ为系统与外界交换的热量；DU和PDV分别为系统进行热交换时，在系统内所引起的流体内能的变化和由于流体体积改变DV后克服外部压力所做的功。对于像我们这里所研究的这种不可逆过程来讲RDDQUPV式中DQR摩擦产生的热量。若以DLW表示摩擦消耗的功，则由上式可得RWL或（32）DDQPVLDWQPVL改写（31）式，可得到两个流动断面之间的能量平衡方程（31A）2SIN0MVUGZP将（31A）式写成微分形式（31B）DSIDDVZVQ将（32）式代入（31B）式，并简化后得（33）SIN0WVPMGL积分上式我们就可得到压力为P1和P2两个流动断面的能量平衡方程（33A）21WDSINVMGZL取单位质量的流体M1，将代入（33）式后得V（33B）W1DSIND0PVGZL式中流体密度，。3公斤米用压力梯度表示，则可写为（34）WDDSIN0IPVGZZ由此可得II式中单位管长上的总压力损失（总压力降）；DPZ由于动能变化而损失的压力或称加速度引起的压力损失；V克服流体重力所消耗的压力；SING克服各种摩擦阻力而消耗的压力。WDIZ令DSINPGZ举高DV加速度WDIPZ摩擦则DPZ摩擦举高加速度根据流体力学管流计算公式2DPVF摩擦式中F摩擦阻力系数；D管径，米。在Z的方向为由下而上的坐标系中为负值，如果我们取为正值，则DPZDPZ（35）2SINVGF（35）式是适合于各种管流的通用压力梯度方程。对于水平管流，因0，。若用X表示水平流动方向的坐标，则D0PZ举高（36）2DVFXD对于垂直管流，,SIN1，若以H表示高度，则90（37）2DDPVGF为了强调多相混合物流动，将方程中的各项流动参数加下角标“M”,则2MMSINDDVPGFZZ式中M多相混合物的密度；VM多相混合物的流速；FM多相混合物流动时的摩擦阻力系数。单相垂直管液流的；单相水平管液流的及均为零。对DPZ加速度DPZ举高加速度于气液多相管流，如果流速不大，则很小，可以忽略不计。DPZ加速度只要求得M、VM及FM就可计算出压力梯度。但是，如前所述，多相管流中这些参数沿程是变化的，而且在不同流动型态下的变化规律也各不相同。所以，研究这些参数在流动过程中的变化规律及计算方法是多相管流研究的中心问题。不同研究者通过实验研究提出了各自计算这些参数的方法。32多相垂直管流压力分布计算步骤按气液两相管流的压力梯度公式计算沿程压力分布时，影响流体流动规律的各相物理参数密度、粘度等及混合物的密度、流速都随压力和温度而变，而沿程压力梯并不是常数，因此气液两相管流要分段计算以提高计算精度。同时计算压力分布时要先给出相应管段的流体物性参数，而这些参数又是压力和温度的函数，压力却又是计算中要求的未知数。因此，通常每一管段的压力梯度均需采用迭代法进行。有两种迭代方法用压差分段、按长度增量迭代和用长度分段、按压力增量迭代。用压差分段、按长度增量迭代的步骤是1已知任一点井口或井底的压力作为起点，任选一个合适的压力降作为计0PP算的压力间隔；2估计一个对应的长度增量,以便根据温度梯度估算该段下端的温度；PL1T3计算该管段的平均温度及平均压力，并确定在该和下的全部流体性质TPTP参数；4计算该管段的压力梯度LPD/5计算对应于的该段管长P；/LPD6将第5步计算得的与第2步估计的进行比较，两者之差超过允许范围，L则以计算的作为估计值，重复25的计算，直至两者之差在允许范围内为止；L07计算该管段下端对应的长度及压力ILIP,I1,2,3,NINI1I08以处的压力为起点，重复第27步，计算下一管段的长度和压力，IL1IL1IP直到各段的累加长度等于或大于管长L时为止。NL第4章气液多相垂直管流压力梯度的摩擦损失系数法41基本压力方程摩擦损失系数法计算压力梯度的基本方程2MVPGHD式中计算段的混合物平均密度，KG/M3；M计算段的混合物平均流速，M/S；V计算段的摩擦损失系数，无因次；D管径，M；G重力加速度，9807M/S2；计算管段的平均压力梯度，PA/M。PH如果用混合物流量表示流速，则上式可写成20TM5M13QWPGHD式中Q0地面脱气原油的产量，M3/S；WT随1M3地面脱气原油同时产出的油、水、气混合物的总质量，KG/M3；其他符号及单位同前。42平均密度、平均流速的确定方法自喷井沿井筒自下而上各个流过断面处油、气、水混合物质量是始终不变的，而体积流量QMT和平均流速逐渐增大，所以油、气、水混合物的重度逐MV渐减小。但是，油井稳定生产时，单位时间内生产的地面（脱气）原油体积是不变的，并从生产日报表中直接查到。为了找出油、气、水混合物的体积流量QMT沿井筒的变化规律，取固定值1M3地面脱气原油的体积，作为研究混合液流的参考值，则有0TTQQV式中Q0产油量，M3/S；VT在某压力和温度下，伴随每生产1M3地面脱气原油的油、水、气总体积M3，/M3,即TOWGV地面每生产1M3脱气原油，在压力P和温度T下油应具有的体积，等于1M3脱气原油乘以该压力、温度条件下的体积系数。BO随压力P和温度T的变化关系，可由高压物性资料得出。当地面每生产1M3脱气原油时，在压力P和温度T下，水应具有的体积可通过生产油水比VW来表示。生产油水比等于产油量比产油量，单位为M3/M3。由于水压缩性很小，可以认为井筒内各流过断面处水的体积是不变的。当地面每生产1M3脱气原油时，在压力P和温度T下天然气应具有的体积VG可作如下分析，并通过气体状态方程式求得。设RP生产油气比，等于产气量比产油量，M3/M3；RS溶解油气比，M3/M3。即在压力P和温度T下，溶解在相当于1M3地面脱气原油中的天然气量。所以原来在压力P和温度T时，伴随每生产1M3地面脱气原油的天然在标准状态下所占有的体积为G0PSVR但是，还需要气休状态方程式将VG0换算到某压力P与温度T下的体积。根据气体状态方程式，知G0PZT式中P0标况压力（绝对），100KPA；T0标况温度，293K；VG0在标准压力和温度下的天然气的体积，M3；P压力（绝对），PA；T温度，K；VG在压力P和温度T下的天然气的体积，M3；Z0、Z气体在标准状态与某压力、温度下的压缩因子、无因次。所以，当Z01时，由上式可整理得0GGZPVT由上面两式可得出当地面每生产1M3脱气原油时，在压力P和温度T下，天然气（自由气）应具有的体积为00GGGSZPTVRA综合以上的分析，当地面每生产1M3脱气时，在某流过断面处油、气、水混合物在压力P和温度T下的体积为0TOGSWZPTVBRV当P1和P2相差不大时，可以用上式来计算某压力（P1和P2）和温度（T1和T2）范围内VT的平均值。只是上式中P应该采用P1和P2的平均值TPAVG，T应该采用该温度范围的平均值TAVG。其它随压力和温度而变化的各值如BO、Z、RS等也应该采用PAVG和TAVG下的值。于是得0AVGTOSWZVBRV平均密度为TMWV式中在平均压力和平均温度下，油、气、水混合物的密度，KG/M3；MWT与1M3地面脱气原油同时产出的油、气、水混合物总质量，KG/M3；与1M3地面脱气原油同时产出的油、气、水混合物在平均压力和TV平均温度下的总体积，M3/M3。总质量WT为TOGPWWRV式中O地面脱气原油密度，KG/M3；G天然气密度（标准条件下），KG/M3；W水的密度，KG/M3；VW水油比，M3/M3。混合物平均流速VTM24OQVVD式中在平均压力和温度下（即计算管段）油、气、水混合物的平均流MV速，M/S。43摩擦损失系数的确定摩擦损失系数是利用由矿场资料相关的关系确定。2ER两相雷诺数与单相雷诺数的关系为2ER2ERBLEAG式中气相雷诺数（纯气体流动的雷诺数）GEGPSOGEAQDR液相雷诺数（纯液体流动的雷诺数）LERLPWSGOLEAVQR平均温度和平均压力下气体的粘度，PAS；G平均温度和平均压力下液体的粘度，PAS；L与气液质量比有关的常数。,AB油、水、气三相混合物的液相粘度在未发生乳化的情况下，可根据相应条件下油的粘度。和水的粘度按体积加权平均求得，WWOLF1式中含水（体积比），小数；WF常数A和B的选取，应当使两相流在任一极端情况下，即只单相流时，两相雷诺数也应随着成为相应的单相雷诺数，一般取A、B为KEKA1；式中K气、液质量比，无因次；WSGOPVRTP下，液相的质量流量及温度在压力下，气相的质量流量及温度在压力根据矿场资料绘制相关曲线时选定的常数。，2ER利用大庆油田自喷井资料做的曲线，在取10、1时有较好的相关性。则两相雷诺数为KELGEERR102再利用下式求取摩擦损失系数23210LGLGEERAA其中2501020854361278019KKKA295712418193176252105490258830400KKKA44油气水高压物性参数的计算方法441溶解油气比需要先计算天然气在6895KPA表压下的相对密度514359201832LG0165OGSPTP式中，6895KPA表压下的天然气相对密度，无因次；S压力（绝对）和温度下的天然气相对密度，无因次；GPT温度，；T压力（绝对），KPA；标准状态下，原油的相对密度，无因次。O求得天然气的在6895KPA表压下的相对密度后，再利用下式即可求得溶解油气比213145078450EXP892COSGSOST式中C1、C2、C3系数，其值见表；P压力（绝对），KPA。系数08762O08762OC1C2C300362109372572400017811870239310442原油体积系数1、当时BP1231453145356856182OOOSSOGSOGSBCSTCST系数C1、C2、C3的值如下表所示系数08762O08762OC1C2C34677104175110518111084670104110010513371092、当时BP0EXPOBBBC其中，01234561356895182/OSGSCASATAAP式中泡点压力下的原油体积系数，M3/M3；OBBBPA114330；A250；A3172；A411800；A51261；A6105。443天然气压缩系数当天然气的压力低于35MPA时，它的压缩系数可以按下式计算2331046758061251553RRRRRRZTT其中，027RCRRRCTPZ式中Z天然气的压缩系数，无因次；对比温度，无因次；RTT温度，K；天然气的假临界温度，KC天然气的对比度，无因次；R天然气的对比压力，无因次；P天然气的假临界压力，KPA。C天然气的假临界温度和假临界压力，可以根据不同情况按下列公式计算CTCP1、富气当天然气的相对密度（空气为1）时07NG3267CNGT5894P当时07NG1062CNGT478P2、贫气当时07NG92176CNGT483P当时07NG92176CNGT48P按上式计算Z值时，需要使用迭代法。一般从设Z1开始，迭代五次即可。444原油粘度1、地面脱气原油的粘度10XON其中16382045134ZOXYT2、饱和原油的粘度BONA其中05138107564SASB445天然气的粘度310EXPYGGC其中1526781305492NGNGTXYX式中，管道条件下天然气的黏度，MPAS；G管道条件下天然气的密度，103KG/M3。446水的粘度252EXP1034791083198203WTT式中水的粘度，MPAS45井温分布计算方法由地面到油层温度是按地温梯度逐渐增加的。所谓地温梯度，即深度每增加100M地层温度的升高值。而在井筒中，由于地层流体不断地向上流动，地层流体便作为热载体将热量也不断地携带上来。通过套管、水泥环向地层传导。因此，井温总是比地温要高。流体的物性参数随温度变化，因此，计算应采用井温来进行流体参数计算。计算常规采油和井筒加热时沿井深温度分布的基本方程为EXP101MLTWKQM对于常规采油来说，可取01式中油管中L位置处原油的温度，；总传热系数，W/M；1K井底原油温度，；0T地层温度梯度，/M；M重力加速度，M/S2；G内热源，W/M；1Q计算段起点高度井底为0，M；L水当量，W/。W水当量可如下计算GWWCMC0式中原油的质量流量，KG/S；OM水的质量流量，KG/S；W井筒中气体质量流量，KG/S；G产出原油的比热，J/KG；OC产出水的比热，J/KG；W产出天然气的比热，J/KG。G在同一口油井，地温梯度M和井底温度都是不变的，传热系数则受地0T1K层物性和地层热阻、油管环形空间介质及其物性和油井的产量等多种因素的影响，而产量对的影响较小。故在一定的地层条件及井筒状况下，也可近似地1K认为为一常数。这样，整个井筒的温度分布就只受与油井产量有关的水当量1W和距井底的距离L的影响。46实例计算某含水自喷井产油量,产气量，油压O405T/DQGQ3247M/DWHP，内径D62MM，油井深1082M，试求井底压力。064MPA解用深度增量迭代方法计算。1选取压力间隔500KPA，假设对应的深度增量50M,P1H2从井口起计算第一段的平均压力及温度平均压力WH120641059KPA绝对）根据井口温度、地温梯度及假定的H1算得的平均温度30249KT即2934。3确定下的流体性质参数PT、溶解油气比RS3607M3/M3，天然气粘度G0010984MPAS,气体压缩因子Z09733，原油体积系数BO10256，原油粘度O20945MPAS,原油密度O841KG/M3,G0929KG/M3水的粘度W08968MPAS4计算混合物平均密度计总质量WTM（1）下的气体体积VGPT、生产油气比M3/M3P24750RP0100KPA,T020GPS03V12947159346075MZ（2）下的混合物总体积VTPT、31025601627MTOGWB（3）混合物的总质量WT84391098160KGTOPWRV4计算混合物的平均密度M3981607259KG/MTMWV5计算摩擦损失系数（1）气相雷诺数34053607924138070785782184OPSGREGQND（2）液相雷诺数RE34050924590785786451860OGSWLLQRV（3）气、液质量比430927281510PSGOWKV（4）两相雷诺数04721010472124359682KEEEEGLR20123LLG2658LG107LG168EEARA489其中20521368025410210765407AKKK2149571281793810832250680425810154907172AKKK6计算压力梯度及深度增量2522521340981667598473109KPA/M74OTMMQWPGHD/9073154MPH7比较深度增量的假设值和计算值0H假设计算如果取001M，则0050134851所以，将13541作为新的假设值，从第2步重新开始计算，即第二次迭代，直到满足要求后再开始计算第二段。第5章设计框图及结果51设计框图开始输入基础数据确定起点压力及计算浓度Z和分段数K1假设温度梯度初设计算段压降P，并计算下端压力确定平均温度及平均压力根据相关公式求出平均温度与平均压力下流体物性参数根据流体物性参数求出气液混合物压力梯度|P22P11|0以下端为起点KN结束NOKK1YESYESP21P22NO52运行界面521程序进入主界面522程序运行主界面523数据计算结果显示界面53设计结果表531计算结果数据表井深压力温度平均密度原油粘度溶解油气比原油体积系数天然气压缩系数天然气粘度水的粘度MMPAKG/M3MPASM3/M3M3/M3MPAS（MPAS007428651284522312561022009800010909120090292115402214132410240097600110090401042976175152067383102700972001100896011730321935920034381029009690011008880129308721020194648910310096600111087100140314222548189553710330096300111086120150319623974184758410350096000112085140161325125322180363010360095800112084160171330426605176167510380095600113083180180335827835172271910400095300113082续表井深压力温度平均密度原油粘度溶解油气比原油体积系数天然气压缩系数天然气粘度水的粘度MMPAKG/M3MPASM3/M3M3/M3MPAS（MPAS200190341129021168476310420095100113081220200346330170164980610440094900114081240209351531286161485010450094700114080260219356732375158189310470094500115079280228361833438154993710490094300115078300237366934480151898010510094100115077320247371935502148810241052009380011607734025637693650614591068105400936001160763602663818374941430111310560093400117075380275386638467140311581058009320011707440028539143942613761203105900930001170744202953961403731349124910610092800118073440305400841307132412951063009270011807246031440544223012991342106500925001190724803244099431431274138910660092300119071500334414344045125014371068009210012007152034441874493712271486107000919001200705403554229458191204153510720091700121069560365427146692118115851074009150012106958037543124755611591636107500913001210686003864353484111138168710770091100122068620397439249257111717391079009090012206764040744305009410961792108100906001230676604184467509231075184510830090400123066680429450351743105619001085009020012406670044045385255510361955108700900001240667204514572533581017201110880089800125065740463460454152998206810900089600125065760474463554938980212510920089400126064780486466555716961218410940089200126064800497469456486944224310960089000127064820509472157247926230310980088800127063840521474657999909236511000088500128063860533477058744892242711020088300128063880545479359479876249011040088100129063900558481360207860255411060087900129062920570483260927844262011080087600130062续表井深压力温度平均密度原油粘度溶解油气比原油体积系数天然气压缩系数天然气粘度水的粘度MMPAKG/M3MPASM3/M3M3/M3MPAS（MPAS94058348496163882826861110008740013006296059548646234081327531112008720013106298060848776303579828221114008690013206210006214888637227832891111500867001320611020634489664474769296211160086500133061104064749036531575530341114008620013306110606604907661587413108111200860001340611080674490867003727318311100085700135061108267549086708872531911110008570013506110406474903653157553034111400862001330611060660490766158741310811120086000134061108067449086700372731831110008570013506154结果图表（1）压力井深曲线图541压力井深曲线压力井深曲线表现为一段初始阶段向上凹其最终接近于一条直线的曲线段，这说明总体上压力随井深的增加而增加，但是压力梯度是不一样的，在井深小于800M时，压力梯度随井深的增加而增加，而当井深大于800M时，压力梯度几乎不随井深的变化而变化。这是由于混合物平均密井深小于800M时，平均密度逐随井深的增加而增大，而油管内压力梯度与混合物的密度成正比，所以在井深小于800M时，曲线的斜率随井深的增大而增大；当井深大于800M时，曲线的斜率几乎不变，曲线接近于一条直线。（2）温度井深曲线图542温度井深曲线由井温曲线分布图可看出随着井深增加井温变化斜率逐渐减小，并在井底处斜率趋近于零。由地面到油层温度是按地温梯度逐渐增加的。所谓地温梯度，即深度每增加100M地层温度的升高值。而在井筒中，由于地层流体不断地向上流动，地层流体便作为热载体将热量也不断地携带上来。通过套管、水泥环向地层传导。因此，井温总是比地温要高。由于井底压力等于油层压力，而井口油管流体温度大于地面温度，所以油管内温度梯度小于地层地温梯度。所以油管内液体随着流体的不断向上流动温度降低速率小于地层温度降低速率，管内温度与地层温度差在井底处为0，并随着井深的减小不断增加，管内流体与地层的传热量逐渐增加，所以管内流体温度降低的速率逐渐增大，即表现为井温分布曲线图中随着井深增加，曲线斜率逐渐减小，并在接近井底处曲线斜率趋近于零。（3）混合物平均密度井深曲线图543密度井深曲线由密度井深曲线可知，当井深小于1000M时，密度随着井深的增加而增加，当井深大于1000M时，混合物的平均密度变化幅度比较小。由压力与深度曲线可知，随着井深增加管内压力增大。井深小于1000M时，随着井深增加，溶解油气比逐渐增加，直至等于生产油气比，即天然气全部溶于原油中。而气液混合物平均密度（其中O为原油密度，OPGWMSRVBG为天然气在标况下的密度，W为水的密度，RP为生产油气比，RS为溶解油气比，BO为原油体积系数，VW为水油比），可知随着井深增大，（RPRS）逐渐减小，混合物平均密度随之增大，即表现为密度曲线图中斜率逐渐增大。井深达到1000M时，随着井深增加，溶解油气比不再发生变化，恒等于生产油气比。这时由气液混合物平均密度公式可知气液混合物平均密度只与生产油气比有关，而生产油气比随井深增加变化不大，所以混合物平均密度随井深增加变化不大，即表现为混合物平均密度曲线中曲线斜率变化较小。结束语进入卓越工们程师班，我觉得收获颇多。在这三周里，虽然比较辛苦，但是我学到了许许多多课堂上学不到的东西。首先，我在C语言的基础上，有学习了一种新的语言，经过使用这两种语言，发现它们各有特色，C语言由于指针的存在而特别灵活，但是VB却有比较好的界面，更方便我们编制小软件，通俗易懂。另外，VB看似简单，但是有许多技巧可言，在这几周里，在张文老师的指导下，学到了非常多的知识，实际在编写代码和运行程序的过程中会出现各种各样的问题，需要把各个窗体、模块联系起来，稍有差错，就不会得到正确的结果。此外，在张老师的指导下，学会了导出数据、画图表，虽然说EXCEL可能也会做成这些事，但是，它又有自己独特的功能。通过这次课程设计，对我来说是一次能力的提升，综合的挑战。与此同时，我也掌握了相关的主要内容，如油气物性参数的计算、能量方程的推导、按深度增量迭代的步骤方法等等。切实体会到了把课堂所学的知识应用于实际资料来解决实际问题，从而把理论与实践有机结合起来。在这一段时间里，同学之间也是不断探讨、搜集相关信息、上网查阅资料，经历一翻波折终于完成了这次课程设计。但是在本次课程设计中我也暴露出来了一些问题，比如说VB使用还不够灵活，这也许是自己刚刚接触VB，练得还是不够的缘故。最后，真诚地感谢张文老师以及所有对同学们的知道，在今后的学习生活中，我一定会更加注重专业素养的提高，为祖国石油事业贡献自己的力量。参考文献1王鸿勋，张琪采油工艺原理M石油工业出版社，北京，2000，74852衣治安，吴雅娟主编实用计算机基础教程M石油工业出版社，20023陈涛平等石油工程M石油工业出版社，20004蒋加伏，张林峰VISUALBASIC程序设计教程北京邮电大学出版社，2009附录程序内容（1）模块中的程序PUBLICJINGHAOASSTRING井号PUBLICHASSINGLE井深PUBLICDASSINGLE油管内径PUBLICPWHASSINGLE油压PUBLICQOASSINGLE日产油量PUBLICQGASSINGLE日产气量PUBLICMIDUOASSINGLE地面脱气原油密度PUBLICMIDUGASSINGLE天然气密度PUBLICMIDUWASSINGLE水密度PUBLICMIDUORASSINGLE地面脱气原油相对密度PUBLICMIDUGRASSINGLE天然气相对密度PUBLICMIDUWRASSINGLE水相对密度PUBLICWOASSINGLE水油比PUBLICGTASSINGLE井口温度PUBLICGRASSINGLE地温梯度PUBLICGCASSINGLE传热系数PUBLICPBASSINGLE饱和压力PUBLICCOASSINGLE原油比热PUBLICCWASSINGLE地层水比热PUBLICCGASSINGLE天然气比热PUBLICFLAGASINTEGER天然气分类富气赋值为1贫气赋值为0PUBLICFWASSINGLE含水率不变量在调用时需赋值PUBLICWTASSINGLE1M3地面脱气原油同时产出的油、水、气混合物的总质量KG/M3PUBLICVTASSINGLE井筒条件下产出WT油气水混合物所对应的体积PUBLICDGR689ASSINGLE6895KPA表压下的天然气相对密度PUBLICTPBASSINGLE饱合压力对应的温度PUBLICP0ASSINGLE标准状态下压力KPA常量在调用时需赋值PUBLICGD0ASSINGLE空气密度KG/M3常量在调用时需赋值PUBLICT0ASSINGLE标准状态下温度常量在调用时需赋值PUBLICGASSINGLE重力加速度M/S2常量在调用时需赋值PUBLICRPASSINGLE生产油气比不变量在调用时需赋值PUBLICAPASSINGLE油管截面积M2不变量在调用时需赋值PUBLICNNASSINGLE计算点的个数数组的定义PUBLICSHUZUH2000ASSINGLE计算点深度MPUBLICSHUZUP2000ASSINGLE计算点压力MPAPUBLICSHUZUT2000ASSINGLE计算点温度PUBLICSHUZUD2000ASSINGLE计算点混合物平均密度KG/M3PUBLICSHUZUUO2000ASSINGLE计算点原油粘度MPASPUBLICSHUZURS2000ASSINGLE计算溶解油气比PUBLICSHUZUBO2000ASSINGLE计算点原油体积系数PUBLICSHUZUZ2000ASSINGLE计算点天然气压缩系数PUBLICSHUZUUG2000ASSINGLE计算点天然气粘度MPASPUBLICSHUZUUW2000ASSINGLE水的粘度MPAS求溶解油气比RS（P,T函数PUBLICFUNCTIONRSPASSINGLE,TASSINGLEASSINGLEDIMC1ASSINGLE,C2ASSINGLEIFMIDUOR08762THENC100362C210937C325724ELSEC100178C21187C323931ENDIFRS01781C1DGR6890145PC2EXPC314151315MIDUOR/MIDUOR18T492IFRSRPTHENRSRPENDFUNCTION求原油体积系数BOP,T函数PUBLICFUNCTIONBOPASSINGLE,TASSINGLEASSINGLEDIMBOBASDOUBLE,C1ASSINGLE,C2ASSINGLE,C3ASSINGLEDIMA1ASSINGLE,A2ASSINGLE,A3ASSINGLE,A4ASSINGLE,A5ASSINGLE,A6ASSINGLEDIMC0ASSINGLEIFP08762THENC14677104C21751105C31811108ELSEC1467104C211105C31337109ENDIFBO15615C1RSP,TC218T2814151315MIDUOR/MIDUORDGR6895615C3RSP,T18T2814151315MIDUOR/MIDUORDGR689ELSEBOBBOPB,TPBA11433A25A3172A41180A51261A6105C06895A15615A2RSP,TA318T32A4DGR689A514151315MIDUOR/MIDUOR/A6PBOBOBEXPC0PPBENDIFENDFUNCTION求天然气压缩系数ZP,T函数富气为1,贫气为0PUBLICFUNCTIONZPASSINGLE,TASSINGLEASSINGLEDIMTCASSINGLE,PCASSINGLEDIMTRASSINGLE,PRASSINGLE,DRASSINGLE（DR为天然气的对比度）DIMIASINTEGERIFFLAG1THENIFMIDUGR07THENTC13211667MIDUGRPC510268948MIDUGRELSETC10615222MIDUGRPC477824821MIDUGRENDIFELSEIFMIDUGR07THENTC9217667MIDUGRPC488138611MIDUGRELSETC9217667MIDUGRPC477824921MIDUGRENDIFENDIFZ1TR273T/TCPRP/PCFORI1TO5DR027PR/ZTRZ103150610467/TR05783/TR3DR0535306123/TR06315/TR3DR2NEXTIENDFUNCTION求原油粘度UO函数PUBLICFUNCTIONUOPASSINGLE,TASSINGLEASSINGLEDIMUONASSINGLEDIMXASSINGLE,YASSINGLE,ZZASSINGLE,AASSINGLE,BASSINGLEZZ3032400202314151315MIDUOR/MIDUORY10ZZXY1820321163地面脱气原油温度取15，对结果影响较大UON10X1A107155615RSP,T1000515B5445615RSP,T1500338UOAUONBENDFUNCTION求管道条件下天然气密度DGP,T函数PUBLICFUNCTIONDGPASSINGLE,TASSINGLEASSINGLEDGMIDUGP27315T0/ZP,TP027315T根据气体状态方程推导,地面标准条件下Z1ENDFUNCTION求管道条件下天然气相对密度DGRP,T函数PUBLICFUNCTIONDGRPASSINGLE,TASSINGLEASSINGLEDGRDGP,T/GD0ENDFUNCTION求天然气粘度UGT函数PUBLICFUNCTIONUGPASSINGLE,TASSINGLEASSINGLEDIMXASSINGLE,YASSINGLE,CASSINGLEX35548/T27315029MIDUGRY2402XC1260078MIDUGR273T15/116306MIDUGR273TUGC103EXPXDGP,T103YENDFUNCTION求水粘度UWT函数PUBLICFUNCTIONUWTASSINGLEASSINGLEUWEXP1003147910218T32198210518T322ENDFUNCTION求以10为底的对数LOG10X的函数PUBLICFUNCTIONLOG10XASSINGLEASSINGLELOG10LOGX/LOG10ENDFUNCTION求混合物的密度DMP,T函数PUBLICFUNCTIONDMPASSINGLE,TASSINGLEASSINGLEVTBOP,TWOZP,TP0273T/P273T0ABSRPRSP,TDMWT/VT混合物的平均密度ENDFUNCTION求压力梯度函数PUBLICFUNCTIONRPASSINGLE,TASSINGLEASSINGLEDIMMASINTEGER,NASINTEGERDIMKASSINGLE,AASSINGLE,BASSINGLEDIMA0ASSINGLE,A1ASSINGLE,A2ASSINGLEDIMREGASSINGLE,RELASSINGLE,RE2ASSINGLEDIMFASSINGLE摩擦阻力损失系数DIMULASSINGLE平均温度和平均压力下的液体粘度DIMAPASSINGLE求气相、液相雷诺数ULUOP,T1FWUWTFWAP314D2/4REGDQOABSRPRSP,TMIDUG/AP/UGP,T/1000天然气粘度为MPAS因此要除以1000RELDQOMIDUOGDGRSP,TMIDUWWO/AP/UL/1000求两相雷诺数KABSRPRSP,TMIDUG/MIDUOGDGRSP,TMIDUWWOM10N1AMK/MK1B1/EXPNKRE2REGARELB求摩擦损失系数FA0201919/100K50123783638606285044100K0212100K2A1102862/100K4957120107321715179193051100K012118100K2A2015604/100K5010732188304025857100K001549100K2F10A0A1LOG10RE2A2LOG10RE22RDMP,TGFQO2WT2/1234D5DMP,TENDFUNCTION井温分布计算PUBLICFUNCTIONTLASSINGLEASSINGLEDIMWASSINGLE水当量DIMQ1ASSINGLE内热源，对于常规采油，Q10DIMMOASSINGLE,MWASSINGLE,MGASSINGLEDIMTWFASSINGLE井底温度Q10MOQOMIDUOMWQOWOMIDUWMGQGMIDUGWMOCOMWCWMGCGTWFGTGRH油层温度等于井口温度加上井深与地温梯度的乘积LHL以井底为零点的深度TWGRQ1/GC1EXPGC/WLTWFGRL井底为零点深度L处的温度LHLENDFUNCTION（2）窗体1中的程序PRIVATESUBCOMMAND1_CLICKFORM1HIDEFORM2SHOWENDSUBPRIVATESUBCOMMAND2_CLICKENDENDSUBPRIVATESUBFORM_LOADPICTURE1PICTURELOADPICTURE“CDOCUMENTSANDSETTINGSADMINISTRATOR桌面井筒压力分布PICTUREJPG“ENDSUB（3）窗体2中的程序PRIVATESUBMIDU_CLICKFORM6SHOWENDSUBPRIVATESUBOPEN_CLICKCD1FILTER“数据文件（TXT|TXT“CD1INITDIRAPPPATHCD1SHOWOPENOPENCD1FILENAMEFORINPUTAS1DIMAASSTRINGFORI0TO16INPUT1,AIFA“THENTEXT1ITEXTANEXTIINPUT1,AIFA1THENOPTION1VALUETRUEELSEOPTION2VALUETRUEENDIFCLOSE1ENDSUBPRIVATESUBPRINT_CLICKCD3SHOWOPENCD3ACTION5打开打印机对话框ENDSUBPRIVATESUBQUIT_CLICKENDENDSUBPRIVATESUBSAVEYUANSHI_CLICKDIMAAASSTRINGDIMFNAMEASSTRINGDIMRESPONSEFNAMEREPLACETEXT10TEXT,CHR13CHR10,“号井原始数据TXT“FNAMEAPPPATH“FNAMEOPENFNAMEFOROUTPUTAS1FORI0TO16AATEXT1ITEXTWRITE1,AANEXTIIFOPTION1VALUETRUETHENWRITE1,1ELSEIFOPTION2VALUETRUETHENWRITE1,0ELSEWRITE1,“ENDIFCLOSE1RESPONSEMSGBOX“保存原始数据“,33,“提示“IFRESPONSEVBOKTHENMSGBOX“原始数据已保存在程序所在文件夹，文件名为“TEXT10TEXT“号井原始数据TXT“,VBOKONLY,“提示“ENDIFENDSUBPRIVATESUBSHUJU_CLICKMSGBOX“需要帮助请与赵二猛联系“,VBOKONLY,“提示“ENDSUBPRIVATESUBWENDUQUXIAN_CLICKFORM5SHOWENDSUBPRIVATESUBYALI_CLICKFORM4SHOWENDSUBPRIVATESUBYIJIAN_CLICKMSGBOX“有建议请与赵二猛联系“,