CN104977271B - 一种红花醇沉过程有效成分近红外在线检测方法 - Google Patents

一种红花醇沉过程有效成分近红外在线检测方法 Download PDF

Info

Publication number
CN104977271B
CN104977271B CN201410135941.6A CN201410135941A CN104977271B CN 104977271 B CN104977271 B CN 104977271B CN 201410135941 A CN201410135941 A CN 201410135941A CN 104977271 B CN104977271 B CN 104977271B
Authority
CN
China
Prior art keywords
near infrared
precipitation
precipitation solution
sample
alcohol precipitation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410135941.6A
Other languages
English (en)
Other versions
CN104977271A (zh
Inventor
姚小青
孙长海
张坤
高俊敏
张桂萍
胡小锋
孙明珍
黎先军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin Chase Sun Pharmaceutical Co Ltd
Original Assignee
Tianjin Chase Sun Pharmaceutical Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin Chase Sun Pharmaceutical Co Ltd filed Critical Tianjin Chase Sun Pharmaceutical Co Ltd
Priority to CN201410135941.6A priority Critical patent/CN104977271B/zh
Publication of CN104977271A publication Critical patent/CN104977271A/zh
Application granted granted Critical
Publication of CN104977271B publication Critical patent/CN104977271B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明提供一种红花醇沉过程有效成分近红外在线检测方法,包括①近红外在线检测循环管路;②在线采集红花醇沉液的近红外透射光谱及醇沉液样本;③采用高效液相色谱法和烘干称重法测得醇沉液样本中各质控指标(包括羟基红花黄色素A和可溶性固形物含量)信息;④采用偏最小二乘算法建立各质控指标模型;⑤将所建模型用于在线分析醇沉过程中各质控指标的变化趋势。本发明将近红外光谱技术用于红花醇沉过程,采用流通池结合光纤探头的方式,实现了红花醇沉过程中羟基红花黄色素A和可溶性固形物含量的在线定量分析,为红花醇沉过程的在线控制提供依据及有效的指导。

Description

一种红花醇沉过程有效成分近红外在线检测方法
技术领域
本发明属于近红外在线检测领域,具体涉及一种红花醇沉过程有效成分近红外在线检测方法。
背景技术
红花具有活血通络的功效,有效成分主要集中在水溶性红花黄色素,如羟基红花黄色素A等。醇沉工艺是红花药材中有效成分纯化以及去除杂质的常用工艺手段,直接关系到红花成品的功效和质量稳定性。目前,醇沉工艺的质量控制主要依靠经验和传统质量分析方法(HPLC等),耗时费力,缺乏有效的指标性成分含量的实时监测手段,易造成不同批次醇沉液质量的不稳定,导致药品批次间的质量差异,以及原药材、能源、时间等的浪费。故研究发展红花醇沉过程中关键质控指标的在线检测方法,有助于解决红花醇沉过程中关键控制指标的质量控制问题,对于中药工业技术进步和产品质量升级具有重大现实意义。
近红外(NIR)光谱技术作为一种快速无损的绿色分析技术,具有快速分析、样品处理简单、无需消耗试剂等特点。近年来,近红外光谱技术已经越来越多的被应用于中药研究,包括药材产地鉴别、有效组分含量测定和制药过程的在线检测和监控。从近年来研究进展情况看来,近红外光谱分析技术是最有希望在中药生产过程实现在线检测及质量控制的过程分析技术之一。目前红花醇沉工艺过程一般不检测有效成分,或者手动离线检测有效成分,其检测过程中,容易污染药液,影响药液质量,此外,由于离线数据的提供往往都滞后于生产过程,无法及时反馈生产现状,容易造成生产成本的提高。
发明内容
本发明的目的在于提供一种红花醇沉过程的在线检测方法。该方法的检测目标为实现红花醇沉过程中各质控指标的在线定量分析,为红花醇沉过程质量控制提供方法。
该方法包括以下步骤:
步骤1,安装近红外在线检测循环管路;
步骤2,在线采集红花醇沉液的近红外透射光谱及醇沉液样本;
步骤3,采用传统分析方法(高效液相色谱法和烘干称重法)测得醇沉液样本中各质控指标信息;
步骤4,剔除异常光谱;
步骤5,选择近红外光谱建模波段和预处理方法;
步骤6,用偏最小二乘算法建立各质控指标模型,并采用各模型评价指标考察模型性能;
步骤7,已建模型用于在线分析红花醇沉过程中各质控指标的变化趋势。
本发明具体方案如下:
(1)安装近红外在线检测循环管路:
近红外在线检测循环管路依次连接以下部件:变频离心泵、双联过滤器、手动阀1、流通池(含近红外光纤探头)、取样口、手动阀2。变频离心泵将醇沉罐内的醇沉药液输送至双联过滤器,过滤后经过手动阀1流过流通池,近红外光纤探头对流通池内的醇沉药液进行光谱采集,通过取样口取样流通池内的醇沉药液,流通池通过手动阀2与醇沉罐连通。
使用时,打开变频离心泵,醇沉液经由双联过滤器滤除固体颗粒后到达流通池,连接在流通池左右两边的光纤探头适时在线采集流通池中醇沉液的近红外光谱,最后醇沉液分为两路,一路回醇沉罐,另一路通过取样口,用于收集醇沉液样品。流通池及循环管路中醇沉液流速控制在2L/min以内。双联过滤器用于滤除醇沉液中的大部分固体杂质颗粒,过滤精度为80微米。
(2)在线采集红花醇沉液的近红外透射光谱及醇沉液样本:
采用透射法采集近红外光谱,光谱范围为4000cm-1~12000cm-1,扫描次数为32次,分辨率为8cm-1,以空气为参比。每隔1分钟在线采集通过流通池的醇沉液近红外光谱;每隔5分钟从取样口采集醇沉液样本。采集醇沉液样本的同时采集近红外光谱。收集不同批次醇沉过程中的醇沉液样品,随机选择其中1/4~1/3批次的数据作为验证集,其余样品作为校正集参与建模。
(3)采用传统分析方法(高效液相色谱法和烘干称重法)测得醇沉液样本中各质控指标信息:
所述的醇沉液样本的各质控指标包括羟基红花黄色素A和可溶性固形物含量。采用高效液相色谱法(HPLC)测定醇沉液样本中羟基红花黄色素A浓度;使用烘干称重法测定可溶性固形物含量。
a、HPLC色谱条件:Agilent eclipse C18分析柱(250×4.6mm,5μm);甲醇-乙腈-0.7%磷酸溶液(v/v,26∶2∶72)为流动相;流速1mL/min;检测波长403nm;柱温40℃;进样量5μL。
醇沉液样本于1500r/min高速离心机中离心10分钟,滤过(0.45μm微孔滤膜),取续滤液用于液相分析。
b.烘干称重法:醇沉液样本于1500r/min高速离心机中离心10分钟,取上清液用于可溶性固形物含量分析。烘干至恒重的扁形瓶(两次烘干后重量差小于5mg)称重X0,取样品约10mL至扁形瓶,称重X1,水浴蒸干,105℃烘5h,取出置干燥器内冷却30min,迅速称重X2。可溶性固形物含量按下式计算:
可溶性固形物=(X2-X0)/(X1-X0)×100%
(4)剔除异常光谱
醇沉过程中产生的气泡以及流通池中的固体颗粒都会影响近红外光谱的采集,导致异常光谱的产生。本发明计算光谱的马氏距离,并使用肖维勒(Chauvenet)准则剔除异常光谱。若测量值Xi(l≤i≤n)的残差满足|Vi|>Wnσ则Xi被视为异常数据,予以剔除。其中,Vi为残差,σ为标准差,Wn可查表得到。
(5)选择近红外光谱建模波段和预处理方法;
采用一阶导数法(Savitzky-Golay平滑)和标准正交变换算法预处理近红外光谱数据,分别用于消除基线漂移、噪音及固体颗粒等对光谱的影响。在选择建模波段时需排除以下波段:4500~5450cm-1波段,即“水峰”;7500~12000cm-1波段,存在较大的噪声,且没有显著的特征吸收。然后,通过光谱和质控指标的相关系数确定建模波段。因此,对于可溶性固形物含量模型使用5450~7500cm-1波段,对于羟基红花黄色素A模型则使用5450~6100cm-1波段。
(6)使用偏最小二乘算法建立各质控指标模型,并采用各模型评价指标考察模型性能;
模型评价指标包括:模型评价指标包括相关系数(R)、校正集和验证集预测误差均方根(RMSEC、RMSEP)、校正集和验证集相对偏差(RSEC和RSEP)。当R值接近于1,RMSEC和RMSEP值相互接近且RMSEP小于2倍的RMSEC,RSEC和RSEP相互接近且小于20%时说明所见模型具有较好的稳定性和预测精度,可以用于红花醇沉过程的在线检测。
(7)将已建模型用于在线分析红花醇沉过程中各质控指标的变化趋势。
在线采集红花醇沉液的近红外光谱图,将光谱数据输入到校正模型中,经过计算即可实时得知醇沉液中各质控指标的信息。
本发明将近红外在线分析技术引入到红花醇沉过程,实现对各质控指标(羟基红花黄色素A和可溶性固形物含量)的实时监测,有利于提高红花醇沉过程的质量控制水平,充分保证产品质量稳定、可靠。
附图说明
图1是醇沉过程近红外在线检测系统简图。
图2是偏最小二乘模型在线分析红花醇沉过程中羟基红花黄色素A浓度预测值与实际测定值的趋势对照图。
图3是偏最小二乘模型在线分析红花醇沉过程中可溶性固形物含量预测值与实际测定值的趋势对照图。
具体实施方式
下面结合附图和实施例做进一步的说明。
实施例1:
1.安装近红外在线检测循环管路
近红外在线检测循环管路依次连接以下部件:变频离心泵、双联过滤器、手动阀1、流通池(含近红外光纤探头)、取样口、手动阀2。变频离心泵将醇沉罐内的醇沉药液输送至双联过滤器,过滤后经过手动阀1流过流通池,近红外光纤探头对流通池内的醇沉药液进行光谱采集,通过取样口取样流通池内的醇沉药液,流通池通过手动阀2与醇沉罐连通。
使用时,打开变频离心泵,醇沉液经由双联过滤器滤除固体颗粒后到达流通池,连接在流通池左右两边的光纤探头适时在线采集流通池中醇沉液的近红外光谱,最后醇沉液分为两路,一路回醇沉罐,另一路通过取样口,用于收集醇沉液样品。流通池及循环管路中醇沉液流速控制在2L/min以内。双联过滤器用于滤除醇沉液中的大部分固体杂质颗粒,过滤精度为80微米。
2.近红外光谱和醇沉液样本的在线采集
醇沉罐内开始加醇时作为醇沉过程计时零点。醇沉过程中每隔1分钟在线采集醇沉液近红外光谱图,每隔5分钟从取样口采集醇沉液样本。近红外在线检测系统参见图1。所取样本分别用于羟基红花黄色素A和可溶性固形物含量的测定。重复15批红花的醇沉实验,每批次的实验都以相同方式进行取样和采集光谱。计算所有近红外光谱的马氏距离,并使用肖维勒(Chauvenet)准则剔除异常光谱。若测量值Xi(l≤i≤n)的残差满足|Vi|>Wnσ则Xi被视为异常数据,予以剔除。其中,Vi为残差,σ为标准差,Wn可查表得到。根据计算结果,有11个样本的光谱异常,将这些样本剔除后剩余数据用于建立定量模型。
3.各指标成分浓度和可溶性固形物含量的测定
采用高效液相色谱法(HPLC)测定醇沉液样本中羟基红花黄色素A浓度。使用烘干称重法测定可溶性固形物含量。
4.选择近红外光谱建模波段和预处理方法
采用一阶导数法(Savitzky-Golay平滑)和标准正交变换算法预处理近红外光谱数据,分别用于消除基线漂移、噪音及固体颗粒等对光谱的影响。
在选择建模波段时需排除以下波段:4500~5445cm-1波段,即“水峰”;7500~12000cm-1波段,存在较大的噪声,且没有显著的特征吸收。为了保证光谱选择区域的准确性,本发明对光谱和质控指标的相关系数进行了考察。选择相关系数大于0.8的光谱区域,根据计算结果,可溶性固形物含量模型使用5445~7500cm-1波段,羟基红花黄色素A模型则使用5445~6100cm-1波段。
5.采用偏最小二乘回归算法建立各质控指标模型
在适当波段范围内,采用一阶导数法(Savitzky-Golay平滑)和标准正交变换算法预处理方法,剔除异常光谱后将剩余数据用于建立偏最小二乘回归模型。可溶性固形物含量和羟基红花黄色素A的校正和验证结果见表1。可以看出,偏最小二乘回归模型的校正集和验证集相关系数均大于0.96,RMSEC和RMSEP值较小而且互相接近,RSEP值都能控制在10%以内。各质控指标偏最小二乘回归模型的校正和验证结果相近,泛化能力强,具有较好的模型稳定性和预测能力。
表1各质控指标偏最小二乘回归模型的校正和验证结果比较
6.在线分析醇沉过程中各质控指标变化趋势
将所建偏最小二乘回归模型用于在线分析5批红花醇沉过程中可溶性固形物含量和羟基红花黄色素A浓度,所得数据归为未知样品集,预测结果如表2所示。通过对比校正集、验证集和未知样品集的各模型评价指标值可以看出,偏最小二乘回归模型具有较高的预测准确度。各质控指标模型的RMSEC值和RMSEP值均与校正集和验证集结果相似,并且RMSEP值小于2倍RMSEC,RSEP值也都控制在10.3%以内。
采用偏最小二乘回归模型在线分析某一批次红花醇沉过程,各质控指标的预测趋势与实际测定值参见图2~3。从图中可以看出,各质控指标的预测趋势与实际测定值的变化趋势基本一致,能够满足中药生产过程实时分析的精度要求。
表2各质控指标偏最小二乘回归模型的在线分析结果(未知样品集)

Claims (1)

1.一种红花醇沉过程有效成分近红外在线检测方法,其特征在于,该方法包括以下步骤:
(1)安装近红外在线检测循环管路
近红外在线检测循环管路依次连接以下部件:变频离心泵、双联过滤器、手动阀1、流通池,含近红外光纤探头、取样口、手动阀2,变频离心泵将醇沉罐内的醇沉药液输送至双联过滤器,过滤后经过手动阀1流过流通池,近红外光纤探头对流通池内的醇沉药液进行光谱采集,通过取样口取样流通池内的醇沉药液,流通池通过手动阀2与醇沉罐连通;
使用时,打开变频离心泵,醇沉液经由双联过滤器滤除固体颗粒后到达流通池,连接在流通池左右两边的光纤探头适时在线采集流通池中醇沉液的近红外光谱,最后醇沉液分为两路,一路回醇沉罐,另一路通过取样口,用于收集醇沉液样品,流通池及循环管路中醇沉液流速控制在2L/min以内,双联过滤器用于滤除醇沉液中的大部分固体杂质颗粒,过滤精度为80微米;
(2)在线采集红花醇沉液的近红外透射光谱及醇沉液样本
采用透射法采集近红外光谱,光谱范围为4000cm-1~12000cm-1,扫描次数为32次,分辨率为8cm-1,以空气为参比,每隔1分钟在线采集通过流通池的醇沉液近红外光谱;每隔5分钟从取样口采集醇沉液样本;采集醇沉液样本的同时采集近红外光谱,收集不同批次醇沉过程中的醇沉液样品,随机选择其中1/4~1/3批次的数据作为验证集,其余样品作为校正集参与建模,
(3)采用传统分析方法,高效液相色谱法和烘干称重法,测得醇沉液样本中各质控指标信息
所述的醇沉液样本的各质控指标包括羟基红花黄色素A和可溶性固形物含量,采用高效液相色谱法测定醇沉液样本中羟基红花黄色素A浓度;使用烘干称重法测定可溶性固形物含量,
a、HPLC色谱条件:AgilenteclipseC18分析柱,250×4.6mm,5μm;甲醇-乙腈-0.7%磷酸溶液,v/v,26∶2∶72为流动相;流速1mL/min;检测波长403nm;柱温40℃;进样量5μL,
醇沉液样本于1500r/min高速离心机中离心10分钟,滤过,0.45μm微孔滤膜,取续滤液用于液相分析,
b.烘干称重法:醇沉液样本于1500r/min高速离心机中离心10分钟,取上清液用于可溶性固形物含量分析,烘干至恒重的扁形瓶,两次烘干后重量差小于5mg,称重X0,取样品约10mL至扁形瓶,称重X1,水浴蒸干,105℃烘5h,取出置干燥器内冷却30min,迅速称重X2,可溶性固形物含量按下式计算:
可溶性固形物=(X2-X0)/(X1-X0)×100%
(4)剔除异常光谱
醇沉过程中产生的气泡以及流通池中的固体颗粒都会影响近红外光谱的采集,导致异常光谱的产生,本发明计算光谱的马氏距离,并使用肖维勒准则剔除异常光谱,若测量值Xi(1≤i≤n)的残差满足|Vi|>Wnσ则Xi被视为异常数据,予以剔除,其中,Vi为残差,σ为标准差,Wn可查表得到,
(5)选择近红外光谱建模波段和预处理方法:
采用一阶导数法,Savitzky-Golay平滑和标准正交变换算法预处理近红外光谱数据,分别用于消除基线漂移、噪音及固体颗粒对光谱的影响,在选择建模波段时需排除以下波段:4500~5450cm-1波段,即“水峰”;7500~12000cm-1波段,存在较大的噪声,且没有显著的特征吸收,然后,通过光谱和质控指标的相关系数确定建模波段,因此,对于可溶性固形物含量模型使用5450~7500cm-1波段,对于羟基红花黄色素A模型则使用5450~6100cm-1波段,
(6)使用偏最小二乘算法建立各质控指标模型,并采用各模型评价指标考察模型性能
模型评价指标包括:模型评价指标包括相关系数,R、校正集和验证集预测误差均方根,RMSEC、RMSEP、校正集和验证集相对偏差,RSEC和RSEP,当R值接近于1,RMSEC和RMSEP值相互接近且RMSEP小于2倍的RMSEC,RSEC和RSEP相互接近且小于20%时说明所见模型具有较好的稳定性和预测精度,可以用于红花醇沉过程的在线检测,
(7)将已建模型用于在线分析红花醇沉过程中各质控指标的变化趋势在线采集红花醇沉液的近红外光谱图,将光谱数据输入到校正模型中,经过计算即可实时得知醇沉液中各质控指标的信息。
CN201410135941.6A 2014-04-08 2014-04-08 一种红花醇沉过程有效成分近红外在线检测方法 Active CN104977271B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410135941.6A CN104977271B (zh) 2014-04-08 2014-04-08 一种红花醇沉过程有效成分近红外在线检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410135941.6A CN104977271B (zh) 2014-04-08 2014-04-08 一种红花醇沉过程有效成分近红外在线检测方法

Publications (2)

Publication Number Publication Date
CN104977271A CN104977271A (zh) 2015-10-14
CN104977271B true CN104977271B (zh) 2018-07-03

Family

ID=54273988

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410135941.6A Active CN104977271B (zh) 2014-04-08 2014-04-08 一种红花醇沉过程有效成分近红外在线检测方法

Country Status (1)

Country Link
CN (1) CN104977271B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106226264A (zh) * 2016-05-05 2016-12-14 江苏康缘药业股份有限公司 青蒿金银花醇沉过程在线实时放行标准建立方法与放行方法及应用
CN106053384A (zh) * 2016-07-15 2016-10-26 江苏康缘药业股份有限公司 一种青蒿金银花醇沉浓缩过程快速定量检测方法
CN107356552A (zh) * 2017-06-12 2017-11-17 浙江大学 一种基于近红外光谱技术的黄芪一次醇沉工艺的过程监控方法
CN111077107A (zh) * 2020-01-08 2020-04-28 山东金璋隆祥智能科技有限责任公司 一种甜菊糖提取液中糖甙含量的在线检测方法
CN112414962B (zh) * 2020-12-14 2023-10-13 华润三九(雅安)药业有限公司 一种测定羟基红花黄色素a含量的方法
CN113588590B (zh) * 2021-08-11 2024-04-16 苏州泽达兴邦医药科技有限公司 一种基于数据挖掘的中药提取过程质量控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102252992A (zh) * 2011-04-28 2011-11-23 天津红日药业股份有限公司 一种对中药提取过程进行近红外在线检测方法
CN102313714A (zh) * 2011-09-21 2012-01-11 浙江大学 一种红花提取液测定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102621092B (zh) * 2012-03-17 2014-07-02 浙江大学 一种丹红注射液醇沉过程在线检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102252992A (zh) * 2011-04-28 2011-11-23 天津红日药业股份有限公司 一种对中药提取过程进行近红外在线检测方法
CN102313714A (zh) * 2011-09-21 2012-01-11 浙江大学 一种红花提取液测定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
红花醇沉液浓缩除醇过程中多元质控指标的近红外快速检测;陈雪英等;《药物分析杂志》;20101231;第2086-2092页(参见"3.2异常点剔除") *
近红外光谱法快速测定红花中羟基红花黄色素A;李耿等;《中成药》;20130731;第35卷(第7期);第1579-1582页 *

Also Published As

Publication number Publication date
CN104977271A (zh) 2015-10-14

Similar Documents

Publication Publication Date Title
CN104977271B (zh) 一种红花醇沉过程有效成分近红外在线检测方法
CN103913433B (zh) 一种丹红注射液双效浓缩过程在线检测方法
CN102621092B (zh) 一种丹红注射液醇沉过程在线检测方法
Jin et al. Near infrared spectroscopy in combination with chemometrics as a process analytical technology (PAT) tool for on-line quantitative monitoring of alcohol precipitation
CN102106939B (zh) 一种测定六味地黄丸浓缩丸提取浓缩液比重及马钱苷、丹皮酚含量的方法
CN105548026A (zh) 一种快速的郁金药材质控检测方法
CN104833651A (zh) 金银花浓缩过程在线实时放行检测方法
CN108241033B (zh) 一种快速检测麦冬醇提取液中6个质量指标物质含量的方法及应用
CN101984343A (zh) 一种中药大孔树脂分离纯化过程关键点的判别方法
CN103115892A (zh) 应用近红外光谱分析技术制备银杏叶提取物的方法
CN102028710B (zh) 一种测定华蟾素醇沉液中吲哚类生物碱含量的方法
CN104390926B (zh) 一种穿心莲浓缩脱色过程在线快速检测方法
CN108051396A (zh) 一种心可舒片有效成分含量的快速检测方法
CN104865322A (zh) 一种栀子萃取液浓缩过程快速检测方法
CN102106950B (zh) 一种测定女金中药提取浓缩液中黄芪苷的方法
CN102323236B (zh) 近红外光谱测定苦参提取过程多种成分含量的方法
CN103913434A (zh) 一种丹红注射液水沉过程在线检测方法
CN110346323B (zh) 一种基于近红外光谱技术在线检测华盖散浓缩液的方法
CN104297441B (zh) 一种红外光谱在线质量监测控制体系在蒙药制剂中的应用
CN107655849A (zh) 一种凉茶近红外在线检测方法
CN103335960A (zh) 一种华蟾素提取浓缩过程关键指标的快速检测方法
Yan et al. Rapid detection of cAMP content in red jujube using near-infrared spectroscopy
CN105784951A (zh) 一种六味地黄丸浓缩丸生药粉多指标快速检测方法
CN110308226A (zh) 一种发酵虫草菌粉生产全过程链质量控制快速检测方法
CN107036997A (zh) 利用近红外光谱法快速检测气滞胃痛颗粒的制备过程的方法及应用

Legal Events

Date Code Title Description
DD01 Delivery of document by public notice

Addressee: Wang Ruiqing

Document name: Notification of Passing Examination on Formalities

C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant