化香樹果序總黃酮提取動力學研究

李 冬，李穩宏*，廉媛媛，高 蓉，李學坤，朱驟海

1西北大學化工學院，西安710069;2陜西地礦醫院藥劑科，西安710014;3第四軍醫大學化學系，西安710032

化香樹果序總黃酮提取動力學研究

李 冬1，李穩宏1*，廉媛媛2，高 蓉3，李學坤1，朱驟海1

1西北大學化工學院，西安710069;2陜西地礦醫院藥劑科，西安710014;3第四軍醫大學化學系，西安710032

以化香樹果序為原料，60%的乙醇水溶液為溶劑提取黃酮類化合物，在不同溫度下對總黃酮的傳質動力學進行了研究。采用平板模型，以Fick第二定律為基礎，建立了化香樹果序總黃酮提取的動力學方程，求得了速率常數、活化能、相對萃余率等一系列動力學參數。研究結果可為化香樹果序總黃酮提取工程放大和深入理論研究提供一定的依據。

化香樹果序;總黃酮;提取;動力學

化香樹果序是化香樹的干燥果序，具有清熱解毒、活血化瘀、清腫排膿、通竅止痛的功效［1］。黃酮類化合物具有抗菌、消炎、抗突變、降壓、清熱解毒、鎮靜、利尿、抗氧化、抗癌防癌、抑制脂肪氧化酶等多種生物活性。其顯著的生物活性對于研究化香樹果序治療鼻炎鼻塞有較大的價值［2-5］。目前，尚無化香樹果序黃酮提取動力學方面的報道。本文以Fick第二定律為基礎，建立化香樹果序總黃酮提取的動力學模型，并依此對速率常數和相對萃余率等一系列動力學參數進行求解，以求為化香樹果序提取工業化放大和深入理論研究提供一定的依據。

1 提取動力學模型的建立

天然產物有效成分的提取過程中，往往同時伴隨其它成分的浸出，機理較為復雜。但一般包括三個步驟，即:(1)溶劑向固體植物顆粒內部滲透，浸潤顆粒并溶解有效成分;(2)被溶解的有效成分以分子形態從植物顆粒內部擴散遷移到固液界面上; (3)有效成分從固液界面向溶劑主體擴散［6］。其中內擴散常被認為是整個提取過程的速率控制步驟［7］。

化香樹果序內部的總黃酮濃度隨提取時間的增加不斷降低，為不穩定擴散過程。為便于分析，本研究中采用平板模型，以單果序葉片為計算對象，現假設:(1)顆粒的厚度是均勻的;(2)有效成分的擴散方向是沿顆粒的厚度方向進行的;(3)萃取開始時刻顆粒內各成分是均勻分布的;(4)顆粒表面的傳質阻力忽略不計且在任意取樣時間間隔內，有效成分的擴散系數不變;(5)顆粒與溶劑的溫度是相同的，且溫度是均勻分布的。

由于平板中的擴散在中心面的兩側對稱，因此，任意時刻濃度分布在兩側必然對稱。將坐標原點取在平板中心，研究半個平板。則根據Fick擴散第二定律，有

將其簡化建立控制方程，并確定初始條件和邊界條件為

初始條件c(d，0)=c0

式中:t為提取時間，min;d為距平板中心的距離，m;c為提取過程中t時刻顆粒內距平板中心為d處的有效成分質量濃度，%;D為單片化香樹果序厚度的一半，m;DS為內擴散系數，m2/min。

定義無量綱濃度c*=(c-c1)/(c0-c1)

將式(2)改寫為

初始條件c*(d，0)=1

用分離變量法求解式(3)，令無量綱濃度c*(d，t)=X(d)T(t)，有

式(4)只有兩側均為常數時才能成立，令該常數為－f2，得到兩個常微分方程

它們的通解分別為

代入邊界條件式，得到常數C1=0，特征值為

由解的疊加原理，得到無量綱濃度的解

由初始條件，確定常數

最后，得到無量綱濃度

式(10)表示了忽略表面傳質阻力情況下，大平板內有效成分擴散濃度是時間和位置的函數。提取是內擴散控制，有效成分一經擴散到表面，即進入溶液主體，在d=D處，c=C，可把式(10)變為

式中:C0為提取液中溶質的初始質量濃度，%; C∞為提取平衡時提取液中溶質的初始質量濃度，%。非穩態傳質的濃度分布式為一無窮級數，一般取第一項計算即可，故

式(13)便是根據Fick第二定律建立的動力學模型。該模型可表達化香樹果序顆粒尺寸、提取時間、提取溫度和料液比與有效成分提取液濃度之間的關系。

2 材料與方法

2.1 材料

化香樹果序由陜西香菊藥業集團有限公司提供;蘆丁，中國醫藥(集團)上海化學試劑公司提供;乙醇、亞硝酸鈉、硝酸鋁、氫氧化鈉等試劑均為分析純。

紫外可見分光光度儀，UV-2501PC，日本島津公司;電熱恒溫水浴鍋，DK-98-1型，天津市泰斯特儀器有限公司;精密增力電動攪拌器，JJ-1型，常州國華電器有限公司;旋轉蒸發儀，RE-52AA型，上海嘉鵬科技公司;分析天平，AL204型，瑞士梅特勒-托利多儀器有限公司。

