CN116500121A - 气体分析仪的分析方法、系统及可读存储介质 - Google Patents
气体分析仪的分析方法、系统及可读存储介质 Download PDFInfo
- Publication number
- CN116500121A CN116500121A CN202310772292.XA CN202310772292A CN116500121A CN 116500121 A CN116500121 A CN 116500121A CN 202310772292 A CN202310772292 A CN 202310772292A CN 116500121 A CN116500121 A CN 116500121A
- Authority
- CN
- China
- Prior art keywords
- peak
- fid
- response amplitude
- component
- fid response
- 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
- 238000004458 analytical method Methods 0.000 title claims abstract description 37
- 238000010586 diagram Methods 0.000 claims abstract description 88
- 238000010606 normalization Methods 0.000 claims abstract description 31
- 238000012545 processing Methods 0.000 claims abstract description 20
- 238000001514 detection method Methods 0.000 claims abstract description 16
- 230000003321 amplification Effects 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims abstract description 15
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 15
- 238000005070 sampling Methods 0.000 claims abstract description 15
- 239000000203 mixture Substances 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 19
- 230000001105 regulatory effect Effects 0.000 claims description 3
- 239000007789 gas Substances 0.000 description 41
- 239000003570 air Substances 0.000 description 10
- 238000005259 measurement Methods 0.000 description 9
- 239000012855 volatile organic compound Substances 0.000 description 5
- 239000012080 ambient air Substances 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- RSJKGSCJYJTIGS-UHFFFAOYSA-N undecane Chemical compound CCCCCCCCCCC RSJKGSCJYJTIGS-UHFFFAOYSA-N 0.000 description 4
- UHOVQNZJYSORNB-UHFFFAOYSA-N Benzene Chemical compound C1=CC=CC=C1 UHOVQNZJYSORNB-UHFFFAOYSA-N 0.000 description 3
- YMWUJEATGCHHMB-UHFFFAOYSA-N Dichloromethane Chemical compound ClCCl YMWUJEATGCHHMB-UHFFFAOYSA-N 0.000 description 3
- IMNFDUFMRHMDMM-UHFFFAOYSA-N N-Heptane Chemical compound CCCCCCC IMNFDUFMRHMDMM-UHFFFAOYSA-N 0.000 description 3
- YXFVVABEGXRONW-UHFFFAOYSA-N Toluene Chemical compound CC1=CC=CC=C1 YXFVVABEGXRONW-UHFFFAOYSA-N 0.000 description 3
- 238000007791 dehumidification Methods 0.000 description 3
- 238000012544 monitoring process Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- FYGHSUNMUKGBRK-UHFFFAOYSA-N 1,2,3-trimethylbenzene Chemical compound CC1=CC=CC(C)=C1C FYGHSUNMUKGBRK-UHFFFAOYSA-N 0.000 description 2
- KAKZBPTYRLMSJV-UHFFFAOYSA-N Butadiene Chemical compound C=CC=C KAKZBPTYRLMSJV-UHFFFAOYSA-N 0.000 description 2
- VGGSQFUCUMXWEO-UHFFFAOYSA-N Ethene Chemical compound C=C VGGSQFUCUMXWEO-UHFFFAOYSA-N 0.000 description 2
- YNQLUTRBYVCPMQ-UHFFFAOYSA-N Ethylbenzene Chemical compound CCC1=CC=CC=C1 YNQLUTRBYVCPMQ-UHFFFAOYSA-N 0.000 description 2
- 239000005977 Ethylene Substances 0.000 description 2
- RRHGJUQNOFWUDK-UHFFFAOYSA-N Isoprene Chemical compound CC(=C)C=C RRHGJUQNOFWUDK-UHFFFAOYSA-N 0.000 description 2
- CBENFWSGALASAD-UHFFFAOYSA-N Ozone Chemical compound [O-][O+]=O CBENFWSGALASAD-UHFFFAOYSA-N 0.000 description 2
- OFBQJSOFQDEBGM-UHFFFAOYSA-N Pentane Chemical compound CCCCC OFBQJSOFQDEBGM-UHFFFAOYSA-N 0.000 description 2
- ATUOYWHBWRKTHZ-UHFFFAOYSA-N Propane Chemical compound CCC ATUOYWHBWRKTHZ-UHFFFAOYSA-N 0.000 description 2
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Chemical compound C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 description 2
- RAHZWNYVWXNFOC-UHFFFAOYSA-N Sulphur dioxide Chemical compound O=S=O RAHZWNYVWXNFOC-UHFFFAOYSA-N 0.000 description 2
- MVPPADPHJFYWMZ-UHFFFAOYSA-N chlorobenzene Chemical compound ClC1=CC=CC=C1 MVPPADPHJFYWMZ-UHFFFAOYSA-N 0.000 description 2
- 238000002485 combustion reaction Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- RWGFKTVRMDUZSP-UHFFFAOYSA-N cumene Chemical compound CC(C)C1=CC=CC=C1 RWGFKTVRMDUZSP-UHFFFAOYSA-N 0.000 description 2
- SNRUBQQJIBEYMU-UHFFFAOYSA-N dodecane Chemical compound CCCCCCCCCCCC SNRUBQQJIBEYMU-UHFFFAOYSA-N 0.000 description 2
- 150000008282 halocarbons Chemical class 0.000 description 2
- NNPPMTNAJDCUHE-UHFFFAOYSA-N isobutane Chemical compound CC(C)C NNPPMTNAJDCUHE-UHFFFAOYSA-N 0.000 description 2
- VLKZOEOYAKHREP-UHFFFAOYSA-N n-Hexane Chemical compound CCCCCC VLKZOEOYAKHREP-UHFFFAOYSA-N 0.000 description 2
- 239000002243 precursor Substances 0.000 description 2
- ODLMAHJVESYWTB-UHFFFAOYSA-N propylbenzene Chemical compound CCCC1=CC=CC=C1 ODLMAHJVESYWTB-UHFFFAOYSA-N 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- SCYULBFZEHDVBN-UHFFFAOYSA-N 1,1-Dichloroethane Chemical compound CC(Cl)Cl SCYULBFZEHDVBN-UHFFFAOYSA-N 0.000 description 1
- PBKONEOXTCPAFI-UHFFFAOYSA-N 1,2,4-trichlorobenzene Chemical compound ClC1=CC=C(Cl)C(Cl)=C1 PBKONEOXTCPAFI-UHFFFAOYSA-N 0.000 description 1
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 1
- OTMSDBZUPAUEDD-UHFFFAOYSA-N Ethane Chemical compound CC OTMSDBZUPAUEDD-UHFFFAOYSA-N 0.000 description 1
- 229920000557 Nafion® Polymers 0.000 description 1
- XSTXAVWGXDQKEL-UHFFFAOYSA-N Trichloroethylene Chemical group ClC=C(Cl)Cl XSTXAVWGXDQKEL-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 150000001335 aliphatic alkanes Chemical class 0.000 description 1
- 150000001336 alkenes Chemical class 0.000 description 1
- HSFWRNGVRCDJHI-UHFFFAOYSA-N alpha-acetylene Natural products C#C HSFWRNGVRCDJHI-UHFFFAOYSA-N 0.000 description 1
- 150000004945 aromatic hydrocarbons Chemical class 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- IAQRGUVFOMOMEM-ARJAWSKDSA-N cis-but-2-ene Chemical compound C\C=C/C IAQRGUVFOMOMEM-ARJAWSKDSA-N 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 125000002534 ethynyl group Chemical group [H]C#C* 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 238000000769 gas chromatography-flame ionisation detection Methods 0.000 description 1
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 description 1
- -1 heterocyclic radical compound Chemical class 0.000 description 1
- 239000001282 iso-butane Substances 0.000 description 1
- 230000014759 maintenance of location Effects 0.000 description 1
- IJDNQMDRQITEOD-UHFFFAOYSA-N n-butane Chemical compound CCCC IJDNQMDRQITEOD-UHFFFAOYSA-N 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001294 propane Substances 0.000 description 1
- QQONPFPTGQHPMA-UHFFFAOYSA-N propylene Natural products CC=C QQONPFPTGQHPMA-UHFFFAOYSA-N 0.000 description 1
- 125000004805 propylene group Chemical group [H]C([H])([H])C([H])([*:1])C([H])([H])[*:2] 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- UBOXGVDOUJQMTN-UHFFFAOYSA-N trichloroethylene Natural products ClCC(Cl)Cl UBOXGVDOUJQMTN-UHFFFAOYSA-N 0.000 description 1
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/213—Feature extraction, e.g. by transforming the feature space; Summarisation; Mappings, e.g. subspace methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/22—Matching criteria, e.g. proximity measures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2123/00—Data types
- G06F2123/02—Data types in the time domain, e.g. time-series data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
- G06F2218/10—Feature extraction by analysing the shape of a waveform, e.g. extracting parameters relating to peaks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/12—Classification; Matching
- G06F2218/14—Classification; Matching by matching peak patterns
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
Abstract
本发明涉及一种气体分析仪的分析方法、系统及可读存储介质,其分析方法包括:分别将标样进样至数台气体分析仪,利用气体分析仪的FID检测器检测得到FID响应幅值时序图,并对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;标样包括高响应组分和若干目标组分的混合;将标样进样至目标气体分析仪,得到目标FID响应幅值时序图;根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。本发明有效提升检测精度,实现检测结果的可比性。
Description
技术领域
本发明属于检测分析技术领域,具体涉及一种气体分析仪的分析方法、系统及可读存储介质。
背景技术
挥发性有机物(VOCs)是PM2.5和臭氧污染的重要前体物,VOCs替代二氧化硫成为了大气污染控制的关键。其中,臭氧前体挥发性有机物(PAMS组分)为工业园区、厂界以及交通站等重点区域的监测重点。
目前,PAMS 57/117种组分监测有GC-MS和GC-FID-FID两种方案,因GC-FID-FID方案相较于系统简单、成本低,逐步成为主流方案,但在实际使用中,存在以下不足:
(1)由于不同VOCs物质中C-H键个数、杂环原子团化合物电离效率和所含官能团极性大小差异,FID检测器对不同组分响应的差异性;同时,由于FID检测器的两个电极间距及电子元件精度差异,同一物质在不同仪器间响应也存在差异,这直接导致测量结果的偏差性较大;如监测卤代烃或含氧衍生物,最大偏差达到170%左右;
(2)环境空气中湿度变化较大,对测量影响较大;湿度对不同组分影响量不同,如对高碳组分(十一烷、十二烷等)以及三氯苯等亲水性组分,影响相对误差达到8%以上;常规的Nafion管除水效率有限,同时可能会导致极性组分损失;深冷富集除水效果好,但过低温度导致结冰会造成保留时间漂移、分离效果变差以及浓度或高或低。
因此,本领域亟需解决FID检测器的差异性问题和平衡湿度影响,实现环境空气中挥发性有机物的质控和准确测量。
发明内容
基于现有技术中存在的上述缺点和不足,本发明的目的之一是至少解决现有技术中存在的上述问题之一或多个,换言之,本发明的目的之一是提供满足前述需求之一或多个的一种气体分析仪的分析方法、系统及可读存储介质。
为了达到上述发明目的,本发明采用以下技术方案:
一种气体分析仪的分析方法,包括以下步骤:
S1、分别将标样进样至数台气体分析仪,利用气体分析仪的FID检测器检测得到FID响应幅值时序图,并对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;其中,标样包括高响应组分和若干目标组分的混合;
S2、将标样进样至目标气体分析仪,得到目标FID响应幅值时序图;
S3、根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。
作为优选方案,所述步骤S1中,对所有FID响应幅值时序图进行归一化处理包括初始基线归一化处理和组分峰强度归一化处理。
作为优选方案,所述初始基线归一化处理包括:
提取FID响应幅值时序图自起点至第一个组分峰的起始点之间的所有采样点,构成初始基线的数列;
计算初始基线的数列对应的FID响应幅值的中位数;
计算所有FID响应幅值时序图的初始基线的数列对应的FID响应幅值的中位数的平均值;
根据FID响应幅值的中位数的平均值对FID响应幅值时序图的初始基线的数列对应的FID响应幅值进行线性归一化,得到标准化初始基线。
作为优选方案,所述组分峰强度归一化处理包括:
分别提取各FID响应幅值时序图对应的多个组分峰,得到各组分峰的峰强度;
计算各FID响应幅值时序图对应同一组分的组分峰的峰强度的中位数和正态偏差;
根据组分峰的峰强度的中位数和正态偏差得到标准化峰强度,并根据标准化峰强度计算得到标准化峰面积;
根据标准化初始基线、标准化峰强度和标准化峰面积构建标准化FID响应幅值时序图。
作为优选方案,所述步骤S3,包括:
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线的数列对应的FID响应幅值的差值,调节偏置电压以使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线匹配;
偏置电压ΔV:
;
其中,ε为常数,SNR为放大电路的信噪比,Q为初始基线的采样点的数量,Δy q 为初始基线的第q个采样点对应的FID响应幅值的差值。
作为优选方案,所述步骤S3,还包括:
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的高响应组分的峰面积的差值,调节前级放大电阻,以使峰面积的差值在预设偏差范围内;
峰面积由FID检测器的放大电路输出的电压信号转换得到,FID检测器的放大电路输出的电压信号V out :
;
其中,V s 为FID检测器的输入信号,F和V b 为放大电路对应的常数,μ为前级放大电阻的调节系数。
作为优选方案,所述步骤S3,还包括:
调节目标组分对应的组分拟合修正系数,以使目标FID响应幅值时序图与标准化FID响应幅值时序图对应的目标组分的峰面积的偏差在目标偏差范围内;
其中,组分拟合修正系数包括峰底、峰高、峰宽和拖尾因子。
作为优选方案,气体分析仪的分析方法,还包括以下步骤:
S4、获取不同湿度下各组分在不同标准浓度的峰基线和峰面积,以计算不同湿度下各组分的湿度影响峰基线系数和湿度影响峰面积系数;其中,峰基线为峰型起始位置所处的时间段和终止位置所处的时间段及其对应的FID响应幅值的线性拟合曲线;其中,湿度检测位于FID检测器的出气口;
S5、判断湿度影响峰基线系数是否大于第一预设阈值;若是,则计算对应湿度下各组分在不同标准浓度的峰基线值的平均值,并将其与基准湿度下各组分在不同标准浓度的峰基线值的差值作为湿度对峰基线的影响量;其中,峰基线值为峰型起始位置所处的时间段和终止位置所处的时间段对应的FID响应幅值的平均值;若否,则湿度对峰基线的影响量为零;
判断湿度影响峰面积系数是否大于第二预设阈值;若是,则根据湿度对峰基线的影响量计算对应湿度下各组分在不同标准浓度的峰面积,并根据组分在不同湿度、不同标准浓度下的峰面积与标准浓度建立加湿工作曲线;
S6、采集待测气样对应的FID响应幅值时序图和当前湿度,根据当前湿度匹配加湿工作曲线,得到待测气样中被测组分的峰面积,并根据被测组分的峰面积计算得到被测组分的浓度。
本发明还提供一种气体分析仪的分析系统,应用如上任一项方案所述的分析方法,所述分析系统包括:
归一化处理模块,用于对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;
调节模块,用于根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。
本发明还提供一种可读存储介质,存储有指令,当指令在计算机上运行时,使得计算机执行如上任一项方案所述的分析方法。
本发明与现有技术相比,有益效果是:
(1)本发明构建标准化FID响应幅值时序图,解决了气体分析仪之间FID检测器的差异性问题,有效提升检测精度,还实现检测结果的可比性;
(2)由于湿度通过影响FID检测器燃烧的氢空比,进而抬高初始基线、峰基线和影响被测组分的峰高或峰面积,本发明建立湿度工作曲线,有效解决湿度的影响,进一步提升检测精度;本发明的湿度工作曲线是在构建标准化FID响应幅值时序图的基础上建立的,湿度工作曲线具有普适性;
(3)本发明的湿度以检测FID检测器出气口的湿度检测为准,主要是考虑到富集处的除湿效率随着时间变化变差,避免湿度工作曲线补偿引入偏差。
附图说明
图1是本发明实施例1的气体分析仪的分析方法的流程图;
图2是本发明实施例1的三个FID检测器对同一标样的FID响应幅值时序图;
图3是本发明实施例1的三个FID检测器在不同时间对同一组分的浓度检测分布图;
图4是本发明实施例1的三个FID检测器在不同时间对同一组分的浓度检测经过标准化之后的浓度分布图;
图5是本发明实施例1的气体分析仪的分析系统的构架图。
具体实施方式
为了更清楚地说明本发明实施例,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
实施例1:
如图1所示,本实施例的气体分析仪的分析方法,包括以下步骤:
一、构建标准化FID响应幅值时序图
1、分别将标样进样至N台气体分析仪,利用气体分析仪的FID检测器检测得到FID响应幅值时序图(即FID检测器检测输出的谱图),N为大于1的整数,气体分析仪的数量具体可根据实际检测场景进行确定。其中,标样包括高响应组分和若干目标组分的混合,高响应组分为FID响应幅值最高的组分或者高于预设幅值的组分,一般选择乙烯;目标组分根据检测应用场景确定,一般根据类别分为四类:烷烃类(乙烷、丙烷、正丁烷、异丁烷、正戊烷、正己烷、正庚烷等)、烯烃类(乙烯、乙炔、丙烯、顺-2-丁烯、1,3丁二烯,异戊二烯等)、芳香烃(苯、甲苯、乙苯、苯乙烯、异丙苯、正丙苯、三甲苯等)、卤代烃(二氯甲烷、二氯乙烷、三氯乙烯、氯苯等);
2、对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;
具体地,对所有FID响应幅值时序图进行归一化处理包括初始基线归一化处理和组分峰强度归一化处理,具体过程如下:
(1)初始基线归一化处理
提取FID响应幅值时序图自起点至第一个组分峰的起始点之间的所有采样点,构成初始基线的数列;
分别计算各FID响应幅值时序图的初始基线的数列对应的FID响应幅值的中位数;
计算所有FID响应幅值时序图的初始基线的数列对应的FID响应幅值的中位数的平均值:
;
其中,A i 为第i台气体分析仪检测得到的FID响应幅值时序图的初始基线的数列对应的FID响应幅值的中位数,i∈[1,N];
根据FID响应幅值的中位数的平均值对FID响应幅值时序图的初始基线的数列对应的FID响应幅值进行线性归一化:
;
其中,为归一化后的FID响应幅值,y为初始基线的数列对应的FID响应幅值,min(y)为初始基线的数列对应的FID响应幅值的最小值,max(y)为初始基线的数列对应的FID响应幅值的最大值;
经过上述线性归一化之后,得到标准化初始基线;
(2)组分峰强度归一化处理
分别提取各FID响应幅值时序图对应的多个组分峰,得到各组分峰的峰强度,组成峰强度数列:
;
其中,m为组分峰的个数,y j , i 为第i台气体分析仪检测得到的FID响应幅值时序图的第j个组分峰的峰强度,j∈[1,m];
计算各FID响应幅值时序图对应同一组分j的组分峰的峰强度的中位数和正态偏差/>;
根据组分j的组分峰的峰强度的中位数和正态偏差得到标准化峰强度的范围为:
;
并根据标准化峰强度计算得到标准化峰面积;
(3)根据标准化初始基线、标准化峰强度和标准化峰面积构建标准化FID响应幅值时序图。
上述根据标准化峰强度计算标准化峰面积以及根据标准化初始基线、标准化峰强度和标准化峰面积构建标准化FID响应幅值时序图的具体过程可参考现有技术,在此不赘述。本实施例通过构建标准化FID响应幅值时序图,解决了气体分析仪设备之间指标可比性问题和品控问题。
二、采集目标FID响应幅值时序图
将标样进样至目标气体分析仪,得到目标FID响应幅值时序图。
三、目标FID响应幅值时序图的标准化
根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。具体的调节过程如下:
1、偏置电压的调节
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线的数列对应的FID响应幅值的差值,调节偏置电压以使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线匹配(即保持一致)。本实施例的偏置电压ΔV:
;
其中,ε为常数,根据放大电路频率进行选择,取值范围为0.7-1.3;SNR为放大电路的信噪比,T为初始基线的采样点的数量,Δy q 为目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线的第q个采样点对应的FID响应幅值的差值,q∈[1,Q]。
通过偏置电压的调节,每个无吸收峰段的初始基线均保持在相同等级的幅值范围内。偏置电压的调节具体可参考现有技术,在此不赘述。
2、前级放大电阻的调节
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的高响应组分的峰面积的差值,调节前级放大电阻,以使峰面积的差值在预设偏差范围内,如峰面积的差值在1%以内;
峰面积由FID检测器的放大电路输出的电压信号转换得到(具体转换方法可参考现有技术,在此不赘述),放大电路一般为对数放大电路,FID检测器的放大电路输出的电压信号V out :;
其中,V s 为FID检测器的输入信号;放大电路确定后,F和V b 为常数,μ为前级放大电阻的调节系数,可根据峰面积的差值大小进行调节。前级放大电阻的调节具体可参考现有技术,在此不赘述。
3、组分拟合修正系数的调节
调节目标组分对应的组分拟合修正系数,以使目标FID响应幅值时序图与标准化FID响应幅值时序图对应的目标组分的峰面积的偏差在目标偏差范围内,如峰面积的差值在2%以内;其中,组分拟合修正系数包括峰底、峰高、峰宽和拖尾因子。
由于FID检测器输出的每个目标组分的峰型在标准化时已确定,可进行偏最小二乘法得到每个目标组分的峰型函数f (D,H,W,T),峰型函数由峰底D、峰高H、峰宽W和拖尾因子T表征;
当针对某台气体分析仪进行调节时,采用偏最小二乘法,调节峰底D、峰高H、峰宽W和拖尾因子T四个系数,确保目标组分的峰面积的最大偏差在2%以内。
如图2所示,1号FID检测器(简称1号FID)、2号FID检测器(简称2号FID)、3号FID检测器(简称3号FID)对同一标样的响应存在明显的差异。
如图3所示,1号FID、2号FID和3号FID在不同时间对同一组分浓度的检测,最大绝对偏差达到0.967nmol/mol;经过本实施例的标准化之后,如图4所示,1号FID、2号FID和3号FID在不同时间对同一组分浓度的检测,最大绝对偏差缩小至0.364nmol/mol。
四、建立加湿工作曲线
另外,由于检测环境湿度的改变会影响FID检测器燃烧的氢空比,进而抬高初始基线、峰基线以及影响被测组分的峰高或峰面积。基于此,本实施例通过自动加湿装置加入不同体积浓度的高湿氮气,获得FID响应幅值时序图,并通过设置在FID检测器出气口的湿度传感器监测除湿后的水汽浓度,与被测组分的色谱时序位置对应,建立加湿工作曲线,具体过程如下:
1、通过自动加湿装置进行加湿,获取不同湿度(以5%湿度为步长,以FID检测器出气口的检测结果为准)下的各组分在不同标准浓度的峰基线和峰面积;
其中,峰基线为峰型起始位置所处的时间段和终止位置所处的时间段及其对应的FID响应幅值的线性拟合曲线f baseline ,相应地,峰基线值为峰型起始位置所处的时间段和终止位置所处的时间段对应的FID响应幅值的平均值。
2、计算不同湿度下各组分的湿度影响峰基线系数和湿度影响峰面积系数;
本实施例在FID检测器出气口检测的湿度R k 对第l个标准浓度的组分的湿度影响峰基线系数:
;
其中,R 0为基准湿度(即常温常压环境下在FID检测器出气口检测的湿度),B l (R 0)为基准湿度R 0下的第l个标准浓度的组分的峰基线值,B l (R k )为湿度R k 下的第l个标准浓度的组分的峰基线值;l∈[1,L],L为组分对应选取的标准浓度的数量;k∈[1,K],K为组分对应选取的湿度的数量;
本实施例的湿度R k 对第l个标准浓度的组分的湿度影响峰面积系数:
;
其中,S l (R 0)为基准湿度R 0下的第l个标准浓度的组分的峰面积,S l (R k )为湿度R k 下的第l个标准浓度的组分的峰面积。
3、判断湿度影响峰基线系数是否大于基线扣除阈值;若是,则计算对应湿度下各组分在不同标准浓度的峰基线值的平均值,并将其与基准湿度下各组分在不同标准浓度的峰基线值的差值作为湿度对峰基线的影响量;若否,则湿度对峰基线的影响量为零。
由于湿度对所影响的组分的峰基线值与组分浓度关联小,故取组分在不同标准浓度的峰基线值的平均值;
相应地,湿度对峰基线的影响量,进而建立湿度和组分关联的基线扣除数据库。
另外,还判断湿度影响峰面积系数是否大于峰面积扣除阈值;若是,则根据湿度对峰基线的影响量计算对应湿度下各组分在不同标准浓度的峰面积,并根据组分在不同湿度、不同标准浓度下的峰面积与标准浓度建立加湿工作曲线。在实际测量中,根据湿度测量值和被测组分,搜寻对应的加湿工作曲线,根据峰面积即可反算得到被测组分的浓度值。
具体地,提取以5%湿度为步长,同一组分的不同标准浓度下的峰面积S l (R k ),湿度R k 下的第l个标准浓度的组分的峰面积计算公式为:
;
其中,和/>为第l个标准浓度的组分的组分峰的峰型起始位置和终止位置的时间,/>为第l个标准浓度的组分的组分峰的峰型函数,/>为第l个标准浓度的组分的组分峰的峰基线。
将不同标准浓度的组分对应的峰面积与标准浓度拟合建立工作曲线,拟合一般为线性拟合,最终得到湿度和组分关联的加湿工作曲线库,为后期实测提供湿度补偿。
由于不同气体分析仪设备的FID检测器已建立标准化,故加湿工作曲线具有普适性。
本实施例的湿度检测以FID检测器的出气口的湿度测量为准,主要是考虑到气体分析仪的富集处的除湿效率随着时间变化变差,导致湿度工作曲线补偿引入偏差的问题。
五、实时在线检测
实时采集环境空气中的待测气样,进样至目标气体分析仪,得到待测气样对应的FID响应幅值时序图和当前湿度R(即FID检测器出气口检测的水汽湿度),根据当前湿度匹配加湿工作曲线,得到待测气样中被测组分g的峰面积:
;
其中,t s 和t e 分别为被测组分g的组分峰的峰型起始位置和终止位置的时间,为被测组分g的组分峰的峰型函数,/>为被测组分g的组分峰的峰基线;/>为当前湿度R对被测组分g的组分峰的峰基线的影响量。
根据被测组分g的峰面积计算得到被测组分g的浓度:
;其中,a和b为上述线性拟合的系数,δ为修正因子。
。本实施例建立的加湿工作曲线能够有效提升浓度测量的精度,以待测气样十一烷为例,如表1所示,不同湿度影响相对误差达到8%以上,但通过加湿工作曲线处理,相对误差能控制在1%以内。
表1有无加湿工作曲线的测量效果表
基于本实施例上述的气体分析仪的分析方法,如图5所示,本实施例还提供气体分析仪的分析系统,包括以下功能模块:归一化处理模块、调节模块、加湿工作曲线构建模块和匹配模块。
本实施例的归一化处理模块用于对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图。
将标样进样至目标气体分析仪,得到目标FID响应幅值时序图。
本实施例的调节模块用于根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。
本实施例的加湿工作曲线构建模块用于建立加湿工作曲线。
实时采集环境空气中的待测气样,进样至目标气体分析仪,得到待测气样对应的FID响应幅值时序图和当前湿度;本实施例的匹配模块用于根据当前湿度匹配加湿工作曲线,以得到待测气样中被测组分的峰面积,并根据被测组分的峰面积计算得到被测组分的浓度。
上述功能模块的具有处理过程可参考上述分析方法中的详细描述,在此不赘述。
本实施例还提供可读存储介质,存储有指令,当指令在计算机上运行时,使得计算机执行本实施例的分析方法,实现智能化分析。
本实施例的气体分析仪的分析方法、系统及可读存储介质可应用于双GC-FID系统或GC-FID-FID系统中,应用范围广。
实施例2:
本实施例的气体分析仪的分析方法与实施例1的不同之处在于:
省略加湿工作曲线的建立过程及其相关处理过程,适用于恒定湿度的场景中使用,简化分析流程;
相应地,本实施例的气体分析仪的分析系统的构架作相应的简化,满足不同应用场景的需求;
本实施例的可读存储介质,存储有指令,当指令在计算机上运行时,使得计算机执行本实施例的分析方法,实现智能化分析;
其他分析流程以及功能模块可参考实施例1。
以上所述仅是对本发明的优选实施例及原理进行了详细说明,对本领域的普通技术人员而言,依据本发明提供的思想,在具体实施方式上会有改变之处,而这些改变也应视为本发明的保护范围。
Claims (10)
1.一种气体分析仪的分析方法,其特征在于,包括以下步骤:
S1、分别将标样进样至数台气体分析仪,利用气体分析仪的FID检测器检测得到FID响应幅值时序图,并对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;其中,标样包括高响应组分和若干目标组分的混合;
S2、将标样进样至目标气体分析仪,得到目标FID响应幅值时序图;
S3、根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。
2.根据权利要求1所述的气体分析仪的分析方法,其特征在于,所述步骤S1中,对所有FID响应幅值时序图进行归一化处理包括初始基线归一化处理和组分峰强度归一化处理。
3.根据权利要求2所述的气体分析仪的分析方法,其特征在于,所述初始基线归一化处理包括:
提取FID响应幅值时序图自起点至第一个组分峰的起始点之间的所有采样点,构成初始基线的数列;
计算初始基线的数列对应的FID响应幅值的中位数;
计算所有FID响应幅值时序图的初始基线的数列对应的FID响应幅值的中位数的平均值;
根据FID响应幅值的中位数的平均值对FID响应幅值时序图的初始基线的数列对应的FID响应幅值进行线性归一化,得到标准化初始基线。
4.根据权利要求3所述的气体分析仪的分析方法,其特征在于,所述组分峰强度归一化处理包括:
分别提取各FID响应幅值时序图对应的多个组分峰,得到各组分峰的峰强度;
计算各FID响应幅值时序图对应同一组分的组分峰的峰强度的中位数和正态偏差;
根据组分峰的峰强度的中位数和正态偏差得到标准化峰强度,并根据标准化峰强度计算得到标准化峰面积;
根据标准化初始基线、标准化峰强度和标准化峰面积构建标准化FID响应幅值时序图。
5.根据权利要求4所述的气体分析仪的分析方法,其特征在于,所述步骤S3,包括:
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线的数列对应的FID响应幅值的差值,调节偏置电压以使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线匹配;
偏置电压ΔV:
;
其中,ε为常数,SNR为放大电路的信噪比,Q为初始基线的采样点的数量,Δy q 为初始基线的第q个采样点对应的FID响应幅值的差值。
6.根据权利要求5所述的气体分析仪的分析方法,其特征在于,所述步骤S3,还包括:
根据目标FID响应幅值时序图与标准化FID响应幅值时序图的高响应组分的峰面积的差值,调节前级放大电阻,以使峰面积的差值在预设偏差范围内;
峰面积由FID检测器的放大电路输出的电压信号转换得到,FID检测器的放大电路输出的电压信号V out :
;
其中,V s 为FID检测器的输入信号,F和V b 为放大电路对应的常数,μ为前级放大电阻的调节系数。
7.根据权利要求6所述的气体分析仪的分析方法,其特征在于,所述步骤S3,还包括:
调节目标组分对应的组分拟合修正系数,以使目标FID响应幅值时序图与标准化FID响应幅值时序图对应的目标组分的峰面积的偏差在目标偏差范围内;
其中,组分拟合修正系数包括峰底、峰高、峰宽和拖尾因子。
8.根据权利要求1-7任一项所述的气体分析仪的分析方法,其特征在于,还包括以下步骤:
S4、获取不同湿度下各组分在不同标准浓度的峰基线和峰面积,以计算不同湿度下各组分的湿度影响峰基线系数和湿度影响峰面积系数;其中,峰基线为峰型起始位置所处的时间段和终止位置所处的时间段及其对应的FID响应幅值的线性拟合曲线;其中,湿度检测位于FID检测器的出气口;
S5、判断湿度影响峰基线系数是否大于第一预设阈值;若是,则计算对应湿度下各组分在不同标准浓度的峰基线值的平均值,并将其与基准湿度下各组分在不同标准浓度的峰基线值的差值作为湿度对峰基线的影响量;其中,峰基线值为峰型起始位置所处的时间段和终止位置所处的时间段对应的FID响应幅值的平均值;若否,则湿度对峰基线的影响量为零;
判断湿度影响峰面积系数是否大于第二预设阈值;若是,则根据湿度对峰基线的影响量计算对应湿度下各组分在不同标准浓度的峰面积,并根据组分在不同湿度、不同标准浓度下的峰面积与标准浓度建立加湿工作曲线;
S6、采集待测气样对应的FID响应幅值时序图和当前湿度,根据当前湿度匹配加湿工作曲线,得到待测气样中被测组分的峰面积,并根据被测组分的峰面积计算得到被测组分的浓度。
9.一种气体分析仪的分析系统,应用如权利要求1-7任一项所述的分析方法,其特征在于,所述分析系统包括:
归一化处理模块,用于对所有FID响应幅值时序图进行归一化处理,得到标准化FID响应幅值时序图;
调节模块,用于根据初始基线以及高响应组分和目标组分的峰面积作为调节参考量,通过调节FID检测器的放大电路的偏置电压、前级放大电阻、组分拟合修正系数,使目标FID响应幅值时序图与标准化FID响应幅值时序图的初始基线和峰面积匹配。
10.一种可读存储介质,存储有指令,其特征在于,当指令在计算机上运行时,使得计算机执行如权利要求1-8任一项所述的分析方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310772292.XA CN116500121B (zh) | 2023-06-28 | 2023-06-28 | 气体分析仪的分析方法、系统及可读存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310772292.XA CN116500121B (zh) | 2023-06-28 | 2023-06-28 | 气体分析仪的分析方法、系统及可读存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116500121A true CN116500121A (zh) | 2023-07-28 |
CN116500121B CN116500121B (zh) | 2023-10-24 |
Family
ID=87317029
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310772292.XA Active CN116500121B (zh) | 2023-06-28 | 2023-06-28 | 气体分析仪的分析方法、系统及可读存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116500121B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101010567A (zh) * | 2004-07-27 | 2007-08-01 | 乌拉帝米尔·阿列克桑德罗维奇·祖布科夫 | 用于产生独立多元校准模型的方法 |
US20080201095A1 (en) * | 2007-02-12 | 2008-08-21 | Yip Ping F | Method for Calibrating an Analytical Instrument |
WO2015169686A2 (en) * | 2014-05-06 | 2015-11-12 | Københavns Universitet | A computer assisted method for quantification of total hydrocarbon concentrations and pollution type apportionment in soil samples by use of gc-fid chromatograms |
CN107817302A (zh) * | 2017-09-27 | 2018-03-20 | 中国石油天然气集团公司 | 一种石油烃组分的气相色谱分析方法 |
US20180321201A1 (en) * | 2017-05-04 | 2018-11-08 | Shimadzu Corporation | Chromatography/mass spectrometry data processing device |
CN108982728A (zh) * | 2018-05-30 | 2018-12-11 | 上海市环境监测中心 | 一种环境空气非甲烷总烃在线监测仪器优化与校准方法 |
CN111735870A (zh) * | 2020-07-31 | 2020-10-02 | 暨南大学 | 一种在线实时分析质谱的校正方法及校正装置 |
CN114397395A (zh) * | 2021-12-31 | 2022-04-26 | 杭州春来科技有限公司 | 基于非甲烷总烃测定用fid检测器的氧干扰修正方法及系统 |
JP2022178382A (ja) * | 2021-05-20 | 2022-12-02 | 株式会社島津製作所 | 臭素系難燃化合物の定量装置 |
CN115639168A (zh) * | 2022-12-21 | 2023-01-24 | 杭州泽天春来科技有限公司 | 气体分析仪的气体检测方法、系统及介质 |
CN116337783A (zh) * | 2023-04-11 | 2023-06-27 | 安荣信科技(南京)有限公司 | 一种气体分析仪的多点标定方法和系统 |
-
2023
- 2023-06-28 CN CN202310772292.XA patent/CN116500121B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101010567A (zh) * | 2004-07-27 | 2007-08-01 | 乌拉帝米尔·阿列克桑德罗维奇·祖布科夫 | 用于产生独立多元校准模型的方法 |
US20080201095A1 (en) * | 2007-02-12 | 2008-08-21 | Yip Ping F | Method for Calibrating an Analytical Instrument |
WO2015169686A2 (en) * | 2014-05-06 | 2015-11-12 | Københavns Universitet | A computer assisted method for quantification of total hydrocarbon concentrations and pollution type apportionment in soil samples by use of gc-fid chromatograms |
US20180321201A1 (en) * | 2017-05-04 | 2018-11-08 | Shimadzu Corporation | Chromatography/mass spectrometry data processing device |
CN107817302A (zh) * | 2017-09-27 | 2018-03-20 | 中国石油天然气集团公司 | 一种石油烃组分的气相色谱分析方法 |
CN108982728A (zh) * | 2018-05-30 | 2018-12-11 | 上海市环境监测中心 | 一种环境空气非甲烷总烃在线监测仪器优化与校准方法 |
CN111735870A (zh) * | 2020-07-31 | 2020-10-02 | 暨南大学 | 一种在线实时分析质谱的校正方法及校正装置 |
JP2022178382A (ja) * | 2021-05-20 | 2022-12-02 | 株式会社島津製作所 | 臭素系難燃化合物の定量装置 |
CN114397395A (zh) * | 2021-12-31 | 2022-04-26 | 杭州春来科技有限公司 | 基于非甲烷总烃测定用fid检测器的氧干扰修正方法及系统 |
CN115639168A (zh) * | 2022-12-21 | 2023-01-24 | 杭州泽天春来科技有限公司 | 气体分析仪的气体检测方法、系统及介质 |
CN116337783A (zh) * | 2023-04-11 | 2023-06-27 | 安荣信科技(南京)有限公司 | 一种气体分析仪的多点标定方法和系统 |
Non-Patent Citations (1)
Title |
---|
尤坤等: "基于FID的机动车尾气THC测量模块设计与对比实验", 大气与环境光学学报, vol. 12, no. 01, pages 43 - 49 * |
Also Published As
Publication number | Publication date |
---|---|
CN116500121B (zh) | 2023-10-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US5032721A (en) | Acid gas monitor based on ion mobility spectrometry | |
US10393763B2 (en) | Odor discriminating apparatus | |
CA2134957C (en) | Transformer oil gas extractor | |
Chen et al. | Applications and technology of electronic nose for clinical diagnosis | |
DE69204000D1 (de) | Direkte Abgasprobennahme für momentane Messungen. | |
CN115639168B (zh) | 气体分析仪的气体检测方法、系统及介质 | |
CN108593586B (zh) | 一种针对高湿污染气体的傅里叶红外光谱检测装置 | |
CN111735870A (zh) | 一种在线实时分析质谱的校正方法及校正装置 | |
CA2469099A1 (en) | Spectroscopic method and apparatus for total hemoglobin measurement | |
Kotiaho et al. | Membrane inlet ion mobility spectrometry for on-line measurement of ethanol in beer and in yeast fermentation | |
CN105911106A (zh) | Gis设备的绝缘气体状态检测系统 | |
CN113267462A (zh) | 一种应用紫外光谱对水体有机污染特征进行表征的方法 | |
Gorlo et al. | Application of solid-phase microextraction to monitoring indoor air quality | |
CN116500121B (zh) | 气体分析仪的分析方法、系统及可读存储介质 | |
WO2012056709A1 (ja) | 分析装置 | |
CN112666241B (zh) | 一种区域tvoc监测方法 | |
NO20170996A1 (en) | Chemical analysis method for measurement of tetrafluoromethane, cf4, with improved selectivity | |
CN112683836B (zh) | 基于bp神经网络的二氧化碳传感器的标定方法及其系统 | |
CN106525998A (zh) | 一种测定固定污染源废气中苯系物成分含量的方法 | |
CN101793820B (zh) | 双模联动水质在线监测方法 | |
Li et al. | Performance of an infrared photoacoustic single gas analyzer in measuring ammonia from poultry houses | |
Pungor et al. | Direct monitoring of the liquid and gas phases during a fermentation in a computer-mass-spectrometer-fermentor system | |
CN115856185B (zh) | 分析仪的处理方法、系统及可读存储介质 | |
CN113917042B (zh) | 一种针对加热卷烟用再造烟叶水分快速测定方法 | |
CN112730574B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |