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

    藥品真空冷凍干燥過(guò)程監(jiān)控技術(shù)研究進(jìn)展

    2015-11-13 00:30:41
    化工進(jìn)展 2015年8期
    關(guān)鍵詞:模擬冷凍干燥控制

    ?

    藥品真空冷凍干燥過(guò)程監(jiān)控技術(shù)研究進(jìn)展

    李俊奇,李保國(guó)

    (上海理工大學(xué)能源與動(dòng)力工程學(xué)院,上海200093)

    摘要:對(duì)藥品凍干過(guò)程進(jìn)行優(yōu)化的關(guān)鍵是在保證藥品質(zhì)量不受損害的情況下盡量縮短干燥時(shí)間。因此,對(duì)凍干過(guò)程進(jìn)行準(zhǔn)確的監(jiān)控是十分重要的,既要保證藥品的溫度保持在合理的范圍內(nèi),對(duì)干燥結(jié)束時(shí)間進(jìn)行準(zhǔn)確地判斷,同時(shí)又要對(duì)凍干過(guò)程壓力和溫度進(jìn)行良好的控制以達(dá)到凍干過(guò)程的最優(yōu)化。本文對(duì)近年來(lái)藥品真空冷凍干燥過(guò)程監(jiān)控技術(shù)的研究進(jìn)展進(jìn)行了綜述,主要有基于動(dòng)態(tài)參數(shù)估計(jì)法(DPE)的監(jiān)控系統(tǒng)、基于卡爾曼濾波法的監(jiān)測(cè)系統(tǒng)、露點(diǎn)法判斷一次干燥結(jié)束點(diǎn)、模型預(yù)測(cè)控制法(MPC)。提出藥品真空冷凍干燥監(jiān)控技術(shù)的研究應(yīng)著重于以下幾點(diǎn):考慮輻射、對(duì)流和導(dǎo)熱3種傳熱方式在凍干傳熱過(guò)程中所占的比重,建立二維、三維凍干模型以更加精確地監(jiān)測(cè)藥品凍干過(guò)程的參數(shù),在此基礎(chǔ)上研究對(duì)加熱隔板溫度和凍干室壓力的實(shí)時(shí)最優(yōu)控制策略,以對(duì)藥品凍干過(guò)程進(jìn)行及時(shí)、有效地控制。

    關(guān)鍵詞:冷凍干燥;一次干燥;模擬;監(jiān)測(cè);控制

    第一作者:李俊奇(1990—),男,碩士研究生,主要從事生物藥品冷凍干燥、制冷及低溫方向的研究。聯(lián)系人:李保國(guó),教授,博士生導(dǎo)師。E-mail lbaoguo@126.com。

    凍干是生物藥品保存的重要方法,在制藥領(lǐng)域備受關(guān)注[1]。由于凍干過(guò)程需要先將物料凍結(jié),然后在一定真空條件下升華干燥,因此干燥時(shí)間長(zhǎng),能量消耗大。對(duì)凍干工藝過(guò)程進(jìn)行最優(yōu)化研究,利用盡可能短的時(shí)間得到高品質(zhì)的凍干產(chǎn)品,是冷凍干燥研究的熱點(diǎn)[2]。有關(guān)試驗(yàn)研究分析表明,在凍干過(guò)程中,產(chǎn)品溫度每升高一度就能使一次干燥時(shí)間減少13%以上[3]。因此,提高干燥效率的關(guān)鍵,就是在保證凍干物料質(zhì)量的前提下,最大限度地提高其溫度。

    監(jiān)控真空冷凍干燥過(guò)程要掌握的主要參數(shù)為:產(chǎn)品溫度,升華界面移動(dòng)位置。對(duì)凍干產(chǎn)品溫度的監(jiān)測(cè)尤為重要,因?yàn)橐刂苾龈蛇^(guò)程溫度不超過(guò)一定溫度(對(duì)非晶體溶質(zhì)系統(tǒng)指“崩解溫度”或“玻璃化轉(zhuǎn)變溫度”,對(duì)結(jié)晶溶質(zhì)系統(tǒng)指“共晶溫度”[2]),防止溫度過(guò)高而對(duì)產(chǎn)品品質(zhì)造成不可逆損傷。升華界面是瓶中的藥品在一次干燥過(guò)程中,位于干燥層和凍結(jié)層之間的移動(dòng)界面,隨著一次干燥的進(jìn)行,升華界面由凍干藥品的頂部逐漸移動(dòng)至底部。升華界面的位置直接反映了一次干燥的進(jìn)程,可由此來(lái)判斷冰升華的結(jié)束點(diǎn),確定該界面的位置十分重要。本文主要論述從2009年以來(lái)國(guó)內(nèi)外在藥品凍干過(guò)程監(jiān)測(cè)與控制的研究進(jìn)展以及發(fā)展趨勢(shì)。

    1 凍干過(guò)程監(jiān)控技術(shù)的研究進(jìn)展

    1.1動(dòng)態(tài)參數(shù)估計(jì)法(DPE)的監(jiān)控系統(tǒng)

    前人對(duì)基于壓力升測(cè)試(pressure rise test,PRT)計(jì)算凍干過(guò)程參數(shù)的方法進(jìn)行了研究,主要有氣壓溫度測(cè)量法(barometric temperature measurement,BTM)、壓力溫度測(cè)量法(manometric temperature measurement,MTM)、動(dòng)態(tài)壓力升法(dynamic pressure rise,DPR)、壓力升分析法(pressure rise analysis,PRA)[4]。動(dòng)態(tài)參數(shù)估計(jì)法(dynamic parameters estimation,DPE)是在前人研究基礎(chǔ)上產(chǎn)生的應(yīng)用于凍干過(guò)程監(jiān)測(cè)的新方法[5-6]。該方法通過(guò)一系列偏微分方程建立了小瓶?jī)?nèi)凍干藥品熱質(zhì)傳遞的一維非穩(wěn)態(tài)傳熱模型。采用總傳熱系數(shù)Kv將凍干過(guò)程發(fā)生的對(duì)流、導(dǎo)熱、輻射等因素表現(xiàn)為熱傳導(dǎo)過(guò)程,其計(jì)算過(guò)程如圖1所示。設(shè)定升華界面溫度Ti、有效擴(kuò)散系數(shù)k1和修正系數(shù)γ(用以修正各小瓶?jī)?nèi)冰升華表面積的不均勻性)的初始估計(jì)值,通過(guò)建立的數(shù)學(xué)模型計(jì)算得出凍干室壓力pc的模擬值,與實(shí)際測(cè)得參數(shù)pc,meas擬合,最終可得出凍結(jié)層厚度Lfrozen、凍結(jié)層內(nèi)部溫度分布T(z,t)等參數(shù)。采用MATLAB軟件可解決該非線性最小二乘擬合問(wèn)題。

    圖1 DPE計(jì)算流程圖

    利用該DPE方法模擬測(cè)得的參數(shù)相比之前的研究方法可靠性更高[5]。Barresi等[7-8]開(kāi)發(fā)了基于DPE模型的測(cè)量系統(tǒng):LyoMonitor和相對(duì)應(yīng)的控制軟件LyoDriver。

    LyoMonitor是一個(gè)創(chuàng)新的、模塊化的系統(tǒng),集成了壓力計(jì)、濕度傳感器和質(zhì)譜儀等測(cè)量設(shè)備和DPE模型,用于監(jiān)測(cè)藥瓶?jī)?nèi)藥品的一次干燥過(guò)程:計(jì)算其升華界面溫度和位置、產(chǎn)品內(nèi)部溫度分布和熱質(zhì)傳遞系數(shù)。Barresi等[7]采用多點(diǎn)無(wú)線測(cè)溫儀對(duì)這些參數(shù)進(jìn)行測(cè)量,并對(duì)LyoMonitor和多點(diǎn)無(wú)線測(cè)溫儀用于確定一次干燥結(jié)束點(diǎn)的效果作了對(duì)比,說(shuō)明這些傳感器都可以用于控制和最優(yōu)化凍干過(guò)程。

    LyoDriver是基于DPE模型開(kāi)發(fā)的控制軟件,該軟件可以用來(lái)確定一次干燥過(guò)程凍干機(jī)加熱隔板的最佳設(shè)定溫度,可以在不超過(guò)產(chǎn)品極限溫度的情況下實(shí)現(xiàn)最快的干燥。Barresi等[8]將LyoDriver控制軟件在工業(yè)樣機(jī)上進(jìn)行試驗(yàn),驗(yàn)證了其有效性。

    DPE系統(tǒng)控制原理如圖2所示。在樣品凍干過(guò)程中執(zhí)行壓力升測(cè)試(PRT),將得到的凍干室壓力pc隨時(shí)間t的變化數(shù)據(jù)進(jìn)行DPE計(jì)算,得到產(chǎn)品溫度T、有效傳熱系數(shù)Kv、和傳質(zhì)系數(shù)k1,然后將計(jì)算數(shù)據(jù)和其他凍干過(guò)程變量(導(dǎo)熱流體溫度、凍干室壓力、凍干室加熱或冷卻速率)配合LyoDriver計(jì)算得到最佳設(shè)定流體溫度Tfluid。

    在國(guó)內(nèi),郭彥偉[9]根據(jù)Velardi提出的DPE數(shù)學(xué)模型,利用MATLAB與VB進(jìn)行混合編程,得出可在凍干機(jī)上運(yùn)行計(jì)算的動(dòng)壓測(cè)溫軟件。該軟件可通過(guò)對(duì)壓力升數(shù)據(jù)進(jìn)行模擬,計(jì)算產(chǎn)品溫度,并在樣機(jī)上對(duì)部分藥品的凍干過(guò)程進(jìn)行了實(shí)驗(yàn)驗(yàn)證,但在生產(chǎn)上使用有待進(jìn)一步完善。

    圖2 DPE系統(tǒng)控制流程圖

    1.2基于卡爾曼濾波算法的監(jiān)測(cè)系統(tǒng)

    Velardi等[10]提出將擴(kuò)展卡爾曼濾波算法(extended kalman filter,EKF)用于凍干過(guò)程的優(yōu)化,并開(kāi)發(fā)了相應(yīng)的軟傳感器。Bosca和Fissore[11]開(kāi)發(fā)了卡爾曼濾波型監(jiān)測(cè)器,提出的測(cè)量模型能夠測(cè)量出一次干燥過(guò)程中產(chǎn)品的溫度變化、隔板與樣品之間的總傳熱系數(shù)以及水蒸氣流過(guò)凍干層的傳質(zhì)阻力。對(duì)具有不同干燥層厚度的樣品進(jìn)行測(cè)試,都得到了準(zhǔn)確的結(jié)果,并在中試凍干機(jī)上對(duì)理論結(jié)果進(jìn)行了試驗(yàn)驗(yàn)證。該方法的缺點(diǎn)是采用熱電偶測(cè)量溫度,不利于藥品凍干過(guò)程的熱質(zhì)傳遞和無(wú)菌化處理[12]。

    Bosca等[13]提出了智能軟傳感器(Smart Soft-Sensor,S3)新測(cè)量工具,工作原理如圖3所示。該系統(tǒng)包含基于卡爾曼濾波算法的監(jiān)測(cè)器(Observer)和數(shù)學(xué)模型(Model),用于監(jiān)測(cè)整個(gè)一次干燥過(guò)程。由于升華干燥后期,部分小瓶?jī)?nèi)的冰已經(jīng)升華完畢導(dǎo)致整個(gè)升華界面面積的減小,與升華干燥的前期存在差異,故S3系統(tǒng)分別用Observer和Model來(lái)監(jiān)測(cè)升華干燥的前期和后期。樣品底部溫度TB用熱電偶測(cè)量,根據(jù)隔板與樣品間的總傳熱系數(shù)Kv初始估計(jì)值和TB通過(guò)Observer計(jì)算得到Kv和升華界面溫度Ti,進(jìn)而得出凍干層厚度Ld、升華通量Jw、凍干層的蒸汽流通阻力Rp;在升華干燥后期,將求得的Kv和Rp用于Model求出Ti和Ld。并采用含有不同輔料[蔗糖,聚乙烯吡咯烷酮(PVP)]的水溶液在不同操作條件下進(jìn)行試驗(yàn)驗(yàn)證該系統(tǒng)的準(zhǔn)確性,結(jié)果表明該軟件能夠?qū)崿F(xiàn)模型參數(shù)和產(chǎn)品特性的快速估算,不需要較貴的硬件或耗時(shí)的分析;通過(guò)實(shí)驗(yàn)證明了使用該軟件能在保證產(chǎn)品質(zhì)量的條件下,顯著縮短干燥時(shí)間。

    圖3 S3 系統(tǒng)工作原理圖

    Bosca等[14]將該方法用于計(jì)算干燥層蒸汽流動(dòng)阻力,進(jìn)而在脫機(jī)狀態(tài)下計(jì)算出給定條件下產(chǎn)品的干燥時(shí)間及溫度,從而實(shí)現(xiàn)凍干過(guò)程的最優(yōu)化。若采用控制成核的方法(這樣一來(lái),插有熱電偶測(cè)溫的小瓶?jī)?nèi)樣品和其他未使用熱電偶測(cè)溫的小瓶?jī)?nèi)水分?jǐn)U散動(dòng)態(tài)是相同的),同時(shí)采用無(wú)線傳感器來(lái)測(cè)量樣品溫度,該系統(tǒng)也可用于符合藥品生產(chǎn)管理規(guī)范(GMP標(biāo)準(zhǔn))的工業(yè)生產(chǎn)規(guī)模的凍干機(jī)。

    1.3露點(diǎn)法判斷一次干燥結(jié)束點(diǎn)

    將露點(diǎn)儀用于確定升華干燥結(jié)束點(diǎn)的原理是通過(guò)凍干過(guò)程中干燥室內(nèi)露點(diǎn)的變化反映出水蒸氣分壓力(水蒸氣含量)的變化,而干燥室內(nèi)水蒸氣含量的增加來(lái)自于凍干樣品中水蒸氣的不斷逸出,由此可判斷干燥進(jìn)程。

    張今等[15]研究了用露點(diǎn)法監(jiān)測(cè)凍干過(guò)程物料的水分遷移,實(shí)驗(yàn)檢驗(yàn)了露點(diǎn)法在監(jiān)測(cè)整批樣品凍干過(guò)程的靈敏性,結(jié)果表明:即使在小瓶?jī)?nèi)藥液填充量不均一或溶液高度很小的情況下露點(diǎn)法仍有足夠的靈敏度,可以清楚地指示一次升華干燥的終點(diǎn)。說(shuō)明露點(diǎn)法是一種可靠的監(jiān)測(cè)一次升華結(jié)束點(diǎn)的方法。

    露點(diǎn)法測(cè)量要用到露點(diǎn)儀,露點(diǎn)儀應(yīng)用于凍干過(guò)程的缺點(diǎn)在于它不能用蒸汽殺菌,然而可以用過(guò)濾器和閥件將露點(diǎn)儀與蒸汽隔離進(jìn)行消毒,在蒸汽消毒的時(shí)候?qū)㈤y關(guān)閉,便可用于實(shí)驗(yàn)室和生產(chǎn)型凍干機(jī)[16]。

    1.4模型預(yù)測(cè)控制法(MPC)

    Daraoui等[17]提出了將模型預(yù)測(cè)控制法(model predictive control,MPC)應(yīng)用于凍干工藝過(guò)程優(yōu)化,并用MATLAB開(kāi)發(fā)了控制軟件MPC@CB。該軟件是基于改進(jìn)的Levenberg-Marquardt算法的控制軟件,可以控制開(kāi)環(huán)或閉環(huán)的連續(xù)過(guò)程,并找到最佳

    約束條件[17-18]。

    Pisano等[19]對(duì)MPC的研究進(jìn)行了分析,指出其不足是控制變量只有熱流體的溫度,凍干室壓力不變;并研究了同時(shí)控制隔板溫度和干燥室壓力來(lái)縮短藥品一次干燥時(shí)間,提出了可應(yīng)用于凍干機(jī)的控制方法。

    2 其他主要模型與監(jiān)控方法

    近年來(lái),在凍干監(jiān)控方面也出現(xiàn)了許多新方法,并且在凍干應(yīng)用上都有所突破。

    Smith等[20]采用阻抗譜方法監(jiān)測(cè)冷凍干燥過(guò)程。這是一種微創(chuàng)的新型過(guò)程分析技術(shù):采用平面電極安裝到凍干樣品的西林瓶外壁,并連接到高阻抗分析儀,將玻璃瓶壁與產(chǎn)品界面的復(fù)合阻抗產(chǎn)生的偽弛豫過(guò)程記錄下來(lái),所測(cè)結(jié)果與樣品溫度和冰升華程度強(qiáng)烈耦合,可由此來(lái)測(cè)量?jī)龈蓸悠犯稍锼俾屎鸵淮胃稍锝Y(jié)束點(diǎn)。該方法與熱電偶測(cè)量相比的明顯優(yōu)勢(shì)在于它不僅能用于檢測(cè)一次干燥結(jié)束點(diǎn),同時(shí)可以測(cè)量干燥速率;提出的阻抗測(cè)量系統(tǒng)可作為過(guò)程控制系統(tǒng)的組成部分,對(duì)干燥速率、擱板溫度和凍干結(jié)束點(diǎn)進(jìn)行監(jiān)測(cè)。

    Dragoi等[21]提出了基于自適應(yīng)差分進(jìn)化算法與反向傳播算法相結(jié)合的“黑箱”模型新方法,用于神經(jīng)網(wǎng)絡(luò)模型的結(jié)構(gòu)和參數(shù)同時(shí)優(yōu)化。在已知隔板加熱溫度和凍干室壓力條件下,可在脫機(jī)狀況下算出產(chǎn)品溫度和剩余冰的含量隨時(shí)間的變化,從而得出產(chǎn)品溫度和升華干燥結(jié)束時(shí)間;而且,“黑箱”模型可以用來(lái)監(jiān)控凍干過(guò)程,將測(cè)得的產(chǎn)品溫度作為神經(jīng)網(wǎng)絡(luò)的輸入變量來(lái)實(shí)現(xiàn)產(chǎn)品狀態(tài)(溫度和剩余冰含量)的實(shí)時(shí)測(cè)量。這為凍干設(shè)計(jì)和優(yōu)化提供了一個(gè)有價(jià)值的工具。

    近紅外光譜技術(shù)(near infrared spectroscopy,NIR)可以用來(lái)監(jiān)測(cè)凍干過(guò)程的物理化學(xué)變化,在凍干過(guò)程優(yōu)化和控制方面富有前景。Beer等[22]將NIR與和拉曼光譜儀(Raman spectroscopy)相結(jié)合用于凍干過(guò)程的監(jiān)控,結(jié)果表明在實(shí)驗(yàn)條件下,該方法可以有效監(jiān)測(cè)凍干過(guò)程參數(shù),并可用于凍干過(guò)程的最優(yōu)化。Rosas等[23]將NIR用于多組分藥物制劑的凍干在線監(jiān)測(cè),結(jié)果表明,在低溫環(huán)境下利用近紅外光譜技術(shù)實(shí)現(xiàn)在線監(jiān)控是可行的。

    3 結(jié)語(yǔ)

    為了能在確保生物藥品冷凍干燥質(zhì)量的同時(shí)縮短凍干時(shí)間,提高凍干效率,真空冷凍干燥過(guò)程監(jiān)控技術(shù)一直是研究熱點(diǎn),相關(guān)研究成果不斷發(fā)表,并且已經(jīng)開(kāi)發(fā)了部分有效的監(jiān)控軟件。隨著國(guó)內(nèi)凍干事業(yè)的發(fā)展,凍干過(guò)程監(jiān)控技術(shù)也逐漸得到重視,研究的重點(diǎn)由凍干數(shù)學(xué)模型[24]逐步轉(zhuǎn)向凍干過(guò)程監(jiān)控技術(shù)。

    基于DPE模型的方法能夠有效地計(jì)算產(chǎn)品溫度、升華界面位置和傳遞系數(shù),是凍干監(jiān)測(cè)的主要方法[7]。與多數(shù)基于PRT的測(cè)量方法能夠計(jì)算出系統(tǒng)狀態(tài)的平均值相比,EKF技術(shù)得到的是單一小瓶的狀態(tài)值,可通過(guò)同時(shí)監(jiān)測(cè)不同位置多個(gè)小瓶來(lái)減小產(chǎn)品不均一問(wèn)題;將EKF與DPE方法結(jié)合可使計(jì)算更快地收斂于精確解,從而顯著改善監(jiān)測(cè)軟件的性能,也可通過(guò)該方法將收集到的信號(hào)實(shí)時(shí)傳送到凍干機(jī)外的數(shù)據(jù)收集系統(tǒng),實(shí)現(xiàn)對(duì)凍干小瓶溫度的非接觸測(cè)量,適用于配備了自動(dòng)裝卸載系統(tǒng)的大型工業(yè)化凍干機(jī)[10]。露點(diǎn)儀應(yīng)用于凍干過(guò)程監(jiān)測(cè)靈敏度高且不對(duì)凍干過(guò)程造成影響,并可以監(jiān)測(cè)一次和二次干燥的結(jié)束點(diǎn)[16]。MPC法可實(shí)現(xiàn)對(duì)小規(guī)?;蚬I(yè)規(guī)模凍干機(jī)的控制[19],是富有前景的控制方法。

    在對(duì)真空冷凍干燥監(jiān)控技術(shù)研究不斷完善的同時(shí),對(duì)凍干工藝過(guò)程模擬研究不斷深入:在一維模型的基礎(chǔ)上,考慮瓶壁輻射及熱量平衡[25]、壓力升測(cè)試的中隔閥關(guān)閉時(shí)間、惰性氣體的導(dǎo)熱等多因素,建立二維、三維凍干模型以更加精確地監(jiān)測(cè)藥品凍干過(guò)程參數(shù),在此基礎(chǔ)上研究?jī)龈蛇^(guò)程的實(shí)時(shí)最優(yōu)控制策略是當(dāng)前研究的熱點(diǎn)。總之,對(duì)于凍干過(guò)程的監(jiān)測(cè)與控制,在方法的完善及發(fā)展和實(shí)際應(yīng)用方面還有待繼續(xù)探索。

    符號(hào)說(shuō)明

    Jw——冰升華通量,kg/(m2?s)k1——有效擴(kuò)散系數(shù),m2/s

    Kv——總傳熱系數(shù),J/(m2?s?K)

    Ld——干燥層厚度,m

    Lfrozen——凍結(jié)層厚度,m

    pc——凍干室壓力,Pa

    pc,means——凍干室壓力實(shí)測(cè)值,Pa

    Rp——干燥層水蒸氣流通阻力,m/s

    T——溫度,K

    TB——凍干小瓶?jī)?nèi)樣品底部溫度,K

    Tfluid——隔板加熱或冷卻溫度,K

    Ti——升華界面溫度,K

    Ti0——壓力升測(cè)試開(kāi)始時(shí)的升華界面溫度,K

    T|t=0——壓力升測(cè)試開(kāi)始時(shí)凍結(jié)層溫度分布,K

    T(z,t)——溫度隨時(shí)間和z軸的一維分布

    t——時(shí)間,s

    γ——修正系數(shù)

    參考文獻(xiàn)

    [1]Kasper J C,Winter G,F(xiàn)riess W.Recent advances and further challenges in lyophilization[J].European Journal of Pharmaceutics and Biopharmaceutics,2013,85:162-169.

    [2]Tang X C,Pikal M J.Design of freeze-drying processes for pharmaceuticals:Practical advice[J].Pharmaceutical Research,2004,21:191-200.

    [3]Costantino H R,Pikal M J.Lyophilization of Biopharmaceuticals(Biotechnology:Pharmaceutical Aspects)[M].Virginia:American Association of Pharmaceutical Scientists,2005.

    [4]Fissore D,Pisano R,Velardi S,et al.PAT tools for the optimization of the freeze-drying process[J].Pharmaceutical Engineering,2009,29(5):58-68.

    [5]Velardi S A,Rasetto V,Barresi AA.Dynamic parameters estimation method:Advanced manometric temperature measurement approach for freeze-drying monitoringof pharmaceutical solutions[J].Industrial and Engineering Chemistry Research,2008,47(21):8445-8457.

    [6]Fissore D,Pisano R,Rasetto V,et al.Applying process analytical technology(PAT)tolyophilizationprocesses[J].Analytical Technologies,2009,27(2):7-11.

    [7]Barresi AA,Pisano R,F(xiàn)issore D,et al.Monitoring of the primary drying of a lyophilization process in vials[J].Chemical Engineering and Processing,2009,48:408-423.

    [8]Barresi A A,Velardi S A,Pisano R.In-line control of the lyophilization process.A gentle PAT approach using software sensors[J].International Journal of Refrigeration,2009,32:1003-1014.

    [9]郭彥偉.真空冷凍干燥動(dòng)壓測(cè)溫技術(shù)研究及應(yīng)用軟件開(kāi)發(fā)[D].上海:上海理工大學(xué),2014.

    [10]Velardi S A,Hammouri H,Barresi A A.In-line monitoring of the primary drying phase of the freeze-drying process in vial by means of a Kalman filter based observer[J].Chemical Engineering Research and Design,2009,87:1409-1419.

    [11]Bosca S,F(xiàn)issore D.Design and validation of an innovative soft-sensorforpharmaceuticalsfreeze-dryingmonitoring[J].Chemical Engineering Science,2011,66:5127-5136.

    [12]華澤釗.冷凍干燥新技術(shù)[M].北京:科學(xué)出版社,2005:244-245.

    [13]Bosca S,Barresi AA,F(xiàn)issore D.Fast freeze-drying cycle design and optimization using a PAT based on the measurement of product temperature[J].EuropeanJournalofPharmaceuticsand Biopharmaceutics,2013,85:253-262.

    [14]Bosca S,Barresi A A,F(xiàn)issore D.Use of a soft sensor for the fast estimation of dried cake resistance during a freeze-drying cycle[J].International Journal of Pharmaceutics,2013,451:23-33.

    [15]張今,周?chē)?guó)燕,郭柏松,等.基于露點(diǎn)法判斷凍干過(guò)程的一次升華結(jié)束點(diǎn)[J].真空,2009,46(4):66-71.

    [16]Patel S M,Pikal M.Process analytical technologies(PAT)in freeze-drying of parenteral products[J].Pharmaceutical Development and Technology,2009,14(6):567-587.

    [17]Daraoui N,Dufour P,Hammouri H,et al.Optimal operation of sublimation time of the freeze drying process by predictive control:Application of the MPC@CB software[J].Computer Aided Chemical Engineering,2008,25:453-458.

    [18]Daraoui N,Dufour P,Hammouri H,et al.Model predictive control during the primary drying stage of lyophilisation[J].Control Engineering Practice,2010,18:483-494.

    [19]Pisano R,F(xiàn)issore D,Barresi AA.Freeze-drying cycle optimization using model predictive control techniques[J].Industrial and Engineering Chemistry Research,2011,50:7363-7379.

    [20]Smith G,Polygalov E,Arshad M S,et al.An impedance-based process analytical technology for monitoring the lyophilisation process[J].International Journal of Pharmaceutics,2013,449:72-83.

    [21]Dragoi E,Curteanu S,F(xiàn)issore D.Freeze-drying modeling and monitoring using a new neuro-evolutive technique[J].Chemical Engineering Science,2012,72:195-204.

    [22]Beer D,Wiggenhorn M,Hawe A,et al.Optimization of a pharmaceutical freeze-dried product and its process using an experimental design approach and innovative process analyzers[J].Talanta,2011,83:1623-1633.

    [23]Rosas J G,Waard H,Beer T D,et al.NIR spectroscopy for the in-line monitoring of a multicomponent formulation during the entire freeze-drying process[J].Journal of Pharmaceutical and Biomedical Analysis,2014,97:39-46.

    [24]詹博,周?chē)?guó)燕,曹斌宏,等.食品藥品冷凍干燥數(shù)值模擬的現(xiàn)狀與展望[J].食品工業(yè)科技,2012,33(7):414-417.

    [25]Velardi S A,Barresi AA.Development of simplified models for the freeze-drying process and investigation of the optimal operating conditions[J].Chemical Engineering Research and Design,2008,86:9-22.

    綜述與專(zhuān)論

    Research progress in monitoring and control technology for pharmaceutical freeze-drying process

    LI Junqi,LI Baoguo

    (School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

    Abstract:For the optimization of pharmaceutical freeze-drying process,the key is to shorten the drying time as much as possible and keep the good quality of the product.Therefore,it is very important to monitor and control the freeze-drying process precisely,i.e.,to keep the temperature of product within a reasonable range,judge the end point of drying time accurately,and a good control of chamber pressure and shelf temperature is needed to optimize the process.Research progress in pharmaceutical freeze-drying process monitoring and control technology in recent years is summarized in this paper,including monitoring system based on dynamic parameters estimation(DPE),monitoring system based on Kalman filter method,dew point method to judge the end point of primary drying,and model predictive control(MPC).The following points should be taken into account for further research on pharmaceutical freeze-drying monitoring and control technology.The contribution of three different heat transfer modes(radiation,convection,and conduction)should be investigated,so that a two-dimensional or three-dimensional model of pharmaceutical lyophilization in vials can be built,and more precise monitoring of process parameters can be obtained;Real-time optimal control strategy of shelf temperature and chamber pressure should be studied,with the aim of controllingbook=8,ebook=223pharmaceutical freeze-drying process timely and effectively.

    Key words:freeze-drying;primary drying;modeling;monitoring;control

    基金項(xiàng)目:上海市聯(lián)盟計(jì)劃(LM201364)及上海市教委科研創(chuàng)新重點(diǎn)資助項(xiàng)目(14ZZ133)。

    收稿日期:2014-12-30;修改稿日期:2015-04-07。

    DOI:10.16085/j.issn.1000-6613.2015.08.036

    文章編號(hào):1000–6613(2015)08–3128–05

    文獻(xiàn)標(biāo)志碼:A

    中圖分類(lèi)號(hào):TQ 025.3

    猜你喜歡
    模擬冷凍干燥控制
    疫苗的泡沫冷凍干燥工藝分析
    冷凍干燥法制備稻殼灰基二氧化硅氣凝膠及其改性研究
    能源工程(2021年6期)2022-01-06 02:04:30
    地黃真空冷凍干燥工藝的優(yōu)化
    中成藥(2019年12期)2020-01-04 02:02:28
    真空冷凍干燥技術(shù)在生物制藥方面的應(yīng)用
    化工管理(2017年6期)2017-03-03 14:12:58
    一個(gè)高分子模擬計(jì)算網(wǎng)格的作業(yè)管理
    工業(yè)機(jī)器人模擬仿真技術(shù)在職業(yè)教育中的應(yīng)用淺析
    淺析柔道運(yùn)動(dòng)員的模擬實(shí)戰(zhàn)訓(xùn)練
    虛擬機(jī)局域網(wǎng)組建技術(shù)應(yīng)用初探
    社會(huì)轉(zhuǎn)型期行政權(quán)控制的路徑探索
    科技視界(2016年21期)2016-10-17 17:06:18
    會(huì)計(jì)預(yù)算控制現(xiàn)狀及方法
    我要看黄色一级片免费的| 色精品久久人妻99蜜桃| bbb黄色大片| √禁漫天堂资源中文www| 女人爽到高潮嗷嗷叫在线视频| 日本欧美视频一区| 国产在线一区二区三区精| 亚洲,欧美,日韩| 天天躁夜夜躁狠狠久久av| 99热网站在线观看| 欧美老熟妇乱子伦牲交| 国产福利在线免费观看视频| 老司机午夜十八禁免费视频| av国产精品久久久久影院| 丝袜脚勾引网站| 少妇精品久久久久久久| 天天添夜夜摸| av片东京热男人的天堂| 高清不卡的av网站| 97在线人人人人妻| 久久中文字幕一级| 亚洲精品乱久久久久久| 精品高清国产在线一区| 午夜精品国产一区二区电影| 性高湖久久久久久久久免费观看| 一区二区av电影网| netflix在线观看网站| 亚洲国产欧美日韩在线播放| 国产精品九九99| 青草久久国产| 高潮久久久久久久久久久不卡| 欧美精品亚洲一区二区| 男人添女人高潮全过程视频| 国产熟女欧美一区二区| 又粗又硬又长又爽又黄的视频| 欧美乱码精品一区二区三区| 亚洲av在线观看美女高潮| 国产一区二区三区综合在线观看| 成年人黄色毛片网站| √禁漫天堂资源中文www| 精品少妇内射三级| 国产黄色视频一区二区在线观看| 成年动漫av网址| 宅男免费午夜| 国产成人影院久久av| 高清av免费在线| 91精品国产国语对白视频| 啦啦啦中文免费视频观看日本| 欧美黄色片欧美黄色片| 国产精品麻豆人妻色哟哟久久| 成人午夜精彩视频在线观看| 赤兔流量卡办理| 国产亚洲av片在线观看秒播厂| 黄色视频不卡| 国产熟女欧美一区二区| 中文欧美无线码| 久久久国产一区二区| 午夜福利视频精品| 侵犯人妻中文字幕一二三四区| 亚洲专区中文字幕在线| 精品一品国产午夜福利视频| 91精品三级在线观看| 纵有疾风起免费观看全集完整版| 日本av免费视频播放| 欧美日韩亚洲国产一区二区在线观看 | 亚洲视频免费观看视频| 啦啦啦 在线观看视频| 亚洲七黄色美女视频| 麻豆国产av国片精品| 精品久久蜜臀av无| 亚洲情色 制服丝袜| 人妻人人澡人人爽人人| 日韩中文字幕欧美一区二区 | 十分钟在线观看高清视频www| 国产人伦9x9x在线观看| 天天躁日日躁夜夜躁夜夜| 国产精品免费视频内射| 丝袜喷水一区| 岛国毛片在线播放| 交换朋友夫妻互换小说| 国产欧美日韩综合在线一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 少妇人妻 视频| 国产成人免费观看mmmm| 男女国产视频网站| 午夜91福利影院| 黄网站色视频无遮挡免费观看| 午夜福利视频在线观看免费| 亚洲视频免费观看视频| 99国产综合亚洲精品| 超色免费av| 国产一区亚洲一区在线观看| 色视频在线一区二区三区| 国产老妇伦熟女老妇高清| 悠悠久久av| 国产精品香港三级国产av潘金莲 | 黄片播放在线免费| 中文乱码字字幕精品一区二区三区| 午夜视频精品福利| 深夜精品福利| 亚洲精品日韩在线中文字幕| 国产97色在线日韩免费| 亚洲 欧美一区二区三区| 亚洲精品第二区| 人妻一区二区av| 午夜av观看不卡| 久久青草综合色| 国产欧美日韩一区二区三区在线| 国产成人精品久久久久久| 亚洲综合色网址| 少妇 在线观看| 亚洲,欧美精品.| 久久九九热精品免费| 中文字幕av电影在线播放| 国产在线一区二区三区精| 中文字幕人妻丝袜制服| 亚洲欧洲精品一区二区精品久久久| 精品国产超薄肉色丝袜足j| 九草在线视频观看| 又粗又硬又长又爽又黄的视频| 日韩 亚洲 欧美在线| 人妻一区二区av| 99热国产这里只有精品6| 亚洲,一卡二卡三卡| 亚洲av综合色区一区| 91麻豆精品激情在线观看国产 | 十分钟在线观看高清视频www| 黄色毛片三级朝国网站| avwww免费| 嫁个100分男人电影在线观看 | 亚洲欧美精品自产自拍| 久久影院123| 美女视频免费永久观看网站| 午夜福利免费观看在线| 国产成人啪精品午夜网站| 一级黄色大片毛片| 19禁男女啪啪无遮挡网站| 两人在一起打扑克的视频| 男人爽女人下面视频在线观看| 成年动漫av网址| 在线观看免费日韩欧美大片| 一本大道久久a久久精品| 青春草视频在线免费观看| 电影成人av| 国产精品一区二区在线不卡| av天堂在线播放| 肉色欧美久久久久久久蜜桃| 久久国产精品人妻蜜桃| av线在线观看网站| 国产精品久久久久久人妻精品电影 | 99国产综合亚洲精品| 精品久久蜜臀av无| 老汉色∧v一级毛片| 国产av精品麻豆| av网站在线播放免费| 国产精品99久久99久久久不卡| 99re6热这里在线精品视频| 国产精品 国内视频| 成年女人毛片免费观看观看9 | 两个人免费观看高清视频| 女人久久www免费人成看片| 日韩制服骚丝袜av| 亚洲欧美一区二区三区久久| 午夜免费成人在线视频| 1024视频免费在线观看| 91精品伊人久久大香线蕉| 一个人免费看片子| 汤姆久久久久久久影院中文字幕| 日日爽夜夜爽网站| 欧美变态另类bdsm刘玥| 七月丁香在线播放| 国产精品久久久久久精品古装| 国产男女超爽视频在线观看| 久久人人爽av亚洲精品天堂| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 久久久久久久精品精品| 午夜影院在线不卡| 亚洲精品久久久久久婷婷小说| 精品第一国产精品| 国产xxxxx性猛交| 国产黄色免费在线视频| 欧美人与性动交α欧美软件| 午夜免费男女啪啪视频观看| 日本wwww免费看| 老司机亚洲免费影院| 国产又色又爽无遮挡免| 七月丁香在线播放| 晚上一个人看的免费电影| 国产熟女欧美一区二区| 搡老乐熟女国产| 999久久久国产精品视频| 最近最新中文字幕大全免费视频 | 波野结衣二区三区在线| 两个人看的免费小视频| 美女中出高潮动态图| 久久精品亚洲av国产电影网| 在线 av 中文字幕| 午夜福利在线免费观看网站| 欧美日韩视频高清一区二区三区二| 人人妻,人人澡人人爽秒播 | 97人妻天天添夜夜摸| 99re6热这里在线精品视频| 在线av久久热| 久久久久久久国产电影| 亚洲一区二区三区欧美精品| 又黄又粗又硬又大视频| 女性生殖器流出的白浆| 久久ye,这里只有精品| 久久精品国产a三级三级三级| 天天躁夜夜躁狠狠久久av| 真人做人爱边吃奶动态| 亚洲第一青青草原| 美女扒开内裤让男人捅视频| av线在线观看网站| 日韩大码丰满熟妇| 男人爽女人下面视频在线观看| 亚洲精品美女久久av网站| 国产男人的电影天堂91| 每晚都被弄得嗷嗷叫到高潮| 国产黄色免费在线视频| 国产老妇伦熟女老妇高清| 日本黄色日本黄色录像| 午夜激情久久久久久久| 欧美亚洲日本最大视频资源| 久久天躁狠狠躁夜夜2o2o | 满18在线观看网站| 日日夜夜操网爽| 一区福利在线观看| 爱豆传媒免费全集在线观看| 国产欧美日韩一区二区三区在线| 色精品久久人妻99蜜桃| 久久精品国产综合久久久| 国产精品久久久人人做人人爽| 菩萨蛮人人尽说江南好唐韦庄| 一级黄片播放器| 免费高清在线观看日韩| 日日夜夜操网爽| 久久亚洲精品不卡| 中文字幕色久视频| av天堂在线播放| 欧美日韩成人在线一区二区| 老司机在亚洲福利影院| 丝袜美腿诱惑在线| 美女福利国产在线| 久久久久精品人妻al黑| 成人18禁高潮啪啪吃奶动态图| 午夜福利视频在线观看免费| 亚洲av成人精品一二三区| 国产亚洲av高清不卡| 中文字幕另类日韩欧美亚洲嫩草| 免费一级毛片在线播放高清视频 | 精品人妻1区二区| 又大又爽又粗| 精品福利观看| 肉色欧美久久久久久久蜜桃| 99久久99久久久精品蜜桃| 亚洲,欧美精品.| av又黄又爽大尺度在线免费看| 午夜视频精品福利| 成年人黄色毛片网站| 一边摸一边抽搐一进一出视频| 亚洲av日韩精品久久久久久密 | 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产av新网站| 亚洲专区中文字幕在线| 亚洲天堂av无毛| 久久久久网色| 少妇裸体淫交视频免费看高清 | 国产成人精品久久二区二区免费| 欧美亚洲 丝袜 人妻 在线| 超碰成人久久| 成人亚洲欧美一区二区av| 久久ye,这里只有精品| 少妇精品久久久久久久| 两人在一起打扑克的视频| 老汉色∧v一级毛片| 日韩一本色道免费dvd| 国语对白做爰xxxⅹ性视频网站| 久久久久精品人妻al黑| 久久久久久久久免费视频了| 欧美变态另类bdsm刘玥| 日本猛色少妇xxxxx猛交久久| 亚洲少妇的诱惑av| 日本wwww免费看| 亚洲av综合色区一区| 一级片免费观看大全| 亚洲国产最新在线播放| 亚洲成色77777| 亚洲av电影在线观看一区二区三区| 国产xxxxx性猛交| 精品一区二区三区四区五区乱码 | 久久免费观看电影| 欧美日韩精品网址| 永久免费av网站大全| 精品少妇内射三级| 久久精品国产亚洲av涩爱| 亚洲精品国产一区二区精华液| av在线app专区| 在线观看免费午夜福利视频| 国产高清视频在线播放一区 | 久久ye,这里只有精品| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 999久久久国产精品视频| 久久精品亚洲av国产电影网| 国产精品国产三级国产专区5o| cao死你这个sao货| 91九色精品人成在线观看| 国产精品九九99| 天堂8中文在线网| 成年动漫av网址| 亚洲国产av新网站| 久久精品亚洲熟妇少妇任你| 日韩熟女老妇一区二区性免费视频| 9色porny在线观看| 免费黄频网站在线观看国产| 国产黄色免费在线视频| 亚洲欧美日韩另类电影网站| 亚洲一区二区三区欧美精品| 侵犯人妻中文字幕一二三四区| 国产成人精品久久二区二区91| 亚洲国产欧美网| 国产精品九九99| 校园人妻丝袜中文字幕| 深夜精品福利| 50天的宝宝边吃奶边哭怎么回事| 亚洲,欧美精品.| 久久久精品94久久精品| 另类亚洲欧美激情| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 中文字幕最新亚洲高清| 午夜免费男女啪啪视频观看| 久久ye,这里只有精品| 午夜福利乱码中文字幕| 国产精品欧美亚洲77777| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品人妻久久久影院| 99久久人妻综合| 久久精品久久久久久噜噜老黄| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 精品少妇一区二区三区视频日本电影| 又大又爽又粗| 久久精品久久久久久久性| 国产黄色免费在线视频| 亚洲中文字幕日韩| 两性夫妻黄色片| 一区二区三区激情视频| 免费人妻精品一区二区三区视频| 精品国产超薄肉色丝袜足j| 日韩大码丰满熟妇| 激情视频va一区二区三区| 麻豆av在线久日| 亚洲欧美一区二区三区国产| 亚洲精品日韩在线中文字幕| 亚洲成人免费电影在线观看 | 国产深夜福利视频在线观看| 少妇人妻久久综合中文| 成人国产一区最新在线观看 | 国产真人三级小视频在线观看| 各种免费的搞黄视频| 久久久国产欧美日韩av| 91麻豆av在线| 日韩精品免费视频一区二区三区| 美女扒开内裤让男人捅视频| 国产午夜精品一二区理论片| 伦理电影免费视频| 深夜精品福利| 日本av手机在线免费观看| 高清视频免费观看一区二区| av国产精品久久久久影院| 国产一区二区三区综合在线观看| 久久国产亚洲av麻豆专区| 国产一区二区激情短视频 | 成人国产一区最新在线观看 | 精品一品国产午夜福利视频| 精品一区在线观看国产| 国产91精品成人一区二区三区 | 亚洲综合色网址| 精品国产一区二区三区四区第35| 大码成人一级视频| 久久亚洲精品不卡| 制服人妻中文乱码| 人妻一区二区av| 黄色a级毛片大全视频| 亚洲久久久国产精品| 国产日韩一区二区三区精品不卡| 精品人妻熟女毛片av久久网站| 一级毛片女人18水好多 | 超碰成人久久| 男的添女的下面高潮视频| 女性被躁到高潮视频| cao死你这个sao货| 秋霞在线观看毛片| 久久精品亚洲av国产电影网| 国产高清视频在线播放一区 | 中文字幕高清在线视频| 成在线人永久免费视频| 老司机亚洲免费影院| cao死你这个sao货| 午夜两性在线视频| 一区二区三区精品91| 妹子高潮喷水视频| 国产成人av教育| 精品人妻1区二区| 国产精品亚洲av一区麻豆| 天天影视国产精品| 亚洲国产中文字幕在线视频| a级片在线免费高清观看视频| 国产高清国产精品国产三级| 日本黄色日本黄色录像| 99九九在线精品视频| 免费观看人在逋| 国产老妇伦熟女老妇高清| 老鸭窝网址在线观看| 久久人人爽人人片av| 国产色视频综合| 悠悠久久av| 少妇人妻 视频| 男女午夜视频在线观看| av又黄又爽大尺度在线免费看| 精品视频人人做人人爽| 久久国产精品人妻蜜桃| 国产一区二区激情短视频 | 男女午夜视频在线观看| 日韩av免费高清视频| 午夜福利乱码中文字幕| 人人澡人人妻人| 一二三四社区在线视频社区8| 黄色毛片三级朝国网站| 一级毛片黄色毛片免费观看视频| 精品福利观看| 国产一区二区三区av在线| 老司机靠b影院| 日本wwww免费看| 丝袜喷水一区| 久久久精品94久久精品| 精品人妻熟女毛片av久久网站| 午夜福利视频精品| 天堂中文最新版在线下载| 精品福利观看| 后天国语完整版免费观看| 中文字幕制服av| 国产麻豆69| 搡老岳熟女国产| 激情五月婷婷亚洲| 啦啦啦 在线观看视频| netflix在线观看网站| 日本a在线网址| 午夜激情av网站| 国产精品偷伦视频观看了| 国产成人欧美| 成人国产av品久久久| 亚洲欧美日韩另类电影网站| 婷婷色综合大香蕉| 久久精品国产a三级三级三级| 久久性视频一级片| 亚洲国产欧美日韩在线播放| 久久久久久亚洲精品国产蜜桃av| 日韩电影二区| 国产一区二区 视频在线| 亚洲国产日韩一区二区| 最近手机中文字幕大全| 一级片免费观看大全| 妹子高潮喷水视频| 免费不卡黄色视频| 久久精品熟女亚洲av麻豆精品| 国产片内射在线| 亚洲国产日韩一区二区| 一本大道久久a久久精品| 日韩av不卡免费在线播放| 久久久精品区二区三区| 日本vs欧美在线观看视频| 国产淫语在线视频| 亚洲精品久久午夜乱码| 欧美日韩综合久久久久久| 9热在线视频观看99| 免费人妻精品一区二区三区视频| av线在线观看网站| 国产精品三级大全| 黄色视频不卡| www日本在线高清视频| 免费看十八禁软件| 欧美+亚洲+日韩+国产| 成人国语在线视频| 欧美老熟妇乱子伦牲交| 亚洲国产欧美网| 天天操日日干夜夜撸| 国产无遮挡羞羞视频在线观看| 老司机午夜十八禁免费视频| 80岁老熟妇乱子伦牲交| 黄色毛片三级朝国网站| 大型av网站在线播放| 热re99久久精品国产66热6| 黄色片一级片一级黄色片| 日韩视频在线欧美| 精品少妇久久久久久888优播| 欧美日韩亚洲综合一区二区三区_| 日韩 亚洲 欧美在线| 99香蕉大伊视频| 午夜影院在线不卡| 欧美日韩视频高清一区二区三区二| 国产高清不卡午夜福利| 欧美日韩av久久| 高清黄色对白视频在线免费看| www日本在线高清视频| 色精品久久人妻99蜜桃| 伊人久久大香线蕉亚洲五| 亚洲精品乱久久久久久| 欧美黄色片欧美黄色片| 一区二区三区乱码不卡18| 欧美精品人与动牲交sv欧美| 少妇 在线观看| 久久精品成人免费网站| 国产精品久久久久久精品古装| 操美女的视频在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲精品国产一区二区精华液| 精品一品国产午夜福利视频| 久久久久久久精品精品| 国产免费福利视频在线观看| 韩国高清视频一区二区三区| 成人国产av品久久久| 狠狠精品人妻久久久久久综合| 色综合欧美亚洲国产小说| 在线精品无人区一区二区三| 2021少妇久久久久久久久久久| 久久精品亚洲av国产电影网| 91成人精品电影| 国产片内射在线| 亚洲av电影在线进入| 久久久久国产精品人妻一区二区| 久久久精品国产亚洲av高清涩受| 一区福利在线观看| 99国产精品99久久久久| 久久久亚洲精品成人影院| 少妇粗大呻吟视频| 国产男女内射视频| 国产欧美亚洲国产| 丝袜脚勾引网站| 18禁黄网站禁片午夜丰满| 欧美在线黄色| 亚洲精品av麻豆狂野| 国产女主播在线喷水免费视频网站| 亚洲国产欧美日韩在线播放| 国产亚洲欧美精品永久| 黄色一级大片看看| 水蜜桃什么品种好| 亚洲国产av影院在线观看| 国产主播在线观看一区二区 | 亚洲欧美一区二区三区黑人| 亚洲精品在线美女| 亚洲国产欧美在线一区| 成年人免费黄色播放视频| 免费日韩欧美在线观看| 99精品久久久久人妻精品| 日韩精品免费视频一区二区三区| 国产1区2区3区精品| 欧美在线一区亚洲| 自线自在国产av| √禁漫天堂资源中文www| 成人亚洲精品一区在线观看| 亚洲精品第二区| av片东京热男人的天堂| 亚洲男人天堂网一区| 日韩免费高清中文字幕av| 日本av免费视频播放| 国产免费又黄又爽又色| 黄色毛片三级朝国网站| 男女边摸边吃奶| 新久久久久国产一级毛片| 黑人欧美特级aaaaaa片| 亚洲七黄色美女视频| 纯流量卡能插随身wifi吗| 黑人猛操日本美女一级片| 夫妻性生交免费视频一级片| 亚洲三区欧美一区| xxx大片免费视频| 黄色一级大片看看| 97精品久久久久久久久久精品| 操美女的视频在线观看| 亚洲,欧美,日韩| 男人爽女人下面视频在线观看| 亚洲欧洲国产日韩| a级片在线免费高清观看视频| 老司机深夜福利视频在线观看 | 中文字幕av电影在线播放| 色精品久久人妻99蜜桃| 精品亚洲成a人片在线观看| 人人妻人人添人人爽欧美一区卜| cao死你这个sao货| 亚洲精品一区蜜桃| 国产成人av激情在线播放| cao死你这个sao货| 精品人妻1区二区| 黄色a级毛片大全视频| 精品福利永久在线观看| av天堂在线播放| 黄色a级毛片大全视频| 精品福利永久在线观看| 精品亚洲乱码少妇综合久久| 91老司机精品| 久久99一区二区三区| 国产深夜福利视频在线观看| 黑人欧美特级aaaaaa片| cao死你这个sao货| 精品亚洲乱码少妇综合久久| 欧美日韩视频高清一区二区三区二| 精品福利永久在线观看| av在线app专区| 脱女人内裤的视频| 麻豆av在线久日| 欧美精品高潮呻吟av久久|