天然气工程课程报告.doc

天然气工程课程报告 题目5.1由已知得：气藏压力4505psia测试点1，流量1152mscf/d 井底压力3025psia测试点2，流量1548mscf/d 井底压力1685psia井底流压1050psia 温度180f天然气相对密度0.65（空气为1）氮气摩尔分数0.1二氧化碳摩尔分数0.08硫化氢摩尔分数0.021.使用真实气体压力和forchheimer方程求解测井点1的压力的真实气体拟压力3025psia：588157460/cp测井点2的压力的真实气体拟压力1685psia:198040416/cp预定压力的真实气体拟压力1050psia：79585534/cp288208000根据：得q=1653mscf/d2.使用压力平方法forchheimer方程求解4.055012根据：得：q=1645 mscf/d3.使用真实气体拟压力和回压模型求解0.610.0053q=1656mscf/d4.使用压力平方法和回压模型求解：0.660.00531656 mscf/d问题5.2pwf（psia）q（mscf/d）压力法压力平方法1519941067245190610644751817105770417261044934163510261164154310041394145097616241355943185412599052084116386223141064941254496576027748647003004762635324365856334635534863693446403392333731241532272164383114111461300井底流压与产量的关系曲线7.1平均参数法计算井底流压：1、第一次试算2、计算平均参数利用经验公式：3、查表方法计算偏差系数z：查得偏差因子z=0.93经计算地层天然气的粘度ug=0.013；4、计算re；5、计算摩阻系数f6、计算无因次量7、计算井底流压8、第二次试算经计算得两次计算pwf结果，所以求得井底流压为6.20mpa。examination 11.欢迎界面以及基本参数的输入2.酸性气体偏差系数计算(hy模型/dpk模型/dpr模型与wa校正法)2.酸性气体偏差系数计算(hy模型/dpk模型/dpr模型与wa校正法)3.天然气相平衡计算(露点模型/泡点模型/等温闪蒸模型)4天然气两相管流压力计算(mukherjeebrill模型)5.天然气水合物生成预测及抑制剂用量计算二、程序代码分析1.气体偏差系数计算（hy模型，dpk模型，dpr模型）1.1hy模型+wa校正法1.1.1hy偏差系数计算模型 该法以 starling-carnahan 状态方程为基础 ,通过对 standing-katz 图版进行拟合 ,得到以下关系式:z = 0. 06125* ppr / (*tpr ) *exp - 1. 2 *(1 - 1/ tpr )2r 为拟对比密度 ,可用牛顿迭代法从如下方程求得:0 = (r +r2 +r3 r4)/(1 - r )3 -(14. 76 / tpr -9. 76 / tpr2+ 4. 58/tpr3)*r2+(90. 7/tpr-242. 2/tpr2+ 42. 4/tpr3)* r (2. 18+2. 82/ tpr)-0. 01652( ppr / tpr )* exp - 1.2*(1 - 1/ tpr )2 该法应用范围是:1. 2 tpr 3 ;0. 1 ppr 24. 0。部分代码如下：private sub command1_click()dim n(11) as double -定义各组分dim i as integerfor i = 0 to 11n(i) = val(欢迎界面.摩尔分数(i) / 100next idim p as double, t as double-定义压力、温度p = val(绝对工作压力)t = val(绝对工作温度)dim tc(11) as double, pc(11) as double-定义临界温度、临界压力tc(0) = 373.5833 pc(11) = 295.5315dim tpc as double, ppc as double-计算拟临界压力、拟临界温度tpc = 0ppc = 0for i = 0 to 11 tpc = tpc + n(i) * tc(i) ppc = ppc + n(i) * pc(i)next idim tpr as double, ppr as double-计算拟对比压力、拟对比温度tpr = t / tpcppr = p / ppcdim r as double-计算对比密度(牛顿迭代法)dim x as double, x1 as double, f as double, f1 as doubleabsolution = 1x = 0do while absolution 0.00001f = (-0.06152) * (ppr / tpr) * exp(-1.2) * (1 - (1 / tpr) 2) + (x + x 2 + x 3 - x 4) / (1 - x) 3) - (14.76 / tpr) - (9.76 / (tpr 2) + (4.58 / (tpr 3) * (x 2) + (90.7 / tpr) - (242.2 / (tpr 2) + (42.4 / (tpr 3) * (x (2.18 + (2.82 / tpr)f1 = (1 + 4 * x + 4 * (x 2) - 4 * (x 3) + x 4) / (1 - x) 4) - (29.52 / tpr) - (19.52 / (tpr 2) + (9.16 / (tpr 3) * x + (2.18 + (2.82 / tpr) * (90.7 / tpr) - (242.2 / (tpr 2) + (42.4 / (tpr 3) * (x (1.18 + (2.82 / tpr)x1 = x - f / f1absolution = abs(x1 - x)x = x1loopr = xdim z as double-用hy模型计算偏差系数z = (1 + r + r 2 - r 3) / (1 - r) 3) - (14.76 / tpr - 9.76 / (tpr 2) + 4.58 / (tpr 3) * r + (90.7 / tpr - 242.2 / (tpr 2) + 42.4 / (tpr 3) * (r (1.18 + (2.82 / tpr)hy模型z = format(z, 0.0000)-hy模型求得未校正的z1.1.2wa校正法wa校正法引入参数,主要考虑了一些常见的极性分子( h2 s、 co2 )的影响 ,希望用此参数来弥补常用计算方法的缺陷。参数的关系式如下:= 15 (m - m2) + 4. 167 ( n0. 5- n2)式中:m 为气体混合物中 h2 s与 co2 的摩尔分数之和; n 为气体混合物中 h2 s的摩尔分数。 根据 wichert2aziz 的观点 ,每个组分的临界温度和临界压力都应与参数有关 ,临界参数的校正关系式如下所示:t ci= tci - p ci= pci * t ci/ tci 式中: tci为i 组分的临界温度 ,k; pci为i 组分的临界压力 ,kpa ; tci为 i 组分的校正临界温度 , k; pci为 i组分的校正临界压力 ,kpa。部分代码如下：dim nhc as double, nh as double, as double, tc1(11) as double, pc1(11) as double-wa校正nhc = n(0) + n(2)nh = n(0) = 15 * (nhc - nhc 2) + 4.167 * (nh 0.5 - nh 2)for i = 0 to 11 tc1(i) = tc(i) - pc1(i) = pc(i) * tc1(i) / tc(i)next idim tpc1 as double, ppc1 as double, z1 as doubletpc1 = 0ppc1 = 0for i = 0 to 11 tpc1 = tpc1 + n(i) * tc1(i) ppc1 = ppc1 + n(i) * pc1(i)next idim tpr1 as double, ppr1 as doubleif t = 17.24 then dim t1 as double t1 = t + 1.94 * (p / 2760 - (2.1 * (10 (-8) * (p 2) tpr1 = t1 / tpc1 ppr1 = p / ppc1 z1 = (1 + r + r 2 - r 3) / (1 - r) 3) - (14.76 / tpr1 - 9.76 / (tpr1 2) + 4.58 / (tpr1 3) * r + (90.7 / tpr1 - 242.2 / (tpr1 2) + 42.4 / (tpr1 3) * (r (1.18 + (2.82 / tpr1) hy模型waz = format(z1, 0.0000)else tpr1 = t / tpc1 ppr1 = p / ppc1 z1 = (1 + r + r 2 - r 3) / (1 - r) 3) - (14.76 / tpr1 - 9.76 / (tpr1 2) + 4.58 / (tpr1 3) * r + (90.7 / tpr1 - 242.2 / (tpr1 2) + 42.4 / (tpr1 3) * (r (1.18 + (2.82 / tpr1) hy模型waz = format(z1, 0.0000)-校正后的偏差系数zend ifend sub1.2dak模型+wa校正法1.2.1dak偏差系数计算模型该模型与 dranchuk-purvis-robinsion 计算法相同 ,但相对密度应采用牛顿迭代法从下式求出:0=1 + (a1 + a2/tpr+ a3/tpr3+ a4/tpr4+a5/tpr5)* r +(a6+a7/tpr+ a8/tpr2)*r2- a9*(a7/tpr+ a8/tpr3) r5 +a10/tpr3*r 2*(1 + a11*r2 )* exp ( - a11*r2 ) - 0. 27*ppr/(r *tpr) 式中: ai为给定系数。 该法应用范围是:1. 0 tpr 3 ,0. 2 ppr 30 ;或者0. 7 tpr 1. 0 , ppr 0.00001f = (-0.27) * ppr / tpr + x + (a1 + a2 / tpr + a3 / (tpr 3) + a4 / (tpr 4) + a5 / (tpr 5) * (x 2) + (a6 + a7 / tpr + a8 / (tpr 2) * (x 3) - a9 * (a7 / tpr + a8 / (tpr 2) * (x 6) + a10 * (1 + a11 * (x 2) * (x 3) / (tpr 3) * exp(-a11 * (x 2)f11 = 1 + 2 * (a1 + a2 / tpr + a3 / (tpr 3) + a4 / (tpr 4) + a5 / (tpr 5) * x + 3 * (a6 + a7 / tpr + a8 / (tpr 2) * (x 2) - 6 * a9 * (a7 / tpr + a8 / (tpr 2) * (x 5)f12 = (a10 / (tpr 3) * (3 * (x 2) + a11 * (3 * (x 4) - 2 * a11 * (x 6) * exp(-a11 * (x 2)f1 = f11 + f12x1 = x - f / f1absolution = abs(x1 - x)x = x1loopr = xdim z as double-用dak模型计算偏差系数z = 1 + (a1 + a2 / tpr + a3 / (tpr 3) + a4 / (tpr 4) + a5 / (tpr 5) * r + (a6 + a7 / tpr + a8 / (tpr 2) * (r 2) - a9 * (a7 / tpr + a8 / (tpr 2) * (r 5) + a10 * (1 + a11 * (x 2) * (x 2) / (tpr 3) * exp(-a11) * (x 2)dak模型z = format(z, 0.0000)-dpk模型求得未校正的z1.2.2wa校正法同上，略。1.3dpr模型+wa校正法1.3.1dpr偏差系数计算模型dranchuk、purvis 和 robinsion 根据 benedict-webb-rubin 状态方程,将偏差系数转换为对比压力和对比温度的函数 ,于 1974 年推导出了带 8 个常数的经验公式 ,其形式为:z = 1 +( a1 + a2/ tpr+ a3/tpr3)* r+ (a4 + a5/tpr)* r2+(a5 * a6/tpr)* r 5 - a7/tpr3* r2*(1 + a8* r2 )* exp ( - a8* r2)式中: ai 为给定系数; ppr为拟对比压力 ,无因次; tpr为拟对比温度 ,无因次。 dpr法使用 newton-raphson迭代法解非线性问题可得到偏差系数的值。这种方法的使用范围是:1. 05 tpr 3 ;0. 2 ppr 30。部分代码如下：-计算对比密度dim a1 as double, a2 as double, a3 as double, a4 as double, a5 as double, a6 as double, a7 as double, a8 as doublea1 = 0.31506237a2 = -1.0467099a3 = -0.57832729a4 = 0.53530771a5 = -0.61232032a6 = -0.10488813a7 = 0.68157001a8 = 0.68446549dim x, x1 as double, f as double, f1 as doubleabsolution = 1x = 0do while absolution 0.00001f = -0.27 * ppr / tpr + x + (a1 + a2 / tpr + a3 / (tpr 3) * (x 2) + (a4 + a5 / tpr) * (x 3) + a5 * a6 / tpr * (x 6) + a7 / (tpr 3) * (1 + a8 * (x 2) * (x 3) * exp(-a8 * (x 2)f1 = 1 + 2 * (a1 + a2 / tpr + a3 / (tpr 3) * x + 3 * (a4 + a5 / tpr) * (x 2) + 6 * a5 * a6 / tpr * (x 5) + a7 / (tpr 3) * (3 * (x 2) + a8 * (3 * (x 4) - 2 * a8 * (x 6) * exp(-a8 * (x 2)x1 = x - f / f1absolution = abs(x1 - x)x = x1loopr = xdim z as double-用dpr模型计算偏差系数z = 1 + (a1 + a2 / tpr + a3 / (tpr 3) * x + (a4 + a5 / tpr) * (x 2) + a5 * a6 / tpr * (x 5) + a7 / (tpr 3) * (1 + a8 * (x 2) * (x 2) * exp(-a8) * (x 2)dpr模型z = format(z, 0.0000)-dpr型求得未校正的z1.3.2wa校正法同上，略。examination 2 节点分析取井底为解节点的节点分析（1） 取井底为解节点，则流入部分包括从地层外边界到井底，流出部分包括从井底到井口（2） 计算流入动态曲线。假设一系列产量，对每一产量，根据地层压力和产能方程（）计算井底压力，该压力就是流入节点压力，见表中第二栏计算流出动态曲线。假设一系列产量，对每一产量，由井口压力和单相气体垂直管流计算方法（） 取井底为节点时的流入和流出动态数据表如下：产量(104m3/d)流入节点压力mpa流出节点压力mpa0302.480.929.165342.6371311.928.063802312.637130992.926.75699912.637130994.124.872012.6371315.222.771212.6371316.320.199832.6371317.416.944972.6371318.512.482992.6371319.55.2368882.6371319.70.8443932.637131取井底为节点时的流入和流出动态数据表如下：求解协调点，协调点处的压力为2.64mpa，产量9.62*104 m3/d,所以该井在井口压力等于2mpa下的产量为9.62*104 m3/d。examination 3 临界流量： g=3.4844*103*0.6*3.21/(0.93*295)=24.46 ucr=3.1*60*10-3*9.8*(1074-g)/ g20.25=3.12 a=3.14*622/4=0.00302(m2) qcr=25000*0.00302*2*3.12/(50*0.93)=10.13examination 41）波诺马列夫法求水合物生成压力 t273.1k时 lgp=-1.0055+0.0541(b+t-273.1) t273.1k，所以lgp=-1.0055+0.0541(12.63+291.5-273.1)= 0.6786故在井口管道温度为18.5度的状态下，水合物生成压力为 p=4.77mpa（2）在此压力下形成温度点左移4度的乙二醇浓度运用hammerschmidt公式作近似计算乙二醇浓度式中 wr乙二醇浓度，% m乙二醇分子量，m=62； t露点降，0c。计算得出乙二醇浓度为16.88%examination 5井底产生积液，如何排液:井底积液雾化携液, 当高速气体到达积液井井底时会产生剧烈的鼓泡效应,利用气液两相流湍流夹带理论推导出气体所能带走的液体涡流的最大半径,利用湍流的能谱分析确定最小涡流尺寸,从而建立起预测井底积液雾化的最小采气产量模型。同时结合现有的描述井筒中液滴运动的携液模型,建立了关于井底积液雾化和井筒携液全过程的携液新模型。气体如何脱水: 1、溶剂吸收法利用适当的液体吸收剂以除去气体混合物中的一部分水份,对吸收后的贫溶剂进行脱吸,使溶剂再生循环使用。常用的脱水剂有二甘醇、三甘醇等。2 、固体干燥剂吸附法利用气体在固体表面上积聚的特性,使某些气体组分吸附在固体吸附剂表面,进行脱除。气体组分不同,在固体吸附剂上的吸附能力也有差异,因而可用吸附方法对气体混合物进行净化。工业上常用的固体吸附剂有硅胶、活性氧化铝和分子筛。吸附是在固体表面张力作用下进行的,根据表面张力的性质可将吸附过程分为物理吸附和化学吸附。物理吸附是可逆过程,可用改变温度和压力的方法改变平衡方向,达到吸附剂的再生。目前广泛采用的用分子筛作吸附剂脱除天然气中水分的过程就是物理吸附过程。3、冷冻分离法将一定温度的混合气体在一定压力下通过干燥的、最低温度可达- 20 以下的冷凝器,使混合气体中的水气变成液滴后分离。常用的设备有冷冻干燥器。4、脱水剂 4.1、三甘醇( teg) 脱水剂甘醇类化合物具有很强的吸水性,其溶液水点较低,沸点高,毒性小,常温下基本不挥发,所以广泛应用于天然气脱水。最先用于天然气脱水的是二甘醇,50 年代后teg 以良好的性能逐步取代了二甘醇成为最主要的脱水溶剂。teg 热稳定性好,易于再生,蒸汽压低,携带损失量小,露点降通常为33 - 47 。 4.2、分子筛吸附剂分子筛具有均一微孔结构,能将不同大小的分子分离的一种高效、高选择性的固体吸附剂。分子筛是一种人工合成的无机吸附剂,天然气脱水常用的是4a 和5a ,它是具有骨架结构的碱金属或碱土金属的硅铝酸盐晶体,其分子式为:m2/ noal2o3xsio2yh2o式中:m 某些碱金属或碱土金属离子, 如li ,na ,mg ,ca 等;n m 的价数;x sio2 的分子数;y h2o 的分子数。在脱水过程中,分子筛作为吸附剂的显著优点是:具有很好的选择吸附性,也就是说分子筛能按照物质的分子大小进行选择吸附。由于一定型号的分子筛其孔径大小一样,所以一般说来只有比分子筛孔径小的分子才能被分子筛吸附在晶体内部的孔腔内,大于孔径的分子就被筛去。经分子筛干燥后的气体, 一般含水量可达到0.1 10 （- 6）10 10（ - 6） 。具有高效的吸附性。分子筛在低水汽分压、高温、高气体线速度等较苛刻的条件下仍保持较高的湿容量,这是因为分子筛的表面积远大于一般吸附剂,达700900m2/ g。如何处理掉酸性气体:按操作特点、脱酸原理，多种脱酸工艺可分类为：间歇法、化学吸收法、物理吸收法、混合溶剂吸收法、直接氧化法和膜分离法。1、 间歇法 特点：脱酸设备只能批量生产，不能连续生产。 按脱酸气原理可分为：（1）化学反应法：有明显的化学反应，与酸气的反应物不能再生，作为废弃物处置，故仅用于气量小、含酸气浓度低的场合。比如海绵铁法脱除酸气。（2）物理吸附法：没有明显的化学反应，只是单纯的吸附，而且吸附过程可逆。比如用分子筛脱除酸气。2、化学吸收法 在一个塔器内以弱碱性溶液作为吸收剂与酸气反应，生成某种化合物。在另一塔器内，改变工艺条件(加热、降压、汽提等)使化学反应逆向进行，碱性溶液得到再生，恢复对酸气的吸收能力，使天然气脱酸气过程循环连续进行。各种醇胺溶液是化学吸收法中使用最广泛的吸收剂，有一乙醇胺(mea)、二乙醇胺(dea)、二甘醇胺(dga)、甲基二乙醇胺(mdea)和二异丙醇胺(dipa)等。3、物理吸收法 以有机化合物作为溶剂，在高压、低温下使酸气组分溶解于溶剂内。吸收了酸气的溶剂又在低压、高温下释放酸气，使溶剂恢复对酸气的吸收能力，使脱酸过程循环持续进行。物理溶剂再生时所需的加热量较少，适用于天然气内酸气负荷高，要求同时进行天然气脱水的场合。如冷甲醇法。4、混合溶剂吸收法 以物理溶剂和化学溶剂配制的混合溶剂作为吸收剂，兼有物理吸收和化学吸收剂作用。如：砜胺法。5、直接氧化法 对h2s直接氧化使其转化成元素硫，如：claus(克劳斯)法。 在天然气工业中常用于天然气脱出酸气的处理，适合于处理流量小、酸气浓度很高的原料气 。6、膜分离法 是利用气体中各组分通过薄膜渗透性能的区别，将某种气体组分从气流中分离和提浓，从而达到天然气脱酸性气的目的。适用于从天然气内分出大量co2的场合。 h2s的腐蚀机理:温度影响硫化铁膜的成分。通常，在室温下的湿h2s气体中，钢铁表面生成的是无保护性的fe9s8。在100含水蒸气的h2s中，生成的也是无保护性的fe9s8和少量fes。在饱和h2s水溶液中，碳钢在50%下生成的是无保护性的fe9s8和少量的fes；当温度升高到100150时，生成的是保护性较好的fes和fes2。1暴露时间在硫化氢水溶液中，碳钢和低合金钢的初始腐蚀速率很大，约为0.7mm/a，但随着时间的增长，腐蚀速率会逐渐下降2碳钢和低合金钢在含h2s流体中的腐蚀速率，通常是随着时间的增长而逐渐下降，平衡后的腐蚀速率均很低，这是相对于流体在某特定的流速下而言的。3氯离子在酸性油气田水中，带负电荷的氯离子，基于电