CN105806825B - 一种天然气组分在线气体拉曼分析方法 - Google Patents

一种天然气组分在线气体拉曼分析方法 Download PDF

Info

Publication number
CN105806825B
CN105806825B CN201610329999.3A CN201610329999A CN105806825B CN 105806825 B CN105806825 B CN 105806825B CN 201610329999 A CN201610329999 A CN 201610329999A CN 105806825 B CN105806825 B CN 105806825B
Authority
CN
China
Prior art keywords
spectrum
gas
sample
component
raman
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610329999.3A
Other languages
English (en)
Other versions
CN105806825A (zh
Inventor
戴连奎
高颖
陈昀亮
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201610329999.3A priority Critical patent/CN105806825B/zh
Publication of CN105806825A publication Critical patent/CN105806825A/zh
Application granted granted Critical
Publication of CN105806825B publication Critical patent/CN105806825B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/65Raman scattering

Landscapes

  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

本发明公开了一种基于拉曼光谱叠加原理的天然气组分在线气体拉曼分析方法,包括以下步骤:配制纯甲烷或纯氮气等基础样本,并配制基础样本与其他气体按一定比例两两混合的建模样本;测量基础样本与建模样本的拉曼光谱,并基于上述已知其组成的双组分气体建模样本,建立一个包含各种纯物质光谱及其参数的光谱数据库;对于组成未知的天然气,先测量其拉曼光谱,并进行光谱预处理,再利用上述建立的数据库和最小二乘法求得各组分的摩尔百分比。本发明方法含有以下特点:检测速度快、准确度较高、重复性好、与样气非接触、现场免维护等,可以对天然气进行在线实时检测。该方法不仅操作简便,检测过程更对样品无损耗,清洁环保。

Description

一种天然气组分在线气体拉曼分析方法
技术领域
本发明属于化学计量学领域,涉及一种气体混合物组成与含量的光谱检测方法,尤其涉及一种天然气各组分含量的在线检测方法。
背景技术
我国既是能源的消费大国,也是能源的生产大国。在目前我国的能源结构中煤炭占据了主要地位,但是煤炭的利用会对环境产生严重的污染。因此新能源的开发迫在眉睫。
天然气作为一种清洁和使用方便的能源,在我国的储量十分可观。应用天然气发电具有低排放、低污染等优点,因此备受各国政府关注,也是我国今后能源发展的一个重要方向。在气体燃料发电中,混合气体成分的实时监测十分重要。原因在于天然气发电厂的发电量取决于蒸汽量,而蒸汽量受气体燃料的流量和组成的影响。
目前多组分气体在线监测的方法主要有气相色谱法、傅立叶红外光谱法等。
气相色谱法是以气体作为流动相的一种色谱分析方法,它利用色谱柱对混合气体进行组分分离,再利用专门的检测器得到组分的色谱信号,最后实现对气体组成的定性和定量分析。它在天然气成分检测方面已发展得十分成熟。在我国天然气成分分析领域也主要采用气相色谱方法,并发布了两个国家标准GB/T 13610‐2014《天然气的组成分析—气相色谱法》和GB/T 10410‐2008《人工煤气和液化石油气常量组分气相色谱分析法》。但由于在线气相色谱法分析周期较长,现场维护工作量大且测量过程复杂,因而限制了其广泛应用。
红外光谱应用广泛,但因其无法检测同核双原子分子,如氮气和氢气等,并且水对其有强大的吸收干扰等,无法满足对天然气等混合气体的在线实时检测要求。
相较于前两种检测方法,用拉曼光谱法来检测混合气体成分有着得天独厚的优势。只需要单个激光源就可以在同一时间检测出样品气体中全部组分,包括不具有红外活性的双原子分子,如氢气、氮气等。该方法检测速度快,且可以获得较高的灵敏度;因其非接触式测量原理,现场免维护。上述优势使得拉曼光谱法尤其适用于天然气成分在线实时检测。
发明内容
本发明提供了一种天然气组分在线气体拉曼分析方法,目的是实现快速、简便、准确地测定天然气的组分和含量。
该气体拉曼分析方法是基于拉曼光谱叠加原理来建立混合气体组分定量分析模型,如
式中,I(v)表示混合气体的拉曼光谱,ki表示在单位摩尔百分比下第i种组分的拉曼光谱的强度系数,ci表示第i种组分的摩尔百分比,Ii(v)表示第i种组分的归一化后的纯物质拉曼光谱,N表示组分数。
在应用该模型求解未知混合气体的组分和含量ci时,ki和Ii(v)为已知量,但由于ki受到激光强度的影响,无法作为一个常数使用,故将(1)式作如下处理,将ki除以k1,得到:
式中,k1表示在单位摩尔百分比下甲烷的拉曼光谱的强度系数,c1表示甲烷的摩尔百分比,I1(v)表示甲烷归一化后的纯物质拉曼光谱。令表示在单位摩尔百分比下i组分相对甲烷的拉曼光谱强度系数,该参数与激光强度等测量条件无关,可作为常数使用,得到:
根据式(1),利用拉曼光谱仪测得组成未知天然气的光谱I(v)和已知组分的纯物质拉曼光谱Ii(v),使用最小二乘回归可求得参数kici,令bi=kici,则再结合各组分的摩尔百分比之和为1,可求得各组分的摩尔百分比ci
为了表示方便,将(4)式简化,并令K1=1,可得
为了应用上述模型,需要根据组成已知的建模样本求取Ii(v)(i=1,……,N)和Ki(i=2,……,N)。首先配制纯甲烷气体样本,测得其归一化拉曼光谱I1(v);再配制甲烷和其他物质m(如乙烷、丙烷等)的双组分气体样本,确保在配制过程中无液化现象,并测量双组分气体样本的光谱I(v);然后利用光谱I(v)扣除光谱I1(v),间接得到纯物质m的归一化拉曼光谱Im(v);最后利用光谱I(v)、I1(v)和Im(v)通过最小二乘法求得回归系数b1和bm,由求得参数Km
本发明涉及的天然气组成在线拉曼分析方法,包括以下步骤:
(1)将一种纯物质作为基础样本;将基础样本与其他纯物质按一定比例两两混合配制若干双组分气体样本,作为建模样本。
(2)使用拉曼光谱仪测量基础样本的拉曼光谱;首先对得到的光谱进行预处理;再对基础样本的光谱进行特征峰归一化,得到光谱I1(v)并将其存入数据库。
(3)基于上述已知其组成的双组分气体建模样本,建立一个包含各种纯物质光谱Ii(v)及其参数Ki的光谱数据库,建立方法如下:
(3.1)使用拉曼光谱仪测量一种双组分气体的拉曼光谱,并对得到的光谱进行预处理。设预处理后的双组分气体光谱为I(v);对光谱I(v)用基础样本光谱I1(v)在其特征谱区进行最小二乘拟合,得到回归系数B;最后从光谱I(v)中扣除BI1(v),并进行特征峰归一化,得到该种双组分气体的另一组分归一化光谱I2(v)。
(3.2)为求取在单位摩尔百分比下另一组分相对于基础样本的拉曼光谱强度系数K2,对光谱I(v)用光谱I1(v)、I2(v)进行最小二乘拟合,得到回归系数b1、b2;再利用双组分气体中基础样本的摩尔百分比c1和另一组分的摩尔百分比c2,得到将光谱I2(v)以及参数K2存入光谱数据库。
(3.3)根据步骤3.1、3.2中方法求得其他种类的双组分气体建模样本中另一组分的光谱Ii(v)及其参数Ki
(4)对于组成未知的天然气样本,首先使用拉曼光谱仪测量其拉曼光谱,并进行预处理,得到光谱I′(v);然后,基于光谱叠加原理,利用所建数据库和最小二乘法得到回归系数bi(i=1,……,N);最后,采用可求得各组分的摩尔百分比。
本发明的有益效果是,相对于气相色谱法,天然气成分在线拉曼光谱分析方法具有检测速度快、重复性好、与样气非接触、现场免维护等优势,可以对天然气进行在线实时检测;相较于红外光谱法,拉曼光谱法具有更广的物质检测范围。该方法不仅操作简便,检测过程更对样品无损耗,清洁环保。
附图说明
图1本发明使用的拉曼光谱测量系统;
图2甲烷预处理归一化后的拉曼光谱;
图3甲烷乙烷双组分气体和甲烷二氧化碳双组分气体的预处理后的拉曼光谱;
图4乙烷归一化后的拉曼光谱;
图5二氧化碳归一化后的拉曼光谱;
图6测试样本预处理后的拉曼光谱。
具体实施方式
以下结合实施例和附图,进一步说明本发明方法。
实施例拉曼光谱测试平台如图1所示,主要光学部件为:激光器、光纤、拉曼探头和光谱仪。上述所说的激光器是中心波长为532nm的激光器;上述所说的拉曼探头是532nm光纤探头;上述所说的光谱仪是Ocean Optics生产的TEC制冷光纤光谱仪。
该系统的测量过程如下:将天然气样气通入密闭的采样管内;激光器发出激光,经过激发光纤传导至拉曼探头,经拉曼探头聚焦后照射样品,产生的拉曼散射光由拉曼探头收集,再经过收集光纤传回光谱仪;拉曼散射光经过光谱仪光栅分光、CCD阵列检测、A/D转换后,将光谱数字信号传输至PC机。
【实施例】
(1)将纯物质甲烷作为基础样本;将甲烷与乙烷、丙烷、异丁烷、正丁烷、氮气、二氧化碳、氢气按一定比例两两混合配制7个双组分气体样本,作为建模样本,建模样本的摩尔百分比如表1所示;
表1:7组建模样本中各成分的摩尔百分比
注:表中数据为各成分的摩尔百分数,单位为%。
(2)使用拉曼光谱仪测量纯甲烷的拉曼光谱;首先对得到的光谱进行预处理,再对基础样本的光谱在波数1501~1600cm-1的范围内进行特征峰归一化,得到光谱I1(v),并将其存入数据库,光谱I1(v)如图2所示。
(3)基于上述已知其组成的双组分气体建模样本,建立一个包含各种纯物质光谱Ii(v)及其参数Ki的光谱数据库,最终求得参数Ki的值如表2所示。下面以组成已知的甲烷乙烷双组分气体和甲烷二氧化碳双组分气体为例,说明光谱数据库的具体建立方法:
(3.1)使用拉曼光谱仪测量甲烷乙烷双组分气体和甲烷二氧化碳双组分气体的拉曼光谱,并对光谱进行预处理,预处理后的双组分气体光谱分别为I(v)和I′(v),如图3所示。然后,对光谱I(v)、I′(v)用I1(v)在波数为甲烷特征谱区1701~1900cm-1的范围内进行最小二乘拟合,得到回归系数B、B′。最后,从光谱I(v)中扣除BI1(v),并在波数为900~1100cm-1范围内进行特征峰归一化,得到乙烷归一化光谱I2(v),如图4所示;从光谱I′(v)中扣除B′I1(v),在波数为1350~1450cm-1范围内进行特征峰归一化,得到二氧化碳归一化光谱I7(v),如图5所示。
(3.2)求取在单位摩尔百分比下乙烷相对甲烷的拉曼光谱强度系数K2:对光谱I(v)在波数为550~2500cm-1的范围内用光谱I1(v)、I2(v)进行最小二乘拟合,得到回归系数b1为0.7532、b2为0.9938;再利用双组分气体中甲烷的摩尔百分比c1和乙烷的摩尔百分比c2,得到将光谱I2(v)以及参数K2存入光谱数据库。
同理,可求取在单位摩尔百分比下二氧化碳相对甲烷的拉曼光谱强度系数K7:对光谱I′(v)在波数为550~2500cm-1的范围内用光谱I1(v)、I7(v)进行最小二乘拟合,得到回归系数b′1为0.6715、b′2为0.9234;再利用该双组分气体中甲烷和二氧化碳的摩尔百分比,得到K7=12.70;将光谱I7(v)以及参数K7存入光谱数据库。
表2:各物质组分所对应的相对拉曼强度参数K的值
组分 CH<sub>4</sub> C<sub>2</sub>H<sub>6</sub> C<sub>3</sub>H<sub>8</sub> iC<sub>4</sub>H<sub>10</sub> nC<sub>4</sub>H<sub>10</sub> N<sub>2</sub> CO<sub>2</sub> H<sub>2</sub>
K<sub>i</sub> K<sub>1</sub> K<sub>2</sub> K<sub>3</sub> K<sub>4</sub> K<sub>5</sub> K<sub>6</sub> K<sub>7</sub> K<sub>8</sub>
数值 1 11.93 15.94 22.50 11.02 11.23 12.70 20.43
(4)将3个组分未知的天然气样本作为测试样本,先采用本发明所提出的拉曼光谱法进行预测,再用标准的色谱分析法进行检验。
(4.1)使用拉曼光谱仪测量测试样本的拉曼光谱,对其进行预处理,得到光谱如图6所示。
(4.2)根据光谱叠加原理,并利用所建数据库和最小二乘法得到回归系数bi(i=1,……,8),将结合各组分的摩尔百分比之和为1,可求得各组分的摩尔百分比测试样本各组分含量的拉曼分析值如表3所示。
表3:用拉曼光谱法测得3组测试样本中各成分的摩尔百分比
注:表中数据为各成分的摩尔百分数,单位为%。
(4.3)为对本发明方法的准确度进行检验,用气相色谱法检测上述3个测试样本各物质的摩尔百分比,所得结果如表4所示,并对两种方法的所得结果做比较,如表5所示。
表4:用气相色谱法测得3组测试样本中各成分的摩尔百分比
注:表中数据为各成分的摩尔百分数,单位为%。
表5:针对测试样本的拉曼光谱法与气相色谱法的分析误差
注:表中数据为各成分的摩尔百分数,单位为%。
比较结果表明:本发明所公开的气体拉曼分析法与气相色谱法测量结果的分析误差均在±0.5%以内,准确度较高。

Claims (1)

1.一种天然气组分在线气体拉曼分析方法,其特征在于,包括以下步骤:
(1)将一种纯物质作为基础样本;将基础样本与其他纯物质按一定比例两两混合配制若干双组分气体样本,作为建模样本;
(2)使用拉曼光谱仪测量基础样本的拉曼光谱;首先对得到的光谱进行预处理;再对基础样本的光谱进行特征峰归一化,得到光谱I1(v)并将其存入数据库;
(3)基于上述已知其组成的双组分气体建模样本,建立一个包含各种纯物质光谱Ii(v)及其参数Ki的光谱数据库,建立方法如下:
(3.1)使用拉曼光谱仪测量一种双组分气体的拉曼光谱,并对得到的光谱进行预处理;设预处理后的双组分气体光谱为I(v);对光谱I(v)用基础样本光谱I1(v)在其特征谱区进行最小二乘拟合,得到回归系数B;最后从光谱I(v)中扣除BI1(v),并进行特征峰归一化,得到该种双组分气体的另一组分归一化光谱I2(v);
(3.2)为求取在单位摩尔百分比下另一组分相对于基础样本的拉曼光谱强度系数K2,对光谱I(v)用光谱I1(v)、I2(v)进行最小二乘拟合,得到回归系数b1、b2;再利用双组分气体中基础样本的摩尔百分比c1和另一组分的摩尔百分比c2,得到将光谱I2(v)以及参数K2存入光谱数据库;
(3.3)根据步骤(3.1)、(3.2)中方法求得其他种类的双组分气体建模样本中另一组分的光谱Ii(v)及其参数Ki
(4)对于组成未知的天然气样本,首先使用拉曼光谱仪测量其拉曼光谱,并进行预处理,得到光谱I′(v);然后,基于光谱叠加原理,N表示组分数,利用所建数据库和最小二乘法得到回归系数bi;最后,采用求得各组分的摩尔百分比。
CN201610329999.3A 2016-05-17 2016-05-17 一种天然气组分在线气体拉曼分析方法 Active CN105806825B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610329999.3A CN105806825B (zh) 2016-05-17 2016-05-17 一种天然气组分在线气体拉曼分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610329999.3A CN105806825B (zh) 2016-05-17 2016-05-17 一种天然气组分在线气体拉曼分析方法

Publications (2)

Publication Number Publication Date
CN105806825A CN105806825A (zh) 2016-07-27
CN105806825B true CN105806825B (zh) 2019-07-23

Family

ID=56451528

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610329999.3A Active CN105806825B (zh) 2016-05-17 2016-05-17 一种天然气组分在线气体拉曼分析方法

Country Status (1)

Country Link
CN (1) CN105806825B (zh)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3060123B1 (fr) * 2016-12-08 2021-12-03 Ifp Energies Now Methode amelioree de mesure en ligne sur les unites en lit mobile simule ou les unites de separation hybrides par lit mobile simule et cristallisation et application au controle et regulation desdites unites.
CN106932378A (zh) * 2017-03-29 2017-07-07 浙江大学 一种基于拉曼光谱的炼厂酸性气成分的在线检测系统和方法
DE102017117709A1 (de) * 2017-08-04 2019-02-21 Carl Zeiss Jena Gmbh Musterbasierte spektroskopische Gasqualitätsüberwachung
US10768115B2 (en) 2017-12-04 2020-09-08 Kaiser Optical Systems Inc. Augmented Raman analysis of a gas mixture
WO2019127352A1 (zh) * 2017-12-29 2019-07-04 深圳达闼科技控股有限公司 基于拉曼光谱的物质识别方法及云端系统
CN108444976B (zh) * 2018-04-26 2019-10-22 浙江大学 一种基于拉曼光谱的天然气热值测量方法
CN108982468B (zh) * 2018-07-13 2020-02-28 浙江大学 一种对氯甲苯中微量杂质的拉曼分析方法
CN111781164B (zh) * 2019-04-04 2023-04-25 武汉米字能源科技有限公司 一种同时分析天然气中酸性组分含量的激光检测方法
CN113466207B (zh) * 2021-06-30 2023-04-14 天津大学 气化焦油监测系统、监测方法及装置
CN113552110A (zh) * 2021-07-16 2021-10-26 中国民航大学 基于拉曼光谱的锂离子电池热失控动态预警系统及方法
CN116008252B (zh) * 2023-02-13 2023-06-06 北京卓立汉光分析仪器有限公司 一种拉曼光谱下的混合物定量分析方法及装置
CN117037935B (zh) * 2023-08-17 2024-02-09 河北敦诚新能源科技有限公司 材料成分确定方法、装置、设备及存储介质

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5754288A (en) * 1995-08-09 1998-05-19 Kyoto Dai-Ichi Kagaku Co., Ltd. Optical measuring method of expiration components
CN101846628A (zh) * 2009-03-27 2010-09-29 上海神开石油化工装备股份有限公司 钻井液的在线拉曼光谱分析方法
CN102590175A (zh) * 2012-02-21 2012-07-18 浙江大学 基于拉曼光谱叠加的甲醇汽油甲醇含量快速测定方法
CN103119448A (zh) * 2010-09-17 2013-05-22 Abbvie公司 用于生物过程操作的拉曼光谱

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5754288A (en) * 1995-08-09 1998-05-19 Kyoto Dai-Ichi Kagaku Co., Ltd. Optical measuring method of expiration components
CN101846628A (zh) * 2009-03-27 2010-09-29 上海神开石油化工装备股份有限公司 钻井液的在线拉曼光谱分析方法
CN103119448A (zh) * 2010-09-17 2013-05-22 Abbvie公司 用于生物过程操作的拉曼光谱
CN102590175A (zh) * 2012-02-21 2012-07-18 浙江大学 基于拉曼光谱叠加的甲醇汽油甲醇含量快速测定方法

Also Published As

Publication number Publication date
CN105806825A (zh) 2016-07-27

Similar Documents

Publication Publication Date Title
CN105806825B (zh) 一种天然气组分在线气体拉曼分析方法
Ying et al. Quantitative analysis of main components of natural gas based on Raman spectroscopy
CN109983338A (zh) 烃燃料气体组成和性质的测量
US9291610B2 (en) Method and apparatus for real-time measurement of fuel gas compositions and heating values
CN103398970B (zh) 一种定性和定量分析食用油并进而检测地沟油的方法
CN102539374B (zh) 一种用于测量煤气成分和热值的方法
CN101210873A (zh) 一种利用太赫兹时域光谱快速检测植物油纯度的方法及设备
CN102778442B (zh) 一种快速鉴别烟用香液料液种类的方法
CN101339150B (zh) 基于介电谱技术测定汽油辛烷值的方法
CN106932378A (zh) 一种基于拉曼光谱的炼厂酸性气成分的在线检测系统和方法
CN102566533B (zh) 一种烟用香精香料调配的在线监控装置及方法
CN102788778B (zh) 一种醇类汽油中甲醇乙醇含量快速测定方法
CN108444976B (zh) 一种基于拉曼光谱的天然气热值测量方法
Jiang et al. Tracing methane dissolved in transformer oil by tunable diode laser absorption spectrum
CN107860743A (zh) 利用反射式近红外光纤探头构建快速预测原油性质的模型的方法及其应用
CN102221533A (zh) 一种基于紫外光谱花生油掺伪定量检测的方法
CN109324014A (zh) 一种自适应的原油性质近红外快速预测方法
CN102914530A (zh) 一种拉曼光谱气体检测系统及其检测方法和应用
CN110414169B (zh) 一种傅立叶红外气测录井方法及其装置
Kinsey et al. Comparison of measurement methods for the characterization of the black carbon emissions from a T63 turboshaft engine burning conventional and Fischer-Tropsch fuels
CN201000429Y (zh) 一种利用太赫兹时域光谱快速检测植物油纯度的设备
CN103134764A (zh) 由透射红外光谱预测原油实沸点蒸馏曲线的方法
CN102928356A (zh) 一种快速测定香精溶剂含量的方法
CN202394096U (zh) 一种烟用香精香料调配的在线监控装置
Li et al. A hard modeling approach to determine methanol concentration in methanol gasoline by Raman spectroscopy

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant