一种用于塑料热裂解过程中模型拟合法的分析方法

1.本发明属于废旧塑料热裂解绿色资源化再利用领域，尤其涉及一种用于塑料热裂解过程中模型拟合法的分析方法。


背景技术：

2.中国拥有着很高的塑料生产、消耗和回收量并且位居世界前列，当下使用广泛的塑料回收手段主要有焚烧、填埋、堆肥化和回收再生等手段。这其中多数地区主要以填埋和焚烧为主，使用回收再生等可循环的方法只占少数。使用焚烧的方法处理塑料会产生大量有害气体从而对人类的健康造成危害，而将其填埋也会对土壤环境造成破坏因而对环境保护产生影响。目前在废弃塑料的处理问题中，如何提高其回收再利用的效率并且做到绿色无污染，成为了一个促进循环经济及节能减排的关键问题。在塑料的回收再生方法中，热裂解是一种效率高，环保性强的方法，而且它具有将废旧塑料转化成不可再生资源的能力，如燃料油、燃料气等，如果广泛普及，还可解决能源短缺等问题，是一种具有经济性和环保性的高效塑料绿色回收方式，实现了对废弃塑料的减量化、资源化和无害化处理，是一种拥有巨大潜力的回收方法。
3.在动力学分析中活化能、机理函数和指前因子是热分析动力学研究的三要素。在完成对动力学三要素的确定工作后即可对塑料裂解过程进行模型函数的假设，在处理实验数据时往往会加入动力学模型与之相配合来分析整个热裂解过程，这种方法被称为模型配合法。而在进行模型的配合过程中往往会出现多个模型都可以与同一个反应过程曲线相配合的情况。
4.做出此项研究的目的在于在进行热裂解实验研究废弃塑料动力学参数的同时，找到最合适的动力学反应机理与之相匹配。在使用ozawa-flynn-wall(o-f-w)、kissinger-akahira-sunose(k-a-s)、starink三种无模型法对废弃塑料热裂解的活化能进行计算后，在使用coats-redfern模型拟合方法对实验过程的反应动力学机理进行了预测时，能够找到与实际过程高度匹配的模型。模型的选择对于机理函数和指前因子的计算具有重要意义。
5.论文thermal degradation of typical plastics under high heating rate conditions 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等模型拟合方法对反应动力学机理进行了预测。将低密度聚乙烯的热解反应机理归纳为收缩球模型，将聚乙烯的热解反应归纳为收缩圆柱模型，而聚氯乙
烯的热解反应机理归纳为第一阶段为二维成核，第二阶段为三维扩散(简化)。但本论文没有揭示含pvc混合塑料热裂解过程活化能与转化率的关系。并且现有技术中计算指前因子的过程中因为机理模型选择不准确会导致指前因子在计算过程中产生的误差较大的问题，在模型拟合法对机理模型的计算中因为机理函数的选择的不准确从而对动力学参数产生的影响。


技术实现要素：

