CN113495096A - 一种基于fid-tvoc检测仪的校准方法 - Google Patents
一种基于fid-tvoc检测仪的校准方法 Download PDFInfo
- Publication number
- CN113495096A CN113495096A CN202010563543.XA CN202010563543A CN113495096A CN 113495096 A CN113495096 A CN 113495096A CN 202010563543 A CN202010563543 A CN 202010563543A CN 113495096 A CN113495096 A CN 113495096A
- Authority
- CN
- China
- Prior art keywords
- concentration
- calibration
- value
- detector
- data
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 230000004044 response Effects 0.000 claims abstract description 41
- 230000009467 reduction Effects 0.000 claims abstract description 11
- 238000011088 calibration curve Methods 0.000 claims abstract description 6
- 238000005070 sampling Methods 0.000 claims description 36
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 32
- 230000003321 amplification Effects 0.000 claims description 26
- DPQUFPIZKSPOIF-UHFFFAOYSA-N methane propane Chemical compound C.CCC.CCC DPQUFPIZKSPOIF-UHFFFAOYSA-N 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000006641 stabilisation Effects 0.000 claims description 3
- 238000011105 stabilization Methods 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 8
- 239000007789 gas Substances 0.000 description 37
- 238000005516 engineering process Methods 0.000 description 5
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000001257 hydrogen Substances 0.000 description 2
- 229910052739 hydrogen Inorganic materials 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000012855 volatile organic compound Substances 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 150000002894 organic compounds Chemical class 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 238000007670 refining Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
- G01N27/626—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode using heat to ionise a gas
Landscapes
- Chemical & Material Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明提供了一种基于FID‑TVOC检测仪的校准方法,包括如下步骤:S1:采集标准气的响应信号数据,进行降噪处理;S2:使用原标准曲线对标准气进行响应预估;S3:计算响应信号数据的校准系数;S4:用校准系数校准原标准曲线的数据;S5:使用校准后的标曲数据重新进行拟合。采用本发明的校准方法校准,仪器的检测精度较高,运行稳定性好。且采用的校准曲线可以经常进行校正,对于不同的FID检测器以及同一检测器不同的工作环境,可以经常使用标气进行校正,确保标准曲线是适用于当前的FID检测器工作环境,保证准确度。
Description
技术领域
本发明涉及FID技术,具体涉及一种基于FID-TVOC检测仪的校准方法。
背景技术
FID技术即Flame Ionization Detector技术,全称火焰离子化检测器技术。是以氢气和空气燃烧产生的火焰为能源,当样品进入火焰时,在高温下发生电离产生大量离子,在高压定向电场作用下形成离子流,经过高阻放大形成与样品浓度成正比的电信号的检测技术。
TVOC即Total Volatile Organic Compound,全称总挥发性有机物。在《炼油与石油化学工业大气污染物排放标准》(DB11/447-2007)中定义为:在20℃条件下蒸气压大于或等于0.01kPa,或者特定适用条件下具有相应挥发性的全部有机化合物的统称。
FID检测器的线性特征为分段线性,即在一定量程范围内具有线性特征。FID检测器的理论线性响应范围为1E7,但由于空气流量控制、氢气流量控制、进样控制、检测器温度控制等控制环节存在***误差,因此在大量程范围上,响应与浓度之间的线性关系往往不够好,导致检测结果的精度不理想。
目前,常见的TVOC-FID检测仪器在检测时,往往只使用单一浓度标气进行标定,只能保证检测值在标定标气的浓度附近有较好精确度,在样品浓度与标定标气浓度差距较大时,检测数据精确性不佳。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种基于FID-TVOC检测仪的校准方法,此方法使用多点拟合建立校准曲线,通过信号的分段采样、校准,实现高精度和高稳定性。本发明的技术方案如下:
一种基于FID-TVOC检测仪的校准方法,包括如下步骤:
S1:采集标准气的响应信号数据,进行降噪处理;
S2:使用原标准曲线对标准气进行响应预估;
S3:计算响应信号数据的校准系数;
S4:用校准系数校准原标准曲线的数据;
S5:使用校准后的标曲数据重新进行拟合。
可选地,所述步骤S1进一步包括:
S11:定义FID-TVOC检测仪的性能系数C;
S12:选取一定浓度的甲烷丙烷1∶1标准气通入FID-TVOC检测仪,用于校准;
S13:所述检测仪以100Hz的频率对信号进行采样,所述检测仪设置有放大倍数不同的两个采样电路:高放大倍数采样电路、低放大倍数采样电路;
前述两个采样电路每秒各采集总计w个信号数据,分别得到各自的数据集合:sig(A)和sig(B);所述低放大倍数采样电路的量程大于所述高放大倍数采样电路的量程;所述低放大倍数采样电路的量程为U1;其中,w为不小于10的正整数;
S14:每秒对sig(A)和sig(B)分别进行相同的降噪处理,分别得到sig(A)和sig(B)各自对应的算术平均值:ksigA,ksigB;
S15:记所述检测仪输出的信号响应为Kout,当ksigA<(U1/2)时,Kout=ksigA;否则Kout=ksigB。
可选地,所述步骤S11中,C的计算公式如下:
C=E*E/U,其中,U为所述检测仪的最大量程信号值,单位为pA;E为所述检测仪的最大示值误差,单位为%;其中,U=U1。
可选地,步骤S14的“降噪处理”包括:
去掉sig(A)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigA;
去掉sig(B)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigB。
可选地,所述步骤S2进一步包括:
S21:所述检测仪为已被标定过的检测仪,其自身预先对应有原标准曲线,所述原标准曲线记录有m个标曲数据点,按浓度已被分为两组;所述m个标曲数据点的第s个标曲数据点的浓度值为Is,s=1,2,…,m;m为正整数;
分组标准为出厂时规定的浓度值Id,浓度值Is小于等于Id的标曲数据点分入低浓度组记为Plow,浓度值Is大于Id的标曲数据点分入高浓度组Phigh;
将所述m个标曲数据点按浓度值Is从低到高排列,记为P=[p1,p2,…,pm],其中ps=[Is,Ms];Ms为第s个标曲数据点的信号值,Is为第s个标曲数据点的浓度值;
对Plow和Phigh的标曲数据点各自使用最小二乘法进行线性方程拟合,得到两条曲线:
记高浓度组的标曲线性方程为LineA:
y=aLA*(x-bLA);
记低浓度组的标曲线性方程为LineB:
y=aLB*(x-bLB);
以下值的含义为:
aLA——FID-TVOC检测仪的高放大倍数采样电路响应系数,单位为ppm/pA;
bLA——FID-TVOC检测仪的高放大倍数采样电路零点响应,单位为pA;
aLB——FID-TVOC检测仪的低放大倍数采样电路响应系数,单位为ppm/pA;
bLB——FID-TVOC检测仪的低放大倍数采样电路零点响应,单位为pA;
x——FID-TVOC检测仪检测到的样气响应,单位为pA;
y——FID-TVOC检测仪输出的样气浓度,单位为ppm;
其中,aLA、bLA、aLB、bLB,为通过最小二乘法拟合线性次方程得到的值;
对拟合得到的两条曲线LineA和LineB,计算其交点z作为分段点,记为pz=[Iz,Mz];Iz,Mz分别为分段点z的浓度值和信号值;
S22:FID-TVOC检测仪存储步骤S1中所输出的Kout,连续取最新的30个Kout,记为[k0 1,k0 30],计算[k0 1,k0 30]的均方根值(RMSkout)和算数术平均值(AVRkout),当RMSkout≤|AVRkout|*C时,认为信号值已稳定;
S23:记步骤S12通入的甲烷丙烷1∶1标准气浓度值为I0,将步骤S22中信号稳定时的AVRkout作为信号值,记为M0,将[I0,M0]记为标曲数据点P0;
S24:记本次校准的标曲为:
y=aLi*(x-bLi);
当I0大于Iz时,认为本次校准使用的标准气为高浓度标气,校准原标准曲线[pz,pm]部分数据,aLi=aLA,bLi=bLA;否则,认为本次校准使用的标准气为低浓度标气,校准原标准曲线[p1,pz]部分数据,aLi=ALB,bLi=bLB。
可选地,步骤S3进一步包括:
S31:计算得标气的实际去零点响应O0:
O0=M0-bLi;
计算得原标准曲线下的预计去零点响应O1:
O1=(I0/aLi+bLi);
S32:计算校准数据用的校准系数J:
J=O0/O1。
可选地,步骤S4进一步包括:
S41:第n个被校准的数据点记为Pn,其对应有初始浓度值In和初始信号值Mn;当输入的标气浓度I0低于Iz时,1≤n≤z;当输入的标准气I0浓度高于Iz时,z<n≤m;
S42:计算校准后原标准曲线上各标曲数据点的响应部分M’n,计算方法如下:
当被校准数据点为[P1,Pz]部分时,
M’n=(Mn-bLi)*J+bLi,1≤n≤z;
M’n=Mn,z<n≤m;
当被校准数据点为[Pz,Pm]部分时,
M’n=Mn,1≤n≤z;
M’n=(Mn-bLi)*J+bLi,z<n≤m;
记校准后原标准曲线上标曲数据点的集合为P’,经过校准后的标曲数据点P’n为(In,M’n)。
可选地,步骤S5进一步包括:
拟合方法同步骤S21,以出厂时规定的浓度值Id为界,校准后的m个标曲数据点按浓度分组后,分别得到高浓度组和低浓度组,对应高、低浓度组分别使用最小二乘法进行线性拟合,获得高低浓度两条曲线。
与现有技术相比,本发明具有如下的有益效果:
(1)采用此数据处理方法校准,仪器的检测精度较高,运行稳定性好
(2)采用的校准曲线可以经常进行校正,对于不同的FID检测器以及同一检测器不同的工作环境,可以经常使用标准气进行校正,确保标准曲线是适用于当前的FID检测器工作环境,保证准确度。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1是本发明具体实施例的一种基于FID-TVOC检测仪的校准方法的流程图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1,本实施例公开了一种基于FID-TVOC检测仪的校准方法,包括如下步骤:
S1:采集标准气的响应信号数据,进行降噪处理;
S2:使用原标准曲线对标准气进行响应预估;
S3:计算响应信号数据的校准系数;
S4:用校准系数校准原标准曲线的数据;
S5:使用校准后的标曲数据重新进行拟合。
其中,所述步骤S1进一步包括:
S11:定义FID-TVOC检测仪的性能系数C;
S12:选取一定浓度的甲烷丙烷1∶1标准气通入FID-TVOC检测仪,用于校准;需要说明的是,选取的甲烷丙烷1∶1标准气的碳计浓度,应当与实际被测气的碳计浓度相近,以获得最好的校准效果。
S13:所述检测仪以100Hz的频率对信号进行采样,所述检测仪设置有放大倍数不同的两个采样电路:高放大倍数采样电路、低放大倍数采样电路;前述两个采样电路每秒各采集总计w个信号数据,分别得到各自的数据集合:sig(A)和sig(B);
所述低放大倍数采样电路的量程大于所述高放大倍数采样电路的量程;所述低放大倍数采样电路的量程为U1;其中,w为不小于10的正整数;本实施例中,w=100。
本实施例中,涉及的采样电路均包括信号采集电路、降噪电路、信号放大电路、模数转换电路。具体应用时按照选取的模数转换芯片、所设计的电路,w可作调整,以获得较好的采集效果。FID-TVOC检测仪内部自设置有两个采样电路,属于现有技术范畴,故这里不再对采样电路的具体结构做进一步的赘述。
S14:每秒对sig(A)和sig(B)分别进行相同的降噪处理,分别得到sig(A)和sig(B)各自对应的算术平均值:ksigA,ksigB;
S15:记所述检测仪输出的信号响应为Kout,当ksigA<(U1/2)时,Kout=ksigA;否则Kout=ksigB。
其中,所述步骤S11中,C的计算公式如下:
C=E*E/U,其中,U为所述检测仪的最大量程信号值,单位为pA;E为所述检测仪的最大示值误差,单位为%;其中,U=U1。
其中,步骤S14的“降噪处理”包括:
去掉sig(A)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigA;
去掉sig(B)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigB。
其中,本实施例中,ksigA,ksigB均为一组频率为1Hz的高低浓度信号数据。
其中,所述步骤S2进一步包括:
S21:所述检测仪为已被标定过的检测仪,其自身预先对应有原标准曲线,所述原标准曲线记录有m个标曲数据点,按浓度已被分为两组;所述m个标曲数据点的第s个标曲数据点的浓度值为Is,s=1,2,…,m;m为正整数;
分组标准为出厂时规定的浓度值Id,浓度值Is小于等于Id的标曲数据点分入低浓度组记为Plow,浓度值Is大于Id的标曲数据点分入高浓度组Phigh;
将所述m个标曲数据点按浓度值Is从低到高排列,记为P=[p1,p2,…,pm],其中ps=[Is,Ms];Ms为第s个标曲数据点的信号值,Is为第s个标曲数据点的浓度值;
对Plow和Phigh的标曲数据点各自使用最小二乘法进行线性方程拟合,得到两条曲线:
记高浓度组的标曲线性方程为LineA:
y=aLA*(x-bLA);
记低浓度组的标曲线性方程为LineB:
y=aLB*(x-bLB);
以下值的含义为:
aLA——FID-TVOC检测仪的高放大倍数采样电路响应系数,单位为ppm/pA;
bLA——FID-TVOC检测仪的高放大倍数采样电路零点响应,单位为pA;
aLB——FID-TVOC检测仪的低放大倍数采样电路响应系数,单位为ppm/pA;
bLB——FID-TVOC检测仪的低放大倍数采样电路零点响应,单位为pA;
x——FID-TVOC检测仪检测到的样气响应,单位为pA;
y——FID-TVOC检测仪输出的样气浓度,单位为ppm;
其中,aLA、bLA、aLB、bLB,为通过最小二乘法拟合线性次方程得到的值;
对拟合得到的两条曲线LineA和LineB,计算其交点z作为分段点,记为pz=[Iz,Mz];
Iz,Mz分别为分段点z的浓度值和信号值;
S22:FID-TVOC检测仪存储步骤S1中所输出的Kout,连续取最新的30个Kout,记为[k0 1,k0 30],计算[k0 1,k0 30]的均方根值(RMSkout)和算数术平均值(AVRkout),当RMSkout≤|AVRkout|*C时,认为信号值已稳定;
S23:记步骤S12通入的甲烷丙烷1∶1标准气浓度值为I0(为已知),将步骤S22中信号稳定时的AVRkout作为信号值,记为M0,将[I0,M0]记为标曲数据点P0;
S24:记本次校准的标曲为:
y=aLi*(x-bLi);
当I0大于Iz时,认为本次校准使用的标准气为高浓度标气,校准原标准曲线[pz,pm]部分数据,aLi=aLA,bLi=bLA;否则,认为本次校准使用的标准气为低浓度标气,校准原标准曲线[p1,pz]部分数据,aLi=ALB,bLi=bLB。本次校准使用的标准气,即为步骤S12通入的甲烷丙烷1∶1标准气;
其中,步骤S3进一步包括:
S31:计算得标气的实际去零点响应O0:
O0=M0-bLi;
计算得原标准曲线下的预计去零点响应O1:
O1=(I0/aLi+bLi);
S32:计算校准数据用的校准系数J:
J=O0/O1。
其中,步骤S4进一步包括:
S41:第n个被校准的数据点记为Pn,其对应有初始浓度值In和初始信号值Mn;当输入的标气浓度I0低于Iz时,1≤n≤z;当输入的标准气I0浓度高于Iz时,z<n≤m;
S42:计算校准后原标准曲线上各标曲数据点的响应部分M’n,计算方法如下:
当被校准数据点为[P1,Pz]部分时,
M’n=(Mn-bLi)*J+bLi,1≤n≤z;
M’n=Mn,z<n≤m;
当被校准数据点为[Pz,Pm]部分时,
M’n=Mn,1≤n≤z;
M’n=(Mn-bLi)*J+bLi,z<n≤m;
记校准后原标准曲线上标曲数据点的集合为P’,经过校准后的标曲数据点P’n为(In,M’n)。
其中,步骤S5进一步包括:
拟合方法同步骤S21,以出厂时规定的浓度值Id为界,校准后的m个标曲数据点按浓度分组后,分别得到高浓度组和低浓度组,对应高、低浓度组分别使用最小二乘法进行线性拟合,获得高低浓度两条曲线。
与现有技术相比,本发明具有如下的有益效果:
(1)采用此数据处理方法校准,仪器的检测精度较高,运行稳定性好
(2)采用的校准曲线可以经常进行校正,对于不同的FID检测器以及同一检测器不同的工作环境,可以经常使用标准气进行校正,确保标准曲线是适用于当前的FID检测器工作环境,保证准确度。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (8)
1.一种基于FID-TVOC检测仪的校准方法,其特征在于,包括如下步骤:
S1:采集标准气的响应信号数据,进行降噪处理;
S2:使用原标准曲线对标准气进行响应预估;
S3:计算响应信号数据的校准系数;
S4:用校准系数校准原标准曲线的数据;
S5:使用校准后的标曲数据重新进行拟合。
2.如权利要求1所述的方法,其特征在于,所述步骤S1进一步包括:
S11:定义FID-TVOC检测仪的性能系数C;
S12:选取一定浓度的甲烷丙烷1∶1标准气通入FID-TVOC检测仪,用于校准;
S13:所述检测仪以100Hz的频率对信号进行采样,所述检测仪设置有放大倍数不同的两个采样电路:高放大倍数采样电路、低放大倍数采样电路;
前述两个采样电路每秒各采集总计w个信号数据,分别得到各自的数据集合:sig(A)和sig(B);所述低放大倍数采样电路的量程大于所述高放大倍数采样电路的量程;所述低放大倍数采样电路的量程为U1;其中,w为不小于10的正整数;
S14:每秒对sig(A)和sig(B)分别进行相同的降噪处理,分别得到sig(A)和sig(B)各自对应的算术平均值:ksigA,ksigB;
S15:记所述检测仪输出的信号响应为Kout,当ksigA<(U1/2)时,Kout=ksigA;否则Kout=ksigB。
3.如权利要求2所述的方法,其特征在于,所述步骤S11中,C的计算公式如下:
C=E*E/U,其中,U为所述检测仪的最大量程信号值,单位为pA;E为所述检测仪的最大示值误差,单位为%;其中,U=U1。
4.如权利要求2所述的方法,其特征在于,步骤S14的“降噪处理”包括:
去掉sig(A)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigA;
去掉sig(B)中w个信号数据中的最大值和最小值后,计算剩下的信号数据的算术平均值,记为ksigB。
5.如权利要求2所述的方法,其特征在于,所述步骤S2进一步包括:
S21:所述检测仪为已被标定过的检测仪,其自身预先对应有原标准曲线,所述原标准曲线记录有m个标曲数据点,按浓度已被分为两组;所述m个标曲数据点的第s个标曲数据点的浓度值为Is,s=1,2,…,m;m为正整数;
分组标准为出厂时规定的浓度值Id,浓度值Is小于等于Id的标曲数据点分入低浓度组记为Plow,浓度值Is大于Id的标曲数据点分入高浓度组Phigh;
将所述m个标曲数据点按浓度值Is从低到高排列,记为P=[p1,p2,…,pm],其中ps=[Is,Ms];Ms为第s个标曲数据点的信号值,Is为第s个标曲数据点的浓度值;
对Plow和Phigh的标曲数据点各自使用最小二乘法进行线性方程拟合,得到两条曲线:
记高浓度组的标曲线性方程为LineA:
y=aLA*(x-bLA);
记低浓度组的标曲线性方程为LineB:
y=aLB*(x-bLB);
以下值的含义为:
aLA——FID-TVOC检测仪的高放大倍数采样电路响应系数,单位为ppm/pA;
bLA——FID-TVOC检测仪的高放大倍数采样电路零点响应,单位为pA;
aLB——FID-TVOC检测仪的低放大倍数采样电路响应系数,单位为ppm/pA;
bLB——FID-TVOC检测仪的低放大倍数采样电路零点响应,单位为pA;
x——FID-TVOC检测仪检测到的样气响应,单位为pA;
y——FID-TVOC检测仪输出的样气浓度,单位为ppm;
其中,aLA、bLA、aLB、bLB,为通过最小二乘法拟合线性次方程得到的值;
对拟合得到的两条曲线LineA和LineB,计算其交点z作为分段点,记为pz=[Iz,Mz];Iz,Mz分别为分段点z的浓度值和信号值;
S22:FID-TVOC检测仪存储步骤S1中所输出的Kout,连续取最新的30个Kout,记为[k0 1,k0 30],计算[k0 1,k0 30]的均方根值(RMSkout)和算数术平均值(AVRkout),当RMSkout≤|AVRkout|*C时,认为信号值已稳定;
S23:记步骤S12通入的甲烷丙烷1∶1标准气浓度值为I0,将步骤S22中信号稳定时的AVRkout作为信号值,记为M0,将[I0,M0]记为标曲数据点P0;
S24:记本次校准的标曲为:
y=aLi*(x-bLi);
当I0大于Iz时,认为本次校准使用的标准气为高浓度标气,校准原标准曲线[pz,pm]部分数据,aLi=aLA,bLi=bLA;否则,认为本次校准使用的标准气为低浓度标气,校准原标准曲线[p1,Pz]部分数据,aLi=ALB,bLi=bLB。
6.如权利要求5所述的方法,其特征在于,步骤S3进一步包括:
S31:计算得标气的实际去零点响应O0:
O0=M0-bLi;
计算得原标准曲线下的预计去零点响应O1:
O1=(I0/aLi+bLi);
S32:计算校准数据用的校准系数J:
J=O0/01。
7.如权利要求6所述的方法,其特征在于,步骤S4进一步包括:
S41:第n个被校准的数据点记为Pn,其对应有初始浓度值In和初始信号值Mn;当输入的标气浓度I0低于Iz时,1≤n≤z;当输入的标准气I0浓度高于Iz时,z<n≤m;
S42:计算校准后原标准曲线上各标曲数据点的响应部分M’n,计算方法如下:
当被校准数据点为[P1,Pz]部分时,
M’n=(Mn-bLi)*J+bLi,1≤n≤z;
M’n=Mn,z<n≤m;
当被校准数据点为[Pz,Pm]部分时,
M’n=Mn,1≤n≤z;
M’n=(Mn-bLi)*J+bLi,z<n≤m;
记校准后原标准曲线上标曲数据点的集合为P’,经过校准后的标曲数据点P’n为(In,M’n)。
8.如权利要求7所述的方法,其特征在于,步骤S5进一步包括:
拟合方法同步骤S21,以出厂时规定的浓度值Id为界,校准后的m个标曲数据点按浓度分组后,分别得到高浓度组和低浓度组,对应高、低浓度组分别使用最小二乘法进行线性拟合,获得高低浓度两条曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010563543.XA CN113495096B (zh) | 2020-06-18 | 2020-06-18 | 一种基于fid-tvoc检测仪的校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010563543.XA CN113495096B (zh) | 2020-06-18 | 2020-06-18 | 一种基于fid-tvoc检测仪的校准方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113495096A true CN113495096A (zh) | 2021-10-12 |
CN113495096B CN113495096B (zh) | 2024-04-09 |
Family
ID=77994867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010563543.XA Active CN113495096B (zh) | 2020-06-18 | 2020-06-18 | 一种基于fid-tvoc检测仪的校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113495096B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114019114A (zh) * | 2022-01-06 | 2022-02-08 | 深圳普门科技股份有限公司 | 标准曲线生成方法、装置、分析设备和可读存储介质 |
CN114397396A (zh) * | 2021-12-29 | 2022-04-26 | 杭州春来科技有限公司 | 一种VOCs检测分析用自动校准方法及*** |
Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030066803A1 (en) * | 2001-10-05 | 2003-04-10 | Wright Adrian Geoffrey | Analytical instrument calibration and data correction |
US20040110306A1 (en) * | 1998-07-08 | 2004-06-10 | Beimesch Wayne Edward | Method for measuring volatile organic compounds and a kit for same |
US20050026294A1 (en) * | 2003-07-29 | 2005-02-03 | Barber Alan C. | Devices and methods for determining and predicting breakthrough times and steady state permeation rates of organics |
JP2007205968A (ja) * | 2006-02-03 | 2007-08-16 | Dkk Toa Corp | 水素炎イオン化形分析計 |
JP2014002082A (ja) * | 2012-06-20 | 2014-01-09 | Chikuno Shokuhin Kogyo Kk | 揮発性有機化合物の定量方法 |
JP2014025854A (ja) * | 2012-07-27 | 2014-02-06 | Apurikusu:Kk | 化学分析装置の校正方法及び化学分析装置 |
WO2016120762A1 (en) * | 2015-01-26 | 2016-08-04 | Dh Technologies Development Pte. Ltd. | Automatic quantitative regression |
CN106018223A (zh) * | 2016-05-24 | 2016-10-12 | 深圳市蜂联科技有限公司 | 一种分段最优标定提高空气质量检测设备测量精度的方法 |
CN107014785A (zh) * | 2017-05-15 | 2017-08-04 | 浙江全世科技有限公司 | 一种发射光谱背景校正的改进方法 |
CN107064281A (zh) * | 2017-04-14 | 2017-08-18 | 中国石油化工股份有限公司 | 废气中总烃在线稀释监测方法 |
CN108982728A (zh) * | 2018-05-30 | 2018-12-11 | 上海市环境监测中心 | 一种环境空气非甲烷总烃在线监测仪器优化与校准方法 |
CN109283511A (zh) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | 一种宽覆盖多波束接收基阵校准方法 |
CN109725275A (zh) * | 2018-12-26 | 2019-05-07 | 刘科宏 | 基于近场多探头天线测量***的探头校准方法 |
CN109870239A (zh) * | 2019-03-12 | 2019-06-11 | 北京环境特性研究所 | 非制冷红外焦平面探测器自适应定标方法 |
CN110007015A (zh) * | 2019-02-13 | 2019-07-12 | 华电智控(北京)技术有限公司 | 一种总挥发性有机物监测装置 |
CN110376324A (zh) * | 2019-07-12 | 2019-10-25 | 北京雪迪龙科技股份有限公司 | 利用氢火焰离子化检测器测定氧浓度的方法及气相色谱仪 |
US20200057024A1 (en) * | 2018-08-17 | 2020-02-20 | Horiba, Ltd. | Gas analysis apparatus and gas analysis method |
CN111141812A (zh) * | 2019-12-31 | 2020-05-12 | 中国地质调查局西安地质调查中心 | 用于榍石U-Pb定年数据校正方法、存储介质及校正设备 |
-
2020
- 2020-06-18 CN CN202010563543.XA patent/CN113495096B/zh active Active
Patent Citations (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040110306A1 (en) * | 1998-07-08 | 2004-06-10 | Beimesch Wayne Edward | Method for measuring volatile organic compounds and a kit for same |
US20030066803A1 (en) * | 2001-10-05 | 2003-04-10 | Wright Adrian Geoffrey | Analytical instrument calibration and data correction |
US20050026294A1 (en) * | 2003-07-29 | 2005-02-03 | Barber Alan C. | Devices and methods for determining and predicting breakthrough times and steady state permeation rates of organics |
JP2007205968A (ja) * | 2006-02-03 | 2007-08-16 | Dkk Toa Corp | 水素炎イオン化形分析計 |
JP2014002082A (ja) * | 2012-06-20 | 2014-01-09 | Chikuno Shokuhin Kogyo Kk | 揮発性有機化合物の定量方法 |
JP2014025854A (ja) * | 2012-07-27 | 2014-02-06 | Apurikusu:Kk | 化学分析装置の校正方法及び化学分析装置 |
WO2016120762A1 (en) * | 2015-01-26 | 2016-08-04 | Dh Technologies Development Pte. Ltd. | Automatic quantitative regression |
CN106018223A (zh) * | 2016-05-24 | 2016-10-12 | 深圳市蜂联科技有限公司 | 一种分段最优标定提高空气质量检测设备测量精度的方法 |
CN107064281A (zh) * | 2017-04-14 | 2017-08-18 | 中国石油化工股份有限公司 | 废气中总烃在线稀释监测方法 |
CN107014785A (zh) * | 2017-05-15 | 2017-08-04 | 浙江全世科技有限公司 | 一种发射光谱背景校正的改进方法 |
CN108982728A (zh) * | 2018-05-30 | 2018-12-11 | 上海市环境监测中心 | 一种环境空气非甲烷总烃在线监测仪器优化与校准方法 |
US20200057024A1 (en) * | 2018-08-17 | 2020-02-20 | Horiba, Ltd. | Gas analysis apparatus and gas analysis method |
CN109283511A (zh) * | 2018-09-01 | 2019-01-29 | 哈尔滨工程大学 | 一种宽覆盖多波束接收基阵校准方法 |
CN109725275A (zh) * | 2018-12-26 | 2019-05-07 | 刘科宏 | 基于近场多探头天线测量***的探头校准方法 |
CN110007015A (zh) * | 2019-02-13 | 2019-07-12 | 华电智控(北京)技术有限公司 | 一种总挥发性有机物监测装置 |
CN109870239A (zh) * | 2019-03-12 | 2019-06-11 | 北京环境特性研究所 | 非制冷红外焦平面探测器自适应定标方法 |
CN110376324A (zh) * | 2019-07-12 | 2019-10-25 | 北京雪迪龙科技股份有限公司 | 利用氢火焰离子化检测器测定氧浓度的方法及气相色谱仪 |
CN111141812A (zh) * | 2019-12-31 | 2020-05-12 | 中国地质调查局西安地质调查中心 | 用于榍石U-Pb定年数据校正方法、存储介质及校正设备 |
Non-Patent Citations (2)
Title |
---|
李相贤;罗桂山;徐亮;高闽光;童晶晶;魏秀丽;刘娜;: "δ~(13)CO_2值傅里叶变换红外光谱检测标准尺度校准方法", 红外与激光工程, no. 10, 25 October 2015 (2015-10-25), pages 2959 - 2964 * |
郭旭东 等: ""六种塑料制品中挥发性有机物分析研究"", 《当代化工》, vol. 47, no. 4, pages 809 - 814 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114397396A (zh) * | 2021-12-29 | 2022-04-26 | 杭州春来科技有限公司 | 一种VOCs检测分析用自动校准方法及*** |
CN114397396B (zh) * | 2021-12-29 | 2023-08-01 | 杭州春来科技有限公司 | 一种VOCs检测分析用自动校准方法及*** |
CN114019114A (zh) * | 2022-01-06 | 2022-02-08 | 深圳普门科技股份有限公司 | 标准曲线生成方法、装置、分析设备和可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN113495096B (zh) | 2024-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113495096B (zh) | 一种基于fid-tvoc检测仪的校准方法 | |
US8197133B2 (en) | System and method for sensor thermal drift offset compensation | |
CN103412202B (zh) | 一种电流分区间校准的方法 | |
US7426848B1 (en) | Gas composition sensing using carbon nanotube arrays | |
TWI442041B (zh) | 氣體分析器、氣體分析器修正方法、及氣體分析器用修正程式之記憶媒體 | |
CN105572191B (zh) | 一种电化学气体传感器的压力补偿方法 | |
CN103475369B (zh) | 基于信号源误差一次性校准识别的高精度adc测试方法 | |
Snowden-Ifft et al. | High precision measurements of carbon disulfide negative ion mobility and diffusion | |
CN111737912A (zh) | 一种基于深度神经网络的mwhts模拟亮温计算方法 | |
US9806685B1 (en) | Trans-impedance amplifier with increased dynamic range | |
CN104950037A (zh) | 一种针对挥发性有机物的在线污染源识别监测方法及*** | |
KR102007241B1 (ko) | 표준 절연유의 가스 분석값 보정 시스템 및 방법 | |
Manard et al. | Differential mobility spectrometry/mass spectrometry: The design of a new mass spectrometer for real-time chemical analysis in the field | |
CN112666241B (zh) | 一种区域tvoc监测方法 | |
CN102445608A (zh) | 一种电能质量监测装置及校准方法 | |
US11402275B2 (en) | Sensor apparatus | |
CN104458630A (zh) | 一种紫外差分气体分析仪的数据处理方法及*** | |
CN114002378B (zh) | 气体浓度传感器的浓度检测方法 | |
EP3250933A1 (en) | Correction circuits for successive-approximation-register analog-to-digital converters | |
CN112836166B (zh) | 基于响应峰分析的等间采样实验监测数据一阶微分算法 | |
CN115469249A (zh) | 一种考虑激励电流影响的霍尔传感器温漂修正方法 | |
CN110542535B (zh) | 脉冲风洞多传感器负压测试***归一化方法 | |
JP2010085339A (ja) | 接触燃焼式ガス検出素子を使用したガス検出装置のゼロ点調整方法 | |
CN112782263A (zh) | 一种六氟化硫气体组分质谱信号的校正方法 | |
CN114544536A (zh) | 一种氮氧化物浓度校准模型的生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 310053 Room 301, floor 3, building 3, No. 611, Dongguan Road, Binjiang District, Hangzhou, Zhejiang Province Applicant after: Zhongkong Quanshi Technology (Hangzhou) Co.,Ltd. Address before: 310053 Room B3155, 3rd floor, 368 Liuhe Road, Binjiang District, Hangzhou City, Zhejiang Province Applicant before: ZHEJIANG TRACETECH TECHNOLOGY CO.,LTD. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |