張永超，徐芳芳，張 欣，徐 冰，趙月梅，朱文燈，吳 云，王振中，肖 偉,

腰痹通膠囊4種中間體粒徑的近紅外光譜通用定量模型研究

張永超1，徐芳芳2, 3, 4*，張 欣2, 3, 4，徐 冰5，趙月梅2, 3, 4，朱文燈2, 3, 4，吳 云2, 3, 4，王振中2, 3, 4，肖 偉1, 2, 3, 4*

1. 南京中醫(yī)藥大學(xué)，江蘇 南京 210000 2. 江蘇康緣藥業(yè)股份有限公司，江蘇 連云港 222001 3. 中藥制藥過程新技術(shù)國家重點(diǎn)實(shí)驗(yàn)室，江蘇 連云港 222001 4. 中成藥智能制造國家地方聯(lián)合工程研究中心，江蘇 連云港 222001 5. 北京中醫(yī)藥大學(xué)，北京 100029

旨在通過近紅外光譜（NIRS）信息與腰痹通膠囊（Yaobitong Capsules，YC）中間體中值粒徑（50）的關(guān)系分析，探討對(duì)該品種生產(chǎn)過程中4種中間體建立50近紅外通用定量模型的可行性。采集YC生產(chǎn)過程中的原料細(xì)粉、干燥顆粒、整粒顆粒和總混顆粒4種中間體的NIRS，考察不同預(yù)處理方法對(duì)模型的影響，并采用間隔偏最小二乘法（iPLS）、組合間隔偏最小二乘法（siPLS）和移動(dòng)窗口偏最小二乘法（mwPLS）優(yōu)選NIRS波段，采用偏最小二乘法（PLS）對(duì)4種中間體建立1個(gè)50通用定量模型。通用模型的交叉驗(yàn)證均方根誤差（RMSECV）為3.918 μm，預(yù)測(cè)均方根誤差（RMSEP）為2.832 μm，預(yù)測(cè)相對(duì)偏差（RSEP）為2.26%，小于5%，性能偏差比（RPD）為4.60，大于3，該模型可以用于定量測(cè)定，且與獨(dú)立模型比，預(yù)測(cè)準(zhǔn)確性相差不大。NIRS通用定量模型可用于YC 4種中間體50的測(cè)定。

腰痹通膠囊；近紅外光譜技術(shù)；偏最小二乘法；粒徑；中值粒徑；通用模型

腰痹通膠囊（Yaobitong Capsules，YC）是由三七、川芎、延胡索、白芍、牛膝、狗脊、熟大黃、獨(dú)活8味中藥組成，具有活血化瘀、祛風(fēng)除濕的功效，主要用于治療腰椎間盤突出癥[1]。顆粒粒徑是評(píng)價(jià)制粒所得顆粒是否合格的關(guān)鍵指標(biāo)[2-3]，并對(duì)多個(gè)關(guān)鍵質(zhì)量屬性有影響，例如吸濕穩(wěn)定性[4]、流動(dòng)性[5]、可壓縮性[6]、溶出度[7]，它還影響制粒工藝參數(shù)的選擇[8]及膠囊的填充[9]。在實(shí)際生產(chǎn)應(yīng)用中，中值粒徑（50）是評(píng)價(jià)顆粒粒徑的主要指標(biāo)，50是指粒徑分布累積達(dá)到50%時(shí)對(duì)應(yīng)的粒徑值[10-11]，它能較準(zhǔn)確的反映樣本的粒徑情況。目前，YC正在進(jìn)行工藝優(yōu)化提升質(zhì)量標(biāo)準(zhǔn)，50擬納入質(zhì)量?jī)?nèi)控標(biāo)準(zhǔn)。