2.2 試驗方法

2.2.1 化香樹果序飽和吸溶劑率的測定

精確稱取不同質量化香樹果序，用滴定管將60%乙醇水溶液滴加，每次滴加后充分振蕩并浸潤數分鐘，使果序充分吸收溶劑，直至充分吸濕且無多余溶劑，記錄消耗的溶液體積。則P=(V1-V2)/M，式中:P為飽和吸溶劑率，mL/g;V1為滴定管初讀數，mL;V2為滴定管終讀數，mL;M為化香樹果序質量，g。

2.2.2 提取方法

將厚度為0.25 mm的化香樹果序12 g加入500 mL的三口燒瓶中，根據果序的飽和吸溶劑率先用適量的60%乙醇水溶液室溫浸潤10 min，再立刻向三口燒瓶中加入己加熱到設定溫度的60%乙醇水溶液至總量200 mL，在此溫度下進行常壓恒溫水浴回流提取，同時在100 rpm下進行攪拌并計時，按試驗設計在不同時間間隔定時取樣200 μL。

提取時間為:t1=10 min，t2=30 min，t3=60 min，t4=90 min，t5=120 min，t6=150 min，t7=180 min，t8=210 min，t9=240 min，t10=270 min，t11=300 min，t12=360 min;

提取溫度為:T1=313 K，T2=323 K，T3=333 K，T4=343 K，T5=353 K。

2.2.3 標準曲線的繪制［8］

精密吸取0.1 mg/mL的蘆丁標準品溶液0、1、2、3、4、5 mL置于10 mL容量瓶中，各加入30%乙醇補齊至5 mL，依次分別加入5%NaNO2溶液0.3 mL，靜置6 min;10%Al(NO3)3溶液0.3 mL，靜置6 min;4%NaOH溶液4 mL;最后用蒸餾水稀釋到10 mL并靜置20 min。用紫外可見分光光度儀在波長510 nm處測定吸光度，繪制標準曲線，得到回歸方程為:C=0.0927A－0.000333，R2=0.9997。

2.2.4 樣品濃度測定方法［9］

每次取樣后浸提液總體積將減少200 μL，多次取樣將導致提取液總體積的減少，同時也造成提取液中有效成分質量的減少。因而，必須對浸提稀釋液所測得的濃度進行修正，修正式為:

式中:Cn為第n次取樣的濃度修正值，μg/mL; C'n為第n次取樣的濃度測定值，μg/mL;n為取樣次數，n=l，2，3，……。C'n的測定方法采用2.2.3得到的標準曲線進行換算得到。

3 結果和討論

3.1 化香樹果序飽和吸溶劑率

根據2.2.1的方法得到化香樹果序的吸溶劑率P=0.68 mL/g。

3.2 速率常數的求解

不同溫度下提取液中總黃酮濃度實驗結果見表1。

表1 不同溫度下提取液中總黃酮濃度(%)Table 1 Concentration of total flavonoids at different temperature(%)

由表1可以看出，不同浸提溫度條件下，總黃酮濃度差異較大，提取溫度越高，提取速度越快，提取液濃度越大。當溫度為333 K、343 K、353 K時，提取時間達到210 min時，基本達到平衡，溫度越低，達到平衡時所需的時間越長。提取時間小于210 min時，達到相同濃度所需時間不同，溫度越高，所需時間越短。在提取過程中，提取液濃度在較長時間內基本保持不變，可認為此時己基本達到浸提平衡，該濃度可作為平衡濃度。

利用表1的數據得到ln［C∞/(C∞-C)和提取時間t的關系結果見表2。

表2 不同溫度下ln［C∞/(C∞-C)］與t的關系Table 2 Relationship between ln［C∞/(C∞-C)］and t at different temperature

323 ln［C∞/(C∞-C)=0.0125t+1.23390.9941 4.02 2.08 2.85 333 ln［C∞/(C∞-C)=0.015t+1.1431 0.9879 4.26 2.50 2.90 343 ln［C∞/(C∞-C)=0.0173t+1.109 0.9960 4.49 2.88 3.01 353 ln［C∞/(C∞-C)=0.0203t+1.124 0.9846 4.58 3.38 3.09

從表2可以看出，所建立的動力學方程式與實驗結果能較好的吻合，ln［C∞/(C∞-C)與t之間線性關系良好，R2均在0.98以上。表觀速率常數隨著溫度的升高較顯著地增大，有效成分的提取速度越快，同時可以看出C0隨溫度的升高也略有升高。

3.3 相對萃余率的求解

設相對萃余率y=(C∞-C)/(C∞-C0)，則式(12)變為y=(4/π)exp(-λt)，分別利用實驗數據作y對提取時間t的關系圖，結果見圖1。

圖1 不同溫度下y與t的關系Fig.1 Relationship between y and t at different temperature

