摘  要：目的  尋找中藥水提液濃縮過程瞬時能耗計算簡便方法。方法  以玉屏風散為模型藥物，測定其水提液濃縮過程中多種溶液環(huán)境特征參數的動態(tài)變化數據，構建中藥溶液環(huán)境特征參數白利度（簡便、易測、準確）與濃縮能耗計算需要的相關參數的相關性模型。結果  白利度與固含率、密度、比熱容、導熱率等具有優(yōu)越的線性相關，相關系數（R2）均高達0.99左右。結論  可通過測定白利度及構建與其他物理參數如比熱容（C_p）、導熱率（k_f）等關系模型，計算濃縮的動態(tài)能耗；為中藥制藥過程工程原理研究開辟新思路。

傳統(tǒng)中藥制藥過程中，提取、濃縮、干燥等工段都需要使用大量蒸汽^[1]。其中，中藥濃縮工段的能耗驚人，濃縮工段所用蒸汽量占總蒸汽量60%以上^[2]。膜濃縮（反滲透、膜蒸餾等）比一般熱法濃縮如單效、雙效蒸發(fā)節(jié)能50%^[3]。而近期Weaver等^[4]發(fā)現，使用反滲透預濃縮技術可比單獨用熱效濃縮技術節(jié)省能耗85%以上。物料濃縮的原理是實現水（或其他溶劑）與可溶性物質（溶質）的分離，目前常見的產業(yè)化濃縮技術有熱法蒸發(fā)（單效、多效蒸發(fā)），膜技術（納濾、反滲透、膜蒸餾等），均屬于物理分離過程。熱法、膜法濃縮是海水淡化的常見手段，已有比較成熟的能耗計算與評估方法。但中藥制藥領域因缺乏相關決策研究，無法采用這些現有的方法。如研究表明中藥與相關天然產物的濃縮過程遠比海水淡化過程復雜，其主要特征是物料的溶液環(huán)境特征參數發(fā)生了明顯的動態(tài)變化，因此海水淡化能耗計算方法對于中藥水提液體系存在較大誤差^[5-6]。因此，亟需提出一種將中藥溶液環(huán)境特征參數變化融入到濃縮能耗計算的新方法。溶液環(huán)境是指溶液體系所具有的電導率、折光指數（實際應用中常換算成白利度）、黏度、pH值、離子強度等特征性質。在中藥應用領域，精制、濃縮等制藥工藝過程所面對的溶液環(huán)境體系即為中藥液體物料。上述溶液環(huán)境特征參數，雖然彼此間存在大量非線性、高噪聲、多因子的復雜關系，但借助理論化學對簡單物質研究的成果，可從中抽提出若干參數和概念，進而運用人工智能技術，從已知數據和實驗事實中抽提規(guī)律性，用于簡化中藥制藥過程所采用膜分離等先進技術的工藝優(yōu)化設計。如張劉紅等^[7]已建立滲透壓與電導率相關性模型，借助中藥物料濃縮過程中電導率的變化推測其他特征參數的變化，可實現中藥溶液環(huán)境各相關特征參數的在線實時監(jiān)測，為中藥制藥的智能監(jiān)控提供更簡便、經濟的手段，并為中藥綠色節(jié)能膜濃縮技術的推廣應用提供有力的計算依據與理論支持^[7-10]。中藥物料的濃縮過程還涉及到其濃度變化所帶來的比熱容、固體含量變化以及密度等溶液環(huán)境表征參數的變化。其中，折光指數是國內外評判天然產物及中藥濃縮物料常用的檢測指標^[11-13]，也是中藥制藥工業(yè)檢測物料比重所用傳感器不可或缺的重要參數^[14-15]。利用白利度的變化推測其他物理參數的變化，對中藥制藥在線檢測和濃縮所需能耗計算可達到事半功倍的作用。但過去的研究中并未對特征參數如何影響中藥物料濃縮工段的能耗變化做出詳細說明（尤其是中藥物料的膜濃縮能耗計算）。多年來中藥制藥工程設計中，因缺乏基本的中藥物性數據（如不同中藥物料的密度、黏度、表面張力、導熱系數、擴散系數等）^[10]，往往憑經驗或采用經驗方式估算，甚至用相近物質的物性代替，其結果必然導致工藝技術選擇或設計的失真甚至失敗。從而使中藥生產的規(guī)范化、現代化難以付諸現實。本實驗擬以玉屏風散為模型藥物，探討濃縮過程對中藥溶液環(huán)境特征參數的影響及各特征參數的相關性，指出相關特征參數變化與中藥物料濃縮工段能耗的關系，為中藥濃縮工藝的瞬時能耗計算提供一種簡便的模式，同時也為探討中藥制藥過程工程原理開辟新視野。

