《优秀毕业设计精品天然气输送过程中水合物的生成防治研究和讨论》由会员分享，可在线阅读，更多相关《优秀毕业设计精品天然气输送过程中水合物的生成防治研究和讨论（63页珍藏版）》请在金锄头文库上搜索。
1、1 绪 论1.1研究的意义和目的随着石油天然气工业的不断发展，在处理和输送天然气过程中发现了气体水合物。水合物是目前科学领域中的热门课题，不仅与石油天然气开采、储存和运输密切相关，而且与环境保护、气候变迁，特别是人类未来赖以生存的能源有关。天然气水合物为白色结晶固体，是在一定温度、压力条件下、天然气中的烃分子与其中的游离水结合而形成的，其中水分子靠氢键形成一种带有大、小孔穴的结晶晶格体，这些孔穴被小的气体分子所充填。在天然气管道输送过程中，水合物在输气干线或输气站某些管段（弯头）阀门、节流装置等处形成后，其流通面积减少从而形成局部堵塞，其上游的压力增大，流量减少，下游的压力降低，因为会影响管道
2、输配气的正常运行。天然气水合物是威胁输气管道安全运行的一个重要因素。天然气水合物一旦形成后，它与金属结合牢固，会减少管道的流通面积，产生节流，加速水合物的进一步形成，进而造成管道、阀门和一些设备的堵塞，严重影响天然气的开采、集输和加工的正常运行。因此，研究和讨论天然气输送过程中水合物的生成和防治，对保障天然气管道的安全运行具有十分重要的实际意义。要形成天然气水合物需要几个必要的条件，一是气体处于水汽的饱和或过饱和状态并存在游离水；二是有足够高的压力和足够地的温度。在具备上述条件时，水合物的形成，还要求有一些辅助条件，如天然气压力的波动，气体因流向的突变而产生的搅动，以及晶种的存在等。因此总结出
3、一些防治天然气水合物生成的方法。通常，在输送天然气过程中清除水合物的方法是用热水或热蒸汽对管道进行加热，在水合物和金属接触点上，将温度提高到3040，使水合物很快分解。据统计防止水合物生成的费用约占生产总成本的58%。在工程上对抑制剂用量不能准确计算，抑制剂的用量往往大于实际需求量，这样一方面不利于节约成本，另一方面导致不必要的环境污染针对上述问题，需要用科学的实验方法，准确测定天然气水合物的生成条件，并筛选和评价抑制剂的抑制效果，从而为天然气集输管道水合物防治工作提供科学依据。1.2国内外发展概况及分析国内外预测水合物形成压力和温度的方法大致可分为经验公式法、图解法、相平衡计算法和统计热力学
4、模型法四大类。经验公式法是通过室外取样室内实验或现场测试水合物形成压力和温度,在此基础上,拟合得到的半经验或经验模型。经验、半经验模型具有简单、计算方便等优点,但使用范围狭窄,计算精度不高。经验图解法是根据甲烷及不同相对密度天然气形成水合物的平衡曲线，大致确定出天然气形成水合物的温度和压力。但对含有硫化氢的天然气误差较大，若相对密度在两条曲线之间,需采用内插法进行近似计算。为了便于计算机运算，图解法被回归成几个公式，若已知天然气的相对密度和温度,可选择其中合适的公式计算水合物形成压力。若已知相对密度和压力可选择其中合适的公式进行迭代求得水合物形成温度。但是这种方法也存在适用范围狭窄，计算精度不
5、高的缺点。相平衡计算法是1940年Katz根据气固平衡常数,提出了一种估算天然气水合物生成条件的方法,可用于计算含有典型烷烃组成的无硫天然气,而对非烃含量多的气体及压力高于6.9MPa时误差较大。水合物的相平衡理论模型大多是在Van der Waals-Platteeuw模型的基础上发展起来的,具有很强的理论基础,计算精度高,但计算较复杂,难掌握。目前预测水合物生成条件的热力学模型几乎都是以经典统计热力学为基础的。Vander Waals和Platteeuw 在1959年提出的 Langmuir气体吸附模型，没考虑到客体分子间的相互作用，而假定洞穴为球形对称等。1964年，Saito等人提出了
6、一种预测水合物生成条件的方法。1972年，Parrish和Pransnitz考虑到Langmuir参数仅为温度的函数，提出了简单的经验关系式计算Langmuir常数，大大简化了Vander Waals模型，并首次将之推广到多元体系。1977年Ng&Robinson模型(简称NR)，考虑不同客体间的相互作用，修正Van Der Waals模型。1983年HolderJohn模型（简称HJ）考虑到洞穴并非球形，水分子与洞穴中心距离不等的实际情况，采用三层球模型来描述水分子和客体分子间的相互作用，并引入扰动因子来矫正球形分子的Langmuir常数；1988年杜亚和和郭天民（简称Du&Guo88模型）
7、利用十点Gauss数值积分法计算出不同温度下Langmuir常数；1996年，Chen- Guo基于水合物生成动力学，提出了一个完全不同于VdW- P模型的新水合物预测模型。他们认为，在生成水合物时，体系存在准化学平衡和物理吸附平衡。左有祥、郭天民基于水合物两步生成动力学机理提出了新的水合物热力学模型，与传统vdW-P模型相比，表达式得到简化且更接近普通的溶液热力学模型，该预测模型都己成功应用于纯水体系、含醇或含盐2 天然气状态方程2.1理想气体状态方程理想气体是为了研究问题和工程计算的方便而设想的一个理想化模型，这种气体的分子本身不占有体积，分子之间没有内聚力，其状态方程式为： (2.1)式
8、中 p气体压力，N/m2; V理想气体体积，m3;v1mol理想气体占有的体积，m3/mol；n气体物质的量，mol；R理想气体常数，8.3143kJ/(kmolK)；T气体温度，K。2.2实际气体状态方程当压力升高或密度增大时，气体分子本身占据的体积大得越来越严重，而分子间的相互作用力也变得越来越明显。未考虑这些效应，范德瓦尔在1837年提出了另一个状态方程，即范德瓦尔状态方程。在范德瓦尔状态方程提出之后，为了考虑实际气体的性质，又提出了大量的状态方程。这里就介绍RK方程、SRK方程和PR方程。 RK方程RK（Redlitch-Kwong）方程是1949年提出的二参数状态方程，基本形式如下：
9、 (2.2)式中a和b是常数，对单组份有：混合气体的a、b常数由各组分的a、b常数及摩尔分数按混合规则求解。 SRK方程SRK是在改进Redlitch-Kwong状态方程基础上，提出的一种状态方程，其形式为： (2.3)式中a和b是常数：上式中与气体对比温度Tr和偏心因子有关，按下式计算：偏心因子根据有关手册的基础物性数据表差得。 PR方程PR是另外一种状态方程，其形式为： (2.4)式中a和b是常数：3 输气管道的水力计算、热力计算3.1输气管道的水力计算根据稳定流动的假设，运动方程整理可得： (3.1)3.1.1输气管的基本公式所谓水平输气管道，是指地形起伏高差小于200m的管道。这种管道
10、阴有雨高差所引起的压降在上式中所占的比例很小，并不足以影响输气管道水力计算的准确性，故可忽略不计，于是方程可写为： （3.2）公式中的和w为随管长而变化的气体的密度和流速；为消去这两个变数，我们可利用连续性方程和气体状态方程。由气体状态方程有： （3.3）由连续性方程有： （3.4）将式（3.3）和（3.4）代入式（3.2）中得：若管路全线长度为L，其起终点压力分别为和，则对式积分可求得输气管流量和压降的关系式：积分得： （3.5）式中 输气管质量流量，kg/s；输气管内径，m；，输气管起、终点压力，MPa；水力摩擦系数，无因次；气体压缩因子，无因次；气体常数，；输气温度，K；管长，m。式（3
11、.5）中，一项表示随着压力下降，流速增加对流量的影响。对长距离管道而言，与相比很小，可以忽略。相反距离短、压降大的管道则必须考虑。故对长输管道而言，上式可简化为： (3.6)这是平坦地区输气管道的质量流量公式。但在工程设计和生产上通常采用的是在标准状况下的体积流量。因此，必须把质量流量M换算成标况下的体积流量。 (3.7)式中为标况下的气体密度。由气体状态方程可知：标况下的气体压缩系数Z0=1，故： （3.8）为计算方便，用空气的气体常数Ra来表示天然气的气体常数R：将代入上式得： （3.9）将（3.9）代入（3.8）得 （3.10）再将（3.7）、（3.8）和（3.10）代入（3.6），得：
12、 (3.11)设，于是 （3.12）式中 Q天然气在标况下的体积流量，m3/s; C常数； Ra空气的气体常数，m2/(s2K)； 天然气的相对密度，无因次。3.1.2输气管道压力分布和平均压力沿线压力分布设有一段输气管道AC长为L。以x表示管段上任一点B至起点A的距离，见图3-1，图3-1 输气管沿线任一点压力则公式（3.12）改写为： (3.13)根据上式，对AB段可写为：对BC段可写为：同一条管路流量相同，因此以上两式相等整理后得： （3.14）平均压力当输气管道停输时，管道内的压力并不像输油管道那样立刻消失，而是仍处于压力状态下，高压端的气体逐渐流向低压端。这样起点压力PQ逐渐下降，而
13、低压端因由高压气体流入，终点压力PZ逐渐上升，最后两段压力达到某个平均值PCP,即平均压力。利用公式（3.14），按管道的全长积分，即可求得输气管道的平均压力：积分并整理的： (3.15)3.2输气管道的热力计算3.2.1温降基本公式由能量方程推导温降公式。稳定流动的能量方程: (3.16)忽略高差和速度变化的影响，则： (3.17)式中 h气体的焓，J/kg; T气体的温度，K或； P气体的压力，MPa; Q单位质量气体通过管长向外所放出的热量，J/kg; x沿管长变化量，m。另外由热力学知识可知：因此： (3.18)由于：式中 CP气体质量定压热容，J/kgK； Di焦耳汤姆逊系数，/MPa。则： (3.19)表示单位质量气体在单位管长上的热力损失，由传热学关系可知： (3.20)式中 K管道的总传热系数，W/(m2K)； D管道的内径，m； T0管道埋深处地温，K； M气体质量流量，kg/s。管道队环境放出的热量应等于气体所放出的热量，因此：令 ，则上式变为： (3.21)非齐次线性微分方程，其通解为：由于x=0时，T=TQ，所以： (3.22)