由圖1可知，擬合方程的復相關系數R2均在0.98以上，曲線的擬合精度較高，所以化香樹果序中總黃酮的提取符合指數模型。由相對萃余率的定義可知，其值越小說明提取率越大，圖中可以看出隨著時間的增加或溫度的升高相對萃余率不斷的減小，提取率在不斷增大，但是增長的幅度下降。所以，實際提取過程中并非時間越長越好，還需考慮提取的效率以及設備的利用率等因素。

3.4 活化能的求解

圖2 –lnk與1000/T關系圖Fig.2 Relationship between–lnk and 1000/T

由化學反應動力學可知，在低溫下表觀擴散速率常數是溫度的函數。那么，在總黃酮的提取過程中其表觀擴散速率常數與提取溫度的關系應遵循Arrhenius公式，對-lnk和1/T作圖，結果見圖2。

由圖2可知，lnk和1/T線性關系良好，其回歸方程為:lnk=1776.4/T+2.9659。據此可以算出化香樹果序提取總酚過程中的活化能，Ea=1367× 8.314=11365 J/mol。

3.5 有效擴散系數的求解

提取過程的內擴散系數可以通過速率常數計算。由λ=(π2Ds)/(4D2)可知，速率常數是內擴散系數Ds和顆粒尺寸的函數。由Ds對T作圖得圖3。得到回歸方程為:104Ds=0.0056exp(0.0153T)，由此可以計算出不同溫度下化香樹果序中總黃酮和總酚提取的內擴散系數。

圖3 Ds和T的關系圖Fig.3 Relationship between Ds and T

4 結論

采用平板模型，根據Fick第二定律為基礎，通過一些假定和簡化，建立了化香樹果序總黃酮提取的動力學模型，模型計算值與實驗數據吻合良好。通過對化香樹果序總黃酮的提取試驗研究，結合所建立的動力學模型，求得了速率常數、相對萃余率、活化能以及有效擴散系數等一系列有價值的動力學參數。研究結果可為化香樹果序有效成分的提取工藝條件優化以及工藝設計提供相應的理論依據和參考。

1 Northwesters Institute of Botany of Academia Sinica.Flora of Qin Ling，1st(秦嶺植物志).Beijing:Science Press，1970.121-125.

2 Pei LP(裴凌鵬)，Hui BL(惠伯棣)，Jin ZL(金宗濂).Review on health function processing technology and market prospects of flavonoids.Food Sci，2005，25:203-207.

3 Peng F(彭芳)，Chen ZH(陳植和).Biological Effects of Flavones.J Da Li Med Coll(大理醫學院學報)，1998，7 (4):52-54.

4 Huang HS(黃河勝)，Ma CG(馬傳庚)，et al.The advance of research on pharmacological effects of flavonoids.Chin J Chin Mat Med(中國中藥雜志)，2000，25:589-592.

5 Zhang AL(張安靈)，Gao JM(高錦明)，et al.The distribution，development and utilization of the flavonoids.J Northwest Forest Univ(西北林學院學報)，2000，15:69-74.

6Chu MQ(儲茂泉)，Gu HC(古宏晨)，Liu GJ(劉國杰).Kinetic model on medicinal herb extraction process.Chin Tradit Herb Drugs，2000，31:504-506.

7 Spiro M，Page CM.The kinetics and mechanism of caffeine infusion from coffee:hydrodynamic aspects.J Sci Food Agric，1984，35:925-930.

8 GaoR(高蓉)，Li WH(李穩宏)，Liu MX(劉明霞).Research on extraction of flavonoid infructescence of Platycarya strobilacea Sieb.et Zucc with diferent distilling methods.Food Sci，2007，28:230-232.

9 Liang Y(梁英).Study on Extraction Kinetics and Extraction Techniques of Flavonnoids from Scutellaria radix.Beijing: China Agricultural University，2005.24-26.

Kinetics on the Extraction of Total Flavonoids from Infructescence of Platycarya Strobilacea Sieb.et Zucc

LI Dong1，LI Wen-hong1*，LIAN Yuan-yuan2，GAO Rong3，LI Xue-kun1，ZHU Zhou-hai11School of Chemical Engineering，Northwest University，Xi’an 710069，China;2Department of Pharmacy，Shaanxi Dikuang Hospital，Xi’an 710014，China; 3Department of Chemistry，Fourth Military Medical University，Xi’an 710032，China

In this paper，60%ethanol was used as solvent to extract flavonoids from Platycarya Strobilacea Sieb.et Zucc.，and the mass transfer kinetics of polyphenol at different temperature was investigated.According to Fick’s second law of diffusion，the kinetic equation for Platycarya Strobilacea Sieb.et Zucc.extraction process was established with plate model，and the parameters，such as K(rate constant)，Ea(activation energy)，y(yield)were gained.The results could provide the valuable theory basis for the technical design and further research of total flavonoids extraction process.

Platycarya Strobilacea Sieb.et Zucc.;total flavonoids;extraction;dynamics

1001-6880(2011)04-0689-05

2009-10-10 接受日期:2010-02-04

陜西省教育廳重大產業化項目(2007JC18);西安市科技創新支撐計劃(F07071)

*通訊作者 E-mail:liwenhong@nwu.edu.cn

TQ028.96;Q946.91

A