[發明專利]中藥綜合品質評價方法

有效

• 申請號: 202010510016.2
• 申請日: 2020-06-08
• 公開（公告）號: CN111398538B
• 公開（公告）日: 2020-10-09
• 發明（設計）人: 耿炤;周美娟;鈕犇
• 申請（專利權）人: 江西匯仁藥業股份有限公司
• 主分類號: G01N33/15
• 分類號: G01N33/15
• 代理公司: 北京華科聯合專利事務所(普通合伙) 11130
• 代理人: 王為

1.一種中藥綜合品質的評價方法，其特征在于，包括以下步驟：

步驟1、收集該中藥生產過程歷史及待評價批次的物料、過程、檢驗指標的相應實際數值，包括但不限于：

步驟1.1 同一種中藥的原料、飲片、中間產品、制劑成品及其留樣的標準檢驗結果，其指標名稱分別為Aa，Ab，Ac …… 其值分別為Aa001，Ab001，Ac001 …… Aa002，Ab002，Ac002 …… 其中的數字代表生產批號；

步驟1.2 同一種中藥的生產過程中提取、分離、濃縮、干燥、制劑及包裝的工藝過程參數，其指標名稱分別為Ba，Bb，Bc ……其值分別為Ba001，Bb001，Bc001 …… Ba002，Bb002，Bc002 …… 其中的數字代表生產批號；

步驟1.3 同一種中藥產品非臨床、臨床研究及上市后研究數據，其指標名稱分別為Ca，Cb，Cc ……其值分別為Ca001，Cb001，Cc001 …… Ca002，Cb002，Cc002 …… 其中的數字代表生產批號；

步驟2、數據處理：

步驟2.1 計算除待評價批次x外，各歷史批次步驟1指標數據的歷史均值和標準偏差，即：

均值Aa=∑（Aa）/N，

標準偏差sAa=sqrt(∑（(Aa-均值Aa)² )/（N-1）)

均值Ab=∑（Ab）/N，

標準偏差sAb=sqrt(∑（(Ab-均值Ab)² )/（N-1）)

……

步驟2.2 計算待評價批次x各指標數據與其歷史均值的差，除以標準偏差，以倍數n計，即：

n_Aa=（Aa00x-均值Aa）/標準偏差sAa

n_Ab=（Ab00x-均值Ab）/標準偏差sAb

……

步驟2.3 將以上所有指標在100%內進行權重分解，作為初始權重E0，即：

E0Aa+E0Ab+……E0Ba+E0Bb+……E0Ca+E0Cb+……=1

步驟3、權重調整：

根據步驟2.2所述n對相應指標的權重進行乘法計算，得到該指標的動態權重E1，即：

E1Aa=E0Aa×n_Aa

E1Ab=E0Ab×n_Ab

……

步驟4、計算及結果判斷：

步驟4.1 將以上各批次指標數據按其初始權重累計求和，以批次、計算結果分別為橫、縱坐標作圖，即可得到該中藥的歷史綜合品質Q趨勢，其逐點坐標為“批號N，綜合品質Q”， N=001時，Q001=Aa001×E0Aa+Ab001×E0Ab+Ac001×E0Ac+……Ba001×E0Ba+Bb001×E0Bb+Bc001×E0Bc+……

步驟4.2 計算步驟4.1各批次值的均值和標準偏差，即得到該中藥總體品質狀態，即:該中藥的綜合品質是以Q均值為中心，具有相應標準偏差分布的一種狀態；

步驟4.3 將待評價批次產品各指標數據按照步驟3權重調整后進行累計同步驟4.1方法求和，所得數值在步驟4.1圖上標記，即可得到本批中藥品質與歷史品質趨勢的偏離狀態，即與趨勢線的距離大小；

步驟4.4 計算各批次動態權重之和，當與數值1偏差較大，說明本批次各指標數值發生了顯著變化；

步驟4.5 計算步驟4.3數值與步驟4.2平均值的算術偏差，再除以步驟4.2標準偏差，得到批次x的品質Qx與其Q歷史分布，即Q（x-1）的關系數值，根據正態分布概率判斷：

當該數值＞3時，說明該中藥批次x與歷史批次的綜合品質比較，具有99%以上的可能性發生了變異，為異常；

當該數值＞2，且≤3時，說明該中藥批次x與歷史批次的綜合品質比較，具有95%以上的可能性發生了變異，為異常；

當該數值＞1，且≤2時，說明該中藥批次x與歷史批次的綜合品質比較，具有68%以上的可能性發生了變異，為異常；

當該數值≤1時，說明該中藥批次x與歷史批次的綜合品質比較，具有68%以下的可能性發生了變異，在此條件下可認為正常；

步驟5、進一步處理：

根據步驟4.4結果，當變異的可能性大于68%時，根據幅度依次列出批次x各指標的動態權重與初始權重狀態差異，即得影響批次x品質Q的指標及其影響程度，確認是否有異常原輔材料、生產工藝、設備、環境要素影響，確認有則本批各指標數據另行保存，不再進行步驟6；變異的可能性低于68%時，直接進入步驟6；

步驟6、將未發現有異常原輔材料、生產工藝、設備、環境要素的批次x各指標，返回步驟1納入歷史數值；

步驟7、按照步驟1-步驟6計算批次x+1。