1  儀器與材料

1.1  藥材

防風購于安徽道源堂中藥飲片有限公司，批號200501，白術購于毫州市永剛飲片廠有限公司，批號200203，蜜炙黃芪購于赤峰榮興堂藥業(yè)有限責任公司蒙中藥飲片廠，批號2005046，上述3味藥材經本課題組郭立瑋教授鑒定，分別為傘形科防風屬植物防風Saposhnikovia divaricata (Turcz.) Schischk.的干燥根；菊科蒼術屬植物白術Atractylodes macrocephala Koidz.的干燥根莖；豆科黃芪屬植物膜莢黃芪Astragalus membranaceus (Fisch.) Bge.干燥根的炮制（蜜炙）加工品。

1.2  儀器

BSA224S-CW型電子分析天平，賽多利斯科學儀器（北京）有限公司；LC-DRT-948數顯折射儀、NDJ-5S數顯黏度計，上海力辰邦西儀器科技有限公司；TP101筆式探針數顯溫度計，測量范圍?50～300 ℃，蘇泊爾C22-IH66E8電磁爐，2200 W，浙江蘇泊爾生活電器有限公司；簡悠便攜式炆火電陶爐1200 W，中山市諾潔仕電器有限公司；DHG-9145A電熱鼓風干燥箱，上海恒勤儀器設備有限公司；OM819C滲透壓儀（摩爾濃度滲透壓儀，摩爾滲透壓儀），北京雅森博科科學儀器有限公司。

2  方法

2.1  提取方法

參照“十二五”普通高等教育本科國家級規(guī)劃教材《方劑學》中玉屏風散方各藥味的劑量比^[16-17]，準確稱取藥材如下：防風50.0 g、白術100.0 g、蜜炙黃芪100.0 g。用2500 mL純凈水（經過實驗室反滲透濾膜制備的純凈水）煎煮1 h，分離藥渣與提取液。再在藥渣中再加入2500 mL純凈水煎煮1 h，分離藥渣與提取液。將2次提取液均勻混合后，150目紗布濾過，收集濾液備用。

2.2  中藥物料濃縮實驗

采用敞口鍋（直徑為24 cm，容積為6800 mL的實驗室用小型敞口鍋）對上述濾液加熱蒸發(fā)進行濃縮。采集不同時間段濃縮物料，進行標記，按照“2.3”項所述方法檢測相應的溶液環(huán)境特征參數。

2.3  溶液環(huán)境特征參數測試

2.3.1  白利度  量筒量取100 mL樣品，其中1～2 mL樣品用來測量白利度，分2組樣品測試讀數，由機器分別自動讀取15次數據，直到數據穩(wěn)定不變再記錄數值。取得的2組數值記錄下來，并且算得平均值。由于實驗室采用的是離線折光儀，需要取樣讀數，而樣品數量僅需1～2滴，因此樣品溫度很快就被降至室溫。為了統(tǒng)一樣品折光度與白利度讀數溫度，此研究的白利度測試均在室溫（30±5）℃下進行。

2.3.2  黏度  待同一批次取得的物料溫度降至室溫（30±5）℃，用黏度計讀取物料的電導率和黏度，讀取3組數據取得平均值。

2.3.3  密度  將干燥好的燒杯放置天平上清零。從上述同一批次取的物料中，用量取精準量取一定體積的物料并記錄體積并倒入干燥的燒杯中。密度測量時物料溫度為（45±5）℃。將稱有物料的燒杯放置天平上，準確稱量其質量并記錄。同一批次、濃度物料的密度可用傳統(tǒng)密度定義公式算得。

2.3.4  固含率  烘干玻璃培養(yǎng)皿并進行準確稱定質量，記為m₁。將同一批次取得的物料放入烘干的玻璃培養(yǎng)皿，并準確稱定質量，記為m₂。將已稱定質量的物料與玻璃培養(yǎng)皿放入電熱鼓風干燥機中，設置50 ℃干燥4 h以上確保水份完全蒸發(fā)。待干燥的物料與玻璃培養(yǎng)皿放置室溫后，準確稱定質量，記為m₃。此樣品的固含率可用公式（1）計算。

固含率＝(m₃－m₂)/(m₂－m₁)                （1）

2.3.5  物料比熱容（C_p）與導熱率（k_f）  濃縮中的物料比熱容與導熱性變化由此部分實驗操作得到。比熱容是由熱量（Q）、質量（M）和溫度變化（ΔT）的比值而定義的^[18]。

C_p＝Q/MΔT                             （2）

在同一組測量實驗中準確稱量相同質量的水和物料。水與物料在相同時間內得到的熱量一樣（均被同樣加熱設備加熱），根據已知純水的比熱容［C_pw，4.18 kJ/(kg·K)］，利用同樣質量的水和物料在相同時間內的溫度變化比，用以下推導得到的公式（6）可以算得采集到的物料的C_p。

Q_p＝Q_pw                                （3）

Q_p＝C_p(T₃－T₁)                          （4）

量取物料被加熱前的溫度（T₁），被加熱一段時間（t）后量取其即時溫度記錄為T₃；水被加熱前的溫度為T₂，同理水被加熱一段時間t后量取其即時溫度記為T₄，則有：

Q_pw＝C_pw(T₃－T₁)                       （5）

聯立（3）、（4）和（5）得到公式（6）。

C_p＝C_pw(T₄－T₂)/(T₃－T₁)                 （6）

為控制水和物料都得到相同的熱量，實驗設計水和物料均被電陶爐（選擇70 W）加熱240 s（t＝240 s）。

另外，根據傳統(tǒng)傳熱學中導熱率的定義，一個體系的導熱率（k）如下^[18]。

k＝QΔx/AdT                             （7）

同理可推得，同一組實驗已經稱得的相同質量的水和物料都得到相同的輸入熱量（Q）、加熱接觸面積（A）和熱量傳遞距離（Δx），物料的k_f可用已知的水的導熱率（k_w）根據公式（8）計算得到。

k_f＝k_w(T₄－T₂)/(T₃－T₁)                   （8）

k_w選擇0.640 6W/(m·K)，水在50 ℃的導熱率

2.3.6  滲透壓（π_f）  按照YASN OM819型滲透壓儀（摩爾濃度滲透壓儀，摩爾滲透壓儀）操作手冊執(zhí)行，由揚子江藥業(yè)集團南京海陵中藥制藥工藝技術國家工程研究中心檢測。每次測量樣品需10 mL，因此，每組樣品收集30mL，共測量3組數據，取平均值，得到滲透濃度b₀（mOsm/kg）。根據公式（9）算得滲透壓（摩爾濃度滲透壓，摩爾滲透壓）（kPa）。

π_f＝b₀RT                                （9）

R為摩爾氣體常數8.314 J/(K?mol)，T為開氏溫度，所有數據在25 ℃（298.15K）時測量得到

3  結果與分析

3.1  白利度與固含率相關性研究

白利度^[19-20]最早被定義為溶液中蔗糖的質量分數，后被廣義定義為與蔗糖質量分數相等的溶液中可溶性固形物質量分數，被廣泛運用在濃縮行業(yè)（如果汁濃縮和涼茶濃縮）的物料濃度。白利度一般用折射儀測得折射指數，并換算成°Bx或%。其中，可溶性固形物也包括溶解性物質與膠態(tài)物質。溶液中固含率定義一般為固體成分質量分數，包括了可溶性與不可溶固形物。因此一個系統(tǒng)里的白利度與固含率不一定相同；一般可認為固含率要比白利度高。本實驗中，所得到的同一時間點的白利度與固含率相近（表1），而測得的白利度基本比固含率低，并且呈現出較好的線性關系，回歸方程為y＝0.011 5 x＋0.001 7，R²＝0.998 8，y代表固含率，x代表白利度。因此，在玉屏風散的濃縮過程研究中，可以利用白利度去預測玉屏風散的固含率。除了經典復方玉屏風散在濃縮過程中具有白利度與固含率的上述關系外，此前Criscuoli等^[11]也考察了棗汁的膜蒸餾濃縮過程中白利度與固含率之間的關系，得出類似結果。

需要特別指出的是，中藥制藥工藝往往涉及澄清、精制工段，淀粉、蛋白質等大分子雜質在濃縮工段之前已被去除。因此中藥物料中一些大分子物質可能因為加熱等問題引起的析出影響可被忽略。

3.2  濃縮過程中中藥物料溶液環(huán)境表征參數的相關性研究

3.2.1  白利度與密度、黏度的關系  濃縮過程中同一時間點玉屏風散物料的溶液環(huán)境特征參數白利度與密度、黏度關系的實測數據分別如表2、3所示。

根據密度的基本定義，物質的密度可以用質量除以體積算得。中藥水提液物料的濃縮過程，可認為是一個去除水，留下固體溶質的過程。一般認為固體溶質的密度比水密度大，因此，除水過程中物料的密度應呈上升趨勢。玉屏風散濃縮過程中，物料的白利度與密度呈現優(yōu)越的線性相關，回歸方程為y＝4.266 1 x＋1 072.3，R²＝0.993 2，y代表密度，x代表白利度。

黏度是流體的一種特性，客觀地表達了流體對流動的阻力，在中藥水提液濃縮體系中則反映了水分子在運動過程中的受阻程度。白利度與黏度是2個不同的物理概念，沒有直接的關聯。本實驗所得到的白利度與黏度的關系方程為y＝3.361 9 e^15.593 x，R²＝0.933 1，y代表黏度，x代表白利度。根據二者的關系方程，雖不能由白利度準確計算物料體系的黏度，但是能粗略地估算玉屏風散物料的黏度。中藥液體物料黏度的數值主要是用來計算流體體系的雷諾數，以預測流體的分類（湍流、層流），為膜蒸餾濃縮中藥物料的即時能耗動態(tài)計算提供一個重要的選取依據。

3.2.2  白利度與比熱容、導熱率的關系  從物質比熱容定義及公式（2）可知，某體系的比熱容是指在該體系中單位質量的物質改變單位溫度（1 K）時需要吸收或釋放的熱量。而在本研究中的比熱容，特指定壓比熱容（在同一壓強下）。

