CN111257265A - 一种钢桥面铺装用沥青针入度指标的快速检测方法 - Google Patents
一种钢桥面铺装用沥青针入度指标的快速检测方法 Download PDFInfo
- Publication number
- CN111257265A CN111257265A CN202010097162.7A CN202010097162A CN111257265A CN 111257265 A CN111257265 A CN 111257265A CN 202010097162 A CN202010097162 A CN 202010097162A CN 111257265 A CN111257265 A CN 111257265A
- Authority
- CN
- China
- Prior art keywords
- asphalt
- penetration
- sample
- infrared spectrum
- model view
- 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.)
- Pending
Links
- 239000010426 asphalt Substances 0.000 title claims abstract description 83
- 230000035515 penetration Effects 0.000 title claims abstract description 46
- 229910000831 Steel Inorganic materials 0.000 title claims abstract description 16
- 239000010959 steel Substances 0.000 title claims abstract description 16
- 238000001514 detection method Methods 0.000 title claims description 16
- 238000012937 correction Methods 0.000 claims abstract description 38
- 238000002329 infrared spectrum Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 25
- 230000002159 abnormal effect Effects 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 18
- 238000013499 data model Methods 0.000 claims abstract description 18
- 238000004458 analytical method Methods 0.000 claims abstract description 15
- 238000001228 spectrum Methods 0.000 claims abstract description 10
- 238000010521 absorption reaction Methods 0.000 claims abstract description 5
- 238000005102 attenuated total reflection Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims abstract description 4
- 238000012216 screening Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 33
- 238000012795 verification Methods 0.000 claims description 29
- 230000003595 spectral effect Effects 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 9
- 239000013598 vector Substances 0.000 claims description 9
- 238000011156 evaluation Methods 0.000 claims description 7
- 238000009826 distribution Methods 0.000 claims description 6
- 238000009499 grossing Methods 0.000 claims description 6
- 230000009467 reduction Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 4
- 230000000737 periodic effect Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 230000003252 repetitive effect Effects 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 238000004566 IR spectroscopy Methods 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 5
- 239000003208 petroleum Substances 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 230000008030 elimination Effects 0.000 description 4
- 238000003379 elimination reaction Methods 0.000 description 4
- 238000010200 validation analysis Methods 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 101100202447 Drosophila melanogaster sav gene Proteins 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000012417 linear regression Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003908 quality control method Methods 0.000 description 2
- 238000011897 real-time detection Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 241000238097 Callinectes sapidus Species 0.000 description 1
- 101000596041 Homo sapiens Plastin-1 Proteins 0.000 description 1
- 102100035181 Plastin-1 Human genes 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000003381 stabilizer Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N11/10—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material
- G01N11/12—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by moving a body within the material by measuring rising or falling speed of the body; by measuring penetration of wedged gauges
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3563—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (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 Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种钢桥面铺装用沥青针入度指标的快速检测方法,包括步骤:S1、收集多个沥青样品,分别制成沥青分析试样;S2、采用衰减全反射红外光谱采集方式对沥青分析试样进行采样获取红外光谱;S3、根据红外光谱的红外特征吸收峰,筛选对应的最优波长作为样本;S4、根据样本建立定量数据模型,将校正集光谱数据导入数据模型进行计算,生成原始模型视图;S5、判断是否存在异常值,若是,则剔除异常值后重新计算新的数据模型并生成新的模型视图后进入步骤S6;若否,则将原始模型视图作为新的模型视图并进入步骤S6;S6、获取待测沥青的红外光谱,并导入新的模型视图中,计算待测沥青的针入度。红外光谱结合化学计量学的分析方法能够更快速准确的检测沥青的针入度。
Description
技术领域
本发明涉及石油化工技术领域,尤其涉及一种钢桥面铺装用沥青针入度指标的快速检测方法。
背景技术
钢桥面铺装是历来桥梁建设的难点与重点。对于钢桥面铺装来说,其受力和变形远较道路路面或机场道面复杂,因而对其强度、变形特性、温度稳定性、疲劳耐久性等均有更高要求。同时由于铺装所处的特殊位置,在使用性能上存在重量轻、粘结性高、不透水等特殊要求。作为桥梁行车系的重要组成部分,桥面铺装的好坏直接影响到行车的安全性、舒适性、桥梁耐久性及投资效益和社会效益。
沥青产品作为钢桥面铺装建设中的重要原材料,其质量对桥面铺装起到决定性的影响。所以,对钢桥面铺装用沥青质量的评价是多年以来各国研究的重点。
据调查我国沥青市场仍存在不合格的沥青产品,假冒伪劣现象时有发生。例如沥青产品个别指标不合格或者指标合格但不同批次产品的波动范围较大,这就导致路面质量的不均匀性。识别冒牌沥青,传统的简单物性检测方法不仅耗时耗力,且易受稳定剂等添加剂或生产工艺的影响,试验结果容易失真。由于沥青化学性质的复杂性,沥青规范开发了针对沥青物理特性的试验,诸如针入度,针入度,延度等,这些物理特性试验在标准测试温度下进行,测试结果被用来确定材料是否满足规范的标准,三大指标是评判沥青性能的传统检测方法,其简单的操作是其他方法不可取代的,至今也是公路行业检测沥青的主要方法。
现有的技术中,沥青传统检测方法耗时较长,只能反映特定温度、特定荷载作用条件下的情况,对于其他条件下沥青的性能差异很难体现出来,检测方法不够全面;SHRP检测指标能够更全面的反映沥青的路用性能,但费时费力,且设备昂贵,不能够及时有效的判定沥青质量。
发明内容
本发明的目的在于提供一种钢桥面铺装用沥青针入度指标的快速检测方法,旨在解决现有技术中传统检测方法耗时较长,只能反映特定温度、特定荷载作用条件下的情况;而SHRP检测费时费力,且设备昂贵,不能够及时有效的判定沥青质量的问题。
本发明实施例提供了一种钢桥面铺装用沥青针入度指标的快速检测方法,包括步骤:S1、收集多个沥青样品,分别制成沥青分析试样;S2、采用衰减全反射红外光谱采集方式对所述沥青分析试样进行采样获取红外光谱;S3、根据所述红外光谱的红外特征吸收峰,筛选对应的最优波长作为样本;S4、根据所述样本建立定量数据模型,将校正集光谱数据导入数据模型进行计算,生成原始模型视图;S5、判断是否存在异常值,若是,则剔除异常值后重新计算新的数据模型并生成新的模型视图后进入步骤S6;若否,则将原始模型视图作为新的模型视图并进入步骤S6;S6、获取待测沥青的红外光谱,并导入新的模型视图中,计算待测沥青的针入度。
在本发明一实施例中,所述步骤S5包括步骤:选取多个建模区间;针对多个建模区间分别建立多个校正模型;针对每个校正模型获取对应校正相关系数及校正标准偏差;对照预设关系剔除异常值后重新优化计算对应的多个校正模型;分别采用验证集验证多个校正模型的验证结果;选取最优建模区间并根据所述样本建立定量数据模型。
在本发明一实施例中,所述步骤“对照预设关系剔除异常值”包括步骤:根据三维空间分布值和线性相关性并参考光谱杠杆值确定异常值后将异常值剔除。
在本发明一实施例中,所述步骤S3包括步骤:对所述样本包含的数据进行平滑、降噪处理;所述平滑、降噪处理采用Savitzky-Golay、百分位滤波器、FFT滤波器、LOVESS、LOESS或二项式法;边界条件选择对称、重复、周期性或外推法。
在本发明一实施例中,所述步骤“剔除异常值后重新计算新的数据模型并生成新的模型视图”包括步骤:生成新的模型视图的相关系数并采用验证集对新的模型视图进行验证。
在本发明一实施例中,所述步骤S6包括步骤:S61、根据待测沥青的红外光谱在新的模型视图中的位置对待测沥青的主成分进行分析,获取全部主成分数及每种主成分的贡献率;S62、根据贡献率选出特征主成分,并对特征主成分进行综合评价,根据综合评价获取待测沥青的针入度。
在本发明一实施例中,所述步骤S62中包括步骤:S621、标准化采集p维随机向量和n个样品构造样本阵,并对样本阵进行标准化变换获取标准化矩阵;S622、根据所述标准化矩阵计算相关系数矩阵;S623、根据所述相关系数矩阵计算预设特征方程的特征根和特征向量,并根据所述特征根确定主成分;S624、根据所述标准化矩阵和所述特征向量获取主成分矩阵;S625、根据所述主成分矩阵对主成分进行综合评价。
在本发明一实施例中,所述步骤S621中:随机向量x=(x1,x2,...,xp)T;样品xij=(xi1,xi2,...,xip)T;标准化矩阵的计算公式为:其中,i=1,2,…,n;j=1,2,…,p;所述步骤S622中:相关系数矩阵的计算公式为:其中,p;所述步骤S623中:所述相关系数矩阵R的特征方程为:|R-λIp|=0;所述步骤S624中:主成分矩阵的计算公式为:其中,j=1,2,…,m。
基于上述技术方案,与现有技术相比,本发明实施例提出的一种钢桥面铺装用沥青针入度指标的快速检测方法,实现了道路石油沥青针入度指标的快速准确检测,为沥青生产质量控制提供了强大的支持,能够实现沥青在线实时检测;通过红外光谱识别道路石油沥青的来源,便于在实际工程建设中进行监控。红外光谱结合化学计量学的分析方法能够更准确的检测沥青的针入度;引入化学计量学方法作为红外数据库建立的新的理念,提出了红外光谱结合化学计量学的分析方法,具有广阔的应用前景和优势。
附图说明
图1显示为本发明中沥青针入度的快速检测方法的流程方框示意图。
图2显示为本发明中校正集光谱计算视图。
图3显示为3138-2797cm-1作为针入度建模区间建立的校正模型视图。
图4显示为3138-2797cm-1作为建模区间验证集的预测结果图。
图5显示为1650-1100cm-1作为针入度建模区间建立的校正模型视图。
图6显示为1650-1100cm-1作为建模区间验证集的预测结果图。
图7显示为965-721cm-1作为针入度建模区间建立的校正模型视图。
图8显示为965-721cm-1作为建模区间验证集的预测结果图。
图9显示为同时选取3138-2797cm-1、1650-1100cm-1和965-721cm-1作为建模区间建立的校正模型视图。
图10显示为同时选取3138-2797cm-1、1650-1100cm-1和965-721cm-1作为建模区间验证集的预测结果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
需要说明的是,当元件被称为“固定于”或“设置于”另一个元件,它可以直接在另一个元件上或可能同时存在居中元件。当一个元件被称为是“连接于”另一个元件,它可以是直接连接到另一个元件或者可能同时存在居中元件。
另外,还需要说明的是,本发明实施例中的左、右、上、下等方位用语,仅是互为相对概念或是以产品的正常使用状态为参考的,而不应该认为是具有限制性的。以下结合具体实施例对本发明的实现进行详细的描述。
如图1所示,图1显示为本发明中沥青针入度的快速检测方法的流程方框示意图。本发明提供了一种钢桥面铺装用沥青针入度指标的快速检测方法,包括步骤:
S1、收集多个沥青样品,分别制成沥青分析试样;
S2、采用衰减全反射红外光谱采集方式对所述沥青分析试样进行采样获取红外光谱;
S3、根据所述红外光谱的红外特征吸收峰,筛选对应的最优波长作为样本;通常,所述步骤S3包括步骤:对所述样本包含的数据进行平滑、降噪处理;所述平滑、降噪处理采用Savitzky-Golay、百分位滤波器、FFT滤波器、LOVESS、LOESS或二项式法;边界条件选择对称、重复、周期性或外推法。红外光谱技术操作简单、分析测试时间短,重复性、再现性好,可显著提高沥青的分析效率。
S4、根据所述样本建立定量数据模型,将校正集光谱数据导入数据模型进行计算,生成原始模型视图;通常,依据道路石油沥青的红外特征吸收峰,采用CARS筛选最优波长建立PLS定量模型。
S5、判断是否存在异常值,若是,则剔除异常值后重新计算新的数据模型并生成新的模型视图后进入步骤S6;若否,则将原始模型视图作为新的模型视图并进入步骤S6;在一实施例中,所述步骤“剔除异常值后重新计算新的数据模型并生成新的模型视图”包括步骤:生成新的模型视图的相关系数并采用验证集对新的模型视图进行验证。
S6、获取待测沥青的红外光谱,并导入新的模型视图中,计算待测沥青的针入度。进一步地,所述步骤S6包括步骤:S61、根据待测沥青的红外光谱在新的模型视图中的位置对待测沥青的主成分进行分析,获取全部主成分数及每种主成分的贡献率;S62、根据贡献率选出特征主成分,并对特征主成分进行综合评价,根据综合评价获取待测沥青的针入度。优选地,所述步骤S62中包括步骤:S621、标准化采集p维随机向量和n个样品构造样本阵,并对样本阵进行标准化变换获取标准化矩阵;S622、根据所述标准化矩阵计算相关系数矩阵;S623、根据所述相关系数矩阵计算预设特征方程的特征根和特征向量,并根据所述特征根确定主成分;S624、根据所述标准化矩阵和所述特征向量获取主成分矩阵;S625、根据所述主成分矩阵对主成分进行综合评价。其中,所述步骤S621中:随机向量x=(x1,x2,...,xp)T;样品xij=(xi1,xi2,...,xip)T;标准化矩阵的计算公式为:其中,i=1,2,…,n;j=1,2,…,p; 所述步骤S622中:相关系数矩阵的计算公式为:其中,所述步骤S623中:所述相关系数矩阵R的特征方程为:|R-λIp|=0;所述步骤S624中:主成分矩阵的计算公式为:其中,j=1,2,…,m。其中,U1称为第一主成分,U2称为第二主成分,…,Up称为第p主成分,以此类推。最后对m个主成分进行加权求和,既可以计算出最后的针入度,权数为每个主成分的方差贡献率。
在一优选实施例中,所述步骤S5包括步骤:选取多个建模区间;针对多个建模区间分别建立多个校正模型;针对每个校正模型获取对应校正相关系数及校正标准偏差;对照预设关系剔除异常值后重新优化计算对应的多个校正模型;分别采用验证集验证多个校正模型的验证结果;选取最优建模区间并根据所述样本建立定量数据模型。在一实施例中,所述步骤“对照预设关系剔除异常值”包括步骤:根据三维空间分布值和线性相关性并参考光谱杠杆值确定异常值后将异常值剔除。
以图2为例,图2左上角视图为三维视图,即样品三维空间分布图。由于做定量分析的样品集为同类样品,它们在主成分空间的分布会因光谱信息的相似性而分布在一个相对集中的空间中,所以可以看到大部分的样品会聚集在一起,少数几个点离散在比较远的区域。这少数的几个离散点就是异常值,模型优化时要将这些点标记出排除掉。
图2右上角视图表示的是光谱杠杆值,即是样品的杠杆值。杠杆值表示样品对模型影响的重要程度。通常,位于被测组分浓度或性质范围两端的样品具有较大的杠杆值,位于被测组分浓度或性质的均值处的样品杠杆值较小。如果某个样品的杠杆值远大于其他值,它有可能不代表被测样品的实际情况,这样的样品可能是异常值,需要给予认真研究。这个指标在模型优化操作中的作用又不是太重要,只要杠杆值不是特别的高,一般不会考虑剔除,该指标模型优化时作为辅助参考指标。
图2左下角视图是PRESS值-主成分数目的关系图。纵坐标是PRESS值,即预测残差平方和。PRESS值是计算方式如下:使用一定数目的主成分建立模型,用这个模型对参加建模的每个样品进行预测,每个样品的预测值和已知值的差的平方和即为该主成分下的PRESS值。PRESS值越小,说明模型的预测能力越好。在计算的多个主成分中,第一主成分最重要,随着主成分数的增加,重要程度依次降低,以致到后来的许多主成分反映的是噪音信息。前面已有论述,如果建立模型时使用的主成分数过少,就不能反应未知样品被测组分产生的光谱数据变化,其模型的预测准确度就会降低,这种情况称之为不充分拟合;如果使用过多的主成分建立模型,就会将一些代表噪音的主成分加到模型中,使模型的预测能力下降,这种情况称为过度拟合。因此,合理确定参加建立模型的主成分数是充分利用光谱信息和滤除噪音的有效方法之一。
图2右下角是线性回归视图,即线性相关性。视图中的横坐标是输入的针入度数据,纵坐标是PLS1回归计算过程中交互验证得出的预测值。图中红点是不包括该点的模型对该点预测得出的验证值;蓝点是所有校正集样品建立的模型对校正集中每一个样品预测得出的校正值。两者都是交互验证过程对模型实际预测能力的模拟,一般以红点的相关性为依据进行异常点的排除。当所有样品的验证值与测试值完全相同时,线性回归视图中的样品点应该都分布在通过原点的斜率为1的直线上,此时的线性相关性为1,这是最理想的状态。我们优化模型的目的就是要剔除线性不好的样品点,使验证值的相关性尽量的接近理想状态的1。那些离回归线比较远的离散点就是异常值,模型优化时要将它们标记出来剔除掉。
模型优化操作中异常值剔除分别采用三维空间分布值和线性相关性这两个统计量来对照检验剔除,同时参考光谱杠杆值,经过异常值的剔除对模型进行优化。
以下以一实施例为例进行说明。
选择3138-2797cm-1、1650-1100cm-1、965-721cm-1作为建模区间,通过PRESS值及变量贡献率选择针入度模型的主因子数。
当选择3138-2797cm-1作为建模区间时,将10组校正集光谱数据导入后,计算,出现计算后的视图,如图2所示。将10组沥青样品的针入度数据和光谱数据一起导入后,采用PLS建立校正模型,得出4个模型视图,如图3所示。由图3可以看出,选择3138-2797cm-1作为建模区间时,建立PLS定量模型,模型校正相关系数为0.8913,校正标准偏差为2.38。此时,需要对对模型性能进行验证,验证结果预测图如图4所示。同时,采用3个验证样本对模型性能进行验证,验证结果如表1所示。
样品编号 | 预测值 | 模型偏差 | 测试值 | 绝对偏差 | 相对偏差(%) |
011016376-1 | 59.176 | 6.088 | 71 | -11.824 | 16.65 |
011016376-2 | 58.817 | 5.918 | 71 | -12.183 | 17.16 |
011016376-3 | 61.024 | 6.505 | 71 | -9.976 | 14.05 |
011016376-4 | 60.792 | 6.307 | 71 | -10.208 | 14.38 |
011016376-5 | 59.735 | 6.68 | 71 | -11.265 | 15.87 |
平均值 | 59.909 | 6.300 | 71 | -11.091 | 15.62 |
标准差 | 0.8698 | 0.2747 | 0 | 0.8698 | 1.2251 |
变异系数 | 0.0145 | 0.0436 | 0 | -0.0784 | 0.0784 |
011016399-1 | 65.346 | 3.652 | 70 | -4.654 | 6.65 |
011016399-2 | 64.361 | 3.787 | 70 | -5.639 | 8.06 |
011016399-3 | 65.877 | 3.977 | 70 | -4.123 | 5.89 |
011016399-4 | 65.579 | 4.1 | 70 | -4.421 | 6.32 |
011016399-5 | 66.074 | 3.995 | 70 | -3.926 | 5.61 |
平均值 | 65.447 | 3.902 | 70 | -4.553 | 6.51 |
标准差 | 0.5975 | 0.1608 | 0 | 0.5975 | 0.8536 |
变异系数 | 0.0091 | 0.0411 | 0 | -0.1312 | 0.1312 |
011016412-1 | 62.076 | 5.501 | 70 | -7.924 | 11.32 |
011016412-2 | 60.846 | 5.343 | 70 | -9.154 | 13.08 |
011016412-3 | 60.036 | 5.492 | 70 | -9.964 | 14.23 |
011016412-4 | 59.65 | 5.418 | 70 | -10.35 | 14.79 |
011016412-5 | 58.607 | 5.187 | 70 | -11.393 | 16.28 |
平均值 | 60.243 | 5.388 | 70 | -9.757 | 13.94 |
标准差 | 1.1657 | 0.1157 | 0 | 1.1657 | 1.6653 |
变异系数 | 0.0194 | 0.0215 | 0 | -0.1194 | 0.1195 |
表1:选取3138-2797cm-1作为建模区间预测验证集样品的结果
当选择1650-1100cm-1作为建模区间时,将10组校正集光谱数据导入后,计算,出现计算后的视图,如图2所示。将10组沥青样品的针入度数据和光谱数据一起导入后,采用PLS建立校正模型,得出4个模型视图,如图5所示。由图5可以看出,选择1650-1100cm-1作为建模区间时,建立PLS定量模型,模型校正相关系数为0.9548,校正标准偏差为1.53。此时,需要对对模型性能进行验证,验证结果预测图如图6所示。同时,采用3个验证样本对模型性能进行验证,验证结果如表2所示。
表2:选取1650-1100cm-1作为建模区间预测验证集样品的结果
当选择965-721cm-1作为建模区间时,将10组校正集光谱数据导入后,计算,出现计算后的视图,如图2所示。将10组沥青样品的针入度数据和光谱数据一起导入后,采用PLS建立校正模型,得出4个模型视图,如图7所示。由图7可以看出,选择965-721cm-1作为建模区间时,建立PLS定量模型,模型校正相关系数为0.8818,校正标准偏差为2.48。此时,需要对对模型性能进行验证,验证结果预测图如图8所示。同时,采用3个验证样本对模型性能进行验证,验证结果如表3所示。
表3:选取965-721cm-1作为建模区间预测验证集样品的结果
当同时选择3138-2797cm-1、1650-1100cm-1、965-721cm-1作为建模区间时,将10组校正集光谱数据导入后,计算,出现计算后的视图,如图2所示。将10组沥青样品的针入度数据和光谱数据一起导入后,采用PLS建立校正模型,得出4个模型视图,如图9所示。由图9可以看出,同时选择3138-2797cm-1、1650-1100cm-1、965-721cm-1作为建模区间时,建立PLS定量模型,模型校正相关系数为0.9720,校正标准偏差为1.21。此时,需要对对模型性能进行验证,验证结果预测图如图10所示。同时,采用3个验证样本对模型性能进行验证,验证结果如表4所示。
样品编号 | 预测值 | 模型偏差 | 测试值 | 绝对偏差 | 相对偏差(%) |
011016376-1 | 63.896 | 1.7 | 71 | -7.104 | 10.00563 |
011016376-2 | 63.429 | 1.578 | 71 | -7.571 | 10.66338 |
011016376-3 | 64.776 | 1.478 | 71 | -6.224 | 8.766197 |
011016376-4 | 64.04 | 1.604 | 71 | -6.96 | 9.802817 |
011016376-5 | 64.796 | 1.486 | 71 | -6.204 | 8.738028 |
平均值 | 64.187 | 1.569 | 71 | -6.813 | 9.595 |
标准差 | 0.5289 | 0.0820 | 0 | 0.5289 | 0.7449 |
变异系数 | 0.0082 | 0.0523 | 0 | 0.0082 | 0.0776 |
011016399-1 | 63.612 | 1.058 | 70 | -6.388 | 9.125714 |
011016399-2 | 63.813 | 1.216 | 70 | -6.187 | 8.838571 |
011016399-3 | 62.93 | 1.438 | 70 | -7.07 | 10.1 |
011016399-4 | 63.264 | 1.639 | 70 | -6.736 | 9.622857 |
011016399-5 | 63.467 | 1.767 | 70 | -6.533 | 9.332857 |
平均值 | 63.417 | 1.424 | 70 | -6.583 | 9.404 |
标准差 | 0.3026 | 0.2613 | 0 | 0.3026 | 0.4323 |
变异系数 | 0.0048 | 0.1835 | 0 | -0.0460 | 0.0460 |
011016412-1 | 59.204 | 1.322 | 70 | -0.15191 | 15.42286 |
011016412-2 | 59.519 | 1.483 | 70 | 3.304386 | 14.97286 |
011016412-3 | 59.835 | 1.599 | 70 | -21.7521 | 14.52143 |
011016412-4 | 59.444 | 1.598 | 70 | -6.5828 | 15.08 |
011016412-5 | 58.758 | 1.491 | 70 | 0.302628 | 16.06 |
平均值 | 59.352 | 1.499 | 70 | -0.0460 | 15.211 |
标准差 | 0.3590 | 0.1014 | 0 | 8.9858 | 0.5128 |
变异系数 | 0.0060 | 0.0677 | 0 | -195.4590 | 0.0337 |
表4:同时选取3138-2797cm-1、1650-1100cm-1、965-721cm-1作为建模区间预测验证集样品的结果
对比发现:选取3138-2797cm-1、1650-1100cm-1、965-721cm-1作为针入度建模区间来预测验证集样品的针入度结果与标准方法结果相比误差较小,从其标准差和变异系数相对比较小可以表明采用红外光谱法测定沥青针入度性质的重复性较好。
同理,在选定建模区间后,采用相同手法来测定具体样品的针入度,此处不再一一赘述。而本发明提出的一种钢桥面铺装用沥青针入度指标的快速检测方法,实现了道路石油沥青针入度指标的快速准确检测,为沥青生产质量控制提供了强大的支持,能够实现沥青在线实时检测;通过红外光谱识别道路石油沥青的来源,便于在实际工程建设中进行监控。红外光谱结合化学计量学的分析方法能够更准确的检测沥青的针入度;引入化学计量学方法作为红外数据库建立的新的理念,提出了红外光谱结合化学计量学的分析方法,具有广阔的应用前景和优势。
以上所述实施例,仅为本发明具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到各种等效的修改、替换和改进等等,这些修改、替换和改进都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (8)
1.一种钢桥面铺装用沥青针入度指标的快速检测方法,其特征在于,包括步骤:
S1、收集多个沥青样品,分别制成沥青分析试样;
S2、采用衰减全反射红外光谱采集方式对所述沥青分析试样进行采样获取红外光谱;
S3、根据所述红外光谱的红外特征吸收峰,筛选对应的最优波长作为样本;
S4、根据所述样本建立定量数据模型,将校正集光谱数据导入数据模型进行计算,生成原始模型视图;
S5、判断是否存在异常值,若是,则剔除异常值后重新计算新的数据模型并生成新的模型视图后进入步骤S6;若否,则将原始模型视图作为新的模型视图并进入步骤S6;
S6、获取待测沥青的红外光谱,并导入新的模型视图中,计算待测沥青的针入度。
2.根据权利要求1所述的沥青针入度的快速检测方法,其特征在于,所述步骤S5包括步骤:
选取多个建模区间;
针对多个建模区间分别建立多个校正模型;
针对每个校正模型获取对应校正相关系数及校正标准偏差;
对照预设关系剔除异常值后重新优化计算对应的多个校正模型;
分别采用验证集验证多个校正模型的验证结果;
选取最优建模区间并根据所述样本建立定量数据模型。
3.根据权利要求2所述的沥青针入度的快速检测方法,其特征在于,所述步骤“对照预设关系剔除异常值”包括步骤:
根据三维空间分布值和线性相关性并参考光谱杠杆值确定异常值后将异常值剔除。
4.根据权利要求1所述的沥青针入度的快速检测方法,其特征在于,所述步骤S3包括步骤:
对所述样本包含的数据进行平滑、降噪处理;所述平滑、降噪处理采用Savitzky-Golay、百分位滤波器、FFT滤波器、LOVESS、LOESS或二项式法;边界条件选择对称、重复、周期性或外推法。
5.根据权利要求1所述的沥青针入度的快速检测方法,其特征在于,所述步骤“剔除异常值后重新计算新的数据模型并生成新的模型视图”包括步骤:
生成新的模型视图的相关系数并采用验证集对新的模型视图进行验证。
6.根据权利要求1所述的沥青针入度的快速检测方法,其特征在于,所述步骤S6包括步骤:
S61、根据待测沥青的红外光谱在新的模型视图中的位置对待测沥青的主成分进行分析,获取全部主成分数及每种主成分的贡献率;
S62、根据贡献率选出特征主成分,并对特征主成分进行综合评价,根据综合评价获取待测沥青的针入度。
7.根据权利要求6所述的沥青针入度的快速检测方法,其特征在于,所述步骤S62中包括步骤:
S621、标准化采集p维随机向量和n个样品构造样本阵,并对样本阵进行标准化变换获取标准化矩阵;
S622、根据所述标准化矩阵计算相关系数矩阵;
S623、根据所述相关系数矩阵计算预设特征方程的特征根和特征向量,并根据所述特征根确定主成分;
S624、根据所述标准化矩阵和所述特征向量获取主成分矩阵;
S625、根据所述主成分矩阵对主成分进行综合评价。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010097162.7A CN111257265A (zh) | 2020-02-17 | 2020-02-17 | 一种钢桥面铺装用沥青针入度指标的快速检测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010097162.7A CN111257265A (zh) | 2020-02-17 | 2020-02-17 | 一种钢桥面铺装用沥青针入度指标的快速检测方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111257265A true CN111257265A (zh) | 2020-06-09 |
Family
ID=70947363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010097162.7A Pending CN111257265A (zh) | 2020-02-17 | 2020-02-17 | 一种钢桥面铺装用沥青针入度指标的快速检测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111257265A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117093841A (zh) * | 2023-10-18 | 2023-11-21 | 中国科学院合肥物质科学研究院 | 小麦透射光谱的异常光谱筛选模型确定方法、装置及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101936882A (zh) * | 2010-08-10 | 2011-01-05 | 江苏大学 | 一种作物氮素和水分无损检测方法及装置 |
CN102081039A (zh) * | 2010-08-17 | 2011-06-01 | 江苏大学 | 一种环境可控的作物营养水分高光谱图像检测装置 |
CN107389608A (zh) * | 2017-07-13 | 2017-11-24 | 山东华唐环保科技有限公司 | Libs激光诱导光谱分析在煤质检测中的在线应用 |
CN107782693A (zh) * | 2017-10-25 | 2018-03-09 | 中石油燃料油有限责任公司研究院 | 一种沥青针入度的红外光谱分析方法 |
CN108593582A (zh) * | 2018-04-12 | 2018-09-28 | 山东建筑大学 | 一种红外光谱快速判定沥青油源的方法 |
CN109374560A (zh) * | 2018-12-14 | 2019-02-22 | 洛阳理工学院 | 一种沥青中蜡含量的测定方法 |
CN109615145A (zh) * | 2018-12-27 | 2019-04-12 | 中国海洋石油集团有限公司 | 一种快速预测不同老化程度基质沥青的物理性能的方法 |
-
2020
- 2020-02-17 CN CN202010097162.7A patent/CN111257265A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101936882A (zh) * | 2010-08-10 | 2011-01-05 | 江苏大学 | 一种作物氮素和水分无损检测方法及装置 |
CN102081039A (zh) * | 2010-08-17 | 2011-06-01 | 江苏大学 | 一种环境可控的作物营养水分高光谱图像检测装置 |
CN107389608A (zh) * | 2017-07-13 | 2017-11-24 | 山东华唐环保科技有限公司 | Libs激光诱导光谱分析在煤质检测中的在线应用 |
CN107782693A (zh) * | 2017-10-25 | 2018-03-09 | 中石油燃料油有限责任公司研究院 | 一种沥青针入度的红外光谱分析方法 |
CN108593582A (zh) * | 2018-04-12 | 2018-09-28 | 山东建筑大学 | 一种红外光谱快速判定沥青油源的方法 |
CN109374560A (zh) * | 2018-12-14 | 2019-02-22 | 洛阳理工学院 | 一种沥青中蜡含量的测定方法 |
CN109615145A (zh) * | 2018-12-27 | 2019-04-12 | 中国海洋石油集团有限公司 | 一种快速预测不同老化程度基质沥青的物理性能的方法 |
Non-Patent Citations (2)
Title |
---|
王策: "光谱反射率预测甘蔗叶片氮磷钾含量及传感器开发", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
袁也: "常用沥青品牌识别及技术性质快速检测研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117093841A (zh) * | 2023-10-18 | 2023-11-21 | 中国科学院合肥物质科学研究院 | 小麦透射光谱的异常光谱筛选模型确定方法、装置及介质 |
CN117093841B (zh) * | 2023-10-18 | 2024-02-09 | 中国科学院合肥物质科学研究院 | 小麦透射光谱的异常光谱筛选模型确定方法、装置及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110687072B (zh) | 一种基于光谱相似度的校正集和验证集的选择及建模方法 | |
JP4856993B2 (ja) | 自己診断型自動分析装置 | |
Botelho et al. | Development and analytical validation of robust near-infrared multivariate calibration models for the quality inspection control of mozzarella cheese | |
CN107703097B (zh) | 利用近红外光谱仪构建快速预测原油性质的模型的方法 | |
US20210247367A1 (en) | Workflow-based model optimization method for vibrational spectral analysis | |
Westerhaus et al. | Quantitative analysis | |
CN117195135B (zh) | 一种水污染异常溯源检测方法及系统 | |
US6477516B1 (en) | System and method for predicting parameter of hydrocarbon with spectroscopy and neural networks | |
CN112179871A (zh) | 一种酱类食品中己内酰胺含量无损检测的方法 | |
CN104316492A (zh) | 近红外光谱测定马铃薯块茎中蛋白质含量的方法 | |
van der Voet et al. | Inter-laboratory, time, and fitness-for-purpose aspects of effective validation | |
CN108663334B (zh) | 基于多分类器融合寻找土壤养分光谱特征波长的方法 | |
CN111257265A (zh) | 一种钢桥面铺装用沥青针入度指标的快速检测方法 | |
de Souza et al. | Validation of the near infrared spectroscopy method for determining soil organic carbon by employing a proficiency assay for fertility laboratories | |
CN110210005A (zh) | 一种无参考值的光谱波数选择方法 | |
CN102998276A (zh) | 由红外光谱测定原油实沸点蒸馏曲线的方法 | |
CN111896497B (zh) | 一种基于预测值的光谱数据修正方法 | |
Lodder et al. | Quantile analysis: a method for characterizing data distributions | |
CN109709060B (zh) | 一种沥青软化点、针入度和质量损失的测定方法 | |
CN111401794A (zh) | 一种基于近红外光谱的饲料品质控制方法 | |
CN111257269A (zh) | 一种沥青软化点的快速检测方法 | |
Bellocchi et al. | Expanding horizons in the validation of GMO analytical methods: fuzzy-based expert systems | |
Christensen et al. | Method evaluation, quality control, and external quality assurance systems of analytical procedures | |
Conti et al. | Uncertainty in environmental analysis: theory and laboratory studies | |
Fajgelj et al. | A critique on available in-house method validation documentation |
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 |