发明内容
基于此,有必要针对上述技术问题,提供一种同时分析天然气中酸性组分含量的激光检测方法,本方法选取天然气中酸性气体吸收光谱重合的光谱波段,通过化学计量学算法实现酸性气体组分的同时测量,解决了传统的激光吸收光谱技术往往选取没有干扰的吸收波段只能测量一种气体的问题,也解决了使用多台分析仪表占地面积更大而且增加了分析仪表的投入和维护成本的问题。
为了解决上述技术问题,本发明提供一种同时分析天然气中酸性组分含量的激光检测方法,采用了如下所述的技术方案:
一种同时分析天然气中酸性组分含量的激光检测方法,其包括以下步骤:
(1)建立i组酸性气体的自变量组,每一组自变量组对应一种酸性气体;每组自变量组建立过程如下:向检测池中通入预设浓度的一酸性气体,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组Xi=(xi1,xi2,...,xin),i代表第i种酸性气体,i≥2,Xi代表第i种酸性气体的自变量组;
(2)建立j组混合气的因变量组,每一组因变量组对应一种由i种酸性气体按配比混合的混合气,每种混合气的第i种酸性气体的浓度Ci是不同的;每组因变量组建立过程如下:向检测池中通入混合气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Yj=(yj1,yj2,...,yjn),j代表代表第j种混合气,j≥2,Yj代表第j种混合气的因变量组;
(3)建立线性回归方程:
Yj=∑ikjiXi+kj0 (1),其中,
kji代表在第j组混合气的因变量组中第i种酸性气体吸收曲线的影响系数,kj0代表在第j组混合气的因变量组中一个偏移量;
根据上述每一组线性回归方程,通过多元线性回归算法得到影响系数kj1、kj2…kji和常数kj0的最优解;
(4)影响系数kji与第j种混合气中的第i种酸性气体浓度Cji呈正比,即第j种混合气中第i种酸性气体浓度Cji与影响系数kji为线性关系,通过一元线性回归算法拟合求解得到回归系数ai和常数bi,以建立浓度Ci与影响系数ki的线性函数:
Ci=aiki+bi (2),其中,
Ci代表混合气中第i种酸性气体的浓度,ki代表对第i种酸性气体吸收曲线的影响系数,ai代表回归系数,bi为常数;
(5)根据待测天然气建立一组待测天然气的因变量组Y’,建立过程如下:向检测池中通入待测天然气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Y’=(y’1,y’2,...,y’n);根据方程(1)建立Y’=∑ik’iXi+k’0,通过多元线性回归算法得到影响系数k’i和k’0;将得到的影响系数k’i代入方程(2),则可求得待测天然气中第i种酸性气体的浓度含量。
一种同时分析天然气中酸性组分含量的激光检测方法,当所述待测天然气中具有在预设波段内对酸性气体的吸收曲线存在干扰的干扰气,则检测方法包括以下步骤:
(1)建立i组酸性气体的自变量组,每一组自变量组对应一种酸性气体;每组自变量组建立过程如下:向检测池中通入预设浓度的一酸性气体,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组Xi=(xi1,xi2,...,xin),i代表第i种酸性气体,i≥2,Xi代表第i种酸性气体的自变量组;
(2)建立m组干扰气的自变量组,过程如下:向检测池中通入预设浓度的一干扰气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X’m=(x’m1,x’m2,...,x’mn),m代表不同的干扰气,m≥1,X’m代表第m种干扰气的自变量组;
(3)建立j组混合气的因变量组,每一组因变量组对应一种由i种酸性气体和m种干扰气按配比混合的混合气,每种混合气的第i种酸性气体的浓度Ci和第m种干扰气的浓度Cm均是不同的;每组混合气的因变量组建立过程如下:向检测池中通入混合气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Yj=(yj1,yj2,...,yjn),j代表代表第j种混合气,j≥3,Yj代表第j种混合气的因变量组;
(4)建立线性回归方程:
Yj=∑ikjiXi+∑mKjmX’m+kj0 (3),其中,
kji、Kjm分别代表在第j组混合气的因变量组中第i种酸性气体吸收曲线和第m种干扰气的影响系数,kj0代表在第j组混合气的因变量组中一个偏移量;
根据上述每一组线性回归方程,通过多元线性回归算法得到影响系数kj1、kj2…kji、Kj1、Kj2…Kjm和常数kj0的最优解;
(5)影响系数kji、Kjm分别与第j种混合气中的第i种酸性气体浓度Cji和第m种干扰气浓度C’m呈正比,即:
(5.1)第j种混合气中第i种酸性气体浓度Cji与影响系数kji为线性关系,通过一元线性回归算法拟合求解得到回归系数ai和常数bi,以建立浓度Ci与影响系数ki的线性函数:
Ci=aiki+bi (2),其中,
Ci代表混合气中第i种酸性气体的浓度,ki代表对第i种酸性气体吸收曲线的影响系数,ai代表回归系数,bi为常数;
(5.2)第j种混合气中第m种干扰气浓度C’jm与影响系数Kjm为线性关系,通过一元线性回归算法拟合求解得到回归系数a’m和常数b’m,以建立浓度C’m与影响系数Km的线性函数:
C’m=a’mKm+b’m (4),其中,
C’m代表混合气中第m种干扰气的浓度,k’i代表对第m种干扰气吸收曲线的影响系数,a’m代表回归系数,b’m为常数;
(6)根据待测天然气建立一组待测天然气的因变量组Y’,建立过程如下:向检测池中通入待测天然气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Y’=(y’1,y’2,...,y’n);根据方程(3)建立Y’=∑ik’iXi+∑mK’mXm+k’0,通过多元线性回归算法得到影响系数k’i、K’m和k’0;将得到的影响系数k’i代入步骤(5.1)中的方程(2),则可求得待测天然气中第i种酸性气体的浓度含量。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述酸性气体至少包括硫化氢和二氧化碳。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述预设波段内的硫化氢和二氧化碳吸收特征峰相隔不超过0.1nm。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述干扰气为甲烷、乙烷的至少一种。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述预设波段为激光器输出波长控制在1569-1571nm,1571-1573nm,1572-1574nm,1574-1575nm,1576-1577nm和1577-1579nm中的任一个波段范围。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述激光器为窄波段可调谐激光光源。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述窄波段可调谐激光光源为分布式反馈激光器、垂直腔面发射激光器或量子级联激光器。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,所述检测池中的背景气为与待测天然气的背景气相同或为对待测酸性气体吸收光谱无干扰的气体。
作为本发明提供的同时分析天然气中酸性组分含量的激光检测方法的一种改进,建立自变量组(Xi,X’m)和建立因变量组(Yj,Y’)步骤中的检测池的压力相同。
与现有技术相比,本发明有以下有益效果:
1、本发明选取天然气中H2S和CO2等酸性气体吸收光谱重合的光谱波段,通过化学计量学算法实现酸性气体组分的同时测量,解决了传统的激光吸收光谱技术往往选取没有干扰的吸收波段只能测量一种气体的问题,也解决了使用多台分析仪表占地面积更大而且增加了分析仪表的投入和维护成本的问题。。
2、本发明的检测方法具有检测速度快、准确度较高、非接触、免维护、寿命长等优点。
3、本方法能够使用一套分析技术、一个激光光源解决两种组分的在线测量,减少了气体在线分析仪器数量,降低了投入成本。
4、相较于传统激光吸收光谱技术采用的差分算法或二次谐波高度算法,利用多元线性回归算法能够分解复杂背景气条件下混叠的吸收峰,抗干扰能力强;此外,由于个别位置上的波动对线型整体影响较小,可以减小随机误差对浓度测量结果的影响,具有更高的稳定性。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”、“顺时针”、“逆时针”、“轴向”、“径向”、“周向”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是至少两个,例如两个,三个等,除非另有明确具体的限定。
在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接或彼此可通讯;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通或两个元件的相互作用关系,除非另有明确的限定。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。
如背景技术所描述的,传统的激光吸收光谱技术往往选取没有干扰的吸收波段,故只能测量一种气体。这就造成了H2S和CO2分析仪无法合并测量的现状,使得天然气净化厂及门站中往往需要购置两台甚至多台分析仪,不仅占地面积更大,而且增加了分析仪表的投入和维护成本。
为了解决上述问题,本发明提出了基于激光技术同时对天然气国标要求的H2S和CO2等酸性气体进行测量,具体地,选取酸性气体吸收光谱重合的光谱波段,通过化学计量学算法实现酸性气体组分的同时测量。
更优选地,是选取H2S和CO2吸收光谱重合的光谱波段,通过化学计量学算法实现H2S和CO2组分的同时测量。下面以此为例进一步说明,即i=2(H2S和CO2这2种酸性气体)的情况。
当背景气对H2S和CO2吸收光谱无干扰时,将H2S和CO2重叠的吸收光谱分解为已知单组份,即H2S和CO2独立吸收光谱分量的加权和形式。
在HITRAN数据库中查阅H2S和CO2吸收光谱,主要集中在1570-1580nm之间,如图1所示。根据混合气体组分的浓度配比、测量量程要求,考虑可能存在的背景气干扰等因素,选择合适的激光器波长。在本发明中H2S和CO2同时测量激光波长选取范围包括:1569-1571nm,1571-1573nm,1572-1574nm,1574-1575nm,1576-1577nm,1577-1579nm等的任一个,但不局限于此。其中,上述激光器应为窄波段可调谐激光光源,如分布式反馈激光器(DFB,Distributed Feedback Laser)、垂直腔面发射激光器(VCSEL,Vertical Cavity SurfaceEmitting Laser)或量子级联激光器(QCL,Quantum Cascade Laser)等。
由于半导体激光器属于窄线宽光源,通过控制激光器的驱动电流,可以改变激光器的发光频率,一般而言驱动电流为锯齿波,最大电流与最小电流之差为10mA,可以将激光器发光频率改变0.1nm。控制激光器的温度,通过调整激光器锯齿波驱动电流的高低位电流,使得激光器发出的激光波长正好完整扫过选取的吸收光谱波段,为了能够同时覆盖到H2S和CO2的吸收光谱,H2S和CO2吸收特征峰相隔最好不超过0.1nm。
确定激光器的波段后,则可对天然气中H2S和CO2进行同步分析。
本发明提供了一种同时分析天然气中H2S和CO2含量的激光检测方法,其包括以下步骤:
(1.1)建立2组酸性气体的自变量组,分别对应H2S和CO2;每组酸性气体的自变量组建立过程如下:
(1.1.1)向检测池中通入参考气a:浓度为α的H2S,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X1=(x11,x12,...,x1n),n为大于1的自然数;
(1.1.2)向检测池中通入参考气b:浓度为β的CO2,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X2=(x21,x22,...,x2n),n为大于1的自然数;
(1.2)建立j组混合气的因变量组,每一组因变量组对应一种由H2S和CO2按配比混合的混合气,每种混合气的H2S和CO2的浓度配比是不同的;每组混合气的因变量组建立过程如下:向检测池中通入混合气,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下所述预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Yj=(yj1,yj2,...,yjn),j代表第j种混合气,每种混合气具有不同浓度配比,j为大于1的自然数,Yj代表第j种混合气的因变量组;
(1.3)建立线性回归方程:
Yj=kj1X1+kj2X2+kj0 (5),其中,
kj1、kj2分别代表在第j组混合气的因变量组中H2S和CO2气体吸收曲线的影响系数,kj0代表在第j组混合气的因变量组中一个偏移量;
根据上述每一组线性回归方程,通过多元线性回归算法得到影响系数kj1、kj2和常数kj0的最优解;
(1.4)影响系数kj1、kj2分别与第j种混合气中的H2S和CO2气体浓度Cj1和Cj2呈正比,即第j种混合气中H2S和CO2气体浓度Cj1和Cj2分别与影响系数kj1、kj2为线性关系,通过一元线性回归算法拟合求解得到回归系数a1、a2和常数b1、b2,以建立浓度与影响系数的线性函数:
C1=a1k1+b1 (6),
C2=a2k2+b2 (7);其中,
C1代表混合气中H2S气体的浓度,k1代表对H2S气体吸收曲线的影响系数,C2代表混合气中CO2气体的浓度,k2代表对CO2气体吸收曲线的影响系数;
(1.5)根据待测天然气建立一组待测天然气的因变量组Y’,建立过程如下:向检测池中通入待测天然气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Y’=(y’1,y’2,...,y’n);根据线性回归方程(5)建立Y’=k’1X1+k’2X2+k’0,通过多元线性回归算法得到影响系数k’1、k’2和k’0;将得到的影响系数k’1、k’2分别代入方程(6)、(7),则可求得待测天然气中H2S和CO2的实际浓度含量。
值得注意的是,本发明涉及到的多元线性回归算法、一元线性回归算法可采用现有的标准算法,也可采用其他回归算法,在此不再赘述。
优选地,本发明中,吸收光谱曲线优选进行归一化处理后再分割建立自变量组,其中,归一化处理可采用现有常规的归一处理方法,在此不再详述。
在具体实现中,如果背景气体中某些组分在所选波段内对H2S和CO2等酸性气体的吸收曲线存在干扰,即所述待测气体中具有在预设波段内对H2S和CO2等酸性气体的吸收曲线存在干扰的干扰气,如甲烷、乙烷等。为了避免背景气中干扰气波动带来的误差,在检测方法中增设干扰气的自变量组X’,再建立因变量组与自变量组的函数关系Y=F(X,X’),通过多元线性回归算法求得影响系数和常数的最优解。具体地,以m=1(1种干扰气如甲烷)为例进一步说明,所述检测方法包括以下步骤:
(2.1)建立2组酸性气体的自变量组,分别对应H2S和CO2;每组酸性气体的自变量组建立过程如下:
(2.1.1)向检测池中通入参考气a:浓度为α的H2S,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X1=(x11,x12,...,x1n),n为大于1的自然数;
(2.1.2)向检测池中通入参考气b:浓度为β的CO2,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X2=(x21,x22,...,x2n),n为大于1的自然数;
(2.2)建立1组干扰气的自变量组,过程如下:向检测池中通入参考气c:浓度为ν的干扰气,背景气应与待测天然气的背景气相同或为N2、H2等对待测酸性气体气体光谱无干扰的气体,记录下所述预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立自变量组X’1=(x’11,x’12,...,x’1n),n为大于1的自然数;
(2.3)建立j组混合气的因变量组,每一组混合气的因变量组对应一种由H2S、CO2和干扰气按配比混合的混合气,每种混合气中H2S、CO2和干扰气的浓度配比是不同的;每组混合气的因变量组建立过程如下:向检测池中通入混合气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Yj=(yj1,yj2,...,yjn),j代表具有不同配比的混合气,j为大于2的自然数,Yj代表第j种混合气的因变量组;
(2.4)建立线性回归方程:
Yj=kj1X1+kj2X2+Kj1X’1+kj0 (8),其中,
kj1、kj2和Kj1分别代表在第j组混合气的因变量组中H2S、CO2和干扰气的影响系数,kj0代表在第j组混合气的因变量组中一个偏移量;
根据上述每一组线性回归方程,通过多元线性回归算法得到影响系数kj1、kj2、Kj1和常数kj0的最优解;
(2.5)影响系数kj1、kj2和Kj1分别与第j种混合气中的H2S浓度Cj1、CO2浓度Cj2和干扰气浓度C’j1呈正比;具体地,
(2.5.1)第j种混合气中H2S和CO2气体浓度Cj1和Cj2分别与影响系数kj1、kj2为线性关系,通过一元线性回归算法拟合求解得到回归系数a1、a2和常数b1、b2,以建立浓度与影响系数的线性函数:
C1=a1k1+b1 (9),
C2=a2k2+b2 (10);其中,
C1代表混合气中H2S气体的浓度,k1代表对H2S气体吸收曲线的影响系数,C2代表混合气中CO2气体的浓度,k2代表对CO2气体吸收曲线的影响系数;
(2.5.2)第j种混合气中干扰气气体浓度C’j1与影响系数Kj1为线性关系,通过一元线性回归算法拟合求解得到回归系数a’1和常数b’1,以建立浓度与影响系数的线性函数:
C’1=a’1K1+b’1 (11);
C’1代表混合气中干扰气气体的浓度,K1代表对干扰气吸收曲线的影响系数;
(2.6)根据待测天然气建立一组待测天然气的因变量组Y’,建立过程如下:向检测池中通入待测天然气,所述检测池中的背景气与待测天然气的背景气相同或为对待测酸性气体的吸收光谱无干扰的气体,记录下预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Y’=(y’1,y’2,...,y’n);根据方程(8)建立Y’=k’1X1+k’2X2+K’1X’1+k’0,通过多元线性回归算法得到影响系数k’1、k’2、K’1和k’0;将得到的影响系数k’1、k’2分别代入方程(9)、(10),则可求得待测天然气中H2S和CO2的实际浓度含量。
需要说明的是,改变检测池的压力,吸收光谱也会随之发生变化,可以记录不同压力下的自变量组曲线,即,在建立自变量组时,事先记录多种压力下的自变量组曲线,建立不同压力下的线性回归方程。当实际测量过程中,检测池压力改变时,选取同种压力下的自变量曲线进行拟合运算。
为了使本技术领域的人员更好地理解本申请方案,下面将结合附图,对本申请实施例中的技术方案进行清楚、完整地描述。
实施例1
以天然气中H2S和CO2同时分析测量为例,假设二类管道天然气内H2S含量为20ppmv(同μL/L、mg/m3),CO2含量为1%。
预设波段选取1574.5-1574.7nm范围,即记录此波段范围的吸收光谱。
以下是在相同的检测池压力下进行测量。
(1)向检测池中通入20ppmv的H2S气体,背景气为N2,记录下归一后的吸收光谱曲线如图2所示,按照图2的吸收光谱曲线,将光谱强度分割成80个独立的波长与相对吸收光谱强度的关系,建立自变量组X1=(x11,x12,...,x180),x11,x12,...,x180分别与表1中编号1到80的相对吸收光谱强度相对应。
表1 20ppmv的H2S吸收强度与波长的对应关系
(2)向检测池中通入1%的CO2气体,背景气为N2,记录下归一后的吸收光谱曲线如图3所示;按照图3的吸收光谱曲线,将光谱强度分割成80个独立的波长与相对吸收光谱强度的关系,建立自变量组X2=(x21,x22,...,x280),x21,x22,...,x280分别与
表2中编号1到80的相对吸收光谱强度相对应。
表2 1%的CO2吸收强度与波长的对应关系
(3)以j=3为例,调配三种同时含有不同浓度配比的H2S和CO2的混合气,三种混合气中H2S和CO2浓度配比分别为20ppmvH2S和1%CO2、15ppmvH2S和1.5%CO2以及10ppmvH2S和2%CO2。
分别向检测池中通入三种混合气,记录下归一后的吸收光谱曲线如图4所示;按照图4的吸收光谱曲线,将每一组光谱强度分割成80个独立的波长与相对吸收光谱强度的关系,建立因变量组Y1=(y11,y12,...,y180)、Y2=(y21,y22,...,y280)、Y3=(y31,y32,...,y380)。
(4)建立线性回归方程:
Y1=k11X1+k12X2+k10,通过多元线性回归算法得到影响系数k11、k12和常数k10的最优解;
Y2=k21X1+k22X2+k20,通过多元线性回归算法得到影响系数k21、k22和常数k20的最优解;
Y3=k31X1+k32X2+k30,通过多元线性回归算法得到影响系数k31、k32和常数k30的最优解。
(5)影响系数kj1(k11、k12和k13)分别与三种混合气中的H2S气体浓度Cj1(20ppmv、15ppmv和10ppmv)呈正比,即混合气中H2S气体浓度Cj1与影响系数kj1为线性关系,通过一元线性回归算法拟合求解得到回归系数a1和常数b1,以建立H2S气体浓度C1与影响系数k1的线性函数:
C1=a1k1+b1 (12),其中,
C1代表混合气中H2S气体的浓度,k1代表对H2S气体吸收曲线的影响系数,a1代表回归系数,b1为常数;
(6)影响系数k21(k21、k22和k23)分别与三种混合气中的CO2气体浓度Cj2(1%、1.5%和2%)呈正比,即混合气中CO2气体浓度Cj2与影响系数kj2为线性关系,通过一元线性回归算法拟合求解得到回归系数a2和常数b2,以建立CO2气体浓度C2与影响系数k2的线性函数:
C2=a2k2+b2 (13),其中,
C2代表混合气中CO2气体的浓度,k2代表对CO2气体吸收曲线的影响系数,a2代表回归系数,b2为常数;
(7)在实际分析待测天然气时,根据待测天然气建立一组待测天然气的因变量组Y’,建立过程如下:向检测池中通入待测天然气,背景气为N2,记录下归一后预设波段内的吸收光谱,将预设波段分割成n个波长,每个波长对应的光谱强度建立因变量组Y’=(y’1,y’2,...,y’n);建立方程Y’=k’1X1+k’2X2+k’0,通过多元线性回归算法得到影响系数k’1、k’2和k’0;将得到的影响系数k’1、k’2分别代入方程(12)、(13),则可求得待测天然气中H2S和CO2的实际浓度含量C1和C2。
需要注意的是,本实施例1中的步骤(5)(6)属于标定过程,其只是建立影响系数ki与浓度的线性Ci关系,在量程范围内至少2组配气就能够完成标定,但为了在量程内保证更高的线性度和准确性,标定时配越多不同浓度配比的混合气,越能够消除偶然误差的影响,标定效果更好。
显然,以上所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例,附图中给出了本申请的较佳实施例,但并不限制本申请的专利范围。本申请可以以许多不同的形式来实现,相反地,提供这些实施例的目的是使对本申请的公开内容的理解更加透彻全面。尽管参照前述实施例对本申请进行了详细的说明,对于本领域的技术人员来而言,其依然可以对前述各具体实施方式所记载的技术方案进行修改,或者对其中部分技术特征进行等效替换。凡是利用本申请说明书及附图内容所做的等效结构,直接或间接运用在其他相关的技术领域,均同理在本申请专利保护范围之内。