基于蟻群算法對有機物熱容的預測研究

白帆　時靜潔（江蘇富淼科技股份有限公司，江蘇　張家港　215613）

基于蟻群算法對有機物熱容的預測研究

白帆時靜潔（江蘇富淼科技股份有限公司，江蘇張家港215613）

根據定量構效關系（QSPR）原理，研究有機物熱容與其分子結構間的內在定量關系。以650種有機化合物作為樣本集，隨機選擇520種作為訓練集，130種作為測試集，采用蟻群算法（ACO）進行變量選擇，得到5個特征描述符作為模型的輸入變量，結合多元線性回歸（MLR）方法建立了蟻群-多元線性回歸（ACO-MLR）預測模型。研究結果表明：ACO-MLR模型的訓練集和測試集的復相關系數分別為0.9430和0.9344，預測效果令人滿意。該研究提供了一種新的預測有機化合物熱容的方法。

熱容；蟻群算法；多元線性回歸；預測

熱容表示系統升高單位熱力學溫度時所吸收的熱［1］，是化合物的基本熱學性質之一。熱容的變化提供了化合物相變信息以及化合物分子結構的變化，獲悉該方面信息是非常重要的。液體的熱容值在工程應用中是一個非常重要的性質，在能量轉移以及熱動力學的計算中也是必需的參數［2］。熱容性質涉及到很多領域，其中包括石油化工行業、制藥和食品工業等。在這些行業的很多過程都涉及到了熱容，如蒸餾過程、蒸發過程、提取過程和加熱過程［3］。比如在某個溫度范圍內進行熱交換，離子液體會呈現出某些性質而使物質更穩定，其中就包括熱容。熱容的實驗測定有很多方法，但測定比較復雜，例如差示掃描量熱儀法［4］、熱線式方法［5］、溫度振蕩量熱法［6］等?？紤]到獲取熱容的重要性以及實驗測定的復雜性，發展方便可靠的理論預測熱容方法迫在眉睫。目前，關于理論預測有機物熱容的文獻比比皆是，但有機物液體熱容的QSPR預測模型卻為數不多。因此，本章擬從分子結構角度出發，應用QSPR研究體系，系統研究有機物熱容與其分子結構間的內在定量關系，建立相應的QSPR預測模型，并評價其模型的預測效果。

1　實驗樣本說明

可靠的預測模型必須建立在可靠的實驗樣本基礎之上。在定量結構-性質關系研究中，數據的精確性對于最終的結果具有至關重要的影響。本文研究的樣本集為650種有機化合物的熱容，包括碳氫化合物、鹵素化合物、醇、酯、醛、酮等。這些化合物涵蓋了廣泛的化學多樣性空間，為建立健壯、有效的預測模型奠定了基礎。此外，為了避免由于數據來源不統一造成的誤差，所有化合物熱容數據均來自于同一個數據庫——《有機化合物實驗物性數據手冊：含碳、氫、氧、鹵部分》。

2　分子描述符的計算與預篩選

分子描述符的計算首先是在HyperChem 7.5軟件中對分子結構進行輸入與優化，導入Dragon 2.1軟件中計算分子描述符，獲得18類共1481種分子描述符，其計算的描述符種類包括結構描述符，幾何描述符，拓撲描述符等。大量的分子描述符中必然存在冗余信息，因此對常數或者近似常數的描述符以及描述符間相關系數達到0.95以上的描述符進行刪除，減少至660個。如此龐大的分子描述符數量還是無法有效進行QSPR模型的建立，因此，本文采用ACO對分子描述符進行進一步篩選。

3　ACO-MLR熱容預測模型

運用ACO對660個描述符進行進一步篩選，ACO是在VC++6.0中采用C語言編程，螞蟻數量設為500，揮發率設為0.9，進行不斷迭代直至收斂，最終獲得5個特征描述符。隨后，以ACO篩選的5個描述符作為輸入參數，運用SPSS17.0統計軟件中的多元線性回歸模塊，在95%的置信區間內，得到了MLR預測模型，結果如下：

Cp=21.594-10.970*ATS4m+8.598*HATS4u+0.787

*TIC4+20.824*S1K+1.605*C-024（5.2）

n=520，R2=0.943，SD=21.250，F=1697.199，p＜0.001

式中，Cp為熱容值，n為訓練集樣本數目，R2為決定系數，SD為模型標準誤差，F為F

檢驗值，F實際=1697.199＞F理論（5，514，0.05）=2.25，模型的顯著性概率p遠小于0.05，認為該回歸方程及所篩選變量的影響均是顯著的。

為檢驗模型的擬合能力和外部預測能力，應用該模型分別對訓練集和測試集進行預測，其訓練集和測試集的復相關系數分別為0.9430和0.9344。

4　結語

本文運用蟻群算法（ACO）作為分子描述符篩選方法，與多元線性回歸（MLR）方法進行組合建立了ACO-MLR模型，對650種有機化合物的熱容進行了QSPR研究。所建立的QSPR熱容預測模型的預測值與實驗值非常接近，其R2均在0.9以上，預測誤差在可接受范圍之內。因此，本文所建模型可用于對有機物熱容的預測，為其提供了一種新的有效方法，對于化工安全設計和風險評價研究具有重要的意義。

［1］趙敏，強曉明，章韋芳.一種新型比熱容測量儀［J］.大學物理.2016（2）.

［2］蔣華龍，周大偉，劉旭焱等.壓力下MoS2的結構相變以及熱動力學性質的第一原理研究［J］.低溫物理學.2014（5）.

［3］陳永艷，田瑞，王亞輝等.膜蒸餾系統中熱容腔結構參數優化的試驗研究［J］.工程熱物理學報.2015（2）.

［4］羅昭強，譚蓬.應用DSC法測量有機熱載體比熱容［J］.鍋爐制造.2016（1）.

［5］鄭明陽，王天為，張國鋒.理想氣體熱力學過程吸放熱情況的圖像判斷法［J］.物理與工程.2014（3）.

［6］張濤，孫冰.航天器再入全過程軸對稱燒蝕熱防護數值仿真研究［J］.宇航學報.2011（5）.