CN109655542B - 一种木蝴蝶配方颗粒uplc特征图谱的构建方法及其应用 - Google Patents

Abstract

本发明公开了一种木蝴蝶配方颗粒UPLC特征图谱的构建方法及其应用，包括步骤：对照品溶液的制备、供试品溶液的制备、超高相液相色谱测定、特征图谱的建立、相似度评价等。本发明建立了木蝴蝶配方颗粒的UPLC特征图谱，并通过液质联用技术对木蝴蝶配方颗粒6种共有成分进行了快速定性及定量分析。该方法为木蝴蝶配方颗粒的质量评价和控制提供了较为全面、系统、有效的快速评价方法。

Description

一种木蝴蝶配方颗粒UPLC特征图谱的构建方法及其应用

技术领域

本发明属于药物分析及药物质量控制技术领域，具体涉及一种木蝴蝶配方颗粒特征图谱的构建方法及其应用。

背景技术

木蝴蝶(Oroxyli semen)为紫葳科木蝴蝶属木蝴蝶Oroxylumindicum(L)Vent.的干燥成熟种子，又名千层纸、玉蝴蝶等，始载于《滇南本草》，具有几百年的药用历史。中医理论认为木蝴蝶性味苦、甘、凉，入肺、肝、胃经，有清肺利咽、疏肝和胃之功效，用于肺热咳嗽，喉痹，喑哑，肝胃气痛，常用于以咳嗽为主症的上呼吸道感染、急性支气管炎、慢性咽炎、肺炎后期等症。现代医学研究表明，木蝴蝶的主要药效成分为黄酮类物质，具有抗炎、抗氧化、抗癌、降糖、抑制病毒及肿瘤生长等多种药理作用。具有药理活性的化学成分为脂肪油、黄芩苷元-7-葡萄糖苷、黄芩苷、黄芩苷元、千层纸素A、木蝴蝶苷B、白杨素、白杨素-7-葡萄糖醛酸苷等黄酮类化合物。

木蝴蝶配方颗粒是由单味中药饮片经水提、浓缩、干燥、制粒而成，目前并没有建立有效控制木蝴蝶配方颗粒的质量的方法。随着现代分析技术的发展，特征图谱作为一种体现中药化学成分整体特征的质量评价方法正在广泛应用，因此，有必要建立了一种快速鉴别木蝴蝶配方颗粒的UPLC特征图谱方法，为有效控制和较全面评价木蝴蝶配方颗粒的质量提供依据。

发明内容

本发明的目的在于提供一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，该方法为木蝴蝶配方颗粒的质量评价和控制提供了较为全面、系统、有效的快速评价方法。

本发明是通过以下技术方案实现：

一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，包括如下步骤：

(a)混合对照品溶液的制备：

分别精密称取白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B对照品，加甲醇溶解定容至25mL，即得白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷、木蝴蝶苷B浓度分别为12.594μg/mL、69.600μg/mL、32.286μg/mL、128.256μg/mL、31.512μg/mL、60.276μg/mL的混合对照品储备液，即得；

(b)供试品溶液的制备：

取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入60-80％乙醇50-100mL，称定重量，超声处理15-30分钟，取出，放冷，用60-80％乙醇补足减失的重量，摇匀，滤膜滤过，即得；

(c)超高相液相色谱测定：分别精密吸取对照品溶液与供试品溶液各1μl，注入超高效液相色谱仪，测定，记录色谱图，其中色谱条件为：

以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.1％磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

(d)特征图谱的建立：

将多批木蝴蝶配方颗粒的色谱图，采用“中药色谱特征图谱相似度评价系统”进行数据匹配，建立木蝴蝶配方颗粒的特征图谱，确定具有6个共有峰，以1号峰为参照峰S，各特征峰的相对保留时间分别为：峰1：1.00、峰2：1.20、峰3：1.22、峰4：1.35、峰5：1.52、峰6：1.68，相对保留时间在规定值的±10％之内；

(e)相似度评价：计算各批次木蝴蝶配方颗粒的相似度，木蝴蝶配方颗粒的相似度均大于0.90。

本发明的木蝴蝶配方颗粒是由单味中药饮片经水提、浓缩、干燥、制粒而成，优选的，所述木蝴蝶配方颗粒的制备方法为：取木蝴蝶饮片，以水为提取溶媒，加水煎煮两次，每次加投料量5～10倍量的水，煎煮0.5～2小时，合并滤液，滤过；滤液真空减压浓缩至相对密度为为1.02～1.10的清膏；取清膏进行喷雾干燥，得干浸膏粉，干浸膏粉干法制粒得粒度在16～40目的木蝴蝶配方颗粒。

本发明对供试品溶液的制备条件进行了优化，对不同提取溶剂、不同提取方式、不同提取时间及不同料液比进行考察，优选的，步骤(b)供试品溶液的制备：取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入70％乙醇50mL，称定重量，超声处理30分钟，功率为300W，频率为40KHz，取出，放冷，用70％乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

作为本发明进一步优选的技术方案，所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，还包括步骤(f)特征图谱峰指认：采用UPLC-Q-TOF-MS方法确认6个共有峰成分，其中，色谱条件为：

以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.01％甲酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

质谱条件为：

电喷雾离子源(ESI⁺)，喷雾电压为1.55kV；脱溶剂温度为597℃，离子源温度为148℃；脱溶剂气流为997L/hr；氮气流为1L/hr；质谱扫描采用多反应检测MRM扫描模式。监测离子及相关电压参数设定如下表所示：

作为本发明进一步优选的技术方案，所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，其特征在于，还包括步骤(g)含量测定：按照供试品溶液的制备方法及质谱条件测定，计算木蝴蝶配方颗粒中的白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B 6个共有峰成分的含量。

本发明还提供了上述木蝴蝶配方颗粒UPLC特征图谱的构建方法在木蝴蝶配方颗粒质量评价和检测中的应用。

本发明与现有技术相比，具有如下有益效果：

本发明建立了木蝴蝶配方颗粒的UPLC特征图谱，并通过液质联用技术对木蝴蝶配方颗粒6种共有成分进行了快速定性及定量分析。该方法为木蝴蝶配方颗粒的质量评价和控制提供了较为全面、系统、有效的快速评价方法。

附图说明

图1为16批木蝴蝶配方颗粒UPLC特征图谱叠加图；

图2为16批木蝴蝶配方颗粒UPLC特征图谱共有模式图；

具体实施方式

下面通过具体实施方式来进一步说明本发明，以下实施例为本发明具体的实施方式，但本发明的实施方式并不受下述实施例的限制。

1、仪器与试药

1.1仪器

Waters液质联用仪(液相部分ACQUITY UPLC H-Class、质谱部分Waters Xevo TQDMS)；百万分之一天平(METTLER TOLEDO，XP26)、万分之一天平(METTLER TOLEDO，ME204E)，Agilent SB C₁₈色谱柱(2.1×100mm,1.8μm)，电热恒温水浴锅(上海一恒科技有限公司，HWS-28)，数控超声波清洗器(昆山市超声仪器有限公司，KQ500DE)，超纯水系统(默克股份有限公司，Milli-Q-Direct)。

1.2对照品和试剂

乙腈(德国merck公司，色谱纯)；乙醇(西陇科学股份有限公司，分析纯)；甲醇(西陇科学股份有限公司，分析纯)；磷酸(天津科密欧化学试剂有限公司，色谱纯)；甲酸(天津科密欧化学试剂有限公司，色谱纯)；水为实验室自制水。

白杨素(批号：111701-200501)、木蝴蝶苷B(批号：111915-201603，含量：91.9％)、黄芩苷(批号：110715-201720，含量：93.5％)、黄芩素(批号：111595-201607，含量：98.5％)对照品购自中国食品药品检定研究院；白杨素-7-O-β-D葡萄糖醛酸苷(批号：wkq17072702，含量：98％)对照品购自四川省维克奇生物科技有限公司；木蝴蝶苷A(批号：17021703，含量：98％)对照品购自成都普菲德生物技术有限公司。

1.3药材

16批木蝴蝶配方颗粒由广东一方制药有限公司生产，制备方法为：取木蝴蝶饮片，以水为提取溶媒，加水煎煮两次，第一次加投料量10倍量的水，煎煮1.5小时，第二次加投料量8倍量的水，煎煮1小时，合并滤液，滤过；滤液真空减压浓缩至相对密度为为1.06的清膏；取清膏进行喷雾干燥，得干浸膏粉，干浸膏粉干法制粒得粒度在16～40目的木蝴蝶配方颗粒。

2、方法与结果

