CN114609258A - 基于液相色谱-高分辨质谱的石油分子表征数据处理方法 - Google Patents

基于液相色谱-高分辨质谱的石油分子表征数据处理方法 Download PDF

Info

Publication number
CN114609258A
CN114609258A CN202011447562.2A CN202011447562A CN114609258A CN 114609258 A CN114609258 A CN 114609258A CN 202011447562 A CN202011447562 A CN 202011447562A CN 114609258 A CN114609258 A CN 114609258A
Authority
CN
China
Prior art keywords
mass
ion
ions
threshold value
petroleum
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
Application number
CN202011447562.2A
Other languages
English (en)
Other versions
CN114609258B (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.)
Dalian University of Technology
Dalian Institute of Chemical Physics of CAS
Petrochina Petrochemical Research Institute
Original Assignee
Dalian University of Technology
Dalian Institute of Chemical Physics of CAS
Petrochina Petrochemical Research Institute
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 Dalian University of Technology, Dalian Institute of Chemical Physics of CAS, Petrochina Petrochemical Research Institute filed Critical Dalian University of Technology
Priority to CN202011447562.2A priority Critical patent/CN114609258B/zh
Publication of CN114609258A publication Critical patent/CN114609258A/zh
Application granted granted Critical
Publication of CN114609258B publication Critical patent/CN114609258B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/62Detectors specially adapted therefor
    • G01N30/72Mass spectrometers
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/88Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/20Identification of molecular entities, parts thereof or of chemical compositions
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/40Searching chemical structures or physicochemical data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02Column chromatography
    • G01N30/88Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86
    • G01N2030/8809Integrated analysis systems specially adapted therefor, not covered by a single one of the groups G01N30/04 - G01N30/86 analysis specially adapted for the sample

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • General Physics & Mathematics (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Analytical Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

本发明公开了基于液相色谱‑高分辨质谱的石油分子表征数据处理方法,基于液相色谱‑高分辨质谱采集到的原始数据文件,通过构建理论分子库、侯选离子筛选、峰提取与检测、去假阳性等步骤实现石油馏分分子表征。本发明方法充分挖掘了石油样本的液相色谱‑高分辨质谱信息,快速、可靠地获得不同石油馏分中各类型化合物的分子式、保留时间及其强度,可用于基于液相色谱‑高分辨质谱的石油馏分数据处理,分子表征信息提取更充分、更可靠。

Description

基于液相色谱-高分辨质谱的石油分子表征数据处理方法
技术领域
本发明涉及一种基于液相色谱-高分辨质谱的石油分子表征数据处理方法,属于分析化学及石油组学分析技术领域。
背景技术
石油资源的重质化,给石油的炼制带来了挑战。石油组学是从分子水平上认识石油的组成及转化规律,从而实现高效精准石油加工,促进炼油技术的进步。高分辨质谱技术凭借其超高的分辨率,可获得化合物的精确质荷比,根据精确质荷比,人们可获得化合物的分子式信息。但是重质化的石油组成十分复杂,在高分辨质谱分析过程中,存在离子抑制效应,获得的化合物信息不全面。因此在重质油进入高分辨质谱前,可对其进行分离,液相色谱是一种高效的分离手段,目前被普遍用于石油的前处理过程中。在线液相色谱-高分辨质谱技术可用于石油的分子表征,但是石油的在线液相色谱数据进一步加大了石油数据的处理难度,目前很少有文献报道用于石油的在线液相色谱-高分辨质谱数据的处理方法。
液相色谱在一定程度上能降低油样的复杂程度,从而降低离子抑制效应扩大高分辨质谱数据的覆盖度。但是重油的组成复杂,化合物在液相色谱中并不能实现基线分离。重油的液相色谱-高分辨质谱数据中,化合物的扫描点数多,色谱峰宽。用商品化软件处理重油的液相色谱-高分辨质谱数据,由于色谱峰的峰很宽,在峰识别的过程中,会造成严重的化合物色谱峰丢失的情况。在不考虑化合物的峰形的情况下对数据进行处理,会出现重复积分,保留时间不准确等问题,因此得到的数据并不能真实体现以及直观表达实际样本中化合物的分布情况。因此有必要建立用于石油的液相色谱-高分辨质谱数据的处理方法。
发明内容
本发明公开了一种基于液相色谱-高分辨质谱的石油馏分分子表征数据处理新方法,基于液相色谱-高分辨质谱采集到的原始数据文件,通过构建理论分子库、侯选离子筛选、峰提取与检测、去假阳性等步骤实现石油馏分分子表征。本发明方法充分挖掘了石油样本的液相色谱-高分辨质谱信息,快速、可靠地获得不同石油馏分中各类型化合物的分子式、保留时间及其强度,可用于基于液相色谱-高分辨质谱的石油馏分数据处理,分子表征信息提取更充分、更可靠。
本发明为实现上述目的所采用的技术方案是:基于液相色谱-高分辨质谱的石油分子表征数据处理方法,包括以下步骤:
1)通过液相色谱-高分辨质谱仪获取石油样本的质谱数据,构建用于数据处理的空白对照样本和待测石油馏分样本数据;
2)将空白对照样本和待测石油馏分样本数据通过质谱分析软件获取空白对照样本和石油馏分样本数据在整个保留时间范围的平均质谱图至excel文件;
3)获取平均质谱图噪音阈值,去除噪音离子;
4)根据平均质谱图间的质量偏差阈值和质谱响应偏差阈值,对步骤3)去噪后的待测石油馏分平均谱图进行扣除空白;
5)建立待测石油馏分的分子类型、等效双键数范围、碳数范围、质荷比范围,并依据石油馏分碳数和等效双键数的极限关系约束最低碳数;生成等效双键数和碳数均连续分布的石油馏分理论分子库,分子库中每个理论分子均包含注释信息,所述注释信息包括分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型及其理论单同位素离子质荷比;
6)获取平均质谱图中实测离子和理论离子的质荷比偏差阈值,以及实测同位素分布和理论计算同位素分布的匹配阈值;查询石油馏分组成理论分子式库,从步骤4)去噪、去空白后的待测石油馏分平均谱图中提取出与石油馏分理论分子库匹配的,且满足小于质荷比偏差阈值的单同位素离子;根据该离子的分子式得到理论同位素的质荷比以及强度,并与样本的平均谱图进行匹配,匹配出满足小于质荷比偏差阈值以及强度阈值的同位素离子;将同时满足质量偏差和同位素分布阈值的离子质荷比及其对应的分子注释信息输出至文件,作为待测石油馏分的候选单同位素离子注释结果;
7)读取液相色谱-高分辨质谱仪所采集数据的原始数据文件,转换并保存为txt文本格式文件;
8)基于步骤6)的候选离子注释结果,获取原始数据中的实测离子与平均谱图中已注释候选离子的质荷比偏差阈值,以及原始数据中单一离子的连续扫描次数阈值,从原始数据中筛选出满足小于质荷比偏差阈值和连续扫描次数阈值的所有离子;将单一离子连续扫描点的累积强度值作为该离子强度;输出各离子的质荷比,起始扫描点和结束扫描点对应的保留时间,强度以及分子注释信息至文件,作为峰检测结果文件;
9)对峰检测结果文件中的离子,提取每个离子所有扫描点的时间和强度,并对其做高斯分布函数,模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子的保留时间;若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,得到其它流出峰与强度最大的流出峰的分离度;根据分离度阈值,去除其中分离度不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件;
10)根据待测石油馏分的等效双键数和碳数连续分布阈值,在初始假阳性去除结果文件中提取出满足小于等效双键数和碳数连续分布阈值的全部离子,输出至文件,作为最终去假阳性结果文件。
步骤2)中:平均谱图中包括质荷比和强度信息。
步骤4)中:同时满足平均质谱图间的质量偏差和质谱响应偏差阈值的石油馏分平均谱图信号将被去除。
步骤5)中:不同质谱离子化模式对应的离子类型包括,正离子模式下的[M+H]+和[M]+离子,负离子模式下的[M-H]-和[M]-
步骤6)中:待测石油馏分的候选单同位素离子注释结果包括质荷比,强度,分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型信息。
步骤7)中:液相色谱-高分辨质谱的文本格式文件包括扫描点数、保留时间以及每个扫描点采集到的离子质荷比及其强度信息。
步骤8)中:峰检测结果文件中包括质荷比,强度,分子类型,分子式,等效双键数,碳数,不同质谱离子化方式对应的离子类型信息以及起始和结束保留时间信息。
步骤10)中,将满足小于等效双键数和碳数连续分布阈值的全部离子进行可视化,包括杂原子类型分布图,同一杂原子类型下的DBE和碳数分布散点图,DBE分布柱状图,碳数分布柱状图,质量缺陷散点图中的至少一种。
基于液相色谱-高分辨质谱的石油分子表征数据处理装置,包括:
去空白模块,用于获取平均质谱图噪音阈值,去除噪音离子;并根据平均质谱图间的质量偏差阈值和质谱响应偏差阈值,对去噪后的待测石油馏分平均谱图进行扣除空白,得到去噪、去空白的待测石油馏分平均谱图;
构建分子库模块,用于获取待测石油馏分的分子类型、等效双键数范围、碳数范围、质荷比范围,并依据石油馏分碳数和等效双键数的极限关系约束最低碳数;生成等效双键数和碳数均连续分布的石油馏分理论分子库,分子库中每个理论分子均包含注释信息,所述注释信息包括分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型及其理论单同位素离子质荷比;
分子注释模块,用于获取平均质谱图中实测离子和理论离子的质荷比偏差阈值,以及实测同位素分布和理论计算同位素分布的匹配阈值;并查询石油馏分组成理论分子式库,从去噪、去空白后的待测石油馏分平均谱图中提取出与石油馏分理论分子库匹配的,且满足小于质荷比偏差阈值的单同位素离子;根据该离子的分子式得到理论同位素的质荷比以及强度,并与样本的平均谱图进行匹配,匹配出满足小于质荷比偏差阈值以及强度阈值的同位素离子;将同时满足质量偏差和同位素分布阈值的离子质荷比及其对应的分子注释信息输出至文件,作为待测石油馏分的候选单同位素离子注释结果;
峰检测模块,用于基于候选离子注释结果,获取原始数据中的实测离子与平均谱图中已注释候选离子的质荷比偏差阈值,以及原始数据中单一离子的连续扫描次数阈值,从原始数据中筛选出满足小于质荷比偏差阈值和连续扫描次数阈值的所有离子;将单一离子连续扫描点的累积强度值作为该离子强度;输出各离子的质荷比,起始扫描点和结束扫描点对应的保留时间,强度以及分子注释信息至文件,作为峰检测结果文件;
去假阳性模块,用于对峰检测结果文件中的离子,提取每个离子所有扫描点的时间和强度,并对其做高斯分布函数,模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子的保留时间;若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,得到其它流出峰与强度最大的流出峰的分离度;根据分离度阈值,去除其中分离度不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件;
处理结果模块,用于根据待测石油馏分的等效双键数和碳数连续分布阈值,在初始假阳性去除结果文件中提取出满足小于等效双键数和碳数连续分布阈值的全部离子,输出至文件,作为最终去假阳性结果文件。
还包括可视化模块,用于将满足小于等效双键数和碳数连续分布阈值的全部离子进行可视化,包括杂原子类型分布图,同一杂原子类型下的DBE和碳数分布散点图,DBE分布柱状图,碳数分布柱状图,质量缺陷散点图中的至少一种。
本发明具有以下有益效果及优点:
本发明公开了一种基于液相色谱-高分辨质谱的石油馏分分子表征数据处理新方法,基于液相色谱-高分辨质谱采集到的原始数据文件,通过构建理论分子库、侯选离子筛选、峰提取与检测、去假阳性等步骤实现石油馏分分子表征。本发明方法充分挖掘了石油样本的液相色谱-高分辨质谱信息,快速、可靠地获得不同石油馏分中各类型化合物的分子式、保留时间及其强度,可用于基于液相色谱-高分辨质谱的石油馏分数据处理,分子表征信息提取更充分、更可靠。
附图说明
图1.待测蜡油样本和空白对照样本的平均质谱图;
图2.基于质谱同位素分布模式的候选离子筛选及分子注释;
图3.蜡油高斯分布函数模拟保留时间图;
图4.蜡油中杂原子类型分布图;
图5.蜡油中N1类化合物的等效双键数和碳数分布散点图。
具体实施方式
下面结合附图及实施例对本发明做进一步的详细说明。
1)自定义用于液相色谱-高分辨质谱数据处理的空白对照样本和待测石油馏分样本数据;
2)输出空白对照样本和石油馏分样本数据在整个保留时间范围的平均质谱图至excel文件,该文件中包括质荷比和强度信息;
3)自定义平均质谱图噪音阈值,去除其中强度低于噪音阈值的离子;
4)自定义平均质谱图间的质量偏差和质谱响应偏差阈值,对步骤3)结果中的石油馏分平均谱图扣除同时满足质量偏差和质谱响应偏差阈值的空白对照;
5)自定义待测石油馏分的分子类型以及等效双键数范围、碳数范围、质荷比范围,并依据石油馏分碳数和等效双键数的极限关系约束最低碳数后,生成等效双键数和碳数均连续分布的石油馏分理论分子库,数据库中的每个理论分子均包含注释信息:分子式,等效双键数,碳数,分子类型,质荷比以及离子类型;
6)自定义平均质谱图中实测离子和理论分子库中离子的质荷比质量偏差阈值,以及实测同位素分布和理论计算同位素分布的匹配阈值;基于石油馏分组成理论分子式库,从步骤4)扣空白后的待测石油馏分平均谱图中提取出与石油馏分理论分子式库匹配的质荷比,匹配出满足质量阈值的单同位素离子根据该离子的分子式计算理论同位素的质荷比以及强度,并与样本的平均谱图进行匹配,匹配出满足质荷比偏差以及强度阈值的同位素离子;将同时满足质量偏差和同位素分布阈值的离子的质荷比及其对应的分子注释信息输出至文件,作为待测石油馏分的候选单同位素离子注释结果,该结果中包括质荷比,强度,分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型信息;
7)读取液相色谱-高分辨质谱数据的原始数据文件,转换并保存为txt文本格式文件,文件中包括扫描点数、保留时间以及每个扫描点采集到的离子质荷比及其强度信息;
8)基于步骤6)的候选离子注释结果,自定义总离子流图中的实测离子与平均谱图中已注释候选离子的质荷比偏差阈值,再自定义单一离子的连续扫描次数阈值,从总离子流中筛选出满足阈值的所有离子;将单一离子连续扫描点的累积强度值作为该离子强度;输出各离子的质荷比m/z,起始扫描点和结束扫描点对应的保留时间,强度以及分子注释信息至文件,作为峰检测结果文件,该文件中包括质荷比,强度,分子类型,分子式,等效双键数,碳数,不同质谱离子化方式对应的离子类型信息以及起始和结束保留时间信息;
9)对峰检测结果文件中的离子,提取每个离子所有扫描点的时间和强度,并对其做高斯分布函数,模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子的保留时间;若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,计算其它流出峰与其的分离度;自定义分离度阈值,去除其中分离度不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件。
10)自定义待测石油馏分的等效双键数和碳数连续分布阈值,在第一次假阳性去除结果文件中提取出满足阈值的全部离子,输出至文件,作为最终去假阳性结果文件。
11)数据可视化:根据峰检测的结果,对数据进行统计并可视化,包括:杂原子类型分布图,同一杂原子类型下的DBE和碳数分布散点图,DBE分布柱状图,碳数分布柱状图,质量缺陷散点图等。
下面通过实例,结合附图对本发明的实施作详细说明:本实施例在以本发明方案为前提下进行实施,给出了详细的实施方式和具体操作过程,但本发明的保护范围不限于下述实施例。
实施例1蜡油液相色谱-电喷雾电离-高分辨质谱数据的处理方法
液相色谱-高分辨质谱得到的待测蜡油样本和空白对照样本的数据文件分别为dgjer.raw和blank-5.raw。
输出平均质谱图:利用XcaliburTM 2.2(Thermo Fisher Scientific,Rockford,U.S.A.)将待测蜡油样本和空白数据在0-30min的平均谱图,如图1A,1B所示,导出至excel中,其中待测蜡油数据得到共235658个质谱离子,空白数据得到227039个质谱离子。
平均质谱谱图去噪、去空白:设置平均质谱图的噪音阈值为1000,从质谱数据中去除噪声信号,以便去除背景信号,提高数据处理速度,去噪后待测样本和空白中分别保留了18589和6177个质谱离子。设置空白和待测蜡油的平均质谱图间的质荷比偏差阈值为2ppm,提高扣除空白信号的精度,避免空白信号的干扰,质谱响应强度的偏差阈值为50%,对待测蜡油样本平均谱图扣空白,去空白后保留16900个质谱离子。
石油馏分理论分子库构建:蜡油液相色谱-电喷雾电离-高分辨质谱采集的主要是蜡油中的杂原子化合物。因此自定义碳数范围为:1-100,等效双键数范围为:1-30,质荷比m/z范围为:150-2000,杂原子分子类型为:Nx(x=1-3)、Ox(x=1-7),Sx(x=1-4)、S1Ox(x=1-5)、N1Ox(x=1-6)、N2Ox(x=1-4)、S2O3、N1Sx(x=1-2)、N1S1O1,限制理论数据库的大小,避免冗余信息;再依据石油馏分组成分子的碳数和等效双键数的极限关系:
碳数>DBE/0.9-氮数
约束最低碳数,生成理论分子库。得到包含分子类型,分子式,等效双键数,碳数,以及在电喷雾电离正、负离子模式下[M+H]+和[M-H]-对应的理论单同位素离子的质荷比。
筛选候选离子及分子注释:将去除空白对照后的平均谱图数据中的16900个特征与理论分子库进行匹配,匹配条件为:平均质谱图中实测离子和理论分子库中离子质荷比的质量偏差≤2ppm,且实测同位素分布模式与理论同位素分布模式相似性阈值≤30,确保分子注释的准确性;相似性计算由归一化的理论同位素和实测同位素强度分布的均方根误差表示。对满足阈值的质谱单同位素离子进行筛选,并基于前述构建的理论分子库对筛选出的单同位素离子进行分子注释。如图2以m/z 334.25258为例说明筛选和分子注释过程,与理论数据库中[C24H31N+H]+单同位素离子质荷比匹配,质荷比偏差Δm=1.05ppm;设置质量偏差2ppm,在实验数据中搜索其同位素峰,找到了对应的M+1峰m/z 335.25592和M+2峰336.25951,相对于单同位素峰的相对强度分别为22%和2%。进一步计算理论同位素和实测同位素强度分布的均方根误差为1.4,满足阈值条件。因此,该离子被筛选为候选离子,注释其分子式为C24H31N,分子类型为N1,等效双键数DBE=10。
共筛选并注释了4235质谱离子用于下一步液相色谱-高分辨质谱数据数据处理。
液相色谱-高分辨质谱数据处理:首先将待测蜡油的液相色谱-高分辨质谱原始数据用MSConvert软件将转换成mzXML格式,然后再利用R语言的apLCMS包将mzXML格式的数据转换成txt文本格式。将前述平均质谱图筛选出的已分子注释的4235个候选单同位素离子依次在txt文件中的各个扫描点中进行搜索,搜索条件为:实测离子与平均谱图中已注释候选单同位素离子的质荷比偏差≤2ppm,且单一离子连续扫描次数≥60次,确保峰检测结果的准确性,降低假阳性率。输出满足搜索条件的离子,并将连续扫描点的起始和终止时间作为该离子的起始和结束保留时间;计算连续扫描点的累积强度作为该离子强度输出,共从原始数据文件提取和识别出3313个候选单同位素离子。
去假阳性:提取出的3313个单同位素离子进一步进行假阳性去除。利用候选离子的连续扫描点,采用高斯分布函数模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子流出峰的保留时间(如图3),以质荷比516.4556为例(分子注释信息为:分子式C37H57N,等效双键数DBE=10),模拟出2个流出峰,对应的保留时间分别为7.12min和15.02min,两个峰的分离度为1.2。若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,计算其它流出峰与其的分离度,设置分离度阈值≤2,确保去假阳性过程真实可靠,避免假阳性过度去除或去除不完全,去除不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件。
再定义等效双键数和碳数的连续阈值分别为3和2,假阳性过度去除或去除不完全,去除同一分子类型,等效双键数和碳数不连续避免的候选离子,作为最终的假阳性去除结果输出,经二步假阳性去除后,共保留3105个单同位素离子为蜡油的液相色谱-高分辨质谱得到分子表征结果。
利用得到的分子表征结果可进一步可视化处理,如图4为该蜡油样本的杂原子类型分布图,从图中可以看出,N1类化合物在蜡油的所有杂原子类型中响应最高,其次是N1O1类化合物。针对N1类化合物,做等效双键数和碳数分布散点图(图5),横坐标为碳数,纵坐标为等效双键数,点的大小代表了化合物的强度,从图中可以看出N1类化合物主要分布范围为等效双键数9-13,碳数10-28。

Claims (10)

1.基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,包括以下步骤:
1)通过液相色谱-高分辨质谱仪获取石油样本的质谱数据,构建用于数据处理的空白对照样本和待测石油馏分样本数据;
2)将空白对照样本和待测石油馏分样本数据通过质谱分析软件获取空白对照样本和石油馏分样本数据在整个保留时间范围的平均质谱图至excel文件;
3)获取平均质谱图噪音阈值,去除噪音离子;
4)根据平均质谱图间的质量偏差阈值和质谱响应偏差阈值,对步骤3)去噪后的待测石油馏分平均谱图进行扣除空白;
5)建立待测石油馏分的分子类型、等效双键数范围、碳数范围、质荷比范围,并依据石油馏分碳数和等效双键数的极限关系约束最低碳数;生成等效双键数和碳数均连续分布的石油馏分理论分子库,分子库中每个理论分子均包含注释信息,所述注释信息包括分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型及其理论单同位素离子质荷比;
6)获取平均质谱图中实测离子和理论离子的质荷比偏差阈值,以及实测同位素分布和理论计算同位素分布的匹配阈值;查询石油馏分组成理论分子式库,从步骤4)去噪、去空白后的待测石油馏分平均谱图中提取出与石油馏分理论分子库匹配的,且满足小于质荷比偏差阈值的单同位素离子;根据该离子的分子式得到理论同位素的质荷比以及强度,并与样本的平均谱图进行匹配,匹配出满足小于质荷比偏差阈值以及强度阈值的同位素离子;将同时满足质量偏差和同位素分布阈值的离子质荷比及其对应的分子注释信息输出至文件,作为待测石油馏分的候选单同位素离子注释结果;
7)读取液相色谱-高分辨质谱仪所采集数据的原始数据文件,转换并保存为txt文本格式文件;
8)基于步骤6)的候选离子注释结果,获取原始数据中的实测离子与平均谱图中已注释候选离子的质荷比偏差阈值,以及原始数据中单一离子的连续扫描次数阈值,从原始数据中筛选出满足小于质荷比偏差阈值和连续扫描次数阈值的所有离子;将单一离子连续扫描点的累积强度值作为该离子强度;输出各离子的质荷比,起始扫描点和结束扫描点对应的保留时间,强度以及分子注释信息至文件,作为峰检测结果文件;
9)对峰检测结果文件中的离子,提取每个离子所有扫描点的时间和强度,并对其做高斯分布函数,模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子的保留时间;若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,得到其它流出峰与强度最大的流出峰的分离度;根据分离度阈值,去除其中分离度不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件;
10)根据待测石油馏分的等效双键数和碳数连续分布阈值,在初始假阳性去除结果文件中提取出满足小于等效双键数和碳数连续分布阈值的全部离子,输出至文件,作为最终去假阳性结果文件。
2.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤2)中:平均谱图中包括质荷比和强度信息。
3.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤4)中:同时满足平均质谱图间的质量偏差和质谱响应偏差阈值的石油馏分平均谱图信号将被去除。
4.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤5)中:不同质谱离子化模式对应的离子类型包括,正离子模式下的[M+H]+和[M]+离子,负离子模式下的[M-H]-和[M]-
5.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤6)中:待测石油馏分的候选单同位素离子注释结果包括质荷比,强度,分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型信息。
6.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤7)中:液相色谱-高分辨质谱的文本格式文件包括扫描点数、保留时间以及每个扫描点采集到的离子质荷比及其强度信息。
7.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤8)中:峰检测结果文件中包括质荷比,强度,分子类型,分子式,等效双键数,碳数,不同质谱离子化方式对应的离子类型信息以及起始和结束保留时间信息。
8.根据权利要求1所述的基于液相色谱-高分辨质谱的石油分子表征数据处理方法,其特征在于,步骤10)中,将满足小于等效双键数和碳数连续分布阈值的全部离子进行可视化,包括杂原子类型分布图,同一杂原子类型下的DBE和碳数分布散点图,DBE分布柱状图,碳数分布柱状图,质量缺陷散点图中的至少一种。
9.基于液相色谱-高分辨质谱的石油分子表征数据处理装置,其特征在于,包括:
去空白模块,用于获取平均质谱图噪音阈值,去除噪音离子;并根据平均质谱图间的质量偏差阈值和质谱响应偏差阈值,对去噪后的待测石油馏分平均谱图进行扣除空白,得到去噪、去空白的待测石油馏分平均谱图;
构建分子库模块,用于获取待测石油馏分的分子类型、等效双键数范围、碳数范围、质荷比范围,并依据石油馏分碳数和等效双键数的极限关系约束最低碳数;生成等效双键数和碳数均连续分布的石油馏分理论分子库,分子库中每个理论分子均包含注释信息,所述注释信息包括分子类型,分子式,等效双键数,碳数以及不同质谱离子化方式对应的离子类型及其理论单同位素离子质荷比;
分子注释模块,用于获取平均质谱图中实测离子和理论离子的质荷比偏差阈值,以及实测同位素分布和理论计算同位素分布的匹配阈值;并查询石油馏分组成理论分子式库,从去噪、去空白后的待测石油馏分平均谱图中提取出与石油馏分理论分子库匹配的,且满足小于质荷比偏差阈值的单同位素离子;根据该离子的分子式得到理论同位素的质荷比以及强度,并与样本的平均谱图进行匹配,匹配出满足小于质荷比偏差阈值以及强度阈值的同位素离子;将同时满足质量偏差和同位素分布阈值的离子质荷比及其对应的分子注释信息输出至文件,作为待测石油馏分的候选单同位素离子注释结果;
峰检测模块,用于基于候选离子注释结果,获取原始数据中的实测离子与平均谱图中已注释候选离子的质荷比偏差阈值,以及原始数据中单一离子的连续扫描次数阈值,从原始数据中筛选出满足小于质荷比偏差阈值和连续扫描次数阈值的所有离子;将单一离子连续扫描点的累积强度值作为该离子强度;输出各离子的质荷比,起始扫描点和结束扫描点对应的保留时间,强度以及分子注释信息至文件,作为峰检测结果文件;
去假阳性模块,用于对峰检测结果文件中的离子,提取每个离子所有扫描点的时间和强度,并对其做高斯分布函数,模拟每个离子的流出峰,将峰顶点对应的保留时间作为该离子的保留时间;若同一离子模拟出多个流出峰,则以强度最大的流出峰为基准,得到其它流出峰与强度最大的流出峰的分离度;根据分离度阈值,去除其中分离度不满足阈值的流出峰,重新输出至文件,作为初始去假阳性结果文件;
处理结果模块,用于根据待测石油馏分的等效双键数和碳数连续分布阈值,在初始假阳性去除结果文件中提取出满足小于等效双键数和碳数连续分布阈值的全部离子,输出至文件,作为最终去假阳性结果文件。
10.根据权利要求9所述的基于液相色谱-高分辨质谱的石油分子表征数据处理装置,其特征在于,还包括可视化模块,用于将满足小于等效双键数和碳数连续分布阈值的全部离子进行可视化,包括杂原子类型分布图,同一杂原子类型下的DBE和碳数分布散点图,DBE分布柱状图,碳数分布柱状图,质量缺陷散点图中的至少一种。
CN202011447562.2A 2020-12-09 2020-12-09 基于液相色谱-高分辨质谱的石油分子表征数据处理方法 Active CN114609258B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011447562.2A CN114609258B (zh) 2020-12-09 2020-12-09 基于液相色谱-高分辨质谱的石油分子表征数据处理方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011447562.2A CN114609258B (zh) 2020-12-09 2020-12-09 基于液相色谱-高分辨质谱的石油分子表征数据处理方法

Publications (2)

Publication Number Publication Date
CN114609258A true CN114609258A (zh) 2022-06-10
CN114609258B CN114609258B (zh) 2022-11-22

Family

ID=81857048

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011447562.2A Active CN114609258B (zh) 2020-12-09 2020-12-09 基于液相色谱-高分辨质谱的石油分子表征数据处理方法

Country Status (1)

Country Link
CN (1) CN114609258B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116046968A (zh) * 2023-03-02 2023-05-02 北京澳诺科技有限公司 一种液相色谱工作站数据处理方法、系统及可存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110870985A (zh) * 2018-08-31 2020-03-10 中国石油化工股份有限公司 重油芳烃组分的分离、表征方法和二维液相色谱及装置

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110870985A (zh) * 2018-08-31 2020-03-10 中国石油化工股份有限公司 重油芳烃组分的分离、表征方法和二维液相色谱及装置

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CAROLE REYMOND等: "Input of an Off-Line, Comprehensive, Three-Dimensional Method (CPC×SFC/HRMS) to Quantify Polycyclic Aromatic Hydrocarbons in Vacuum Gas Oils", 《ANALYTICAL CHEMISTRY》 *
SAMI LABABIDI等: "Online normal-phase high-performance liquid chromatography/Fourier transform ion cyclotron resonance mass spectrometry:Effects of different ionization methods on the characterization of highly complex crude oil mixtures", 《RAPID COMMUN. MASS SPECTROM》 *
YUEYI XIA等: "Characterization of nitrogen-containing compounds in petroleum fractions by online reversed-phase liquid chromatography-electrospray ionization Orbitrap mass spectrometry", 《FUEL》 *
刘美等: "沙特常压渣油中含硫化合物的分子组成表征", 《化工进展》 *
马苏甜等: "基于质谱技术的重油分子组成表征研究进展", 《石油化工》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116046968A (zh) * 2023-03-02 2023-05-02 北京澳诺科技有限公司 一种液相色谱工作站数据处理方法、系统及可存储介质