近紅外光譜（near infrared spectrum，NIRS）技術(shù)與化學(xué)計(jì)量法結(jié)合，能夠?qū)崿F(xiàn)快速檢測(cè)化學(xué)成分含量及物理性質(zhì)指標(biāo)，在藥材鑒別與化學(xué)成分含量快速測(cè)定中應(yīng)用較多[12?14]，但較少用來快速檢測(cè)物性指標(biāo)。粒徑信息可以通過NIRS基線的變化反映出來[15]，已有研究者將NIRS技術(shù)用于粒徑的快速檢測(cè)。Rosas等[16]以實(shí)驗(yàn)室自制尼美舒利為研究對(duì)象，成功開發(fā)出NIRS快速檢測(cè)顆粒粒徑分布（＜125 μm、125～250 μm、＞250 μm）的模型，模型的預(yù)測(cè)均方根誤差（root mean square error of prediction，RMSEP）為2.5%。Otsuka等[17]采用NIRS技術(shù)實(shí)時(shí)監(jiān)控實(shí)驗(yàn)室規(guī)模的撲熱息痛散流化床制粒過程，建立了顆粒50的快速檢測(cè)方法，模型的校正集相關(guān)系數(shù)達(dá)0.961。Pauli等[18]開發(fā)了3個(gè)在線NIRS定量分析方法，用于測(cè)定雙氯芬酸鈉濕法制粒和流化床干燥過程中的顆粒粒徑分布，預(yù)測(cè)值在允許的誤差范圍內(nèi)。潘曉寧[19]以黨參與玄參藥材為研究對(duì)象，分別構(gòu)建了粒徑分布累積為10%、50%、90%對(duì)應(yīng)的粒徑值（10、50、90）的NIRS定量預(yù)測(cè)模型。上述研究表明，NIRS技術(shù)結(jié)合化學(xué)計(jì)量法是可以測(cè)定藥物顆粒粒徑的，但多數(shù)研究都是針對(duì)某一類樣本進(jìn)行分析，且都集中于化學(xué)藥領(lǐng)域。通用模型是指針對(duì)某個(gè)指標(biāo)建立一個(gè)模型，可以用于分析2種及以上的樣本，通用性強(qiáng)，可以節(jié)約較多成本。建立一個(gè)穩(wěn)定、可靠的NIRS定量預(yù)測(cè)模型是一項(xiàng)較為復(fù)雜的工作，耗時(shí)長(zhǎng)且成本高。如果能針對(duì)兩種及以上的樣本建立一個(gè)通用的模型，將會(huì)提升檢測(cè)效率，節(jié)約成本。本研究以YC生產(chǎn)過程中的4種中間體為研究對(duì)象，通過采集4種中間體的近紅外漫反射光譜，以中間體50為參考值，采用偏最小二乘法（partial least squares，PLS）建立一個(gè)檢測(cè)4種中間體50的通用方法。

1 儀器與材料

1.1 儀器

Antaris II型傅里葉近紅外變換光譜儀，配有積分球漫反射采樣系統(tǒng)、Result光譜采集軟件，美國Thermo公司；BT-2600激光粒度分布儀，丹東百特儀器有限公司。

1.2 材料

YC生產(chǎn)過程中的中間體，包含原料細(xì)粉、干燥顆粒、整粒顆粒和總混顆粒，由江蘇康緣藥業(yè)股份有限公司提供。本研究收集的樣品為2019年8月至2020年5月生產(chǎn)的中間體樣品，收集到的中間體樣品批次為Z190801～Z190820、Z190901～Z190918、Z200301～Z200320、Z200401～Z200420、Z200501～Z200502，包含4種中間體共320個(gè)樣品。

2 方法

2.1 NIRS采集

取4 g左右樣品，置于配備的樣品杯里，輕輕壓實(shí)，采用積分球漫反射方式采集NIRS。掃描范圍為10 000～4000 cm?1，分辨率為8 cm?1，2倍增益，掃描次數(shù)64次，以空氣為背景，每小時(shí)掃描1次背景。每個(gè)樣品掃描3次，平均值用于分析。

2.2 參考值測(cè)定

激光粒度分布儀法測(cè)粒徑代表性強(qiáng)且適用于粒徑小的顆粒[10]，“1.2”項(xiàng)下的樣本粒徑小于550 μm，采用該法測(cè)粒徑結(jié)果較為準(zhǔn)確。取適量樣品置于進(jìn)料斗中，使用BT-2600激光粒度分布儀測(cè)定樣品的粒徑分布，并計(jì)算50（粒徑分布累積為50%時(shí)對(duì)應(yīng)的粒徑值）。參數(shù)設(shè)置：分散介質(zhì)為空氣，空氣壓力為0.25 MPa，遮光率為3%～12%。

2.3 光譜預(yù)處理方法

采集光譜時(shí)，由于受到環(huán)境溫濕度、儀器狀態(tài)以及顆粒狀態(tài)等各種因素的影響，NIRS會(huì)產(chǎn)生噪聲信號(hào)、基線漂移等，獲取的光譜中包含了自身的信息以及其他不必要的信號(hào)。對(duì)NIRS進(jìn)行合適的預(yù)處理，可以減少噪聲，濾過無關(guān)信息，提高模型穩(wěn)健性。常見的預(yù)處理方法有矢量歸一化法、導(dǎo)數(shù)法（一階導(dǎo)數(shù)、二階導(dǎo)數(shù)）、平滑法［Savitzky-Golay（SG）平滑、Norris Derivative（ND）平滑］、多元散射校正（multiplicative scatter correction，MSC）、標(biāo)準(zhǔn)正態(tài)變量變換法（standard normal variate transformation，SNV）、基線校正以及上述幾種方法的結(jié)合。矢量歸一化法能增強(qiáng)光譜差異，可以校正由光程或樣品稀釋等導(dǎo)致的光譜變化；MSC和SNV可以消除顆粒大小、分布不勻帶來的干擾；導(dǎo)數(shù)法可以消除基線漂移；平滑法能夠?yàn)V除噪聲[20-21]。

2.4 數(shù)據(jù)處理與評(píng)價(jià)方法

采用Unscramble X 10.4（Camo software AS，Norway）軟件對(duì)NIRS進(jìn)行預(yù)處理，采用Matlab 2016a（Mathwork Inc.，USA）軟件對(duì)樣本進(jìn)行變量篩選與建模，采用留一交叉驗(yàn)證法以交叉驗(yàn)證均方根誤差（root mean square errors of cross validation，RMSECV）為評(píng)價(jià)指標(biāo)確定主因子數(shù)。以中間體50為因變量，以對(duì)應(yīng)的NIRS值為自變量建立PLS定量模型。本研究以下列參數(shù)來評(píng)價(jià)模型性能，以優(yōu)選最佳模型。校正集相關(guān)系數(shù)（cal）、驗(yàn)證集相關(guān)系數(shù)（pre）分別表示校正模型與驗(yàn)證模型的擬合程度，cal和pre越大，模型擬合效果越好；校正均方根誤差（root mean square errors of calibration，RMSEC）、RMSEP分別指校正模型與驗(yàn)證模型中參考值與預(yù)測(cè)值之間的均方根誤差，越小模型預(yù)測(cè)性能越高；校正集偏差（BIAScal）和驗(yàn)證集偏差（BIASpre）分別表示校正模型與驗(yàn)證模型中參考值與預(yù)測(cè)值之間的偏差，其越小，模型預(yù)測(cè)精度越 高[22-23]。預(yù)測(cè)相對(duì)偏差（relative standard error of prediction，RSEP）是模型中參考值與預(yù)測(cè)值之間的相對(duì)偏差，一般RSEP越小模型預(yù)測(cè)性能越好；性能偏差比（ratio of performance to deviation，RPD）表示模型預(yù)測(cè)性能，是標(biāo)準(zhǔn)差（SD）與RMSEP的比值，當(dāng)RPD＞3時(shí)，表示模型預(yù)測(cè)精度高[24]。各項(xiàng)評(píng)價(jià)參數(shù)不能孤立參考，需要綜合起來評(píng)價(jià)。

3 結(jié)果與分析

3.1 樣本劃分

采用隨機(jī)抽樣（random sampling，RS）法劃分校正集與驗(yàn)證集，并保證驗(yàn)證集中參考值范圍包含于校正集中[12]。每種中間體按照3∶1的比例劃分校正集與驗(yàn)證集，隨機(jī)選取60個(gè)樣品作為校正集，20個(gè)樣品作為驗(yàn)證集；通用模型的校正集為4種中間體校正集的總和，共240個(gè)樣品，驗(yàn)證集為4種中間體驗(yàn)證集的總和，共80個(gè)樣品，劃分結(jié)果見表1。

表1 樣品校正集與驗(yàn)證集參考方法測(cè)定值的統(tǒng)計(jì)結(jié)果

3.2 光譜預(yù)處理方法的選擇

YC中間體的近紅外原始光譜見圖1。本研究考察了以下預(yù)處理方法：移動(dòng)窗口平滑（9點(diǎn)）、SNV、S-G 1st（9點(diǎn)）、基線校正、MSC、歸一化法。使用上述方法對(duì)光譜進(jìn)行預(yù)處理后，建立50的PLS定量模型，不同預(yù)處理方法對(duì)模型性能的影響如表2所示。以RPD與RSEP為評(píng)價(jià)標(biāo)準(zhǔn)，篩選最優(yōu)的預(yù)處理方法。由表2可知，原料細(xì)粉模型中，一階求導(dǎo)結(jié)合S-G平滑對(duì)光譜進(jìn)行預(yù)處理建模效果最優(yōu)，RPD為4.50，RSEP為2.30%；干燥顆粒模型中，采用移動(dòng)窗口平滑（9點(diǎn)）預(yù)處理方法最佳，RPD為4.12，RSEP為1.88%；整粒顆粒模型里，采用原始光譜建模比進(jìn)行預(yù)處理后建模效果好，RPD為4.15，RSEP為1.80%；總混顆粒模型采用移動(dòng)窗口平滑（9點(diǎn)）方法對(duì)光譜進(jìn)行預(yù)處理后建模效果最好，RPD為4.84，RSEP為1.54%；通用模型中，采用原始光譜建模性能最佳，RPD為4.60，RSEP為2.26%。

圖1 YC中間體近紅外原始光譜

表2 不同預(yù)處理方法對(duì)模型的影響

3.3 特征變量篩選

篩選特征變量可以剔除無關(guān)信息，提高模型性能。本研究在上述篩選出的最佳預(yù)處理方法基礎(chǔ)上進(jìn)一步篩選特征變量。

3.3.1 基于間隔偏最小二乘法（interval PLS，iPLS）篩選變量 iPLS是把全光譜劃分成若干個(gè)子區(qū)間，然后在每個(gè)子區(qū)間進(jìn)行建模[25]。本研究將全光譜劃分成20個(gè)子區(qū)間，以RMSECV為評(píng)價(jià)指標(biāo)，篩選最佳建模波段。

3.3.2 基于組合間隔偏最小二乘法（synergy interval PLS，siPLS）篩選變量 siPLS是基于iPLS，它是將全光譜劃分成若干個(gè)子區(qū)間后，再把子區(qū)間任意組合起來建模。本研究是把全光譜劃分成20個(gè)子區(qū)間，再以子區(qū)間組合數(shù)為4建立模型，以RMSECV為評(píng)價(jià)指標(biāo)，篩選最佳建模波段。

3.3.3 基于移動(dòng)窗口偏最小二乘法（moving window PLS，mwPLS）篩選變量 mwPLS是從整個(gè)光譜的第一個(gè)波長(zhǎng)點(diǎn)開始移動(dòng)，沿波長(zhǎng)變化的方向截取選定窗口寬度的區(qū)間，建立一系列的PLS模型[26]。本研究以初始窗口寬度為31，以10為步長(zhǎng)依次增加窗口寬度，建立了窗口寬度為31～311的PLS模型，并根據(jù)RMSECV選取最佳建模波段。