2.1混合对照品溶液的制备：分别精密称取白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B对照品，加甲醇溶解定容至25mL，即得白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷、木蝴蝶苷B浓度分别为12.594μg/mL、69.600μg/mL、32.286μg/mL、128.256μg/mL、31.512μg/mL、60.276μg/mL的混合对照品储备液，即得。

2.2供试品溶液的制备：取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入70％乙醇50mL，称定重量，超声处理30分钟，功率为300W，频率为40KHz，取出，放冷，用70％乙醇补足减失的重量，摇匀，滤膜滤过，即得。

2.3超高相液相色谱测定

分别精密吸取对照品溶液与供试品溶液各1μl，注入超高效液相色谱仪，测定，记录色谱图，其中色谱条件为：

以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.1％磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

3、供试品溶液制备方法考察

3.1提取溶媒考察

本次实验分别考察了不同提取溶剂对木蝴蝶配方颗粒特征图谱的影响，选取甲醇、70％乙醇、乙醇作为提取溶剂，对不同提取溶剂的样品溶液特征图谱进行测定，确定最佳提取溶剂。

取木蝴蝶配方颗粒，研细，取0.1g，共6份，精密称定，置锥形瓶中，精密加入甲醇、70％乙醇、乙醇50ml，称定重量，超声处理30分钟，功率为300W，频率为40KHz，取出，放冷，用甲醇、70％乙醇、乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

通过对比3种不同提取溶剂的色谱图可发现，三种提取溶剂色谱图的色谱峰数目与色谱峰峰形一致，总峰面积大小的顺序为：70％乙醇>甲醇>乙醇；综合考虑各溶剂的提取能力、色谱峰峰形、峰面积/称样量的值和溶剂效应等，最终采用70％乙醇作为提取溶剂。

3.2提取溶剂方式考察

本次实验考察不同提取方式对木蝴蝶配方颗粒特征图谱的影响，选取加热回流与超声处理两种方式。对不同提取方式的样品溶液特征图谱进行测定，确定最佳提取方式。

取木蝴蝶配方颗粒，研细，取0.1g，共2份，精密称定，置锥形瓶中，精密加入70％乙醇50ml，称定重量，分别加热回流、超声处理(功率为300W，频率为40KHz)30分钟，取出，放冷，用70％乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

对比了两种提取方式，不同提取方式的色谱图色谱峰数目与色谱峰峰形一致，考虑到实验操作的简便性，故选择超声处理作为提取方式。

3.3提取时间考察

本次实验考察不同提取时间对木蝴蝶配方颗粒特征图谱的影响，选取超声处理提取时间分别为：15分钟、30分钟、60分钟。对不同提取时间的样品溶液特征图谱进行测定，确定最佳提取时间。

取木蝴蝶配方颗粒，研细，取0.1g，共3份，精密称定，置锥形瓶中，精密加入70％乙醇50ml，称定重量，分别超声处理15分钟、30分钟、60分钟，取出，放冷，用70％乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

通过对比不同提取时间木蝴蝶配方颗粒色谱图，可发现随着超声处理时间的增加，色谱图的色谱峰数目与色谱峰峰形一致，色谱图总面积与称样量比值基本一致，说明超声处理15分钟已经能提取完全，为保证方法的耐用性，提取时间选择30分钟。

3.4提取溶剂用量的考察

本次实验考察不同提取溶剂用量对木蝴蝶配方颗粒特征图谱的影响，选取提取溶剂用量分别为：15ml、25ml、50ml、100ml。对不同提取溶剂用量的样品溶液特征图谱进行测定，确定最佳提取溶剂用量。

取木蝴蝶配方颗粒，研细，取0.1g，共4份，精密称定，置锥形瓶中，分别精密加入70％乙醇15ml、25ml、50ml、100ml，称定重量，分别超声处理(功率为300W，频率为40KHz)30分钟，取出，放冷，用70％乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

通过对比不同提取溶剂用量木蝴蝶配方颗粒色谱图，可发现不同提取溶剂用量，色谱图的色谱峰数目与色谱峰峰型一致，色谱图总面积与称样量比值基本一致，说明提取溶剂用量为15ml已提取完全，为保证方法的耐用性及选择合适的进样浓度，提取溶剂用量选择50ml。