一般認為一個復雜體系中，比熱容與混合物質質量呈一定比例。因此在中藥溶液環(huán)境的復雜體系中，水的去除應該使物料比熱容的大小有所改變。濃縮過程中同一時間點的溶液環(huán)境特征參數白利度與比熱容、導熱率的實測數據分別見表4、5。白利度與玉屏風散溶液環(huán)境中比熱容的關系方程為y＝0.044 x＋4.039 7，R²＝0.987 9，y代表比熱容，x代表白利度；在實驗初始白利度為2.6°Bx時，溶液的比熱容與水的比熱容非常接近［4.18 kJ/(kg?K)］，而比熱容隨著水分的減少而增加。比熱容是能耗測算公式的關鍵參數，但因為中藥水體液物料體系是一種高度多維的復雜系統(tǒng)，其比熱容難于測定。因此利用白利度去預測不同濃度溶液下的比熱容，對濃縮能耗計算具有填補空白的突破性意義。

導熱率（又稱導熱系數）是一個形容系統(tǒng)熱傳導能力的物理量，由其計算公式（7）可知，導熱率即為在相同的熱量傳導條件下，溫差為1 K的單位厚度的系統(tǒng)在傳熱為單位面積單位時間條件下的傳熱量，導熱率與系統(tǒng)的成分組成、狀態(tài)等有關系。根據公式（8）可推知，本研究測試的玉屏風散系統(tǒng)中導熱率與系統(tǒng)物料溫差成正比。利用白利度與導熱率的關系方程（y＝0.006 7x＋0.619 1，R²＝0.9879，y代表導熱率，x代表白利度），可以成功地預測出玉屏風散在不同濃度溶液下的導熱率。

3.2.3  白利度與滲透壓的關系  濃縮過程中同一時間點的溶液環(huán)境特征參數白利度與滲透壓（摩爾濃度滲透壓，摩爾滲透壓）實測數據見表6；白利度與滲透壓的關系方程為y＝160.54   x－366.66，R²＝0.965 0，y代表滲透壓（摩爾濃度滲透壓，摩爾滲透壓），x代表白利度。

滲透壓（摩爾濃度滲透壓，摩爾滲透壓）是與單位體積溶劑中溶質微粒的濃度有關，可以理解一個溶液系統(tǒng)的滲透壓（摩爾濃度滲透壓，摩爾滲透壓）是由溶解性溶質產生。因此在玉屏風散的物料體系中，可以觀察到白利度與滲透壓（摩爾濃度滲透壓，摩爾滲透壓）呈現良好的線性關系，而白利度表示一個溶液體系的溶解性固體物質量分數。溶解性固體物質量分數增加，滲透壓（摩爾濃度滲透壓，摩爾滲透壓）也隨著增加，符合理論期待值。實驗中觀察到在白利度高于10°Bx時，測得的滲透壓（摩爾濃度滲透壓，摩爾滲透壓）比線性方程預測的數值稍微低一些，這可能是因為物料中有些物質經過一段時間（樣品由廣州寄往南京期間）儲存后沉淀析出，導致了測得的滲透壓（摩爾濃度滲透壓，摩爾滲透壓）比樣品收集時的滲透壓（摩爾濃度滲透壓，摩爾滲透壓）稍低。利用液體物料白利度可以較好地推算中藥物料的瞬時滲透壓（摩爾濃度滲透壓，摩爾滲透壓），對反滲透膜濃縮在中藥物料濃縮中的瞬時能耗變化和總能耗計算有重要意義。

3.3  基于本研究結果的熱法濃縮及膜蒸餾濃縮的能耗計算

3.3.1  熱法濃縮能耗計算相關性  如熱法單效蒸發(fā)濃縮中藥物料，需要大量蒸汽或者電轉化成熱能，因此，計算單效蒸發(fā)最低所需熱能（Q_T）的公式如下^[3]。

Q_T＝Q_v＋Q_L                            （10）

Q_v＝M_v(ΔH_v－ΔH_f)                      （11）

Q_L＝M_LC_pΔT                           （12）

式（10）中，總能耗Q_T由Q_L與Q_v組成，分別代表物料升溫需要的能量以及蒸發(fā)需要的能量。M_v代表水揮發(fā)速率，ΔH_v代表水蒸氣在此壓強、溫度下的焓，ΔH_f代表飽和水蒸氣在此壓強、溫度下的焓。物料守恒定律，可以算得濃縮物料的水揮發(fā)速率，見式（13）。