3.3.4 篩選特征變量 采用3種方法建立模型的性能參數(shù)如表3～7所示。分別以RMSEC、RMSECV、RMSEP為評(píng)價(jià)指標(biāo)，綜合評(píng)價(jià)篩選出建模的最佳波段。

原料細(xì)粉模型采用全光譜建模與mwPLS優(yōu)選的波段建模性能相差不大，但是篩選波段后變量減少至115個(gè)，建模時(shí)間會(huì)縮短，因此選擇建模波段3 999.64～4 018.92、5 739.12～6 136.38 cm?1。

干燥顆粒模型采用siPLS和mwPLS篩選變量后，模型性能均有所提高，變量數(shù)分別減少到392、278個(gè)，但mwPLS篩選變量后的波點(diǎn)數(shù)更少且cal、pre較大，RMSEC、RMSEP較小，故采用波段4 030.50～4 219.49、4 443.19～4 504.90、6 957.91～7 767.86 cm?1建模。

表3 原料細(xì)粉優(yōu)選區(qū)間建模與全光譜建模的比較

表4 干燥顆粒優(yōu)選區(qū)間建模與全光譜建模的比較

表5 整粒顆粒優(yōu)選區(qū)間建模與全光譜建模的比較

表6 總混顆粒優(yōu)選區(qū)間建模與全光譜建模的比較

表7 通用模型優(yōu)選區(qū)間建模與全光譜建模的比較

整粒顆粒模型采用siPLS與mwPLS篩選變量模型性能有所提升，變量數(shù)分別減少至312、83個(gè)，但mwPLS的cal、pre更大，RMSEC、RMSEP也較小，且波點(diǎn)數(shù)也較少，故選用波段4 933.02～ 4 971.59、5 982.10～6 213.52、7 359.03～7 397.60 cm?1建模。

總混顆粒模型采用mwPLS法進(jìn)行波段篩選后，cal、pre均變大，RMSEC、RMSEP均變小，模型性能提高，因此選用波段4 015.07～42 342.6、 4 396.90～4 524.18、7 667.58～7 698.44、8 369.54～8 454.40 cm?1建模。

通用模型篩選變量后雖然降低了模型復(fù)雜程度，但是模型性能卻不及全光譜模型，可能是因?yàn)橥ㄓ媚Ｐ褪前?個(gè)中間體信息的共用模型，采用波段篩選可能會(huì)使重要信息的流失，造成模型性能降低，故通用模型采用全波段建模。

3.4 主因子數(shù)的確定

以RMSECV值為評(píng)價(jià)指標(biāo)確定主因子數(shù)，主因子數(shù)選擇不當(dāng)會(huì)導(dǎo)致擬合不充分或擬合過度，一般以RMSECV值最小時(shí)對(duì)應(yīng)的主因子數(shù)為最佳。本研究采用留一交叉驗(yàn)證法，以RMSECV為指標(biāo)，考察主因子數(shù)對(duì)模型的影響。最終確定原料細(xì)粉模型、干燥顆粒模型、整粒顆粒模型、總混顆粒模型和通用模型的最佳主因子數(shù)分別為6、15、12、11、14。

3.5 D50定量預(yù)測(cè)模型的建立

采用篩選出來的最佳預(yù)處理方法及最優(yōu)NIRS波段建立50近紅外定量預(yù)測(cè)模型，見圖2。優(yōu)選出來的原料細(xì)粉、干燥顆粒、整粒顆粒、總混顆粒和通用模型的評(píng)價(jià)參數(shù)見表8。各項(xiàng)模型的cal、pre接近于1，說明參考值與預(yù)測(cè)值相關(guān)性較高；RMSEC、RMSECV、RMSEP、RSEP較小，RPD大于3，說明模型預(yù)測(cè)性能較高，可以用于50的定量預(yù)測(cè)。

3.6 模型驗(yàn)證與評(píng)價(jià)

將驗(yàn)證集樣品的NIRS導(dǎo)入校正模型中，預(yù)測(cè)中間體的50，并與參考值進(jìn)行比較，結(jié)果見表9～10。原料細(xì)粉模型、干燥顆粒模型、整粒顆粒模型、總混顆粒模型及通用模型50的平均相對(duì)預(yù)測(cè)誤差分別為1.94%、1.43%、1.46%、1.21%、1.87%，均小于5%，說明2種模型的預(yù)測(cè)準(zhǔn)確性均較高。由表11可知，獨(dú)立模型與通用模型對(duì)各中間體50的預(yù)測(cè)誤差的差值均小于1%，說明二者的預(yù)測(cè)性能相近。

對(duì)原料細(xì)粉、干燥顆粒、整粒顆粒、總混顆粒和通用模型的驗(yàn)證集中參考值與預(yù)測(cè)值進(jìn)行配對(duì)檢驗(yàn)，結(jié)果值依次為0.22、0.17、?0.42、?0.38、0.75，值依次為0.825、0.866、0.678、0.706、0.454，值均＞0.05，說明參考值與預(yù)測(cè)值之間無顯著性差異。

表8 最佳模型的評(píng)價(jià)參數(shù)

表9 NIRS參考值與預(yù)測(cè)值的比較(獨(dú)立模型)

表10 NIRS參考值與預(yù)測(cè)值的比較(通用模型)

表11 獨(dú)立模型與通用模型相對(duì)預(yù)測(cè)誤差比較

4 討論

本研究以YC為研究對(duì)象，收集其生產(chǎn)過程中的原料細(xì)粉、干燥顆粒、整粒顆粒及總混顆粒共4種中間體，采用NIRS技術(shù)結(jié)合PLS建立了預(yù)測(cè)中間體50的通用定量模型。各項(xiàng)評(píng)價(jià)指標(biāo)表明，通用定量模型預(yù)測(cè)準(zhǔn)確性高、能滿足實(shí)際需求，且與單一模型比預(yù)測(cè)性能差異不大，可以實(shí)現(xiàn)對(duì)4種中間體50的測(cè)定。

本研究中的通用模型是基于一個(gè)品種中差異不大的中間體建立起來的，通用模型預(yù)測(cè)效果較好，可能是因?yàn)橹虚g體的50有一定差異但還在可接受范圍內(nèi)，且由于這種差異使得50的分布范圍變寬，更利于建模，但對(duì)不同品種建立相關(guān)指標(biāo)的通用模型的可行性還有待探索。相比激光粒度分布儀法，NIRS法測(cè)粒徑會(huì)更快，若應(yīng)用于在線研究將無需取樣，提高分析效率。目前，NIRS通用模型在制藥領(lǐng)域研究較少，若將NIRS技術(shù)與通用模型結(jié)合起來，建立的NIRS通用模型將具備快速準(zhǔn)確、無損綠色、節(jié)約成本等優(yōu)點(diǎn)，可以為中藥生產(chǎn)過程實(shí)現(xiàn)在線檢測(cè)相關(guān)指標(biāo)提供更好的方法。

本研究是初步證明了YC 4種中間體近紅外50通用定量模型的可行性，后續(xù)完善該模型還應(yīng)當(dāng)加入更多的樣本，以及后續(xù)生產(chǎn)的持續(xù)驗(yàn)證，并朝著在線監(jiān)測(cè)方向上研究。

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

[1] 中國藥典 [S]. 一部. 2015.

[2] 臧恒昌, 臧立軒, 張惠, 等. 近紅外光譜分析技術(shù)在制藥領(lǐng)域中的應(yīng)用研究進(jìn)展 [J]. 藥學(xué)研究, 2014, 33(3): 125-128.

[3] 姜爽, 楊海龍, 臧恒昌. 近紅外光譜分析技術(shù)在中藥生產(chǎn)過程質(zhì)量控制領(lǐng)域的應(yīng)用 [J]. 食品與藥品, 2013, 15(2): 135-138.

[4] 李瓊, 趙磊, 王權(quán). 不同粒度黃芪三七膠囊粉體學(xué)性質(zhì)及質(zhì)量的初步對(duì)比 [J]. 現(xiàn)代中藥研究與實(shí)踐, 2019, 33(6): 32-35.