根据上述实验结果，木蝴蝶配方颗粒特征图谱样品前处理方法可确定为：取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入70％乙醇50mL，称定重量，超声处理30分钟，取出，放冷，用70％乙醇补足减失的重量，摇匀，滤膜滤过，即得。

4、色谱条件优化

4.1色谱柱的选择

经比较发现，Agilent SB C18柱(2.1×100mm,1.8μm)柱分离效果最好，峰型更好，因此，选择Agilent SB C18柱(2.1×100mm,1.8μm)柱为色谱柱。

4.2检测波长的选择

采用二极管阵列检测器对样品进行全波长扫描，发现在276nm处基线稳定，各峰峰面积较高，因此，检测波长定为276nm。

4.3流动相及洗脱条件的选择

本发明流动相及洗脱条件的选择采用如下方法进行比较：

流动相条件优化1：

以Agilent SB C18柱(2.1×100mm,1.8μm)柱为色谱柱；以乙腈为流动相A，以0.1％磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.4mL/min，进样体积为1μL，柱温为30℃。

流动相条件优化2：

以Agilent SB C18柱(2.1×100mm,1.8μm)柱为色谱柱；以乙腈为流动相A，以0.1％磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.35mL/min，进样体积为1μL，柱温为30℃。

流动相条件优化3：

以Agilent SB C18柱(2.1×100mm,1.8μm)柱为色谱柱；以乙腈为流动相A，以0.1％磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃。

由三种洗脱条件得到的色谱图结果可以看出：优化条件1与优化条件2中峰2包含了其他色谱峰，峰纯度没有达到要求，且没有达到基线分离要求，而优化条件3中个各色谱峰满足基线分离要求，故选择优化条件3进行实验。

5、方法学考察

5.1精密度试验取木蝴蝶配方颗粒供试品溶液，在上述色谱条件下连续进样6次，测定峰面积，计算峰面积相对标准偏差RSD，结果RSD<3，表明仪器精密度良好。

5.2重复性试验取同一批木蝴蝶配方颗粒(批号：8016072)，按照供试品制备方法平行制备6份供试品溶液，按照上述色谱条件进样，测定峰面积，计算峰面积相对标准偏差RSD，结果RSD<3，实验结果显示本方法重现性良好。

5.4稳定性试验精密吸取同一供试品溶液，分别于0、2、4、6、8、12h后进样，测定峰面积，计算峰面积相对标准偏差RSD，结果RSD<3，实验结果显示供试品溶液在12h内稳定性良好。

6、木蝴蝶配方颗粒特征图谱的建立

取16批木蝴蝶配方颗粒样品，按照“2.2”项下的方法制备供试品溶液，进样1μL，测定得色谱图谱。将图谱导入“中药色谱指纹图谱相似度评价系统”进行数据进行数据匹配，建立木蝴蝶配方颗粒的特征图谱，叠加图见图1，生成的共有模式见图2。

以1号峰为参照峰，标示出木蝴蝶配方颗粒UPLC-TUV色谱图中的6个共有峰，共有峰的峰面积之和大于90％。计算各色谱峰与参照峰相对保留时间的RSD值。结果表明，样品各色谱图共有峰的相对保留时间RSD值均小于3.0％；而相对峰面积的RSD值相差较大，见表1、表2。

表1 16批木蝴蝶特征图谱共有峰的相对保留时间

表2 16批木蝴蝶特征图谱共有峰的相对峰面积

分别计算16批木蝴蝶配方颗粒样品的相似度，见表3，16批木蝴蝶配方颗粒相似度均大于0.90，说明木蝴蝶配方颗粒化学成分一致性较好，质量比较稳定。

表3 16批木蝴蝶配方颗粒相似度

7、特征图谱峰指认：

采用UPLC-Q-TOF-MS方法确认6个共有峰成分，其中，色谱条件为：以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.01％甲酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

质谱条件为：电喷雾离子源(ESI⁺)，喷雾电压为1.55kV；脱溶剂温度为597℃，离子源温度为148℃；脱溶剂气流为997L/hr；氮气流为1L/hr；质谱扫描采用多反应检测MRM扫描模式。Cone参数和Collision质谱参数优化通过蠕动针泵注射的方法调谐各标准化合物。所得参数见表4。

表4 LC-MS检测标准化合物的质谱参数

备注：其中带有“*”为定量离子。

8、含量测定：

