在天然气的管输过程中,当管道中存在液态水时,天然气中的H2S和CO2将会对管道及相关设备造成腐蚀,因此,管输过程中必须特别关注水含量/水露点的控制[1-2]。GB/T 37124-2018《进入天然气长输管道的气体质量要求》给出水露点指标为“水露点应比输送条件下最低环境温度低5 ℃”和“在输送条件下,当管道管顶埋地温度为0 ℃时,水露点应不高于-5 ℃”[3-4]。为使水露点指标符合要求,油气生产单位一般在现场安装在线水露点仪,用以监控天然气的水露点是否达到标准[5]。目前,在线水露点仪的检测原理主要有冷镜法、石英晶体振荡法、激光法、光纤法、电解法以及阻容法,除光纤法外,均有相应的国家或行业标准对其进行规定(见表 1)[6-8]。
不同检测原理的在线水露点仪直接测得的数据不同,只有冷镜法和光纤法是直接测定得到温度值,其余原理均是先测得水含量,再通过不同的换算方法得到温度值。目前的换算方法主要有ISO 18453:2004《天然气水含量与水露点的相关性》中的“改进PR方程”和ASTM D1142-95(2021)《露点温度法测定气体燃料水蒸气含量的标准试验方法》中的水含量与水露点的换算关联式方法[9-10]。我国参考这两者制定了GB/T 22634-2008《天然气水含量与水露点之间的换算》,该标准考虑到我国法规和工业的特殊需求,进行了相应的修改[11]。GB/T 22634-2008中水含量与水露点的关联关系也是基于改进的PR状态方程,并给出了方程的工作范围与不确定度(见表 2)。
ASTM D1142-95(2021)根据IGT研究报告8给出的天然气水含量平衡数据关联式进行水露点与水含量的换算,其计算公式见式(1)。该关联式根据Raoult定律修改完善,除在接近气体临界温度测定的露点外,该方式的准确度能够满足气体燃料工作的要求[10]。
式中:A为与温度有关的常数,mg·MPa/m3;p为压力,MPa;B为与温度有关的常数,mg/m3;ρw为水的质量浓度,mg/m3。
前期调研了中国石油西南油气田公司在线水露点仪的配置情况,并针对不同检测原理报出的结果及换算方法进行了统计(见表 4)。
冷镜法作为水露点检测的仲裁方法相比其他方法具有快速响应和高准确性的优点,根据GB/T 17283-2014规定在-25~5 ℃的测量范围,该原理自动检测准确度可达±1 ℃,手动检测准确度可达±2 ℃[6]。因此,本研究采用便携式冷镜仪测定结果作为标准值,通过测量不同水含量甲烷中水气体标准物质在不同压力下的露点值,与上述两种换算方法得到的温度值进行比较。同时,采用数学软件Matlab通过三次样条插值技术建立一种新型水含量-水露点换算方法,并与前文叙述的两种换算方法进行比较,探讨其适用性。最后,计算不同方法在不同水含量及压力下的平均绝对偏差,给出相应的使用建议,以减少不同换算方法造成的误差。
根据GB/T 17283-2014的规定,经处理的管输天然气水露点范围一般为-25~5 ℃,在相应的气体压力下,水的体积分数为(50~200)×10-6[5]。经调研,生产单位的气质压力多为4~5 MPa。因此,配制了6瓶8 L甲烷中水气体标准物质,水的体积分数为(10~210)×10-6,制备压力为4.0~9.0 MPa(见表 5)。
采用俄罗斯VYMPEL公司生产的mini-SL便携式冷镜露点仪,技术规格见表 6。冷镜法作为水露点检测仲裁方法具有在带压条件下不需要进行换算则可直接得到露点的优点,因此,以该装置的测定值作为准确值。
为了验证该仪器的准确度,将其送至中国测试技术研究院进行验证。使用标准水发生器在常压下发生出不同水含量的标准气,分别通入便携式冷镜露点仪和精密露点仪,再比较这两台仪器所得的数值,结果见表 7。
根据表 7的数据,便携式冷镜露点仪与精密露点仪的差值绝对值基本在1 ℃以内,无超过1.5 ℃的情况,平均差值为0.68 ℃,说明便携式冷镜露点仪的测量准确度能够保证,其实验值能够作为标准值进行研究。
将不同水含量的标准气在不同压力下以稳定流速通入便携式冷镜仪,以露点仪在该压力下初露点和消露点的平均值作为1次测定结果,连续测定7次,以7次的平均值作为最后结果(见表 8)。
将表 8的实验数据作图,结果如图 1所示。由图 1可知,在同等测试压力下,水露点与水含量呈正相关,随着压力的升高,相同水含量标准气测得的水露点值不断升高,其增速逐渐放缓。
采集表 8中水露点数据,并利用Matlab软件采用三次样条插值法进行处理,该技术拟合曲面具有平滑性、高精确度的特点,在许多物理系统中,平滑的变化比直线或折线更符合实际,三次样条插值能较好地模拟这些平滑过程[12-14]。三次样条插值法并不能得到1个单一的数学关系式,而是在每相邻数据点区间上都有不同参数的分段多项式。对于一维三次样条插值,存在一系列的数据点,例如(x0,y0),(x1,y1)…(xn,yn),通过这些点拟合一个平滑的曲线时,整个曲线会被分成n段,每一段都是一个独立的三次多项式,对于每个区间[xi, xi+1], Si(x)定义为式(2):
式中:ai通常是直接通过在xi点上的函数值得出;bi、ci和di需要考虑到整体的连续性和边界条件,通常通过解一个线性方程组来得到,该方程组关联了所有的数据点和插值条件。
每个Si(x)只在其对应的区间[xi, xi+1]内有效,不同区间的多项式Si(x)和Si+1(x)在xi+1处的函数值、一阶导数和二阶导数都应相等,确保整个插值曲线的平滑性。这种方法让插值曲线在整个数据范围内具有很好的平滑性和适应性,适用于需要高精度和平滑曲线表达的多种场合。
对于二维数据(压力、水含量),插值函数会在两个方向拟合不同的关系式,可以通过在两个方向上各自应用三次样条插值来扩展这个方法。这意味着对于二维数据点集需要在x方向和y方向上分别进行三次样条插值。最终得到的插值函数可以表示为式(3),拟合曲面见图 2。
式中:S(x,y)为差值曲面函数,℃;x为水的体积分数,10-6;y为压力,MPa;Sij(x, y)表示一个分段三次多项式,用来对变量水含量与压力进行插值。
根据三次样条插值原理,本次建立的模型不能直接得到数学方程式。因此,采用留一交叉验证法计算该模型的平均绝对误差和均方根误差,以对已建立的插值模型进行评估。留一交叉验证法涉及将数据集分为训练集和测试集,测试集只保留单点,每次移除该点后再用训练集的点建立插值曲面,最后预测被移除那个点的值。重复此过程,以便每个点都被预测过一次,最后计算所有预测值和实际值之间的误差。而平均绝对误差和均方根误差是衡量预测模型或拟合曲面准确性的常用统计量,这两个指标的值越小,表示预测的准确性越高。
该模型的平均绝对误差为1.163 1 ℃,均方根误差为2.278 4 ℃,相比于GB/T 22634-2008与GB/T 17283-2014所规定的水露点值准确度(±2 ℃)而言,该模型的误差是可以接受的。
确定3种换算方法的适用温度与压力区间,如表 9所列。
基于3种不同的换算方法,参考表 9确认换算方法适用范围,分别计算6种不同水含量的标气在不同压力下的换算水露点值,并由式(4)计算平均绝对偏差[15-16],结果见表 10。
式中:D为平均绝对偏差,%;i为第i组数据;N为参与计算的组数;W实验值为采用便携式冷镜仪测得的水露点,K;W计算值为通过不同换算方法得到的水露点,K。
由表 10可知,在3种换算方法当中,采用三次样条插值法进行换算的平均误差最小,只有0.46%,其次是采用GB/T 22634-2008的换算方法的平均误差,为0.69%,最后是采用ASTM D1142-95(2021)的换算方法进行换算的平均误差,为0.84%。
ASTM D1142-95(2021)的换算方法没有给出不确定度,GB/T 17283-2014和GB/T 22634-2008采用ISO 18453:2004的PR方程换算方法的不确定度均为2 ℃,两种换算方法合理可接受的最大偏差应为两者不确定度的平方和根,即为2.8 ℃。将3种换算方法得到的露点值与实验值的差值的绝对值按照不同压力进行作图,结果见图 3。
由图 3可知,采用三次样条插值法建立的换算模型在不同水含量、压力段下均基本满足在最大偏差以内的要求,相比其他两种方法,三次样条插值法的偏差都保持在比较稳定的状态。GB/T 22634-2008方法在高水含量、高压段的平均绝对偏差已经超过了可接受范围,但在水含量较低段的表现比较优异。ASTM D1142-95(2021)在低水含量、高压段的偏差较高,但在其余段都表现比较稳定。因此,根据3种换算方法在不同水含量和压力段的平均绝对偏差进一步评估其适用性,具体数据见表 11。
由表 11可知,从整体趋势上来看,采用三次样条插值法的平均绝对偏差都保持在相对较低的水平并维持稳定,说明该方法建立的模型具有较好的适应性,但其在中水含量、低压段的表现还是不如其余两种方法。在低水含量段,GB/T 22634-2008方法相比ASTM D1142-95(2021)具有更好的表现,且压力越高,其换算误差越低,但是随着水含量的逐渐增加,换算误差逐渐上升。ASTM D1142-95(2021)在低水含量段的换算误差比较大,但是随着水含量的不断上升,其误差出现逐渐降低的趋势。根据对表 11所列数据的分析,得出在不同条件下使用的换算方法(见表 12),以供工程人员使用,以减少因不同换算方法所带来的误差。
1) 目前,天然气中水含量与水露点的换算关系主要参考的是GB/T 22634-2008与ASTM D1142-95(2021),这两个标准在目的和方法上侧重点不同,GB/T 22634-2008需要求解天然气气液两相平衡状态方程而计算水露点,而ASTM D1142-95(2021)则提供了一种换算公式求解水露点,该方法虽不知准确度,但计算速度快,永远有解。
2) 基于三次样条插值法建立了一种水含量与水露点的换算方法,并将该方法与标准方法进行对比,三者平均误差从低至高分别为0.46%(三次样条插值法)、0.69%(GB/T 22634-2008)和0.84%(ASTM D1142-95(2021))。相比其他两种标准方法,该方法在一定的温压范围内具有更好的准确度和稳定性。因此,该方法可作为一种新型快速的水含量、水露点换算方法。
3) 通过对比3种换算方法在不同水含量及压力段的平均绝对偏差,给出了3种换算方法的建议使用区间,可供工程人员进行参考,以降低因不同换算方法而造成的误差。