[5] Goh H P, Heng P W S, Liew C V. Comparative evaluation of powder flow parameters with reference to particle size and shape [J]., 2018, 547(1/2): 133-141.

[6] Fu X W, Huck D, Makein L,. Effect of particle shape and size on flow properties of lactose powders [J]., 2012, 10(2): 203-208.

[7] 梁兆昌, 褚洪標(biāo), 肖琳, 等. 杜仲超微粉體理化特性及體外溶出性能研究 [J]. 中草藥, 2015, 46(11): 1609-1614.

[8] 李淞明, 周永全, 王玲, 等. 中藥浸膏粉物理性質(zhì)及其表征方法的研究進(jìn)展 [J]. 中成藥, 2015, 37(11): 2491-2494.

[9] Faulhammer E, Llusa M, Radeke C,. The effects of material attributes on capsule fill weight and weight variability in dosator nozzle machines [J]., 2014, 471(1/2): 332-338.

[10] 許俊男, 涂傳智, 陳穎翀, 等. 激光粒度測(cè)定法在中藥粉體粒徑測(cè)定中的應(yīng)用與思考 [J]. 世界科學(xué)技術(shù)—中醫(yī)藥現(xiàn)代化, 2016, 18(10): 1776-1781.

[11] 宋曉嬌, 夏俊, 陳峰, 等. 粒度檢測(cè)方法研究及其在口服吸入制劑研發(fā)中的應(yīng)用 [J]. 藥物評(píng)價(jià)研究, 2019, 42(12): 2319-2324.

[12] 雷曉晴, 王秀麗, 李耿, 等. 近紅外光譜法快速測(cè)定當(dāng)歸中7種成分的含量 [J]. 中草藥, 2019, 50(16): 3947-3954.

[13] 王晴, 徐芳芳, 張欣, 等. 在線近紅外光譜監(jiān)測(cè)桂枝茯苓膠囊流化床干燥過程水分的方法研究 [J]. 中草藥, 2019, 50(22): 5429-5438.

[14] 高鴻彬, 瞿敏, 劉浩, 等. 不同產(chǎn)地覆盆子的近紅外漫反射快速無損鑒別 [J]. 世界中醫(yī)藥, 2020, 15(10): 1386-1390.

[15] Liu R, Li L, Yin W,. Near-infrared spectroscopy monitoring and control of the fluidized bed granulation and coating processes-A review [J]., 2017, 530(1/2): 308-315.

[16] Rosas J G, Blanco M, González J M,. Real-time determination of critical quality attributes using near-infrared spectroscopy: A contribution for process analytical technology (PAT) [J]., 2012, 97: 163-170.

[17] Otsuka M, Koyama A, Hattori Y. Real-time release monitoring for water content and mean particle size of granules in lab-sized fluid-bed granulator by near-infrared spectroscopy [J]., 2014, 4(34): 17461.

[18] Pauli V, Roggo Y, Kleinebudde P,. Real-time monitoring of particle size distribution in a continuous granulation and drying process by near infrared spectroscopy [J]., 2019, 141:90-99.

[19] 潘曉寧. 基于Kubelka-Munk理論的中藥近紅外漫反射定量研究 [D]. 北京: 北京中醫(yī)藥大學(xué), 2016.

[20] 范劍. 近紅外光譜技術(shù)在干姜和桂枝混合蒸餾提取過程中的應(yīng)用研究 [D]. 濟(jì)南: 山東大學(xué), 2017.

[21] 楊越. 基于近紅外光譜技術(shù)的中藥生產(chǎn)過程質(zhì)量控制方法研究 [D]. 杭州: 浙江大學(xué), 2018.

[22] 李春燕. 近紅外光譜法在葛根活性成分含量快速測(cè)定中的應(yīng)用研究 [D]. 新鄉(xiāng): 新鄉(xiāng)醫(yī)學(xué)院, 2018.