M_v＝M_L(1－x)                          （13）

式（13）中，M_L代表物料進入單效濃縮器中的速率，x代表物料固含率。其中，物料固含率會隨著物料濃縮過程變化。此時，根據測得物料的白利度讀數，可用“3.1”項得到的白利度與固含率關系模型得到某一時刻的固含率，從而算出M_v。

另外，式（12）中C_p（物料的比熱容）會隨著中藥物料的濃縮過程而變化^[3]。因此Q_L也會隨濃縮過程中，中藥物料濃度升高而變化。根據“3.2”項中白利度與熱容率的關系，可計算出熱法單效濃縮的能耗。

3.3.2  膜蒸餾能耗計算相關性  同理，膜蒸餾濃縮中藥物料的能耗也與中藥濃縮物料特征有關系。膜蒸餾MD理論最低所需能耗（Q_f）計算公式如下^[21]。

Q_f＝A_fh_f (T_f,b－T_f,m)                     （14）

式（14）中，A_f代表膜接觸物料端面積，h_f代表中藥物料端熱對流導熱系數，T_f,b代表中藥物料端溫度以及T_f,m代表中藥物料端膜表面溫度。其中，膜蒸餾物料端膜表面的溫度T_f,m不能被直接測出來，因此，需要根據物料衡算和傳熱傳質模型的推導。膜表面溫度T_f,m可用以下公式進行計算^[22]。

T_f,m＝[k_m(T_p,b＋h_fT_f,b/h_p)/δ＋h_fT_f,b－JΔH_latent]/[k_m(1＋h_f/h_p)/δ＋h_f]                                （15）

式（15）中，k_m表示膜的導熱系數，δ表示膜的厚度，T_p,b代表純水透過液端的溫度，T_f,b代表中藥濃縮物料端的溫度，h_f代表物料端的熱對流導熱系數，h_p代表純水透過液端的熱對流導熱系數。J代表膜通量，以及ΔH_latent代表飽和液體水蒸發(fā)成水蒸氣的汽化潛熱。

其中，物料端的膜表面表面溫度也需要得到物料的熱對流導熱系數h_f。然而根據經典流體力學和傳熱理論，計算熱對流的導熱系數h_f需要物料的比熱容C_p、密度ρ、黏度μ和物料導熱率k_f等濃縮過程變化的物理特征，涉及公式（16）～（20）等復雜、繁瑣的較多計算過程如下^[22]。

根據（17），首先利用濃縮物料膜端的水力直徑d_h、物料端流速ν、物料的黏度μ以及物料的密度ρ算出雷諾數Re。再用（18）與物料的比熱容C_p、物料的粘度μ與物料導熱率k_f算出普朗特數Pr。根據Re與Pr數，利用（19）或（20）算出新的努塞爾數Nu_f。利用Nu_f代入（16）中，算出新的物料的熱對流的導熱系數h_f。

h_f＝Nu_fk_f/d_h                            （16）

Re＝d_hνρ/μ                             （17）

Rr＝μC_p/k_f                             （18）

Nu_f＝1.62(RePrd_h/L)^1/3Re＜2300           （19）

Nu_f＝0.023Re^4/5Pr^1/3Re＞2300             （20）

利用“3.2”項得出白利度與其他物理參數如C_p、k_f等關系模型，計算膜蒸餾濃縮中藥的能耗。

3.3.3  反滲透能耗相關性  反滲透濃縮中藥物料的能耗的計算必須根據物料濃縮過程中滲透壓（摩爾濃度滲透壓，摩爾滲透壓）的變化來計算反滲透技術最低所需能耗（W_min）計算公式如下^[22]。

W_min＝V₀ΔπdR                       （21）

V₀代表中藥物料初始體積，R代表濃縮率，Δπ代表物料端與滲透液端滲透壓差