8.1线性关系精密量取“2.1“项下的混合对照品溶液，按稀释倍数为1、1.2、1.5、2、3、6依次逐级稀释，定容，制得6个浓度梯度的系列混合对照品溶液，按照“7”中的色谱条件和质谱条件测定，以峰面积对进样量浓度进行线性回归，绘制标准曲线，求得各成分回归方程、相关系数、线性范围，见表5。

表5木蝴蝶配方颗粒中6个成分的线性关系考察结果

8.2按照“2.2”中供试品溶液的制备方法及“7”中的色谱条件和质谱条件测定，计算16批木蝴蝶配方颗粒中的白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B 6种成分的含量，结果见表6。

表6样品含量测定结果(mg·g^-1，n＝2)

Claims (5)

1.一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，其特征在于，包括如下步骤：

（a）混合对照品溶液的制备：

分别精密称取白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B对照品，加甲醇溶解定容至25mL，即得白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷、木蝴蝶苷B浓度分别为12.594μg/mL、69.600μg/mL、32.286μg/mL、128.256μg/mL、31.512μg/mL、60.276μg/mL的混合对照品储备液，即得；

（b）供试品溶液的制备：

取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入60-80%乙醇50-100mL，称定重量，超声处理15-30分钟，取出，放冷，用60-80%乙醇补足减失的重量，摇匀，滤膜滤过，即得；

所述木蝴蝶配方颗粒的制备方法为：取木蝴蝶饮片，以水为提取溶媒，加水煎煮两次，每次加投料量5~10倍量的水，煎煮0.5~2小时，合并滤液，滤过；滤液真空减压浓缩至相对密度为1.02~1.10的清膏；取清膏进行喷雾干燥，得干浸膏粉，干浸膏粉干法制粒得粒度在16～40目的木蝴蝶配方颗粒；

（c）超高相液相色谱测定：分别精密吸取对照品溶液与供试品溶液各1µl，注入超高效液相色谱仪，测定，记录色谱图，其中色谱条件为：

以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.1%磷酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

（d）特征图谱的建立：

将多批木蝴蝶配方颗粒的色谱图，采用“中药色谱特征图谱相似度评价系统”进行数据匹配，建立木蝴蝶配方颗粒的特征图谱，确定具有6个共有峰，以1号峰为参照峰S，各特征峰的相对保留时间分别为：峰1：1.00、峰2：1.20、峰3：1.22、峰4：1.35、峰5：1.52、峰6：1.68，相对保留时间在规定值的±10%之内；

（e）相似度评价：计算各批次木蝴蝶配方颗粒的相似度，木蝴蝶配方颗粒的相似度均大于0.90；

（f）特征图谱峰指认：采用UPLC-Q-TOF-MS方法确认6个共有峰成分，其中，色谱条件为：

以十八烷基硅烷键合硅胶为填充剂；以乙腈为流动相A，以0.01%甲酸为流动相B，按下表中的规定进行梯度洗脱；检测波长为276nm，流速为0.28mL/min，进样体积为1μL，柱温为30℃；

质谱条件为：

电喷雾离子源ESI⁺，喷雾电压为1.55kV；脱溶剂温度为597℃，离子源温度为148℃；脱溶剂气流为997L/hr；氮气流为1L/hr；质谱扫描采用多反应检测MRM扫描模式。

2.根据权利要求1所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，其特征在于，步骤(b) 供试品溶液的制备：

取木蝴蝶配方颗粒，研细，取0.1g，精密称定，置锥形瓶中，精密加入70%乙醇50mL，称定重量，超声处理30分钟，功率为300W，频率为40KHz，取出，放冷，用70%乙醇补足减失的重量，摇匀，用0.22μm滤膜滤过，即得。

3.根据权利要求1所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，其特征在于，还包括步骤（g）含量测定：按照供试品溶液的制备方法及质谱条件测定，计算木蝴蝶配方颗粒中的白杨素、黄芩素、白杨素-7-O-β-D葡萄糖醛酸苷、木蝴蝶苷A、黄芩苷和木蝴蝶苷B 6个共有峰成分的含量。

4.根据权利要求3所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法，其特征在于，质谱条件中监测离子及相关电压参数设定如下表所示：

5.权利要求1-4任一项所述的一种木蝴蝶配方颗粒UPLC特征图谱的构建方法在木蝴蝶配方颗粒质量评价和检测中的应用。

Images