• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于近、中紅外光譜的熱毒寧注射液制劑過(guò)程投料和二次熱處理工序快速檢測(cè)方法研究

    2022-11-05 06:53:24徐芳芳李執(zhí)棟章晨峰王振中
    中草藥 2022年21期
    關(guān)鍵詞:模型

    童 楓,徐芳芳,張 欣,李執(zhí)棟,吳 云,章晨峰,王振中

    基于近、中紅外光譜的熱毒寧注射液制劑過(guò)程投料和二次熱處理工序快速檢測(cè)方法研究

    童 楓1,徐芳芳2, 3*,張 欣2, 3,李執(zhí)棟2, 3,吳 云2, 3,章晨峰2, 3,王振中2, 3*

    1. 南京中醫(yī)藥大學(xué),江蘇 南京 210023 2. 江蘇康緣藥業(yè)股份有限公司,江蘇 連云港 222001 3. 中藥制藥過(guò)程新技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室,江蘇 連云港 222001

    應(yīng)用近紅外光譜(near-infrared reflectance spectroscopy,NIRS)與中紅外光譜(mid-infrared reflectance spectroscopy,MIRS)技術(shù),對(duì)熱毒寧注射液(Reduning Injection,RI)制劑過(guò)程的投料和二次熱處理工序中6種質(zhì)控指標(biāo)進(jìn)行快速檢測(cè),提高制劑過(guò)程的質(zhì)量控制水平。利用NIRS透射技術(shù)與MIRS衰減全反射技術(shù),結(jié)合偏最小二乘法(partial least squares,PLS),經(jīng)過(guò)光譜預(yù)處理方法的優(yōu)選以及波段篩選,分別建立綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷和固含量的快速預(yù)測(cè)模型,以校正集相關(guān)系數(shù)(cal)、驗(yàn)證集相關(guān)系數(shù)(pre)、校正均方根誤差(root mean square error of calibration,RMSEC)、交叉驗(yàn)證均方根誤差(root mean square error of cross validation,RMSECV)、預(yù)測(cè)相對(duì)誤差(relative standard error of prediction,RSEP)為評(píng)價(jià)指標(biāo),評(píng)價(jià)模型性能。NIRS預(yù)測(cè)模型的6個(gè)質(zhì)控指標(biāo)的RMSEC和RMSECV均小于0.3,RSEP小于4.0%;MIRS預(yù)測(cè)模型的6個(gè)質(zhì)控指標(biāo)的RMSEC和RMSECV均小于0.4,RSEP均小于5.0%,建立的PLS模型具有模型性能好、預(yù)測(cè)精度高的優(yōu)點(diǎn)。NIRS及MIRS分析技術(shù),均可用于RI投料和二次熱處理工序中6種質(zhì)控指標(biāo)的快速檢測(cè),模型RSEP在5%以內(nèi),方法操作簡(jiǎn)單,結(jié)果可靠。

    近紅外光譜;中紅外光譜;熱毒寧注射液;制劑過(guò)程;光譜預(yù)處理;波段篩選;偏最小二乘法;綠原酸;新綠原酸;隱綠原酸;梔子苷;斷氧化馬錢子苷;固含量

    熱毒寧注射液(Reduning Injection,RI)是江蘇康緣藥業(yè)股份有限公司的獨(dú)家品種,由青蒿、金銀花、梔子3味藥材制得,具有清熱、疏風(fēng)、解毒的功效,常用于上呼吸道感染所致的高熱、微惡風(fēng)寒、頭身痛、咳嗽、痰黃等癥狀及抑制流感病毒[1-2]。在中華人民共和國(guó)國(guó)家衛(wèi)生健康委員會(huì)發(fā)布的《新型冠狀病毒感染的肺炎診療方案(試行第九版)》中,針對(duì)重型、危重型新型冠狀病毒肺炎(coronavirus disease 2019,COVID-19)患者推薦的中成藥中,RI也名列其中[3]。

    RI的生產(chǎn)過(guò)程可以分為提取、精制和制劑過(guò)程,目前,利用近紅外光譜(near-infrared reflectance spectroscopy,NIRS)或中紅外光譜(mid-infrared reflectance spectroscopy,MIRS)技術(shù)結(jié)合化學(xué)計(jì)量學(xué),針對(duì)RI提取、精制過(guò)程的快速檢測(cè)已較為成熟,并取得一定的應(yīng)用效果[4-8]。但是,基于NIRS和MIRS的RI制劑過(guò)程的快速檢測(cè)研究尚未見報(bào)道。因此,基于RI的制劑過(guò)程,研究合適的快速檢測(cè)方法,可為RI制劑過(guò)程的質(zhì)量控制水平的提升提供技術(shù)支撐。NIRS以及MIRS作為過(guò)程分析技術(shù)的典型代表,近年來(lái)發(fā)展迅速。與傳統(tǒng)的化學(xué)分析方法相比,這2種光譜分析方法具有以下優(yōu)點(diǎn):(1)可同時(shí)用于定性和定量分析;(2)分析速度快,只需掃描得出光譜即可測(cè)定樣品中的各種成分;(3)不使用化學(xué)試劑,不損害樣品,對(duì)環(huán)境友好;(4)分析成本低,操作簡(jiǎn)單[9-11]。

    投料和二次熱處理工序指的是投料和二次熱處理2個(gè)連續(xù)工序,即在一定的溫濕度以及壓力下,依次將梔子提取物以及青金提取物投入到合并一定比例的注射用水和鹽酸的配液罐中,待冷卻后轉(zhuǎn)移至冷藏罐;冷藏罐內(nèi)藥液通過(guò)板框過(guò)濾器濾過(guò)至熱配間熱配罐中,調(diào)節(jié)蒸氣壓力以及藥液溫度,一段時(shí)間后對(duì)藥液進(jìn)行冷卻。

    投料和二次熱處理工序作為RI制劑過(guò)程的初始工序,在實(shí)際生產(chǎn)中,會(huì)出現(xiàn)提取物溶解不充分等情況,導(dǎo)致有效成分的損失,直接影響最后的成品質(zhì)量,因此,亟需引入1種快速檢測(cè)方法,降低生產(chǎn)過(guò)程控制風(fēng)險(xiǎn)。本實(shí)驗(yàn)將NIRS和MIRS技術(shù)與化學(xué)計(jì)量學(xué)結(jié)合,首次應(yīng)用于RI制劑過(guò)程,建立投料和二次熱處理工序6種關(guān)鍵質(zhì)控指標(biāo)的偏最小二乘(partial least squares,PLS)定量校正模型,實(shí)現(xiàn)對(duì)綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷及固含量的快速測(cè)定。

    1 材料與儀器

    1.1 材料

    59批次的投料工序樣本(批號(hào)Z220123~Z2201240、Z220201~Z220232、Z220301~Z220309)及對(duì)應(yīng)批次的59批次的二次熱處理工序樣本,共118批次樣品,由江蘇康緣藥業(yè)股份有限公司水針車間提供。對(duì)照品綠原酸(批號(hào)110753-202018,質(zhì)量分?jǐn)?shù)96.1%)、梔子苷(批號(hào)110749-201919,質(zhì)量分?jǐn)?shù)97.1%)購(gòu)自中國(guó)食品藥品檢定研究院;對(duì)照品新綠原酸(批號(hào)DSTDX001503,質(zhì)量分?jǐn)?shù)99.58%)、斷氧化馬錢子苷(批號(hào)DST211224-111,質(zhì)量分?jǐn)?shù)98.87%)購(gòu)自成都樂美天醫(yī)藥-德思特生物技術(shù)有限公司;隱綠原酸對(duì)照品(批號(hào)ST07850120,質(zhì)量分?jǐn)?shù)98.04%)購(gòu)自上海詩(shī)丹德標(biāo)準(zhǔn)技術(shù)服務(wù)有限公司。

    1.2 儀器

    UltiMate 3000型高效液相色譜儀、Antaris II型傅立葉變換NIRS儀,配有透射檢測(cè)器、采樣軟件以及Result、CAMO等數(shù)據(jù)處理軟件,(賽默飛世爾科技(中國(guó))有限公司);React IR 702L型MIRS儀、XP6型電子天平,梅特勒-托利多儀器(上海)有限公司;BSA 224S-CW型萬(wàn)分之一電子天平,德國(guó)Sartorius公司;KQ-500DB型數(shù)控超聲波清洗器,昆山市超聲儀器有限公司;Milli-Q Academic型純水機(jī),美國(guó)Millipore公司。

    1.3 數(shù)據(jù)處理軟件

    采用Unscramble X 10.4(挪威Camo Analytics公司)軟件進(jìn)行光譜預(yù)處理;采用GraphPad Prism 8.0(美國(guó)GraphPad Software)繪圖;采用Matlab R2018a(美國(guó)MathWorks公司)軟件進(jìn)行樣本劃分、波段篩選及模型構(gòu)建。

    2 方法與結(jié)果

    2.1 NIRS的采集

    以空氣為掃描背景,采用NIRS技術(shù)于室溫下采集光譜,選擇1 mm比色皿,光譜掃描范圍4000~10 000 cm?1,掃描次數(shù)為64次,分辨率4 cm?1,衰減器選擇B模塊,1倍增益,每個(gè)樣品采集3次,取平均光譜,得到59批投料工序樣本及59批二次熱處理工序樣本,共計(jì)118個(gè)樣本。NIRS結(jié)果見圖1。

    2.2 MIRS的采集

    在室溫條件下,使用光纖探頭掃描樣品,以空氣為掃描背景;光譜掃描范圍3000~650 cm?1;掃描次數(shù)32次;分辨率8 cm?1;增益選擇“l(fā)ow”,每個(gè)樣品掃描2次,取其一作為樣本光譜,得到59批投料工序樣本及59批二次熱處理工序樣本,共計(jì)118個(gè)樣本。MIRS結(jié)果見圖2。

    圖1 118批樣本的NIRS平均光譜圖

    圖2 118批樣本的MIRS圖

    2.3 5種質(zhì)控指標(biāo)含量的測(cè)定

    2.3.1 HPLC分析色譜條件 參照崔小弟等[12]的實(shí)驗(yàn)結(jié)果,色譜條件設(shè)為Kromasil C18(150 mm×4.6 mm,6 μm)色譜柱;流動(dòng)相為乙腈-0.1%磷酸水溶液,梯度洗脫:0~5 min,15%~18%乙腈;5~15 min,18%乙腈;15~38 min,18%~31%乙腈;38~38.5 min,31%~90%乙腈;38.5~40 min,90%乙腈;體積流量0.7 mL/min;進(jìn)樣量10 μL;柱溫20 ℃;檢測(cè)波長(zhǎng)238、327 nm。

    2.3.2 供試品溶液的制備 精密吸取1 mL樣品置于20 mL量瓶中,用50%甲醇溶解并稀釋至刻度,搖勻,0.45 μm微孔濾膜濾過(guò),即得。

    2.3.3 樣品測(cè)定 取118批次樣品,按照“2.3.2”項(xiàng)方法制備供試品溶液,按照“2.3.1”項(xiàng)色譜條件進(jìn)樣測(cè)定,每份樣品平行2次測(cè)定,取平均值,作為建模的質(zhì)控指標(biāo)含量實(shí)測(cè)值。6種質(zhì)控指標(biāo)的含量測(cè)定結(jié)果見表1。

    2.4 固含量的測(cè)定

    參照文獻(xiàn)方法[13-14],稱取約5 g樣品至已烘干至恒定質(zhì)量的稱量瓶(0)中,稱定質(zhì)量(1),置烘箱105 ℃條件下烘干5 h至恒定質(zhì)量,計(jì)為2。

    表1 樣本中6種質(zhì)控指標(biāo)統(tǒng)計(jì)值(n = 2)

    Table 1 Statistical values of six quality control indicators in samples (n = 2)

    質(zhì)控指標(biāo)質(zhì)量濃度/(mg·mL?1)固含量/% 綠原酸新綠原酸隱綠原酸梔子苷斷氧化馬錢子苷 最小值8.768 82.066 82.304 312.285 41.116 27.295 3 下四分位數(shù)9.405 62.349 42.601 212.791 31.178 77.765 5 上四分位數(shù)10.279 12.619 42.898 014.254 11.432 78.221 8 最大值10.606 12.809 43.052 015.311 11.506 98.725 7 均值9.851 82.478 32.740 813.633 51.302 37.995 1

    固含量=(2-0)/(1-0)

    2.5 校正集與驗(yàn)證集的劃分

    采用Kennard-Stone(K-S)劃分法[15],將收集的59批投料工序樣品以及對(duì)應(yīng)批次的59批二次熱處理工序樣品,共計(jì)118批樣品,以4∶1的比例劃分為校正集和驗(yàn)證集,最終得到校正集樣本94個(gè),驗(yàn)證集樣本24個(gè)。

    2.6 光譜預(yù)處理

    采集光譜時(shí),由于受到環(huán)境溫度、濕度、儀器狀態(tài)、雜散光和隨機(jī)噪聲等因素的影響,紅外光譜會(huì)出現(xiàn)基線漂移和噪聲干擾等問題,因此,在建模前對(duì)光譜進(jìn)行適當(dāng)?shù)念A(yù)處理是有必要的,可以有效地減少噪聲,凈化無(wú)用信息,提高模型的精度和預(yù)測(cè)效果。常見的預(yù)處理方式有:移動(dòng)平均法(moving average,MA)、Savitzky-Golay(S-G)平滑法、導(dǎo)數(shù)法(一階導(dǎo)數(shù)、二階導(dǎo)數(shù))、矢量歸一化法、基線校正、標(biāo)準(zhǔn)正態(tài)變換(standard normal transformation,SNV)、多元散射校正(multivariate scattering correction,MSC)以及上述方法的結(jié)合。MA和S-G平滑法可以消除光譜信息中的隨機(jī)噪聲,提高信噪比;一階導(dǎo)數(shù)可以消除背景的常數(shù)平移,二階導(dǎo)數(shù)可以消除背景的線性平移,消除背景漂移造成的影響;矢量歸一化法通過(guò)增強(qiáng)光譜數(shù)據(jù)間的差異,校正由光程變化或者樣品濃度改變引起的光譜變化;基線校正可以消除儀器背景漂移對(duì)信號(hào)的影響;SNV和MSC可以消除由于樣品分布不均勻?qū)庾V造成的影響[16]。

    本研究采用MA、S-G平滑、S-G 1st、基線校正、歸一化法以及SNV分別對(duì)NIRS和MIRS進(jìn)行預(yù)處理,并分別建立投料和二次熱處理工序的6種質(zhì)控指標(biāo)的PLS定量校正模型。根據(jù)模型評(píng)價(jià)指標(biāo)篩選出最佳的光譜預(yù)處理方法,使用留一交叉驗(yàn)證法確定校正模型中的最佳的潛變量數(shù)(latent variables,LVs),以預(yù)測(cè)相對(duì)誤差(relative standard error of prediction,RSEP)為主要評(píng)價(jià)指標(biāo),RSEP值越小越好。一般而言,性能優(yōu)異的預(yù)測(cè)模型應(yīng)具有較高的校正集相關(guān)系數(shù)(cal)及驗(yàn)證集相關(guān)系數(shù)(pre)、較小且接近的校正均方根誤差(root mean square error of calibration,RMSEC)和交叉驗(yàn)證均方根誤差(root mean square error of cross validation,RMSECV)及較低的RSEP[17]。上述評(píng)價(jià)指標(biāo)的相關(guān)公式為式(1)~(4)。

    cal=[1-∑(-y)2/∑(-y)2]1/2(1)

    pre=[1-∑(-y)2/∑(-y)2]1/2(2)

    RMSEC=[∑(-y)2/]1/2(3)

    RSEP=[∑(-y)2/∑2]1/2(4)

    、為校正集、驗(yàn)證集樣本數(shù),為實(shí)測(cè)值,y為預(yù)測(cè)值,y、y分別為校正集和驗(yàn)證集實(shí)測(cè)值的平均值

    NIRS和MIRS預(yù)處理結(jié)果見表2、3,可以發(fā)現(xiàn):基于NIRS,當(dāng)綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷及固含量分別以SNV、MA、MA、無(wú)預(yù)處理、SNV和歸一化法為預(yù)處理方法時(shí),RSEP最小,且都在4.0%以下,此時(shí),cal與pre值較大,RMSEC與RMSECV值較小且接近;基于MIRS,當(dāng)綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷及固含量的預(yù)處理方法分別為無(wú)預(yù)處理、MA、MA、MA、SNV和無(wú)預(yù)處理時(shí),模型具有較低的RSEP,均小于5%,較大的cal與pre值以及較小且接近的RMSEC與RMSECV值,從預(yù)處理方法的比較中可以猜測(cè),在MIRS中,MA適合多數(shù)質(zhì)控指標(biāo),可能是由于其減輕了光譜信息自身攜帶的隨機(jī)誤差,提高信噪比。

    2.7 光譜波段的篩選

    全光譜模型的準(zhǔn)確性和穩(wěn)健性可能會(huì)因?yàn)橐恍o(wú)關(guān)信息的存在而減弱,篩選光譜波段可以剔除無(wú)用信息,提高模型的預(yù)測(cè)精度和穩(wěn)定性。本研究在上述篩選出的最佳預(yù)處理方法的基礎(chǔ)上,對(duì)比全光譜、間隔PLS(interval PLS,iPLS)、組合間隔PLS(synergy interval PLS,siPLS)及移動(dòng)窗口PLS(moving windows PLS,mwPLS)建模的效果。

    2.7.1 iPLS篩選最優(yōu)波段 iPLS是由N?rgaard等[18]提出的一種光譜區(qū)間篩選方法,通過(guò)將全光譜平均分為若干子區(qū)間,在各子區(qū)間內(nèi)以不同的潛變量數(shù)建立獨(dú)立的PLS回歸模型[19]。本研究將光譜區(qū)間均分為20個(gè)子區(qū)間,并以RMSECV作為評(píng)價(jià)指標(biāo),篩選出最優(yōu)建模波段。

    2.7.2 siPLS篩選最優(yōu)波段 siPLS是在iPLS的基礎(chǔ)上,對(duì)不同子區(qū)間進(jìn)行任意組合,建立所有可能的2、3或4個(gè)區(qū)間的PLS回歸模型[20]。本研究將光譜區(qū)間均分為20個(gè)子區(qū)間,以子區(qū)間組合數(shù)為3建立模型,以RMSECV作為評(píng)價(jià)指標(biāo),篩選出最優(yōu)建模波段。

    2.7.3 mwPLS篩選最優(yōu)波段 mwPLS的基本思想是將1個(gè)窗口沿著光譜軸連續(xù)移動(dòng),每移動(dòng)1個(gè)波長(zhǎng)點(diǎn),采用交叉驗(yàn)證法確定最佳潛變量數(shù)并建立相應(yīng)的PLS最優(yōu)模型[21]。本研究以全光譜波數(shù)的10%作為窗口寬度,以RMSECV作為評(píng)價(jià)指標(biāo),篩選出最優(yōu)建模波段。

    表2 NIRS不同預(yù)處理方法建模結(jié)果分析

    Table 2 Analysis of modeling results of different NIRS pretreatment methods

    質(zhì)控指標(biāo)預(yù)處理方法rcalRMSECRMSECVrpreRESP/%LVs質(zhì)控指標(biāo)預(yù)處理方法rcalRMSECRMSECVrpreRESP/%LVs 綠原酸無(wú)預(yù)處理0.979 80.088 10.158 50.933 31.717 515梔子苷無(wú)預(yù)處理0.975 40.187 40.257 10.974 61.100 912 MA0.976 00.096 00.164 50.924 01.828 116 MA0.979 00.203 50.276 30.972 81.168 512 S-G平滑0.975 70.096 60.159 70.926 41.788 715 S-G平滑0.972 50.197 90.269 60.974 41.120 012 S-G 1st0.982 50.082 20.186 00.922 21.818 210 S-G 1st0.998 80.041 40.290 10.948 21.551 116 基線校正0.980 00.087 80.176 40.932 01.709 315 基線校正0.974 70.189 90.274 40.975 01.132 612 歸一化法0.971 30.104 70.153 80.945 51.544 112 歸一化法0.998 90.039 80.177 00.975 51.077 620 SNV0.969 90.107 20.152 50.951 91.473 412 SNV0.999 00.038 70.175 40.974 71.086 620 新綠原酸無(wú)預(yù)處理0.834 30.096 80.127 00.829 93.934 610斷氧化馬無(wú)預(yù)處理0.973 00.030 2 0.049 1 0.960 22.573 814 MA0.811 00.101 50.130 30.846 63.777 010錢子苷MA0.984 70.023 0 0.048 7 0.947 52.860 518 S-G平滑0.818 50.099 60.128 80.843 13.820 910 S-G平滑0.984 80.023 0 0.050 2 0.956 52.708 017 S-G 1st0.957 10.050 20.118 10.686 45.084 410 S-G 1st0.993 80.014 7 0.046 8 0.965 92.546 312 基線校正0.799 00.104 30.134 50.846 93.807 19 基線校正0.993 70.014 8 0.045 5 0.963 12.556 418 歸一化法0.849 10.094 30.123 90.837 43.847 010 歸一化法0.974 40.029 7 0.047 3 0.963 82.533 114 SNV0.838 30.103 10.130 80.840 13.850 69 SNV0.973 00.030 5 0.047 1 0.968 22.294 614 隱綠原酸無(wú)預(yù)處理0.872 40.085 10.111 60.942 42.591 110固含量無(wú)預(yù)處理0.973 40.077 70.131 60.941 21.378 415 MA0.865 50.087 00.111 50.951 12.507 310 MA0.970 50.081 80.133 20.936 11.423 917 S-G平滑0.868 00.086 50.111 60.950 22.527 310 S-G平滑0.967 70.085 50.128 80.943 91.388 215 S-G 1st0.978 00.036 30.098 50.922 23.170 511 S-G 1st0.995 50.031 90.151 70.906 21.704 714 基線校正0.882 10.082 10.110 80.941 22.725 010 基線校正0.975 10.075 30.138 60.927 71.488 515 歸一化法0.881 60.082 10.108 10.949 82.548 210 歸一化法0.978 20.070 40.120 60.952 91.259 115 SNV0.885 40.080 80.106 40.945 32.551 610 SNV0.978 50.069 90.119 90.948 61.309 915

    除固含量的RMSEC與RMSECV的單位為“%”外,其余5種指標(biāo)的RMSEC與RMSECV的單位均為“mg?mL?1”,下表同

    except the unit of RMSEC and RMSECV of solid content is “%”, the units of RMSEC and RMSECV of the other five indicators are all “mg?mL?1”, same as the following table

    表3 MIRS不同預(yù)處理方法建模結(jié)果分析

    Table 3 Analysis of modeling results of different MIRS pretreatment methods

    質(zhì)控指標(biāo)預(yù)處理方法rcalRMSECRMSECVrpreRESP/%LVs質(zhì)控指標(biāo)預(yù)處理方法rcalRMSECRMSECVrpreRESP/%LVs 綠原酸無(wú)預(yù)處理0.924 50.168 20.267 10.876 42.324 46梔子苷無(wú)預(yù)處理0.911 40.350 00.594 30.723 13.883 16 MA0.916 40.176 60.268 20.870 02.374 89 MA0.968 10.212 90.577 80.822 93.093 019 S-G平滑0.938 20.152 60.280 30.866 62.352 69 S-G平滑0.908 10.356 20.602 90.704 34.062 38 S-G 1st0.929 20.163 10.313 20.849 62.529 38 S-G 1st0.963 30.228 10.667 60.745 64.248 716 基線校正0.925 90.166 40.293 30.861 02.390 27 基線校正0.912 70.347 40.608 10.720 03.761 06 歸一化法0.900 10.191 90.308 80.841 82.623 86 歸一化法0.922 00.329 10.657 60.696 94.207 28 SNV0.909 90.182 70.312 60.836 82.639 76 SNV0.897 60.374 70.655 20.652 64.355 56 新綠原酸無(wú)預(yù)處理0.802 10.103 50.150 10.738 54.691 85斷氧化馬無(wú)預(yù)處理0.870 70.0650 0.1133 0.780 15.630 66 MA0.806 40.107 10.160 40.765 14.447 58錢子苷MA0.836 70.0724 0.0972 0.736 46.110 47 S-G平滑0.744 80.115 70.156 40.701 74.935 95 S-G平滑0.869 90.0652 0.1112 0.739 86.084 88 S-G 1st0.798 30.104 40.163 10.774 14.728 16 S-G 1st0.967 30.0335 0.0977 0.727 37.128 516 基線校正0.815 60.100 40.154 90.721 84.877 85 基線校正0.900 60.0576 0.1121 0.802 85.308 27 歸一化法0.812 90.100 90.144 10.750 74.566 55 歸一化法0.883 00.0620 0.1080 0.787 85.455 56 SNV0.808 90.101 90.146 20.743 34.634 85 SNV0.897 00.0584 0.1085 0.806 85.247 96 隱綠原酸無(wú)預(yù)處理0.893 60.087 80.134 20.799 64.848 55固含量無(wú)預(yù)處理0.914 60.138 50.230 20.784 22.547 06 MA0.830 80.096 90.145 00.870 24.172 68 MA0.959 30.095 70.221 90.796 62.822 916 S-G平滑0.863 50.087 70.148 10.802 04.798 07 S-G平滑0.882 00.159 90.239 10.718 92.882 66 S-G 1st0.821 60.099 10.151 40.831 94.441 3 5 S-G 1st0.922 70.121 30.265 90.775 12.748 010 基線校正0.842 70.093 70.145 10.801 74.770 05 基線校正0.901 30.147 00.241 60.765 82.596 66 歸一化法0.852 30.091 00.134 10.813 24.855 35 歸一化法0.908 80.141 40.256 20.725 12.931 17 SNV0.847 00.092 40.137 80.822 94.587 05 SNV0.910 80.139 90.253 80.707 83.010 77

    2.7.4 最優(yōu)波段篩選結(jié)果 由3種波段篩選方法建立的模型與全光譜模型的性能參數(shù)對(duì)比結(jié)果見表4,以RMSEC、RMSECV和RSEP作為評(píng)價(jià)指標(biāo),綜合評(píng)價(jià),篩選出最優(yōu)建模波段。結(jié)果表明,基于NIRS,綠原酸和斷氧化馬錢子苷的預(yù)測(cè)模型經(jīng)3種方法篩選波段后,雖然建模變量數(shù)減少,但RMSEC、RMSECV普遍增大,且RSEP均增大,模型預(yù)測(cè)性能下降,故選擇全光譜區(qū)間(3 999.64~9 999.10 cm?1)建模;新綠原酸、隱綠原酸及梔子苷預(yù)測(cè)模型,采用siPLS篩選波段后,建模變量數(shù)降為465個(gè),cal、pre均增大,RMSEC、RMSECV均減小,且RSEP均減小,模型預(yù)測(cè)性能提升,故新綠原酸、隱綠原酸及梔子苷分別選擇4 601.32~4 900.23、 5 503.84~5 802.76、7 609.73~7 906.71 cm?1, 5 503.84~5 802.76、5 804.69~6 103.60、7 308.89~7 607.80 cm?1及4 300.48~4 599.39、5 203.00~5 501.92、5 804.69~6 103.60 cm?1作為最優(yōu)建模波段;固含量預(yù)測(cè)模型采用siPLS篩選波段后,雖然cal、pre略微減小,RSEP略微增大,但是RMSEC、RMSECV減小且接近,故固含量選擇4 300.48~ 4 599.39、5 503.84~5 802.76、7 008.05~7 306.96 cm?1作為最優(yōu)建模波段。

    表4 不同波段篩選方法對(duì)NIRS及MIRS的影響

    Table 4 Influences of different band screening methods on NIRS and MIRS spectra

    質(zhì)控指標(biāo)光譜篩選方法ν/cm?1rcalRMSECRMSECVrpreRESP/%LVs 綠原酸NIRS全光譜3 999.64~9 999.100.969 90.107 20.152 50.951 91.473 412 siPLS4 300.48~4 599.39, 7 008.05~7 306.96,0.978 40.091 10.138 90.940 81.606 213 7 609.73~7 906.71 iPLS4 300.48~4 599.390.974 10.099 60.162 80.903 72.014 610 mwPLS4 387.26~4 987.020.922 70.169 70.194 30.923 91.865 78 MIRS全光譜3 000.00~648.000.924 50.168 20.267 10.876 42.324 46 siPLS1 572.00~1 460.00, 1 224.00~1 112.00,0.940 60.149 60.194 30.902 42.093 712 1 108.00~996.00 iPLS1 224.00~1 112.000.886 30.205 90.244 90.826 22.743 910 mwPLS1 176.00~944.000.865 50.220 60.296 30.863 32.453 411 新綠原酸NIRS全光譜3 999.64~9 999.100.811 00.101 50.130 30.846 63.777 010 siPLS4 601.32~4 900.23, 5 503.84~5 802.76,0.877 30.083 30.108 90.870 53.431 812 7 609.73~7 906.71 iPLS4 601.32~4 900.230.767 20.111 90.147 50.890 03.192 06 mwPLS6 067.60~6 676.350.673 00.128 30.145 60.762 14.545 66 MIRS全光譜3 000.00~648.000.806 40.107 10.160 40.765 14.447 58 siPLS1 456.00~1 344.00, 1 340.00~1 228.00,0.800 80.103 80.121 90.738 54.661 210 1 108.00~996.00 iPLS1 224.00~1 112.000.645 40.132 70.146 70.717 24.811 18 mwPLS1 401.00~1 172.000.767 60.129 90.147 90.790 34.317 99 隱綠原酸NIRS全光譜3 999.64~9 999.100.865 50.087 00.111 50.951 12.507 310 siPLS5 503.84~5 802.76, 5 804.69~6 103.60,0.900 60.075 70.099 30.952 62.589 510 7 308.89~7 607.80 iPLS4 601.32~4 900.230.844 40.094 40.121 10.921 73.433 07 mwPLS4 697.75~5 297.500.872 40.085 60.169 70.855 74.340 07 MIRS全光譜3 000.00~648.000.830 80.096 90.145 00.870 24.172 68 siPLS1 340.00~1 228.00, 1 108.00~996.00,0.763 80.112 30.135 10.829 44.464 46 992.00~880.00 iPLS1 224.00~1 112.000.777 90.109 50.129 80.834 24.453 09 mwPLS1 412.00~1 180.000.828 70.097 20.115 20.848 44.283 211

    續(xù)表4

    基于MIRS,隱綠原酸預(yù)測(cè)模型經(jīng)3種方法篩選波段后,雖然建模變量數(shù)減少,但cal、pre減小,RMSEC、RMSECV普遍增大,且RSEP增大,模型預(yù)測(cè)性能下降,故選擇全光譜區(qū)間(3 000.00~648.00 cm?1)建模;綠原酸、梔子苷、斷氧化馬錢子苷和固含量的預(yù)測(cè)模型經(jīng)采用siPLS篩選最優(yōu)波段后,建模變量數(shù)降為87個(gè),cal、pre普遍增大,RMSEC、RMSECV較小且接近,RSEP均增大,模型預(yù)測(cè)性能提升,故綠原酸、梔子苷、斷氧化馬錢子苷和固含量分別選擇1 572.00~1 460.00、 1 224.00~1 112.00、1 108.00~996.00 cm?1, 1 572.00~1 460.00、1 340.00~1 228.00、1 108.00~996.00 cm?1,1 688.00~1 576.00、1 224.00~ 1 112.00、1 108.00~996.00 cm?1及1 456.00~ 1 344.00、1 224.00~1 112.00、1 108.00~996.00 cm?1作為最優(yōu)建模波段;新綠原酸預(yù)測(cè)模型采用mwPLS篩選最優(yōu)波段后,建模變量數(shù)降為59個(gè),RMSEC、RMSECV較小且接近,且RSEP略微減小,模型預(yù)測(cè)性能提升,故新綠原酸選擇1 401.00~1 172.00 cm?1作為最優(yōu)建模波段。

    2.8 PLS模型的建立

    模型經(jīng)過(guò)光譜預(yù)處理方法的選擇、最佳的潛變量數(shù)的確定及最優(yōu)建模波段的確定后,運(yùn)用PLS分別建立了6個(gè)質(zhì)控指標(biāo)的NIRS和MIRS最佳定量校正模型,6個(gè)質(zhì)控指標(biāo)的NIRS模型預(yù)測(cè)效果均強(qiáng)于MIRS模型,故選擇6個(gè)質(zhì)控指標(biāo)的NIRS預(yù)測(cè)模型作為最優(yōu)模型,結(jié)果見表5。圖3為對(duì)應(yīng)的6種質(zhì)控指標(biāo)的最佳模型的預(yù)測(cè)值與實(shí)測(cè)值的相關(guān)性。從表5可見,對(duì)于RI制劑過(guò)程的投料和二次熱處理工序的6種質(zhì)控指標(biāo)綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷及固含量,利用NIRS建模,其RSEP分別為1.473 4%、3.431 8%、2.589 5%、1.070 6%、2.294 6%、1.329 8%,均小于4%。而利用MIRS建模,其RSEP分別為2.093 7%、4.317 9%、4.172 6%、2.320 5%、3.816 4%、1.720 3%,均小于5%。鑒于NIRS及MIRS模型的RMSEC及RMSECV也均令人滿意,可以認(rèn)為NIRS及MIRS模型均滿足實(shí)際運(yùn)用的需要,可以快速有效地預(yù)測(cè)投料和二次熱處理工序的6種質(zhì)控指標(biāo)的含量預(yù)測(cè)。

    表5 6種質(zhì)控指標(biāo)的最優(yōu)模型性能參數(shù)

    Table 5 Optimal model performance parameters of six quality control indexes

    質(zhì)控指標(biāo)光譜預(yù)處理方法波段篩選LVsrcalRMSECRMSECVrpreRSEP/% 綠原酸NIRSSNV全光譜120.969 90.107 20.152 50.951 91.473 4 MIRS無(wú)預(yù)處理siPLS60.940 60.149 60.194 30.902 42.093 7 新綠原酸NIRSMAsiPLS120.877 30.083 30.108 90.870 53.431 8 MIRSMAmwPLS90.767 60.129 90.147 90.790 34.317 9 隱綠原酸NIRSMAsiPLS100.900 60.075 70.099 30.952 62.589 5 MIRSMA全光譜80.830 80.096 90.145 00.870 24.172 6 梔子苷NIRS無(wú)預(yù)處理siPLS150.989 70.125 60.252 40.975 71.070 6 MIRSMAsiPLS170.966 50.218 30.324 50.888 32.320 5 斷氧化馬錢子苷NIRSSNV全光譜140.973 00.030 5 0.047 1 0.968 22.294 6 MIRSSNVsiPLS110.930 00.048 6 0.069 1 0.892 73.816 4 固含量NIRS歸一化法siPLS130.970 20.082 20.108 10.948 31.329 8 MIRS無(wú)預(yù)處理siPLS80.903 10.145 60.166 70.899 21.720 3

    圖3 最佳模型中預(yù)測(cè)值與實(shí)測(cè)值的的相關(guān)性

    將驗(yàn)證集的光譜數(shù)據(jù)導(dǎo)入用已建立的最佳模型中,即NIRS預(yù)測(cè)模型。根據(jù)樣品的實(shí)測(cè)值和模型預(yù)測(cè)值,計(jì)算其相對(duì)誤差。綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷和固含量的模型預(yù)測(cè)值與樣本實(shí)測(cè)值的平均絕對(duì)偏差(mean absolute,MAD)及平均相對(duì)偏差(mean relative bais,MRB)見表6,結(jié)果顯示這6種質(zhì)控指標(biāo)的MRB均小于3.00%。

    表6 驗(yàn)證集樣本在NIRS預(yù)測(cè)模型中的預(yù)測(cè)值與實(shí)測(cè)值的對(duì)比

    Table 6 Comparison of predicted and measured values in validation set samples

    質(zhì)控指標(biāo)平均絕對(duì)偏差/(mg·mL?1)平均相對(duì)偏差/% 綠原酸0.114 11.20 新綠原酸0.070 82.90 隱綠原酸0.055 32.11 梔子苷0.118 10.89 斷氧化馬錢子苷0.024 81.85 固含量0.080 21.00

    3 討論

    本研究首次以RI制劑過(guò)程為研究對(duì)象,采用NIRS和MIRS技術(shù),結(jié)合PLS法,經(jīng)過(guò)光譜預(yù)處理方法及建模波段篩選,建立了投料和二次熱處理工序的綠原酸、新綠原酸、隱綠原酸、梔子苷、斷氧化馬錢子苷和固含量的定量預(yù)測(cè)模型。對(duì)比NIRS和MIRS模型對(duì)6種質(zhì)控指標(biāo)的預(yù)測(cè)效果,可以發(fā)現(xiàn),NIRS模型均優(yōu)于MIRS模型,這可能與質(zhì)控指標(biāo)的質(zhì)量濃度有關(guān)系。5種質(zhì)控指標(biāo)成分的平均質(zhì)量濃度分別為9.851 8、2.478 3、2.740 8、13.633 5、1.302 3 mg/mL,均高于1 mg/mL,可以認(rèn)為均為高質(zhì)量濃度分析物,平均固含量7.995 1%亦屬于高質(zhì)量分?jǐn)?shù)。這與Tao等[22]的研究結(jié)果相符合,即對(duì)于高質(zhì)量濃度的分析物,NIRS表現(xiàn)出比MIRS更好的預(yù)測(cè)性能。

    根據(jù)平均質(zhì)量濃度對(duì)6種質(zhì)控指標(biāo)進(jìn)一步分組,將新綠原酸、隱綠原酸及斷氧化馬錢子苷分為高質(zhì)量濃度組,綠原酸、梔子苷及固含量分為超高質(zhì)量濃度組,以預(yù)測(cè)模型的RSEP作為決定評(píng)價(jià)指標(biāo),可以發(fā)現(xiàn),超高質(zhì)量濃度組預(yù)測(cè)模型的RSEP均小于較高質(zhì)量濃度組,可以猜測(cè),NIRS模型對(duì)高質(zhì)量濃度樣品具有優(yōu)異的預(yù)測(cè)效果,并且質(zhì)量濃度越高,預(yù)測(cè)效果越好,后續(xù)可以通過(guò)增大樣本量來(lái)驗(yàn)證。NIRS模型與MIRS模型預(yù)測(cè)效果存在差異,也有可能是由于NIRS相較于MIRS,其波數(shù)更多,建模變量數(shù)更多,數(shù)據(jù)中包含的信息更加豐富,使得其最終模型預(yù)測(cè)效果更佳,這也需要后續(xù)的研究證明。

    本研究表明,NIRS與MIRS模型都適用于這6種質(zhì)控指標(biāo)的預(yù)測(cè),2種光譜技術(shù)在離線或者在線檢測(cè)方面各有所長(zhǎng),NIRS模型預(yù)測(cè)精度高、檢測(cè)速度快、穩(wěn)定性好[23];MIRS建模參數(shù)少、儀器小巧便攜,而且可以適用于混懸體系樣品的檢測(cè),減少了樣品的預(yù)處理過(guò)程,擴(kuò)大了紅外光譜的應(yīng)用范圍[24],但同時(shí)MIRS對(duì)環(huán)境的要求更高,儀器位置的改變、溫濕度及壓力的變化、光纖的彎折程度等都會(huì)對(duì)光譜產(chǎn)生影響。與傳統(tǒng)化學(xué)分析方法相比,這2種光譜分析技術(shù)更環(huán)保,耗時(shí)更少,均可以對(duì)RI投料和二次熱處理工序的質(zhì)控指標(biāo)進(jìn)行快速檢測(cè),也為后續(xù)的RI制劑過(guò)程的全工序快速檢測(cè)研究奠定了基礎(chǔ)。

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

    [1] 孫蘭, 周軍, 王振中. 熱毒寧注射液藥理作用與臨床研究進(jìn)展 [J]. 中國(guó)中醫(yī)急癥, 2014, 23(12): 2247-2249.

    [2] 王高舉, 焦紅軍. 基于網(wǎng)絡(luò)藥理學(xué)和分子對(duì)接技術(shù)的熱毒寧注射液抗SARS、MERS和COVID-19的潛在共性作用機(jī)制與活性成分研究 [J]. 藥物評(píng)價(jià)研究, 2022, 45(5): 842-852.

    [3] 新型冠狀病毒肺炎診療方案(試行第九版) [J]. 中國(guó)醫(yī)藥, 2022, 17(4): 481-487.

    [4] Sinelli N, Cerretani L, Egidio V D,. Application of near (NIR) infrared and mid (MIR) infrared spectroscopy as a rapid tool to classify extra virgin olive oil on the basis of fruity attribute intensity [J]., 2010, 43(1): 369-375.

    [5] Urí?ková V, Sádecká J. Determination of geographical origin of alcoholic beverages using ultraviolet, visible and infrared spectroscopy: A review [J]., 2015, 148: 131-137.

    [6] Nenadis N, Tsimidou M Z. Perspective of vibrational spectroscopy analytical methods in on-field/official control of olives and virgin olive oil [J]., 2017, 119(1): 1600148.

    [7] 耿姝, 陳勇, 金葉, 等. 基于近紅外技術(shù)的金銀花藥材多指標(biāo)成分快速檢測(cè) [J]. 中國(guó)現(xiàn)代應(yīng)用藥學(xué), 2016, 33(2): 144-150.

    [8] 王磊, 楊越, 潘紅燁, 等. 熱毒寧注射液金銀花提取過(guò)程軌跡研究 [J]. 世界科學(xué)技術(shù)—中醫(yī)藥現(xiàn)代化, 2017, 19(8): 1277-1282.

    [9] 徐芳芳, 馮雙雙, 李雪珂, 等. 青蒿濃縮過(guò)程在線近紅外快速檢測(cè)模型的建立 [J]. 中草藥, 2016, 47(10): 1690-1695.

    [10] 吳莎, 劉啟安, 吳亞男, 等. 近紅外光譜對(duì)熱毒寧注射液梔子萃取過(guò)程中的可行性分析 [J]. 中國(guó)實(shí)驗(yàn)方劑學(xué)雜志, 2015, 21(1): 53-56.

    [11] 王璇, 王永香, 畢宇安, 等. 熱毒寧注射液中金銀花和青蒿干燥工藝過(guò)程的批放行標(biāo)準(zhǔn) [J]. 中成藥, 2017, 39(10): 2059-2063.

    [12] 崔小弟, 李執(zhí)棟, 徐芳芳, 等. 熱毒寧注射液多指標(biāo)成分的一測(cè)多評(píng)法同時(shí)測(cè)定 [J]. 時(shí)珍國(guó)醫(yī)國(guó)藥, 2019, 30(5): 1116-1119.

    [13] 金葉, 丁海櫻, 吳永江, 等. 近紅外光譜技術(shù)用于血必凈注射液提取過(guò)程的在線檢測(cè)研究 [J]. 藥物分析雜志, 2012, 32(7): 1214-1221.

    [14] 杜文俊, 劉雪松, 陶玲艷, 等. 熱毒寧注射液金銀花和青蒿(金青)醇沉過(guò)程中多指標(biāo)的近紅外快速檢測(cè) [J]. 中草藥, 2015, 46(1): 61-66.

    [15] 陳奕云, 趙瑞瑛, 齊天賜, 等. 結(jié)合光譜變換和Kennard-Stone算法的水稻土全氮光譜估算模型校正集構(gòu)建策略研究 [J]. 光譜學(xué)與光譜分析, 2017, 37(7): 2133-2139.

    [16] 第五鵬瑤, 卞希慧, 王姿方, 等. 光譜預(yù)處理方法選擇研究 [J]. 光譜學(xué)與光譜分析, 2019, 39(9): 2800-2806.

    [17] 張永超, 徐芳芳, 張欣, 等. 腰痹通膠囊4種中間體粒徑的近紅外光譜通用定量模型研究 [J]. 中草藥, 2021, 52(1): 55-64.

    [18] N?rgaard L, Hahn M T, Knudsen L B,. Multivariate near-infrared and Raman spectroscopic quantifications of the crystallinity of lactose in whey permeate powder [J]., 2004, 15(12): 1261-1270.

    [19] 冼瑞儀, 黃富榮, 黎遠(yuǎn)鵬, 等. 可見和近紅外透射光譜結(jié)合區(qū)間偏最小二乘法(iPLS)用于橄欖油中摻雜煎炸老油的定量分析 [J]. 光譜學(xué)與光譜分析, 2016, 36(8): 2462-2467.

    [20] 王彩虹, 黃林, 劉木華, 等. 基于SiPLS模型的稻殼中重金屬鉻LIBS檢測(cè) [J]. 激光與光電子學(xué)進(jìn)展, 2016, 53(11): 298-305.

    [21] 史智佳, 田寒友, 鄒昊, 等. MWPLS法在羊肉中TVB-N近紅外定量分析中的應(yīng)用 [J]. 食品科學(xué), 2015, 36(20): 218-221.

    [22] Tao L Y, Lin Z L, Chen J S,. Mid-infrared and near-infrared spectroscopy for rapid detection ofby a liquid-liquid extraction process [J]., 2017, 145: 1-9.

    [23] 汪東升. 在線近紅外檢測(cè)系統(tǒng)的研發(fā)及在油菜籽含油率檢測(cè)中的應(yīng)用 [D]. 鎮(zhèn)江: 江蘇大學(xué), 2020.

    [24] 徐芳芳, 杜慧, 張欣, 等. 在線中紅外光譜監(jiān)測(cè)熱毒寧注射液金銀花與青蒿醇沉過(guò)程7種指標(biāo)成分研究 [J]. 中草藥, 2021, 52(10): 2909-2917.

    Research on fast detection method of feeding and secondary heat treatment process of Reduning Injection preparation process by near-infrared and mid-infrared reflectance spectroscopy

    TONG Feng1, XU Fang-fang2, 3, ZHANG Xin2, 3, LI Zhi-dong2, 3, WU Yun2, 3, ZHANG Chen-feng2, 3, WANG Zhen-zhong2, 3

    1. Nanjing University of Chinese Medicine, Nanjing 210023, China 2. Jiangsu Kangyuan Pharmaceutical Co., Ltd., Lianyungang 222001, China 3. State Key Laboratory of New Technology for Pharmaceutical Process of Traditional Chinese Medicine, Lianyungang 222001, China

    The near-infrared reflectance spectroscopy (NIRS) and mid-infrared reflectance spectroscopy (MIRS) techniques were used to detect six quality control indexes in the feeding and secondary heat treatment process of Reduning Injection (熱毒寧注射液, RI) preparation process, so as to improve the quality control level of preparation process.The quantitative correction models of chlorogenic acid, neochlorogenic acid, cryptochlorogenic acid, gardenoside, secoxyloganin and solid content were established by using NIRS technology and attenuated total MIRS technology, combined with partial least squares (PLS) method, after spectral pretreatment method selection and band selection respectively. Calibration set correlation coefficient (cal), validation set correlation coefficient (pre), root mean square error of calibration (RMSEC), root mean square error of cross validation (RMSECV) and relative standard error of prediction (RSEP) were used as evaluation indexes to evaluate the model performance.The RMSEC and RMSECV of the six quality control indexes of the NIRS prediction model were both less than 0.3 and close to each other, and the RSEP value were all less than 4.0%. The RMSEC and RMSECV of the six quality control indexes of the MIRS prediction model were all less than 0.4 and close to each other, and the RSEP value were all less than 5.0%. These PLS models had the advantages of good model performance and high prediction accuracy.Both NIRS and MIRS technology can be used for rapid detection of six quality control indexes in RI feeding-secondary heat treatment process, and the RSEP value is within 5%, and the method is simple and reliable.

    near-infrared reflectance spectroscopy; mid-infrared reflectance spectrum; Reduning Injection; preparation process; spectral pretreatment; band screening; partial least squares; chlorogenic acid; neochlorogenic acid; cryptochlorogenic acid; gardenoside;secoxyloganin; solid content

    R283.6

    A

    0253 - 2670(2022)21 - 6706 - 10

    10.7501/j.issn.0253-2670.2022.21.007

    2022-05-16

    2022年中央財(cái)政轉(zhuǎn)移支付地方項(xiàng)目:基于重點(diǎn)研究室研究領(lǐng)域的中醫(yī)藥多學(xué)科研究能力提升項(xiàng)目——中藥提取精制新技術(shù)

    童 楓(1998—),男,碩士研究生,研究方向?yàn)橹兴幹扑庍^(guò)程新技術(shù)。Tel: 18851193208 E-mail: tongfeng5324@126.com

    王振中,研究員級(jí)高級(jí)工程師,碩士生導(dǎo)師,研究方向?yàn)橹兴幹扑庍^(guò)程新技術(shù)。E-mail: kanionlunwen@163.com

    徐芳芳,女,博士,主管藥師,研究方向?yàn)橹兴幹扑庍^(guò)程新技術(shù)。E-mail: 879164331@qq.com

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

    猜你喜歡
    模型
    一半模型
    一種去中心化的域名服務(wù)本地化模型
    適用于BDS-3 PPP的隨機(jī)模型
    提煉模型 突破難點(diǎn)
    函數(shù)模型及應(yīng)用
    p150Glued在帕金森病模型中的表達(dá)及分布
    函數(shù)模型及應(yīng)用
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    3D打印中的模型分割與打包
    午夜福利网站1000一区二区三区| 校园人妻丝袜中文字幕| 男女免费视频国产| 欧美在线黄色| 天天躁夜夜躁狠狠躁躁| 不卡av一区二区三区| 国产高清国产精品国产三级| 国产男女内射视频| 亚洲国产欧美在线一区| 午夜av观看不卡| 成年美女黄网站色视频大全免费| 国产精品久久久人人做人人爽| 高清在线视频一区二区三区| 欧美亚洲日本最大视频资源| 亚洲欧美中文字幕日韩二区| 中文字幕高清在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 永久免费av网站大全| 黄色视频不卡| 少妇被粗大的猛进出69影院| 91成人精品电影| 久久久久久久国产电影| 国产精品免费视频内射| 老司机深夜福利视频在线观看 | 亚洲自偷自拍图片 自拍| 日本爱情动作片www.在线观看| 男女边吃奶边做爰视频| 中国三级夫妇交换| 精品一区二区三区四区五区乱码 | 国产精品国产三级国产专区5o| 午夜91福利影院| 天美传媒精品一区二区| 精品午夜福利在线看| 亚洲色图 男人天堂 中文字幕| 亚洲av日韩精品久久久久久密 | 在线观看免费高清a一片| 精品久久蜜臀av无| 最黄视频免费看| 亚洲天堂av无毛| 国产熟女午夜一区二区三区| 人成视频在线观看免费观看| 久久国产精品大桥未久av| 伦理电影大哥的女人| 高清视频免费观看一区二区| 亚洲精华国产精华液的使用体验| 久热这里只有精品99| 精品一品国产午夜福利视频| 欧美精品高潮呻吟av久久| 久久精品国产a三级三级三级| 三上悠亚av全集在线观看| 久久久久久免费高清国产稀缺| 超碰97精品在线观看| 国产精品免费视频内射| 丝袜脚勾引网站| 久久久久精品性色| 亚洲自偷自拍图片 自拍| 无限看片的www在线观看| 黄片播放在线免费| 午夜久久久在线观看| 中国三级夫妇交换| 别揉我奶头~嗯~啊~动态视频 | 9191精品国产免费久久| 黑人欧美特级aaaaaa片| 亚洲精品国产av成人精品| 又粗又硬又长又爽又黄的视频| 亚洲av成人精品一二三区| 操美女的视频在线观看| 免费人妻精品一区二区三区视频| 这个男人来自地球电影免费观看 | 一二三四中文在线观看免费高清| 自拍欧美九色日韩亚洲蝌蚪91| 精品一区二区三区四区五区乱码 | 秋霞在线观看毛片| xxxhd国产人妻xxx| 制服人妻中文乱码| 制服人妻中文乱码| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久成人aⅴ小说| av在线老鸭窝| 国产日韩欧美亚洲二区| 国产免费又黄又爽又色| 国产野战对白在线观看| 亚洲成人一二三区av| 亚洲专区中文字幕在线 | 99香蕉大伊视频| 久久精品久久久久久噜噜老黄| 日本色播在线视频| 亚洲熟女精品中文字幕| 精品国产乱码久久久久久小说| 极品人妻少妇av视频| 一级片'在线观看视频| 99久久精品国产亚洲精品| 国产一区二区在线观看av| 亚洲国产av影院在线观看| 精品一区二区免费观看| 色吧在线观看| 午夜久久久在线观看| 亚洲av电影在线观看一区二区三区| 久久综合国产亚洲精品| 黄色怎么调成土黄色| 久久影院123| 最近中文字幕高清免费大全6| 丰满迷人的少妇在线观看| 在线观看免费高清a一片| 99久久综合免费| 午夜日韩欧美国产| 亚洲婷婷狠狠爱综合网| 国语对白做爰xxxⅹ性视频网站| 色婷婷久久久亚洲欧美| 欧美日韩精品网址| 免费看av在线观看网站| 丝瓜视频免费看黄片| 极品少妇高潮喷水抽搐| 美女午夜性视频免费| 午夜免费男女啪啪视频观看| 久久久精品免费免费高清| 久久毛片免费看一区二区三区| 国产伦理片在线播放av一区| 亚洲,欧美精品.| 女人久久www免费人成看片| 欧美日韩亚洲综合一区二区三区_| 国产精品 国内视频| 午夜日本视频在线| 在线观看免费午夜福利视频| 日韩大片免费观看网站| 丰满迷人的少妇在线观看| 咕卡用的链子| 欧美精品亚洲一区二区| 亚洲五月色婷婷综合| 超碰97精品在线观看| 人人妻,人人澡人人爽秒播 | 激情五月婷婷亚洲| 色综合欧美亚洲国产小说| 多毛熟女@视频| 免费观看人在逋| 最近最新中文字幕大全免费视频 | 满18在线观看网站| 成人毛片60女人毛片免费| 亚洲av欧美aⅴ国产| 又大又爽又粗| 波多野结衣一区麻豆| 国产成人欧美在线观看 | 国产亚洲精品第一综合不卡| 在线观看免费视频网站a站| 亚洲国产精品999| 人妻一区二区av| 19禁男女啪啪无遮挡网站| 精品人妻熟女毛片av久久网站| 亚洲成色77777| 国产男女超爽视频在线观看| 国产免费现黄频在线看| 你懂的网址亚洲精品在线观看| 肉色欧美久久久久久久蜜桃| 久久99一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 亚洲欧美精品自产自拍| 欧美亚洲 丝袜 人妻 在线| 别揉我奶头~嗯~啊~动态视频 | 青春草国产在线视频| 久久久久久人人人人人| 日韩 欧美 亚洲 中文字幕| 日韩电影二区| 性高湖久久久久久久久免费观看| 欧美精品一区二区大全| 人人妻人人澡人人看| 男女边吃奶边做爰视频| av国产精品久久久久影院| 国精品久久久久久国模美| av女优亚洲男人天堂| 日韩一区二区三区影片| 国产精品蜜桃在线观看| 18禁观看日本| 成人午夜精彩视频在线观看| 久久久久精品人妻al黑| 免费观看av网站的网址| 亚洲一区二区三区欧美精品| 少妇人妻久久综合中文| 国产成人av激情在线播放| 午夜福利影视在线免费观看| 亚洲欧美精品综合一区二区三区| 亚洲成人一二三区av| 日日啪夜夜爽| 赤兔流量卡办理| 性高湖久久久久久久久免费观看| 免费高清在线观看日韩| 香蕉国产在线看| 欧美少妇被猛烈插入视频| 欧美乱码精品一区二区三区| 国产成人欧美| 久久毛片免费看一区二区三区| 欧美日韩亚洲综合一区二区三区_| 国产精品一国产av| 黑人猛操日本美女一级片| 中文字幕高清在线视频| 老司机靠b影院| 亚洲精品久久成人aⅴ小说| 国产一区二区激情短视频 | 午夜激情av网站| 欧美成人精品欧美一级黄| 成人午夜精彩视频在线观看| 国产高清不卡午夜福利| 欧美人与性动交α欧美精品济南到| 2018国产大陆天天弄谢| 中文字幕av电影在线播放| 在线免费观看不下载黄p国产| 久久久国产欧美日韩av| 亚洲精品美女久久av网站| 日韩大码丰满熟妇| 狂野欧美激情性xxxx| 免费观看性生交大片5| 9热在线视频观看99| 婷婷成人精品国产| 国产伦理片在线播放av一区| 日韩大码丰满熟妇| 欧美少妇被猛烈插入视频| 曰老女人黄片| 国产野战对白在线观看| av有码第一页| 中文字幕另类日韩欧美亚洲嫩草| 一区在线观看完整版| 9色porny在线观看| 一二三四在线观看免费中文在| 午夜免费观看性视频| 大话2 男鬼变身卡| 色婷婷av一区二区三区视频| 国产99久久九九免费精品| av女优亚洲男人天堂| 日本欧美视频一区| 啦啦啦 在线观看视频| 久久99精品国语久久久| 亚洲欧洲日产国产| 婷婷色综合www| 久久久久久人妻| 国产毛片在线视频| 又大又爽又粗| 性少妇av在线| 亚洲情色 制服丝袜| 最近的中文字幕免费完整| 精品视频人人做人人爽| 一级毛片黄色毛片免费观看视频| 午夜老司机福利片| 999精品在线视频| 精品一区二区三卡| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 丰满乱子伦码专区| 一级a爱视频在线免费观看| 久久婷婷青草| 在线看a的网站| 中国三级夫妇交换| 国产免费一区二区三区四区乱码| 女的被弄到高潮叫床怎么办| 一本一本久久a久久精品综合妖精| 婷婷色麻豆天堂久久| 国产精品 欧美亚洲| 热99国产精品久久久久久7| 嫩草影视91久久| 成人亚洲欧美一区二区av| 日韩伦理黄色片| 日本爱情动作片www.在线观看| 黄色 视频免费看| 99re6热这里在线精品视频| 国产 一区精品| 成人漫画全彩无遮挡| 一级a爱视频在线免费观看| 男女无遮挡免费网站观看| 中文字幕另类日韩欧美亚洲嫩草| 精品久久久久久电影网| 啦啦啦在线免费观看视频4| 9热在线视频观看99| 国产极品粉嫩免费观看在线| 国产国语露脸激情在线看| 校园人妻丝袜中文字幕| 999精品在线视频| 亚洲,欧美精品.| 青草久久国产| 在线观看国产h片| 亚洲欧美激情在线| 亚洲久久久国产精品| 日本av免费视频播放| 天堂8中文在线网| 黄色 视频免费看| 亚洲av电影在线进入| 十分钟在线观看高清视频www| 色播在线永久视频| 亚洲伊人色综图| 国产免费现黄频在线看| 国产成人免费观看mmmm| 亚洲欧美清纯卡通| 国产亚洲精品第一综合不卡| 国产成人精品久久二区二区91 | 久久久久久人人人人人| 国产精品女同一区二区软件| 国产老妇伦熟女老妇高清| 少妇精品久久久久久久| 国产成人精品福利久久| 亚洲av成人精品一二三区| 午夜91福利影院| 国产亚洲av片在线观看秒播厂| 最近中文字幕2019免费版| 国产成人欧美在线观看 | 国产成人精品久久久久久| av又黄又爽大尺度在线免费看| 日韩免费高清中文字幕av| 高清在线视频一区二区三区| 你懂的网址亚洲精品在线观看| 免费不卡黄色视频| 毛片一级片免费看久久久久| 少妇 在线观看| 十八禁高潮呻吟视频| 亚洲第一av免费看| 精品免费久久久久久久清纯 | 十八禁网站网址无遮挡| 99精品久久久久人妻精品| 国语对白做爰xxxⅹ性视频网站| 天天影视国产精品| 亚洲国产欧美一区二区综合| 中文字幕最新亚洲高清| 性色av一级| 亚洲av日韩精品久久久久久密 | 精品国产露脸久久av麻豆| 一级片'在线观看视频| 日韩中文字幕欧美一区二区 | 亚洲成人av在线免费| 可以免费在线观看a视频的电影网站 | 婷婷色综合www| 青春草视频在线免费观看| 婷婷色麻豆天堂久久| 午夜免费鲁丝| 亚洲美女搞黄在线观看| 免费女性裸体啪啪无遮挡网站| 中文精品一卡2卡3卡4更新| 男人爽女人下面视频在线观看| 日韩欧美一区视频在线观看| 青草久久国产| 色播在线永久视频| 精品视频人人做人人爽| 欧美精品人与动牲交sv欧美| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区蜜桃| 你懂的网址亚洲精品在线观看| 国产片内射在线| 国产无遮挡羞羞视频在线观看| 亚洲国产日韩一区二区| 久久精品久久久久久噜噜老黄| 国语对白做爰xxxⅹ性视频网站| 秋霞伦理黄片| 成人手机av| 久久国产亚洲av麻豆专区| 波多野结衣av一区二区av| 午夜福利网站1000一区二区三区| 一级毛片电影观看| 如何舔出高潮| 90打野战视频偷拍视频| 一边亲一边摸免费视频| 美女视频免费永久观看网站| 欧美日韩综合久久久久久| 亚洲欧洲日产国产| 黑人巨大精品欧美一区二区蜜桃| 51午夜福利影视在线观看| 一级毛片电影观看| 热re99久久精品国产66热6| 色94色欧美一区二区| 熟女少妇亚洲综合色aaa.| 超碰97精品在线观看| 国产男女超爽视频在线观看| 亚洲精品乱久久久久久| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 精品一区二区三区四区五区乱码 | 国产亚洲欧美精品永久| 国产日韩欧美视频二区| 国产一区二区三区综合在线观看| 国产男人的电影天堂91| 热re99久久精品国产66热6| 美女脱内裤让男人舔精品视频| 2018国产大陆天天弄谢| 女的被弄到高潮叫床怎么办| 免费观看av网站的网址| 黄色怎么调成土黄色| 久久久国产欧美日韩av| 亚洲精品,欧美精品| 欧美人与性动交α欧美软件| 亚洲欧美成人精品一区二区| 国产国语露脸激情在线看| 欧美日韩精品网址| 成人免费观看视频高清| 日本猛色少妇xxxxx猛交久久| 欧美 日韩 精品 国产| 国产免费福利视频在线观看| 亚洲国产精品一区三区| 国产乱人偷精品视频| 18在线观看网站| 一区二区三区乱码不卡18| 蜜桃在线观看..| 精品一区二区免费观看| 一级黄片播放器| 精品酒店卫生间| 午夜久久久在线观看| 欧美日韩亚洲高清精品| 永久免费av网站大全| 久久热在线av| 精品亚洲成a人片在线观看| 国产精品国产三级专区第一集| 亚洲国产最新在线播放| 午夜福利,免费看| 国产一区亚洲一区在线观看| 999精品在线视频| 国产成人精品无人区| 日日摸夜夜添夜夜爱| 精品一区二区三卡| 亚洲欧美成人综合另类久久久| 日韩一区二区三区影片| 国产免费现黄频在线看| 久久国产精品大桥未久av| 日韩不卡一区二区三区视频在线| 亚洲欧美精品综合一区二区三区| 美女脱内裤让男人舔精品视频| 成人影院久久| 极品少妇高潮喷水抽搐| 青草久久国产| 曰老女人黄片| 午夜福利免费观看在线| 成人影院久久| 国产精品麻豆人妻色哟哟久久| 波多野结衣一区麻豆| 9热在线视频观看99| 亚洲av电影在线进入| 国产日韩一区二区三区精品不卡| av在线app专区| 国产精品无大码| www日本在线高清视频| 国产欧美亚洲国产| 国产精品免费大片| 成年av动漫网址| 国产无遮挡羞羞视频在线观看| 少妇被粗大的猛进出69影院| 国产麻豆69| 一边摸一边做爽爽视频免费| 国产伦人伦偷精品视频| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品在线美女| 亚洲av日韩在线播放| 国精品久久久久久国模美| kizo精华| 精品午夜福利在线看| 久久人人97超碰香蕉20202| 婷婷色麻豆天堂久久| 不卡av一区二区三区| 午夜福利免费观看在线| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 别揉我奶头~嗯~啊~动态视频 | 日韩av不卡免费在线播放| 亚洲欧美清纯卡通| 黄色视频不卡| 久久久久久久久免费视频了| 久久久精品免费免费高清| 桃花免费在线播放| 国产精品蜜桃在线观看| 亚洲一区二区三区欧美精品| 国语对白做爰xxxⅹ性视频网站| 十八禁网站网址无遮挡| 免费av中文字幕在线| 女人高潮潮喷娇喘18禁视频| 只有这里有精品99| 老司机在亚洲福利影院| 你懂的网址亚洲精品在线观看| 啦啦啦视频在线资源免费观看| 亚洲视频免费观看视频| 一边亲一边摸免费视频| 老熟女久久久| 亚洲色图 男人天堂 中文字幕| av女优亚洲男人天堂| 天天添夜夜摸| 综合色丁香网| av视频免费观看在线观看| 国产成人欧美在线观看 | 免费观看性生交大片5| 人妻 亚洲 视频| 免费不卡黄色视频| 国产爽快片一区二区三区| 久久性视频一级片| 秋霞在线观看毛片| 国产福利在线免费观看视频| 又黄又粗又硬又大视频| 90打野战视频偷拍视频| 亚洲情色 制服丝袜| 国产精品国产av在线观看| 多毛熟女@视频| 国产精品麻豆人妻色哟哟久久| 亚洲av成人精品一二三区| 欧美亚洲日本最大视频资源| 一级,二级,三级黄色视频| 老司机影院成人| 午夜福利在线免费观看网站| 一区二区三区四区激情视频| 亚洲国产精品999| 久久久国产一区二区| 日本91视频免费播放| bbb黄色大片| 美女视频免费永久观看网站| 色精品久久人妻99蜜桃| 日韩免费高清中文字幕av| 性高湖久久久久久久久免费观看| 精品一区二区三区av网在线观看 | www.自偷自拍.com| av福利片在线| 国产成人欧美| 亚洲美女黄色视频免费看| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| av国产久精品久网站免费入址| 女人久久www免费人成看片| 日韩成人av中文字幕在线观看| 蜜桃在线观看..| 伦理电影免费视频| avwww免费| 日韩av免费高清视频| 97在线人人人人妻| 哪个播放器可以免费观看大片| 黄网站色视频无遮挡免费观看| 午夜日本视频在线| 欧美精品亚洲一区二区| 男女免费视频国产| 亚洲自偷自拍图片 自拍| 叶爱在线成人免费视频播放| 王馨瑶露胸无遮挡在线观看| 性高湖久久久久久久久免费观看| 久久99热这里只频精品6学生| 欧美人与善性xxx| 久久人人爽av亚洲精品天堂| 黄片小视频在线播放| 久久人人爽av亚洲精品天堂| 久久亚洲国产成人精品v| 亚洲国产精品一区三区| 午夜福利,免费看| 精品一区二区三区av网在线观看 | 男女无遮挡免费网站观看| 国产男女内射视频| 午夜免费观看性视频| 欧美另类一区| h视频一区二区三区| 亚洲国产中文字幕在线视频| 在线观看免费视频网站a站| 夫妻午夜视频| 亚洲视频免费观看视频| 成人漫画全彩无遮挡| 国产精品二区激情视频| 日本av手机在线免费观看| 国产在线视频一区二区| 精品少妇久久久久久888优播| 99国产精品免费福利视频| 一区二区三区四区激情视频| 五月开心婷婷网| 99九九在线精品视频| 一级毛片黄色毛片免费观看视频| 欧美 亚洲 国产 日韩一| 国产成人精品福利久久| 久久99热这里只频精品6学生| 亚洲av日韩精品久久久久久密 | 最近最新中文字幕免费大全7| 国产一区有黄有色的免费视频| 成人国语在线视频| 国产日韩欧美在线精品| 男女边吃奶边做爰视频| 久久亚洲国产成人精品v| 久久精品久久精品一区二区三区| 看十八女毛片水多多多| 丰满少妇做爰视频| 大话2 男鬼变身卡| 黄色视频在线播放观看不卡| 亚洲精品视频女| 国产亚洲一区二区精品| 只有这里有精品99| 晚上一个人看的免费电影| 精品一区二区免费观看| 久久久国产欧美日韩av| 国产男人的电影天堂91| 热99久久久久精品小说推荐| 黄色怎么调成土黄色| 国产在线免费精品| 国产1区2区3区精品| 丝袜美腿诱惑在线| a级毛片黄视频| 丝袜人妻中文字幕| 亚洲国产欧美日韩在线播放| 日本wwww免费看| 无遮挡黄片免费观看| 日韩中文字幕视频在线看片| 亚洲精品久久成人aⅴ小说| 久久影院123| 大香蕉久久网| 国语对白做爰xxxⅹ性视频网站| 国产成人精品无人区| 一级黄片播放器| 别揉我奶头~嗯~啊~动态视频 | 哪个播放器可以免费观看大片| 天堂俺去俺来也www色官网| 香蕉国产在线看| 国产一级毛片在线| 最黄视频免费看| 91成人精品电影| 日韩大片免费观看网站| 天天躁夜夜躁狠狠躁躁| av福利片在线| 极品少妇高潮喷水抽搐| 亚洲欧美成人综合另类久久久| 亚洲精品中文字幕在线视频| 国产熟女午夜一区二区三区| 亚洲国产精品国产精品|