Δπ＝π_f－π_p≈π_f                         （22）

π_f為物料滲透壓，π_p為滲透液端（清水）滲透壓

物料的滲透壓（π_f）會隨著物料濃縮的濃度上升而上升。根據白利度與滲透壓（摩爾濃度滲透壓，摩爾滲透壓）的關系方程，可利用測得的白利度得到某時反滲透濃縮中藥物料的瞬時能耗，從而計算出整個濃縮過程的總能耗。

4 討論

4.1  有關溶液環(huán)境特征參數的討論

本實驗所發(fā)現的中藥水提液濃縮過程中若干“溶液環(huán)境”特征參數的相關性，有關文獻亦有過相關報道^[9,17]。這是一種偶然現象，還是具有普遍意義的規(guī)律，尚有待進一步深入系統(tǒng)的研究。中藥溶液環(huán)境的黏度、pH值、離子強度等表征參數，既來源于中藥水提液中各種物質的化學組成，又是水提液體系中各種物質不同表現的綜合反饋，當然也必定與制藥分離過程物料體系中各種物質的熱力學和傳遞性質、多相流、多組分傳質等影響分離的因素密切相關^[23]。本團隊將在進一步開展系統(tǒng)實驗研究的基礎上，借助現代分離科學與物理化學等多學科交叉研究的方法，從理論上加以詮釋。

4.2  化工過程控制原理是探討中藥制藥過程工程原理的重要基石

化工過程控制原理提出：對于復雜的化工過程，不能滿足于在現有裝置上通過測試獲得所需的對象動態(tài)特性知識，更重要的是在對象處于設計階段，就能利用計算方法預估其特性，以改變黑匣無知狀態(tài)，指導工藝設計中的原理性設計。中藥水提液由幾十甚至上百種小分子化合物（包括揮發(fā)油、氨基酸及生物堿、有機酸、黃酮類、皂苷等化學成分）和生物大分子物質（包括肽、蛋白、糖肽及多糖等）組成，是一種典型的高度復雜的大系統(tǒng)客體，因此，十分有必要引入現代自然科學化工過程控制的方法對其進行研究。節(jié)能減排是中藥綠色制造的核心目標，智能制造則是中藥現代化的重要內容。而傳統(tǒng)的濃縮能耗計算公式涉及多達近10個物料理化參數和濃縮過程工藝參數，其中部分參數采集方法繁瑣、缺乏精準性，且無法實現在線檢測，加上繁瑣的推算過程，難以在大生產中推廣應用。

本研究將化工過程控制基本原理引入中藥生產工藝，結合數學建模的基本方法，將原型的某一部分信息簡化、壓縮、提煉而構造成的原型替代物，在開展中藥制劑、分析化學、現代分離科學及計算機技術等多學科聯合攻關的基礎上，對中藥制藥濃縮過程進行深入系統(tǒng)的基礎研究。

本研究提出的能耗方法，僅需對某一常年生產的大品種建立以白利度為核心參數，與物料黏度、滲透壓（摩爾濃度滲透壓，摩爾滲透壓）、比熱容、密度和導熱率等參數的相關性模型，采用有關算法編制相關軟件，即可根據實時采集的白利度數據，得出瞬時能耗及同一時刻的物料溶液環(huán)境各相關特征參數，為膜濃縮這一綠色制造先進技術在大生產中的推廣應用，也為中藥生產過程實現在線監(jiān)測、智能控制提供技術支撐。

利益沖突  所有作者均聲明不存在利益沖突

參考文獻（略） 

來  源：鐘文蔚，黎萬鈺，丁  菲，李除夕，鄭東陽，郭立瑋．基于中藥水提液濃縮過程溶液環(huán)境特征參數相關性的瞬時能耗計算方法探索——以玉屏風散水提液為例  [J]. 中草藥, 2021, 52(7): 1937-1944 .