CN113188942B - 一种用于含pvc混合塑料热裂解过程中的活化能分析方法 - Google Patents

一种用于含pvc混合塑料热裂解过程中的活化能分析方法 Download PDF

Info

Publication number
CN113188942B
CN113188942B CN202110371790.4A CN202110371790A CN113188942B CN 113188942 B CN113188942 B CN 113188942B CN 202110371790 A CN202110371790 A CN 202110371790A CN 113188942 B CN113188942 B CN 113188942B
Authority
CN
China
Prior art keywords
activation energy
pyrolysis stage
pyrolysis
thermal cracking
pvc
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
CN202110371790.4A
Other languages
English (en)
Other versions
CN113188942A (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.)
Qingdao University of Science and Technology
Original Assignee
Qingdao University of Science and Technology
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 Qingdao University of Science and Technology filed Critical Qingdao University of Science and Technology
Priority to CN202110371790.4A priority Critical patent/CN113188942B/zh
Publication of CN113188942A publication Critical patent/CN113188942A/zh
Application granted granted Critical
Publication of CN113188942B publication Critical patent/CN113188942B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N5/00—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid
    • G01N5/04—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder
    • G01N5/045—Analysing materials by weighing, e.g. weighing small particles separated from a gas or liquid by removing a component, e.g. by evaporation, and weighing the remainder for determining moisture content
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00—Sampling; Preparing specimens for investigation
    • G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/40—Concentrating samples
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00—Sampling; Preparing specimens for investigation
    • G01N1/28—Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/44—Sample treatment involving radiation, e.g. heat
    • 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
    • Y02W—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO WASTEWATER TREATMENT OR WASTE MANAGEMENT
    • Y02W30/00—Technologies for solid waste management
    • Y02W30/50—Reuse, recycling or recovery technologies
    • Y02W30/62—Plastics recycling; Rubber recycling

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (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 Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开一种用于含PVC混合塑料热裂解的活化能分析方法，具有在混合塑料热裂解过程中动力学方法分析数据处理部分的应用，并提供了一种将含PVC混合塑料热裂解过程分阶段的分析方法。本发明选择了混合塑料热裂解分析曲线的独立分析方法，数据处理的计算方法中采用了三种无模型配合法计算混合塑料热裂解中活化能，可以在先不涉及动力学模型的情况下，直接求取混合塑料热裂解过程中的活化能，从而避免了因为假设的反应机理函数不同，产生计算结果的误差，揭示了含PVC混合塑料热裂解过程活化能于转化率的关系，具有高效性和准确性。

Description

一种用于含PVC混合塑料热裂解过程中的活化能分析方法
技术领域
本发明涉及一种含PVC混合塑料热裂解过程中的活化能的分析方法，属于废旧塑料热裂解绿色资源化再利用领域。
背景技术
中国是一个塑料生产、消耗和回收的大国，年消费量位居世界前位。塑料的使用寿命较短，在大量塑料生产的背景下，塑料的绿色循环回收便成为一个需要解决的问题。因为塑料制品的特殊性，在自然环境下降解的效率非常低，而且采用堆积填埋的方式会对环境造成巨大的污染,是一种耗资高、效率低的方法。塑料的绿色回收就变成了一个需要重视的问题。
以往塑料制品的回收方式主要是填埋、焚烧、堆肥化和回收再生。其中回收再生所处理的废旧塑料仅仅占少部分，填埋焚烧等方式为塑料制品的主要回收方式。而填埋焚烧会对环境造成严重污染，效益较低，而在焚烧等方式中产生的有害物质会对环境造成二次污染，威胁人的身体健康。将废塑料高效无污染的绿色回收并实现循环再利用，成为了一个循环经济及节能减排领域亟待解决的重要问题；废旧塑料回收再生中，塑料的热裂解作为一种环保、高效的绿色回收方式，可以将废旧塑料转化成燃料油、燃料气等不可再生资源，是一种经济、高效、环保的塑料绿色回收方式，实现了废塑料的无害化、资源化和减量化，具有巨大的节能于回收潜力。
在研究废旧混合塑料热裂解中，需要对混合塑料热裂解的过程进行实验研究。研究废旧塑料在热解过程中的变化就需要对裂解过程进行动力学分析，为混合塑料热裂解等热性能提供理论数据。
热分析动力学研究的三要素分别为活化能、指前因子和机理函数。确定动力学三要素可以选择对塑料热裂解进行动力学模型函数的假设，一般会选择假定为简单的级数反应，采用单升温速率法，在一个升温速率下，求取表观反应级数、活化能和指前因子，而由于混合塑料热裂解实际过程是非常复杂的，通常是采取实验数据和动力学模型相配合的方式描述这个热分解反应过程，此过程称为模型配合法。在模型配合的过程中难免会产生多种动力学模型函数都可以配合同一反应过程曲线的情况发生。
活化能作为动力学三要素之一，对于准确确定含PVC混合塑料热裂解过程的指前因子和机理函数，以及准确解释反应过程具有重要的意义。
论文Thermal degradation of typical plastics under high heating rateconditions by TG-FTIR:Pyrolysis behaviors and kinetic analysis一文中，以低密度聚乙烯(LDPE)、聚丙烯(PP)和聚氯乙烯(PVC)为研究对象，在高加热速率条件下研究了它们的热解行为，并找出最适宜的动力学反应机理。低密度聚乙烯的TG和DTG曲线与聚丙烯相似，但与聚氯乙烯不同。结果表明，高升温速率导致初始温度、终止温度和峰值温度向高温方向移动。红外光谱结果表明，低密度聚乙烯和聚丙烯的主要产物为烷烃和烯烃，聚氯乙烯的主要产物为盐酸、烯烃和少量芳香族化合物。对无模型方法如Ozawa-Flynn-Wall(OFW)、Kissinger-Akahira-Sunose(KAS)和Friedman(FR)方法计算活化能进行了比较研究。利用Coats-Redfern和Criado等模型拟合方法对反应动力学机理进行了预测。将低密度聚乙烯的热解反应机理归纳为收缩球模型，将聚乙烯的热解反应归纳为收缩圆柱模型，而聚氯乙烯的热解反应机理归纳为第一阶段为二维成核，第二阶段为三维扩散(简化)。但是对于PP/PVC/PS混合塑料，因动力学模型假设不同而带来的对整个含PVC混合塑料热裂解过程的动力学参数的误差，并没有揭示含PVC混合塑料热裂解过程活化能于转化率的关系。
发明内容
有鉴于此，本发明针对上述问题，提出一种用于含PVC混合塑料热裂解过程中的活化能分析方法。
为了实现上述目的，本发明采用如下技术方案实现：一种用于含PVC混合塑料热裂解过程中的活化能分析方法，包含以下步骤：
(1)选择塑料种类和混合塑料的配比，对混合塑料试样进行干燥预处理，经过破碎，研磨处理后，用天平称重分组；
(2)采用多升温速率法，设定混合塑料热裂解过程的升温区间、升温速率以及每组实验试样实验用量；
(3)对含PVC混合塑料热裂解实验试样进行热重实验分析，得到多升温速率下的TG曲线和DTG曲线及相关实验数据；
(4)对TG和DTG曲线进行数值模拟分析，结合含PVC混合塑料热裂解的实验过程，以热解速率近乎为零的时间节点或者热解过程中失重率近乎不变的时间节点为热裂解过程分为第一热解阶段和第二热解阶段的分界节点；
(5)选用无模型配合法计算热裂解的活化能，将含PVC混合塑料热裂解的TG和DTG曲线得到的实验试样数据代入到无模型配合法的分析方程，多段升温速率下构成的分析方程用origin进行线性回归拟合，得到拟合直线图像；
(6)将回归拟合方程近似为线性方程，以线性方程参数于热解阶段活化能之间的函数关系，计算出含PVC混合塑料热裂解第一热解阶段和第二热解阶段的活化能。
优选的，为避免模型拟合法的不确定性，选用无模型配合法直接求取动力学参数；
无模型配合法又称为等转化率法，无模型配合法采用的是多升温速率法，即在多条升温速率曲线下进行热分析动力学处理的方法。
优选的，(1)步骤中，首先是选择含PVC混合塑料热裂解的实验试样组成成分，要对实验试样进行干燥、破碎和研磨预处理后，按照比例均匀混合后选取m克放入实验用的坩埚中；
进一步的，实验前对热重分析仪进行预设置。
优选的，(2)步骤中，实验反应的环境采用高纯度氮气进行吹扫气，设置吹扫气流量，设置实验温度上升区间T1～T2，选择三组到四组升温速率，分别为β1 K/min、β2 K/min、β3K/min、β4 K/min。为消除环境因素影响，每组实验采取重复三次取平均值的方式，保证实验得到的曲线真实可靠。
在热重法分析含PVC混合塑料热裂解的过程中，转化率为：
Figure GDA0003995188670000031
式中：m0为热分解前塑料样品的质量，mt为反应过程中t时刻的试样质量，mf为反应结束后试样的质量。
2)采用多升温速率法对混合塑料的热裂解在不同升温速率过程下，进行塑料热裂解动力学分析；
优选的，(3)步骤中，通过热重实验分析得到不同升温速率下的TG和DTG曲线，由TG和DTG曲线获得失重率α和温度T的实验数据及关系图。
优选的，(4)步骤中，在含PVC混合塑料热裂解整个过程中，以连续两个热解速率的高峰中间热解速率几乎为零的时间节点为整个热解过程的第一热解阶段和第二热解阶段的分界节点。
3)优选的，(3)步骤中，取含PVC混合塑料热裂解过程中的失重率区间为0.1～0.9，在β1 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段的温度分别为TAi(i＝1，2，3...9)，第二个热解阶段的温度分别为Tai(i＝1，2，3...9)；在β2 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为TBi(i＝1，2，3...9)，第二个热解阶段温度分别为Tbi(i＝1，2，3...9)；在β3 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为TCi(i＝1，2，3...9)，第二个热解阶段温度分别为Tci(i＝1，2，3...9)；在β4 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段的温度分别为TDi(i＝1，2，3...9)，第二个热解阶段的温度分别为Tdi(i＝1，2，3...9)。三种升温速率下各失重率所对应温度数据如表1所示：
表1各升温速率下失重率对应的温度
Figure GDA0003995188670000041
4)优选的，(5)步骤中，将表格中的实验数据制成失重率α和温度T的关系图，将实验数据带入到以下三种无模型配合法的分析方法：
热分析曲线的数据分析主要分为积分法和微分法，积分法代表为Flynn-Wall-Ozawa(F-W-O)法，微分法代表为Kissinger-Akahira-Sunose(K-A-S)法和Starink法。这三种分析方法方程如下：
(1)Flynn-Wall-Ozawa(F-W-O)法：
Figure GDA0003995188670000042
其中B＝1.052，G(α)为积分形式的机理函数。
(2)Kissinger-Akahira-Sunose(K-A-S)法：
Figure GDA0003995188670000043
其中B＝1，G(α)为积分形式的机理函数。
(3)Starink法：
Figure GDA0003995188670000044
其中B＝1.0037，G(α)为积分形式的机理函数。
以上三种分析方法可以采用一个通式方程进行表示：
Figure GDA0003995188670000045
其中对F-W-O方程，U＝0，B＝1.052，C＝ln[AE/RG(α)]-5.331；K-A-S方程，U＝2，B＝1，C＝ln[AR/EG(α)]；Starink法，U＝1.8，B＝1.0037，C＝ln[AR/EG(α)]。
第一热解阶段：
将失重率α和温度T的实验数据带入到三种无模型配合法的分析方程(2)(3)(4)中，以同一失重率不同升温速率为一组数据，用origin进行曲线回归拟合，分别得到关于lm(β)、ln(β/T2)及lm(β/T1.8)与1/T之间的线性关系曲线。再以每个失重率下不同升温速率下，构建一个关于x、y线性方程。即：
y＝ax+b (6)
其中：假设y＝ln(β/TU)；x＝1/T；将a＝-BE/R；b＝C。分别对应于公式(5)
优选的，(6)步骤中，结合拟合曲线图像的获得斜率a，通过斜率a和E的关系式求得活化能。由失重率α从0.1～0.9对应共9组拟合直线。
第一个热解阶段的每组拟合直线的活化能为：
Eli＝-aR/B(i＝1，2，3...9) (7)
将第一个热解阶段9组活化能求取平均值，得到第一个热解阶段活化能E1st：
Figure GDA0003995188670000051
第二热解阶段：
与第一个热解阶段方法类似，将第二个热解阶段的失重率α和温度T的实验数据带入到三种无模型配合法的分析方程，得到关于ln(β)、ln(β/T2)及ln(β/T1.8)与1/T之间的线性关系曲线，对线性关系曲线进行线性回归拟合得到一组线性方程。即构建一个关于x、y线性方程，对应公式(6)。
5)结合拟合曲线图像的获得斜率a，通过斜率a和E的关系式求得活化能。第二个热解阶段的每组拟合直线的活化能为：
E2i＝-aR/B(i＝1，2，3...9) (9)
将第二个热解阶段9组活化能值求取平均值，得到第二个热解阶段活化能E2nd：
Figure GDA0003995188670000052
6)通过热解曲线的分析，这两段总体上的热解反应速率及反应物的热解变化性质具有一些差异，所以要将这两个阶段分开进行活化能及相关参数的确定，方便后续对于含PVC塑料热裂解过程中对两个阶段进行定性和定量分析。
本发明的技术方案至少具有如下优点和有益效果：
1.本发明不需要提前假设含PVC混合塑料热裂解过程的动力学模型，避免了动力学模型假设不同而带来的对整个含PVC混合塑料热裂解过程的动力学参数的误差，有利于揭示含PVC混合塑料热裂解过程活化能于转化率的关系，又在一定程度上解释了热解反应的复杂性。
2.本发明根据在实际生活中的废旧塑料大多为多种塑料制品混合组成的现状，本次实验试样组分选择含PVC混合塑料热裂解过程进行研究，是塑料制品的典型成分构成，对废旧塑料经过热裂解进行绿色循环回收具有代表作用，能够清晰地表征含PVC混合塑料热裂解过程中的活化能。
3.本发明实现了结合热解反应过程的特征温度，将整个热解过程分为第一热解阶段和第二热解阶段，选择在热解速率几乎为零的时间节点为第一热解阶段和第二热解阶段的分界节点，或者为整个热解阶段中失重率近乎不变的温度阶段为分界节点，从而提出一种用于含PVC混合塑料热裂解过程中的活化能计算方法。
附图说明
图1为本发明PP/PVC/PS混合塑料在非等温条件下的DTG曲线图；
图2本发明PP/PVC/PS混合塑料在非等温条件下的TG曲线图；
图3为本发明PP/PVC/PS混合塑料失重率和温度的关系曲线图；
图4为PP/PVC/PS混合塑料F-W-O法第一热解阶段拟合曲线图；
图5为PP/PVC/PS混合塑料F-W-O法第二热解阶段拟合曲线图；
图6为PP/PVC/PS混合塑料K-A-S法第一热解阶段拟合曲线图；
图7为PP/PVC/PS混合塑料K-A-S法第二热解阶段拟合曲线图；
图8为PP/PVC/PS混合塑料Starink法第一热解阶段拟合曲线图；
图9为PP/PVC/PS混合塑料Starink法第二热解阶段拟合曲线图；
图10为一种用于含PVC混合塑料热裂解过程中的活化能计算方法的流程图；
具体实施方式
下面结合附图和具体实施例对本发明进行说明。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图，对本发明实施例中的技术方案进行清楚、完整的描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。
因此，以下对本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的部分实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。
还需要说明的是，在本发明的描述中，除非另有明确的规定和限定，对于本领域技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
如图1-10所示，
作为本发明一个优选的实施例，本实验实施例采用PP、PVC和PS的三种混合塑料制品进行试验，这三种塑料为中国塑料制品的主要构成。首先将这三种塑料制品放在105℃的环境中静止干燥12个小时，再分组放入破碎机中进行破碎处理，经过研磨混合后，用天平进行称重，按照配比放入离心管作为热失重实验样品组。每次实验取固定的量放入坩埚中。实验混合塑料试样配比及实验用量如下表2所示：
表2含PVC混合塑料给组分配比及实验用量
Figure GDA0003995188670000071
实验前要对热重分析仪进行预调整。实验反应的环境采用99.9999％的高纯度氮气进行吹扫气，从而排出空气，设置吹扫气流量，设置实验温度上升区间T1～T2，选择三组升温速率，分别为β1 K/min、β2 K/min、β3 K/min。为消除环境因素影响，每组实验采取重复三次取平均值的方式，保证实验得到的曲线真实可靠。
图1和图2分别为含PVC的混合塑料热裂解在三组不同升温速率下的DTG曲线、TG曲线，结合含PVC混合塑料热裂解的TG和DTG曲线，以含PVC的混合塑料热裂解过程是整体分为两个阶段，本次混合塑料热裂解的活化能将依照热解第一阶段和第二阶段进行分别求取。由DTG和TG图像曲线得到的实验统计数据如表3所示：
表3在三种升温速率下各失重率对应的温度
Figure GDA0003995188670000072
取含PVC混合塑料热裂解过程中的失重率区间为0.1～0.9，在β1 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段的温度分别为TAi(i＝1，2，3...9)，第二个热解阶段的温度分别为Tai(i＝1，2，3...9)；
在β2 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为TBi(i＝1，2，3...9)，第二个热解阶段温度分别为Tbi(i＝1，2，3...9)；
在β3 K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为Tci(i＝1，2，3...9)，第二个热解阶段温度分别为Tci(i＝1，2，3...9)。
由实验得到的含PVC混合塑料热裂解两个热解阶段的实验数据，做成失重率与实验温度的关系图，如图3所示。
作为本发明一个优选的实施例，将实验DTG和TG图像中的得到的表3的实验数据，分别代入以下三种无模型方法进行活化能的计算：
(1)F-W-O法：
将三种升温速率下β1 K/min、β2 K/min、β3 K/min升温速率下第一个热解阶段的ln(β1)与1/TAi(i＝1，2，3...9)、ln(β2)与1/TBi(i＝1，2，3...9)、ln(β3)与1/Tci(i＝1，2，3...9)的数据代入到建立的以ln(β)为纵坐标，1/T为横坐标的坐标系中，以同一失重率不同升温速率为一组数据进行曲线回归拟合，共得到9组且近乎平行拟合直线，每一组拟合直线对应一个斜率，典型拟合图像如4所示。根据F-W-O法分析方程公式(2)：
其中C＝ln[AE/RG(α)]-5.331，A为指前因子，E为活化能，R为气体常数，G(α)为机理函数，B＝1.052，C在此认为是一个定值；ln(β)为因变量y，1/T为自变量x，-BE/R为该线性方程的斜率，由此可构成线性方程(6)。
根据图4拟合直线图像，失重率为0.1～0.9的9组拟合直线斜率为a1i(i＝1，2，3...9)，将其代入到公式(7)，得到每组拟合直线的活化能E1i。
将9组活化能值代入到公式(8)，求其平均值即为通过F-W-O法求得的含PVC混合塑料热裂解过程的第一个热解阶段的稳定活化能值E1st。
第一个热解阶段的ln(β)和1/T的拟合曲线图如图4所示：
第二个热解阶段数据用以上方法代入到F-W-O法分析方程中，将三种升温速率下β1 K/min、β2 K/min、β3 K/min升温速率下第二个热解阶段的ln(β1)与1/Tai(i＝1，2，3...9)、ln(β2)与1/Tbi(i＝1，2，3...9)、ln(β3)与1/Tci(i＝1，2，3...9)的数据带入F-W-O法分析方程(2)中。
并建立的以ln(β)为纵坐标，1/T为横坐标的坐标系中，R为气体常数，G(α)为机理函数，B＝i.052，C在此认为是一个定值；ln(β)为因变量y，1/T为自变量x，-BE/R为该线性方程的斜率，由此可构成线性方程(6)。
第二个热解阶段的ln(β)和1/T的拟合曲线图如图5所示。根据图5，失重率为0.1～0.9的9组拟合直线斜率为a2i(i＝1，2，3...9)，将9组拟合直线斜率代入到公式(9)，得到每组拟合直线的活化能为E2i。
将9组活化能值代入到公式(10)，求其平均值即为通过F-W-O法求得的含PVC混合塑料热裂解过程的第二个热解阶段的稳定活化能值E2nd。
(2)K-A-S法：
与F-W-O法类似，以ln(β/T2)为纵坐标，1/T为横坐标，在β1 K/min、β2 K/min、β3 K/min三种升温速率下，将ln(β1/T2)与1/TAi(i＝1，2，3...9)、ln(β2/T2)与1/TBi(i＝1，2，3...9)、ln(β3/T2)与1/Tci(i＝1，2，3...9)的第一热解阶段的实验数据带入到上述以ln(β/T2)与1/T的坐标系中，以通以失重率下的不同升温速率为一个组进行曲线的回归拟合，得到9组斜率相近的拟合直线，每一组拟合直线对应的一个斜率，典型拟合图像如图6所示。根据K-A-S法分析方程(3)中。
其中C＝ln[AR/EG(α)]，A为指前因子，E为活化能，R为气体常数，G(α)为机理函数，B＝i，C认为是一个定值；ln(β/T2)为因变量y，1/T为自变量x，-BE/R为该线性方程的斜率，由此构成的线性方程(6)。
根据拟合直线图像，求出9组拟合直线的斜率α1i(i＝1，2，3...9)，将其代入到公式(7)，得到每组拟合直线的活化能E1i。
将9组活化能值代入到公式(8)，求其平均值即为通过K-A-S法求得的稳定活化能值E1st。
第一热解阶段的ln(β/T2)和1/T的拟合曲线图如图6所示：
第二个热解阶段数据用以上方法代入到K-A-S法分析方程中，将三种升温速率下β1K/min、β2 K/min、β3 K/min升温速率下第二个热解阶段的ln(β1/T2)与1/Tai(i＝1，2，3...9)、ln(β2/T2)与1/Tbi(i＝1，2，3...9)、ln(β3/T2)与1/Tci(i＝1，2，3...9)的数据带入K-A-S法分析方程(3)中。
并建立的以ln(β/T2)为纵坐标，1/T为横坐标的坐标系中，其中C＝lm[AR/EG(α)]，A为指前因子，E为活化能，R为气体常数，G(α)为机理函数，B＝1，C认为是一个定值；ln(β/T2)为因变量y，1/T为自变量x，-BE/R为该线性方程的斜率，由此构成的线性方程为(6)。
第二个热解阶段的ln(β/T2)和1/T的拟合曲线图如图7所示。根据拟合直线图像，失重率为0.1～0.9的9组拟合直线斜率为a2i(i＝1，2，3...9)，将9组拟合直线斜率代入到公式(9)，得到每组拟合直线的活化能为E2i。
将9组活化能值代入到公式(10)，求其平均值即为通过K-A-S法求得的含PVC混合塑料热裂解过程的第二个热解阶段的稳定活化能值E2nd：
(3)Starink法：
与上述两种热分析方法类似，建立以ln(β/T1.8)为纵坐标，以1/T为横坐标的坐标系，在建立的坐标中做出在β1 K/min、β2 K/min、β3 K/min三种升温速率下，将ln(β1/T1.8)与1/TAi(i＝1，2，3...9)、ln(β2/T1.8)与1/TBi(i＝1，2，3...9)、ln(β3/T1.8)与1/Tci(i＝1，2，3...9)的第一热解阶段的实验数据带入到上述以ln(β/T1.8)与1/T的坐标系中，以不同升温速率同一失重率为一组进行曲线的回归拟合，共拟合得到9组直线，每一组直线对应在此失重率下的拟合直线的斜率，典型拟合图像如图8所示。根据Starink法分析方程(4)中。
其中C＝ln[AR/EG(α)]，A为指前因子，E为活化能，R为气体常数，G(α)为机理函数，B＝1.0037，C在此认为是一个定值。ln(β/T1.8)为因变量y，1/T为自变量x，-BE/R为该线性方程的斜率，此可构成线性方程(6)。
根据拟合直线图像，失重率为0.1～0.9的9组拟合直线斜率为a1i(i＝1，2，3...9)，将其代入到公式(7)，得到每组拟合直线的活化能E1i。
将9组活化能值代入到公式(8)，求其平均值即为通过K-A-S法求得的稳定活化能值E1st。
第一个热解阶段的ln(β)和1/T的拟合曲线图如图8所示：
第二个热解阶段数据用以上方法代入到Starink法分析方程中，将三种升温速率下β1 K/min、β2 K/min、β3 K/min升温速率下第二个热解阶段的ln(β1/T1.8)与1/Tai(i＝1，2，3...9)、lm(β2/T1.8)与1/Tbi(i＝1，2，3...9)、ln(β3/T1.8)与1/Tci(i＝1，2，3...9)的数据带入Starink法分析方程(4)中。
并建立的以ln(β/T1.8)为纵坐标，1/T为横坐标的坐标系中，以同一失重率不同升温速率为一组数据进行曲线回归拟合，构成线性方程(6)。
第二个热解阶段的ln(β/T1.8)和1/T的拟合曲线图如图9所示。根据拟合直线图像，失重率为0.1～0.9的9组拟合直线斜率为a2i(i＝1，2，3...9)，将9组拟合直线斜率代入到公式(9)，得到每组拟合直线的活化能为E2i。
将9组活化能值代入到公式(10)，求其平均值即为通过Starink法求得的含PVC混合塑料热裂解过程的第二个热解阶段的稳定活化能值E2nd。
以上实施方案仅用于说明而非限制本发明的技术方案。不脱离本发明精神的任何修改或局部替换，均应涵盖在本发明的权利要求范围当中。

Claims (1)

1.一种用于含PVC混合塑料热裂解过程中的活化能分析方法，其特征在于：所述PVC混合塑料为PP/PVC/PS混合塑料，包括以下步骤：
(1)选择塑料种类和混合塑料的配比，对混合塑料试样进行干燥预处理，经过破碎，研磨处理后，用天平称重分组；
(2)采用多升温速率法，设定混合塑料热裂解过程的升温区间、升温速率以及每组实验试样实验用量；
(3)对含PVC混合塑料热裂解实验试样进行热重实验分析，得到多升温速率下的TG曲线和DTG曲线及相关实验数据；
(4)对TG和DTG曲线进行数值模拟分析，结合含PVC混合塑料热裂解的实验过程，以热解速率接近零的时间节点或者热解过程中失重率接近不变的时间节点为热裂解过程分为第一热解阶段和第二热解阶段的分界节点；
(5)选用无模型配合法计算热裂解的活化能，将含PVC混合塑料热裂解的TG和DTG曲线得到的实验试样数据代入到无模型配合法的分析方程，多段升温速率下构成的分析方程用origin进行线性回归拟合，得到拟合直线图像；
(6)将回归拟合方程近似为线性方程，以线性方程参数于热解阶段活化能之间的函数关系，计算出含PVC混合塑料热裂解第一热解阶段和第二热解阶段的活化能；
其中，(1)步骤中，首先是选择含PVC混合塑料热裂解的实验试样组成成分，要对实验试样进行干燥、破碎和研磨预处理后，按照比例均匀混合后选取m克放入实验用的坩埚中，实验前对热重分析仪进行预设置；
(2)步骤中，实验反应的环境采用高纯度氮气进行吹扫气，设置吹扫气流量，设置实验温度上升区间T1～T2，选择三组到四组升温速率，分别为β1K/min、β2K/min、β3K/min、β4K/min，为消除环境因素影响，每组实验采取重复三次取平均值的方式，保证实验得到的曲线真实可靠；
在热重法分析含PVC混合塑料热裂解的过程中，转化率为：
Figure FDA0003995188660000011
式中：m0为热分解前塑料样品的质量，mt为反应过程中t时刻的试样质量，mf为反应结束后试样的质量；
(4)步骤中，在含PVC混合塑料热裂解整个过程中，以连续两个热解速率的高峰中间热解速率接近零的时间节点为整个热解过程的第一热解阶段和第二热解阶段的分界节点；
在(3)步骤中，取含PVC混合塑料热裂解过程中的失重率区间为0.1-0.9，在β1K/min升温速率下，失重率α分别取0.1，0.2，0.3…0.9，对应第一个热解阶段的温度分别为TAi(i＝1，2，3...9)，第二个热解阶段的温度分别为Tai(i＝1，2，3...9)；
在β2K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为TBi(i＝1，2，3...9)，第二个热解阶段温度分别为Tbi(i＝1，2，3...9)；
在β3K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段温度分别为Tci(i＝1，2，3...9)，第二个热解阶段温度分别为Tci(i＝1，2，3...9)；
在β4K/min升温速率下，失重率α分别取0.1，0.2，0.3...0.9，对应第一个热解阶段的温度分别为TDi(i＝1，2，3...9)，第二个热解阶段的温度分别为Tdi(i＝1，2，3...9)；(5)步骤中，将表格中的实验数据制成失重率a和温度T的关系图，将实验数据带入到以下三种无模型配合法的分析方法：Flynn-Wall-Ozawa法，Kissinger-Akahira-Sunose法，Starink法；
在第一热解阶段，将失重率α和温度T的实验数据带入到三种无模型配合法的分析方程中，以同一失重率不同升温速率为一组数据，用origin进行曲线回归拟合，分别得到关于ln(β)、ln(β/T2)及ln(β/T1，8)与1/T之间的线性关系曲线；
再以每个失重率下不同升温速率下，构建一个关于x、y线性方程，即：
y＝ax+b
其中：假设y＝ln(β/TU)；x＝1/T；将a＝-BE/R；b＝C，分别对应于公式
Figure FDA0003995188660000021
Figure FDA0003995188660000022
(6)步骤中，结合拟合曲线图像的获得斜率a，通过斜率a和E的关系式求得活化能；由失重率α从0.1～0.9对应共9组拟合直线；
第一个热解阶段的每组拟合直线的活化能为：
E1i＝-aR/B(i＝1，2，3...9)；
将第一个热解阶段9组活化能求取平均值，得到第一个热解阶段活化能E1st：
Figure FDA0003995188660000023
在第二热解阶段，将第二个热解阶段的失重率α和温度T的实验数据带入到三种无模型配合法的分析方程，得到关于ln(β)、ln(β/T2)及ln(β/T1.8)与1/T之间的线性关系曲线，对线性关系曲线进行线性回归拟合得到一组线性方程；即构建一个关于x、y线性方程，对应公式y＝ax+b；
结合拟合曲线图像的获得斜率a，通过斜率a和E的关系式求得活化能；第二个热解阶段的每组拟合直线的活化能为：
E2i＝-aR/B(i＝1,2,3...9)；
将第二个热解阶段9组活化能值求取平均值，得到第二个热解阶段活化能E2nd：
Figure FDA0003995188660000031
CN202110371790.4A 2021-04-07 2021-04-07 一种用于含pvc混合塑料热裂解过程中的活化能分析方法 Active CN113188942B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110371790.4A CN113188942B (zh) 2021-04-07 2021-04-07 一种用于含pvc混合塑料热裂解过程中的活化能分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110371790.4A CN113188942B (zh) 2021-04-07 2021-04-07 一种用于含pvc混合塑料热裂解过程中的活化能分析方法

Publications (2)

Publication Number Publication Date
CN113188942A CN113188942A (zh) 2021-07-30
CN113188942B true CN113188942B (zh) 2023-02-03

Family

ID=76974770

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110371790.4A Active CN113188942B (zh) 2021-04-07 2021-04-07 一种用于含pvc混合塑料热裂解过程中的活化能分析方法

Country Status (1)

Country Link
CN (1) CN113188942B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110793891A (zh) * 2019-10-18 2020-02-14 江苏大学 一种柴油机分级颗粒热解活化能定量评价方法
CN113945480A (zh) * 2021-11-03 2022-01-18 辽宁工程技术大学 一种基于特征温度划分的煤二次氧化活化能分析方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499511A (zh) * 2013-10-16 2014-01-08 南京林业大学 基于多阶段热分析动力学模型预测沥青燃烧过程的方法
CN111366603A (zh) * 2020-04-13 2020-07-03 国网天津市电力公司电力科学研究院 一种基于kissinger法的干式绝缘设备活化能获取方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS6153556A (ja) * 1984-08-22 1986-03-17 Mitsubishi Electric Corp 熱重量分析装置
JP3712837B2 (ja) * 1997-07-02 2005-11-02 株式会社リガク 固体の熱分解反応の解析方法
JP3856173B2 (ja) * 1997-11-11 2006-12-13 株式会社リガク 固体の熱分解反応の解析方法
CN102331437A (zh) * 2011-07-30 2012-01-25 常州天合光能有限公司 预测光伏组件用复合材料的热降解寿命的方法
RU2469299C1 (ru) * 2011-09-02 2012-12-10 Открытое акционерное общество "Корпорация "Росхимзащита" (ОАО "Корпорация "Росхимзащита") Способ определения энергии активации процесса десорбции
JP6096567B2 (ja) * 2012-07-04 2017-03-15 日本電信電話株式会社 昇温脱離分析方法、昇温脱離分析装置およびプログラム
CN102967620A (zh) * 2012-12-11 2013-03-13 常州大学 一种评定高分子材料使用寿命的方法
RU2627552C1 (ru) * 2016-10-05 2017-08-08 Российская Федерация, от имени которой выступает Государственная корпорация по атомной энергии "Росатом" (Госкорпорация "Росатом") Термоаналитический способ определения энергии активации термодеструкции полимерного материала
CN110793891A (zh) * 2019-10-18 2020-02-14 江苏大学 一种柴油机分级颗粒热解活化能定量评价方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103499511A (zh) * 2013-10-16 2014-01-08 南京林业大学 基于多阶段热分析动力学模型预测沥青燃烧过程的方法
CN111366603A (zh) * 2020-04-13 2020-07-03 国网天津市电力公司电力科学研究院 一种基于kissinger法的干式绝缘设备活化能获取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
PP/PVC共混体系的热反应性研究;何小芳等;《材料导报B：研究篇》;20140228;正文1.3节-2.2节内容和图1-3、表1-3 *
聚氯乙烯等塑料废弃物热解特性及动力学研究;韩斌;《中国优秀博硕士论文全文数据库工程科技I辑》;20130515;正文2.2.1节-2.2.4节内容 *