Also Published As

Publication number Publication date
CN114609258B (zh) 2022-11-22

Similar Documents

Publication Publication Date Title
US10198630B2 (en) Peak detection method
CN107860845B (zh) 自动解析gc-ms重叠峰准确识别化合物的方法
US6449584B1 (en) Measurement signal processing method
CN105334279B (zh) 一种高分辨质谱数据的处理方法
CN108875170B (zh) 一种基于改进变分模态分解的噪声源识别方法
CN110243806B (zh) 拉曼光谱下基于相似度的混合物组分识别方法
US20210088495A1 (en) Edible oil analysis system and method
WO2021174901A1 (zh) 基于数据非依赖采集质谱的分子组学数据结构的实现方法
CN114609258B (zh) 基于液相色谱-高分辨质谱的石油分子表征数据处理方法
JP2020527695A (ja) 質量分析装置におけるデータ取得方法
EP3254126B1 (en) Interference detection and peak of interest deconvolution
CN113554176B (zh) 代谢特征谱推断方法、系统、计算机设备及存储介质
CN113194332B (zh) 基于多策略的新广告发现方法、电子设备和可读存储介质
JP6738816B2 (ja) 曲線減算を介する類似性に基づく質量分析の検出
CN114858958B (zh) 质谱数据在质量评估中的分析方法、装置和存储介质
CN113720952B (zh) 用于储层解释评价的图版生成方法、装置、设备及介质
CN111896609B (zh) 一种基于人工智能分析质谱数据的方法
CN113921086A (zh) 基于质谱分析的蛋白质从头肽测序方法及系统
CN109917245B (zh) 考虑相位差的超声检测局部放电信号模式识别方法
US10236167B1 (en) Peak waveform processing device
CN108108721A (zh) 一种利用高光谱进行道路提取的方法
WO2006125863A1 (en) Analysis techniques for liquid chromatography/mass spectrometry
CN109829513B (zh) 一种顺序式波长色散x射线荧光光谱智能分析方法
WO2005015199A1 (en) Methods and systems for chromatography/mass-spectrometry analysis
CN113889196B (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