6.有鉴于此，本发明针对上述问题，提出一种用于塑料热裂解过程中模型拟合法的分析方法。
7.为了实现上述目的，本发明采用如下技术方案实现：
8.一种用于塑料热裂解过程中模型拟合法的分析方法，用所选塑料在热裂解过程中的数据进行计算，实现对其活化能的计算以及对其机理函数的选择，主要包含以下步骤：
9.(1)对所选废弃塑料试样进行热重实验分析，得到该试样在多种升温速率下的tg曲线和dtg曲线及其相关实验数据；
10.(2)通过热重实验中得到的热裂解的tg和dtg曲线，使用无模型法根据实验数据计算出所选塑料试样在热裂解过程中的平均活化能值；
11.(3)将所选的废弃塑料试样在热重实验中得到的热裂解的tg和dtg曲线进行数值模拟分析，将分析后得到的实验数据代入到c-r方程中，计算出各模型所对应的活化能的值，将这些值与无模型法计算得来的活化能值进行比较；
12.(4)找到与无模型法计算得来的活化能值最为相近的两个模型的活化能值，确定出平均活化能值的最小位置区间；
13.(5)找到此区间的中点，在中点处建立新的机理模型，计算新建立的模型对应的活化能值，由此建立出新的最小位置区间并重复此步骤；
14.(6)规定差值范围最终此区间被不断缩小，达到所要求的精确度后，或设定范围区间，对应的模型即为理想模型，即为所求机理模型的反应级数进行精确优化。
15.优选的，模型拟合法中可能与真实实验过程存在差距，选用无模型配合法直接进行对活化能的计算过程；
16.无模型配合法又被叫做等转化率法，其在实验过程中采用多种升温速率进行试验的方法，即设置多条升温速率曲线并对其进行热分析的动力学处理方式。
17.优选的，在(1)中，第一步首先是要明确所选的废弃塑料试样的组成成分，在对其进行充分的烘干处理后，对塑料试样进行破碎和研磨的预处理，依照实验比例选取m克研磨好的试样放入坩埚中做准备；
18.做好以上准备工作后，还要对实验所要使用的热裂解炉设备进行预设置处理。
19.优选的，在(2)中，对于塑料式样的热裂解环境，选用使用高纯度氮气对其进行吹扫气的方法，并且需要设置合适的吹扫气的流量避免对实验结果造成影响，在试验前设置热裂解温度的上升区间t
1-t2，选择四组升温速率进行试验，分别为β1k/min、β2k/min、β3k/min、β4k/min。为降低由于环境因素所造成的实验误差，每组实验都要重复三到四次并且取其平均值进行计算，以此降低由于环境因素所造成的误差。
20.优选的，1)在(1)中，通过对所选试样进行热重实验分析的方法得到的四种升温速
率下的tg和dtg曲线，通过对tg和dtg曲线的分析得到转化率α和温度t的实验数据及关系图。
21.在(2)中选择β1、β2、β3、β4四种升温速率，通过在热重实验中得到的热裂解的tg和dtg曲线分别找到在四种不同的升温速率下，失重率α从0.1-0.9所对应的不同温度。
22.2)在得到了四种升温速率在9组不同的失重率下所对应的36个温度后，将这36个温度带到三种无模型法方程中对其活化能进行计算：
23.无模型法主要有flynn-wall-ozawa(f-w-o)法、kissinger-akahira-sunose(k-a-s)法和starink法三种方法。三种无模型法可以总结为一个方程：
[0024][0025]
将四组升温速率对应的36个温度分别带入到上述三种无模型法所对应的方程中，以失重率相等但升温速率不相等的条件下为一组数据，分别将1/t作为x值，三种方法对应的u值分别为0、2和1.8，故将ln(β)、ln(β/t2)、ln(β/t
1.8
)作为y值带入origin并对其进行线性拟合，三种方法可分别得到9条不同的线性方程。
[0026]
在进行线性拟合后可以获得斜率a的值，通过方程中活化能e与斜率a的关系可求出e的值。9条线性方程共可以求得9组e的值，将这9组活化能值取平均值即可将该值作为无模型法所求得的活化能值e。
[0027]
3)在模型拟合法中主要是应用c-r方程来对不同模型对应的活化能值进行计算，c-r方程的具体形式为：
[0028][0029]
对塑料的活化能和热裂解温度区间来说，ln[ar/βe(1-2rt/e)]可视为固定常数，因此式(2)可看作线性方程处理，将ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，使用origin进行线性拟合，通过拟合后所得到的斜率的值与活化能时间的关系可计算出不同模型所对应的活化能的值，将无模型法计算出的活化能平均值与各模型所对应的活化能值比较，即可得到与实验较为相近得模型。常用机理函数列在表1中：
[0030]
表1常用的固体热裂解动力学机理函数
[0031][0032]
4)在使用无模型法计算出平均活化能值e后，在拟合法中可能会面临平均活化能值与所选择机理函数模型的活化能值相差较大，从而导致在拟合法中所选择的模型与实际过程模型差距较大的问题，可以在使用拟合法选择机理函数模型时引入二分法，从而使所选择的机理函数模型可以无限接近于实际过程模型。
[0033]
5)当选定某种塑料后，对其进行热裂解试验，随着转化率的提高，三种方法得到混
0.5x5y50.6x6y60.7x7y70.8x8y80.9x9y9[0042]
计算出α从0.1-0.9所对应的九个1/t(
×
103k-1
)和ln[g(α)/t2]后，再以ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，对此线性方程进行拟合，可得到该线性方程的斜率，由c-r方程可以计算出a
(a+b)/2
所对应的活化能e
(a+b)/2
的值。
[0043]
8)此时用无模型法计算得出的平均活化能e若大于e
(a+b)/2
，则e位于e
(a+b)/2
与ea之间，平均活化能所在位置区间缩减了一半，然后再找到(a+b)/2与a的中点即(3a+b)/4，计算出a
(3a+b)/4
所对应的活化能e
(3a+b)/4
的值，则平均活化能所在位置区间将再次缩减一半，此时对e
(3a+b)/4
进行判定其是否达到指定精度，这里规定p为所选模型对应的活化能与平均活化能e之间的最大差值，所选模型对应活化能的精确度取决于p的精确度；
[0044]
如果|e
(3a+b)/4-e|≤p则模型a
(3a+b)/4
达到指定精度，此时便可输出a
(3a+b)/4
作为所选择的模型，若未达到指定精度即|e
(3a+b)/4-e|＞p则继续重复上述步骤到达到指定精度为止；
[0045]
若e等于e
(a+b)/2
，则a
(a+b)/2
即为精准模型；若e小于e
(a+b)/2
，则e位于eb与e
(a+b)/2
之间，此时再找到b与(a+b)/2的中点即(3a+b)/4，计算出a
(a+3b)/4
所对应的活化能e
(a+3b)/4
的值，此时对e
(a+3b)/4
进行判定其是否达到指定精度，若|e
(a+3b)/4-e|≤p则模型a
(a+3b)/4
达到指定精度，此时输出a
(a+3b)/4
作为所选模型，若未达到指定精度即|e
(a+3b)/4-e|＞p则继续重复上述步骤继续循环至到达到指定精度为止。如此循环往复，使平均活化能所在的位置区间的两端点无限接近与平均活化能的值，此时所取的模型即为最贴合真实实验过程的模型，其所对应的机理函数g(α)也更为准确。
[0046]
本发明的技术方案至少具有如下优点和有益效果：
[0047]
1.本发明结合了二分法的方式，在模型拟合法的计算过程中联合无模型法得到的平均活化能值，提高了模型拟合法所求得的机理模型精度，在原有计算指前因子的过程中充分解决了机理模型选择不准确的问题从而导致指前因子在计算过程中产生的误差较大的问题。
[0048]
2.本发明在模型拟合法对机理模型的计算中使得无模型法计算获得的平均活化能和与之相邻的两个模型对应的活化能的区间不断缩小，无模型法计算得来的平均活化能值与所在的位置区间两端点的差值不断缩小，从而使得无模型法计算得来的平均活化能值与c-r法计算得来的活化能值的契合度能够不断提高。
[0049]
3.本发明避免了在模型拟合法对机理模型的计算中因为机理函数的选择的不准确从而对动力学参数产生的影响，使所选择的模型无限接近实际过程模型。
[0050]
4.本发明将二分法与模型拟合法进行了结合，在原有的模型拟合法上进行了升级，使之能够选择出与实际过程模型高度契合的模型，该模型对应的活化能值可以无限接近于无模型法计算出的活化能值。
[0051]
5.本发明使用无模型法获得平均活化能值后，再通过模型拟合法计算动力学模型参数，根据模型拟合法获得的不同模型的活化能值，选出活化能值最接近无模型法计算出的平均活化能值的两个模型，这两组活化能构成的区间即为平均活化能所在的最小位置区
间，使用二分法不断缩小平均活化能所在的最小位置区间，从而使c-r法计算出的活化能值与使用无模型法获得的平均活化能值达到一个极高的贴合度，进而找到更精确的机理函数g(α)，使所选择的模型无限接近实际过程模型。
附图说明
[0052]
图1为本发明二分法用于模型拟合法流程图；
[0053]
图2本发明二分法计算f-w-o法数据流程图；
具体实施方式
[0054]
下面结合附图和具体实施例对本发明进行说明。
[0055]
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合附图，对本发明实施例中的技术方案进行清楚、完整的描述。显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0056]
如图1-2所示，作为本发明一个优选的实施例，这里以pc的热裂解实验为例。经实验，pc在热裂解过程的活化能总结在表5中：
[0057]
表5pc混合物热裂解过程的活化能
[0058][0059]
可以看出随着转化率的提高，三种方法得到混合物的活化能逐渐增大，其中k-a-s、starink和f-w-o法得到的活化能范围分别为87.69-178.14kj/mol、88.25-178.49kj/mol和91.76-178.91kj/mol，对应的活化能平均值分别为114.24kj/mol、114.76kj/mol和117.58kj/mol。
[0060]
由于pc热裂解机理均属于随机成核随后生长，将转化率α从0.1-0.9及其所对应的不同温度带入到c-r方程并进行线性拟合后，计算出不同模型对应的活化能表6所示：
[0061]
表6pc热裂解实验各模型所对应活化能
[0062][0063][0064]
这里以f-w-o法为例，可以看出，用f-w-o法计算得出的平均活化能117.58kj/mol处在a1和a
3/4
之间。此时，可继续计算a
7/8
所对应的e
7/8
的值。则由表1可以得到g(α)的机理函数为[-ln(1-α)]
8/7
，这里先对1/t(
×
103k-1
)进行计算，然后再对ln[g(α)/t2]进行计算，可以得到表7中的数据：
[0065]
表7机理函数取值
[0066]
α1/t(
×
103k-1
)ln[g(α)/t2]0.11.9401691-15.06180870.21.8815242-14.26556340.31.8502038-13.76312510.41.8230012-13.38223030.51.8032151-13.05524000.61.7869716-12.75437630.71.7646612-12.46744820.81.7309034-12.17435510.91.6812007-11.8233142
[0067]
计算出α从0.1-0.9所对应的九个ln[g(α)/t2]后，再以ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，对此线性方程进行拟合，可得到该线性方程的斜率为-13.1271，假设p的精确度为0.0001，由c-r方程可知，a
7/8
所对应的活化能值为：
[0068]e7/8
＝13.1271
×
8.314＝109.1387kj/mol
[0069]
此时用f-w-o法计算得出的平均活化能117.58kj/mol位于a
3/4
与a
7/8
之间。然后再找到3/4与7/8的中点即13/16，计算出模型a
13/16
所对应的活化能e
13/16
的值，判定e
13/16
是否达到指定精度，若|e
13/16-e|≤p，则模型a
13/16
达到指定精度，此时便可输出a
13/16
作为所选模型，若未达到指定精度即|e
13/16-e|＞p，则重复上述步骤继续计算至达到指定精度为止。
[0070]
以上实例若采用黄金分割法分割区间，则用f-w-o法计算得出的平均活化能117.58kj/mol处在a1和a
3/4
之间时，采用黄金分割比0.618，此时，可继续计算a
0.9045
所对应的e
0.9045
的值。则由表1可以得到g(α)的机理函数为[-ln(1-α)]
1/0.9045
，这里先对1/t(
×
103k-1
)进行计算，然后再对ln[g(α)/t2]进行计算，可以得到表8中的数据：
[0071]
表8黄金分割法分割区间机理函数取值
[0072]
α1/t(
×
103k-1
)ln[g(α)/t2]0.11.9401691-14.97792860.21.8815242-14.20965470.31.8502038-13.72469830.41.8230012-13.35719240.51.8032151-13.04157860.61.7869716-12.75111770.71.7646612-12.47436730.81.7309034-12.19209320.91.6812007-11.8544019
[0073]
计算出α从0.1-0.9所对应的九个ln[g(α)/t2]后，再以ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，对此线性方程进行拟合，可得到该线性方程的斜率为-12.66289，假设p的精确度为0.0001，由c-r方程可知，a
0.9045
所对应的活化能值为：
[0074]e0.9045
＝12.66289
×
8.314＝105.2793kj/mol
[0075]
此时用f-w-o法计算得出的平均活化能117.58kj/mol位于a
3/4
与a
0.9045
之间。然后再找到3/4与0.9045的黄金分割点即0.845481，计算出模型a
0.845481
所对应的活化能e
0.845481
的值，判定e
0.845481
是否达到指定精度，若|e
0.845481-e|≤p，则模型a
0.845481
达到指定精度，此时便可输出a
0.845481
作为所选模型，若未达到指定精度即|e
0.845481-e|＞p，则重复上述步骤继续循环计算至达到指定精度为止。
[0076]
以上实施方案仅用于说明而非限制本发明的技术方案。不脱离本发明精神的任何修改或局部替换，均应涵盖在本发明的权利要求范围当中。

技术特征：
1.一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于，用所选塑料在热裂解过程中的数据进行计算，实现对其活化能的计算以及对其机理函数的选择，该方法包含以下步骤：(1)对所选废弃塑料试样进行热重实验分析，得到该试样在多种升温速率下的tg曲线和dtg曲线及其相关实验数据；(2)通过热重实验中得到的热裂解的tg和dtg曲线，使用无模型法根据实验数据计算出所选塑料试样在热裂解过程中的平均活化能值；(3)将所选的废弃塑料试样在热重实验中得到的热裂解的tg和dtg曲线进行数值模拟分析，将分析后得到的实验数据代入到c-r方程中，计算出各模型所对应的活化能的值，将这些值与无模型法计算得来的活化能值进行比较；(4)找到与无模型法计算得来的活化能值最为相近的两个模型的活化能值，确定出平均活化能值的最小位置区间；(5)找到此区间的中点，在中点处建立新的机理模型，计算新建立的模型对应的活化能值，建立出新的最小位置区间并重复此步骤；(6)规定差值范围最终此区间被不断缩小，达到所要求的精确度后或设定范围区间，对应的模型即为理想模型。2.根据权利要求1所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：在(1)中，明确所选的废弃塑料试样的组成成分，在对其进行充分的烘干处理后，对塑料试样进行破碎和研磨的预处理，依照实验比例选取m克研磨好的试样放入坩埚中做准备；对实验所要使用的热裂解炉设备进行预设置处理；在(2)中，对于塑料式样的热裂解环境，选用使用高纯度氮气对其进行吹扫气的方法，，在试验前设置热裂解温度的上升区间t
1-t2，选择四组升温速率进行试验，分别为β1k/min、β2k/min、β3k/min、β4k/min。3.根据权利要求2所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：在(1)中，通过对所选试样进行热重实验分析的方法得到的四种升温速率下的tg和dtg曲线，通过对tg和dtg曲线的分析得到转化率α和温度t关系图；在(2)中选择β1、β2、β3、β4四种升温速率，通过在热重实验中得到的热裂解的tg和dtg曲线，分别找到在四种不同的升温速率下，失重率α从0.1-0.9所对应的不同温度。4.根据权利要求3所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：在得到四种升温速率在9组不同的失重率下所对应的36个温度后，将这36个温度带到三种无模型法方程中对其活化能进行计算；三种无模型法方程用通式表示为，将四组升温速率对应的36个温度分别带入到所述三种无模型法所对应的方程中，以失重率相等但升温速率不相等的条件下为一组数据，分别将1/t作为x值，三种方法对应的u值分别为0、2和1.8；
将ln(β)、ln(β/t2)、ln(β/t
1.8
)作为y值带入origin并对其进行线性拟合，三种方法分别得到9条不同的线性方程；在进行线性拟合后获得斜率a的值，通过方程中活化能e与斜率a的关系求出e的值；9条线性方程共求得9组e的值，将9组活化能值取平均值即将该值作为无模型法所求得的活化能值e。5.根据权利要求4所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：在模型拟合法中应用c-r方程来对不同模型对应的活化能值进行计算，具体形式为：对于塑料的活化能和热裂解温度区间，ln[ar/βe(1-2rt/e)]为固定常数，将ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，使用origin进行线性拟合，通过拟合后所得到的斜率的值与活化能时间的关系计算出不同模型所对应的活化能的值，将无模型法计算出的活化能平均值与各模型所对应的活化能值比较，得到与实验相近的模型。6.根据权利要求5所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：使用拟合法选择机理函数模型时引入二分法，使所选择的机理函数模型无限接近于实际过程模型。7.根据权利要求6所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：当选定某种塑料后，对其进行热裂解试验；热裂解机理均属于随机成核随后生长，将转化率α从0.1-0.9带入到各模型的机理函数中计算出不同转化率下各模型所对应的g(α)，每个模型将得到9组不同的g(α)值，再将这9组g(α)值以及每组所对应的温度t带入ln[g(α)/t2]进行计算，得到9组不同的ln[g(α)/t2]值，每个模型以9组不同的ln[g(α)/t2]作为y，9组不同的1/t(
×
103k-1
)作为x进行线性拟合；每个模型均得到一条拟合出来的直线，取该直线的斜率，运用式c-r方程中斜率与活化能e的关系计算出每个模型所对应的活化能值。8.根据权利要求7所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：设通过无模型法计算得出的平均活化能e处在e
a
和e
b
之间；继续计算a
(a+b)/2
所对应的e
(a+b)/2
的值；先对1/t(
×
103k-1
)进行计算，然后再对ln[g(α)/t2]进行计算，计算出α从0.1-0.9所对应的九个1/t(
×
103k-1
)和ln[g(α)/t2]后，再以ln[g(α)/t2]作为y，1/t(
×
103k-1
)作为x，对此线性方程进行拟合，得到该线性方程的斜率。9.根据权利要求8所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：用无模型法计算得出的平均活化能e若大于e
(a+b)/2
，则e位于e
(a+b)/2
与e
a
之间，平均活化能所在位置区间缩减一半，找到(a+b)/2与a的中点(3a+b)/4；计算出a
(3a+b)/4
所对应的活化能e
(3a+b)/4
的值，平均活化能所在位置区间将再次缩减一半，对e
(3a+b)/4
进行判定其是否达到指定精度，规定p为所选模型对应的活化能与平均活化能
e之间的最大差值，所选模型对应活化能的精确度取决于p的精确度；若|e
(3a+b)/4-e|≤p则模型a
(3a+b)/4
达到指定精度，此时输出a
(3a+b)/4
作为所选择的模型，若未达到指定精度即|e
(3a+b)/4-e|>p，则继续重复上述步骤到达到指定精度。10.根据权利要求8所述的一种用于塑料热裂解过程中模型拟合法的分析方法，其特征在于：若e等于e
(a+b)/2
，则a
(a+b)/2
即为精准模型；若e小于e
(a+b)/2
，则e位于e
b
与e
(a+b)/2
之间，再找到b与(a+b)/2的中点(a+3b)/4，计算出a
(a+3b)/4
所对应的活化能e
(a+3b)/4
的值，对e
(a+3b)/4
进行判定其是否达到指定精度；若|e
(a+3b)/4-e|≤p则模型a
(a+3b)/4
达到指定精度，输出a
(a+3b)/4
作为所选模型；若未达到指定精度即|e
(a+3b)/4-e|>p，则继续重复上述步骤继续循环至到达到指定精度。

技术总结
本发明公开一种用于塑料热裂解过程中模型拟合法的分析方法，步骤包括：1)对废弃塑料试样进行热重实验分析，得到多种升温速率下的TG/DTG曲线及数据；2)使用无模型法计算出塑料试样在热裂解过程中的平均活化能值；3)TG/DTG曲线进行数值模拟分析得到的数据代入到C-R方程中，得到各模型所对应的活化能的值，与无模型法的该值比较；4)找到与无模型法的活化能值最为相近的两个模型的活化能值，确定出平均活化能值的最小位置区间；5)找到此区间的中点并在此新的机理模型，计算对应的活化能值，建立出新的最小位置区间，并重复此步骤；6)区间缩小达到所要求的精确度后对应的模型即为理想模型。本发明有效解决了机理模型选择不准确，指前因子误差较大的问题。指前因子误差较大的问题。指前因子误差较大的问题。


技术研发人员：尹凤福 李亚峰 刘硕 张茗昱 范方硕 田晓龙 汪传生
受保护的技术使用者：青岛科技大学
技术研发日：2023.04.18
技术公布日：2023/7/21