Also Published As

Publication number Publication date
CN113188942A (zh) 2021-07-30

Similar Documents

Publication Publication Date Title
CN113188942B (zh) 一种用于含pvc混合塑料热裂解过程中的活化能分析方法
Egel et al. Premeiotic DNA synthesis in fission yeast
CN105930618A (zh) 一种针对复合材料层合板的混合疲劳可靠性优化方法
Plessis et al. Production of neutral species in Titan’s ionosphere through dissociative recombination of ions
CN111337530A (zh) 一种基于热分析和质谱法分析金属有机骨架材料热解动力学的方法
Chaurasia et al. Modeling & simulation of pyrolysis of biomass: effect of thermal conductivity, reactor temperature and particle size on product concentrations
Katemee et al. Control charts for zero-inflated Poisson models
Newcomb et al. Rate constants for the analysis of alkyl iodide mechanistic probe studies
Babu et al. Optimization of pyrolysis of biomass using differential evolution approach
Bakker ANOSIM
CN113759081B (zh) 新型环保绝缘气体与设备内固体材料相容性的评估方法
CN116451476A (zh) 一种用于塑料热裂解过程中模型拟合法的分析方法
Lech et al. Automatic control of the activated sludge process—I. Development of a simplified dynamic model
CN114707784A (zh) 一种基于组合赋权法的权重确定方法及系统
Alpana et al. Automated char classification using image analysis and artificial intelligence
Hamann Analyses of Wolf-Rayet Spectra
CN114638144A (zh) 一种三元正极材料前驱体分解反应模拟方法及装置
Baltay et al. Search for exotics produced opposite a forward Λ in 15 GeV/c π+ p interactions
CN109358579B (zh) 等离子体裂解煤制乙炔生产过程的全流程协调控制系统及方法
Kunioka et al. BIOBASED CONTENTS OF MODEL RUBBER COMPOUNDS BASED ON ISO 19984 SERIES
CN116453631A (zh) 一种厌氧水解酸化过程中预测和调节挥发性脂肪酸的方法
Cioslowski et al. Symmetry handling in calculations of properties of atoms in molecules
Pohlmann Controlling the Type I Error Rate in Stepwise Regression Analysis.
Hamann et al. Iteration with approximate lambda operators, and its application to the expanding atmospheres of WR stars
Jing et al. A hybrid genetic algorithm for Bin Packing problem based on item sequencing

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