[23] 王永香, 鄭偉然, 米慧娟, 等. 熱毒寧注射液青蒿金銀花濃縮過程近紅外快速定量檢測(cè)方法的建立 [J]. 中草藥, 2017, 48(1): 102-108.

[24] 張娜, 徐冰, 賈帥蕓, 等. 丹參提取過程多源信息融合建模方法研究 [J]. 中草藥, 2018, 49(6): 1304-1310.

[25] 吳靜珠, 石瑞杰, 陳巖, 等. 食用油油酸的近紅外特征譜區(qū)優(yōu)選 [J]. 中國糧油學(xué)報(bào), 2015, 30(2): 118-121.

[26] 陳玥洋, 高志山, 郁曉暉, 等. 基于區(qū)間組合移動(dòng)窗口法篩選近紅外光譜信息[J]. 應(yīng)用光學(xué), 2017, 38(1): 99-105.

Research on near-infrared universal quantitative model of four intermediate particle sizes of Yaobitong Capsules

ZHANG Yong-chao1, XU Fang-fang2, 3, 4, ZHANG Xin2, 3, 4, XU Bing5, ZHAO Yue-mei2, 3, 4, ZHU Wen-deng2, 3, 4, WU Yun2, 3, 4, WANG Zhen-zhong2, 3, 4, XIAO Wei1, 2, 3, 4

1. Nanjing University of Chinese Medicine, Nanjing 210000, China 2. Jiangsu Kanion Pharmaceutical Co., Ltd., Lianyungang 222001, China 3. State Key Laboratory of New-tech for Chinese Medicine Pharmaceutical Process, Lianyungang 222001, China 4. The National & Local Joint Engineering Research Center on the Intelligent manufacturing of TCM, Lianyungang 222001, China 5. Beijing University of Chinese Medicine, Beijing 100029, China

To explore the feasibility of establishing a50near-infrared universal quantitative model for the four intermediates in the productive process of Yaobitong Capsules (YC) by analyzing the relationship between near infrared spectrum (NIRS) information and the median particle size (50).The NIRS of the four intermediates in the production of YC, including raw material fine powder, dry particles, sieving particles and total mixed particles, were collected. Different pretreatment methods were investigated. The interval partial least squares (iPLS), the synergy interval partial least squares (siPLS), and moving window partial least squares (mwPLS) were used to screen the optimal NIRS interval. The partial least squares (PLS) method was employed to build a50universal quantitative model for the four intermediates.In the universal model, the root mean square error of cross-validation (RMSECV) was 3.918 μm, the root mean square error of prediction (RMSEP) was 2.832 μm; The relative standard error of prediction (RSEP) was 2.26%, which was less than 5%; The ratio of predicton to deviation (RPD) was 4.60, which was greater than 3. The universal model could be used for quantitative determination, and its prediction accuracy was close to that of the independent model.The NIRS universal quantitative model can be used to determine the50of four intermediates in YC.

Yaobitong Capsules; near infrared spectroscopy; partial least squares; particle size; median diameter; universal model

R283.6

A

0253 - 2670(2021)01 - 0055 - 10

10.7501/j.issn.0253-2670.2021.01.008

2020-08-24

國家科技重大專項(xiàng)：基于功效成分群的中藥口服固體制劑先進(jìn)制藥與信息化技術(shù)融合示范應(yīng)用（2018ZX09201010-004）；康緣中藥智能化固體制劑工廠項(xiàng)目（KYYY190820）

張永超（1995—），女，碩士研究生，研究方向?yàn)橹兴幹扑幖夹g(shù)與工程。E-mail: zyc020896@163.com

肖 偉，研究員級(jí)高級(jí)工程師，博士，研究方向?yàn)橹兴幮滤幍难芯颗c開發(fā)。Tel: (025)87181851 E-mail: kanionlunwen@163.com

徐芳芳，博士，研究方向?yàn)檫^程分析技術(shù)。E-mail: 879164331@qq.com

[責(zé)任編輯 鄭禮勝]