2.845

2023影响因子

(CJCR)

  • 中文核心
  • EI
  • 中国科技核心
  • Scopus
  • CSCD
  • 英国科学文摘

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

一种双判别器GAN的古彝文字符修复方法

陈善雄 朱世宇 熊海灵 赵富佳 王定旺 刘云

赖轩, 曲延云, 谢源, 裴玉龙. 基于拓扑一致性对抗互学习的知识蒸馏. 自动化学报, 2023, 49(1): 102−110 doi: 10.16383/j.aas.c200665
引用本文: 陈善雄, 朱世宇, 熊海灵, 赵富佳, 王定旺, 刘云. 一种双判别器GAN的古彝文字符修复方法. 自动化学报, 2022, 48(3): 853−864 doi: 10.16383/j.aas.c190752
Lai Xuan, Qu Yan-Yun, Xie Yuan, Pei Yu-Long. Topology-guided adversarial deep mutual learning for knowledge distillation. Acta Automatica Sinica, 2023, 49(1): 102−110 doi: 10.16383/j.aas.c200665
Citation: Chen Shan-Xiong, Zhu Shi-Yu, Xiong Hai-Ling, Zhao Fu-Jia, Wang Ding-Wang, Liu Yun. A method of inpainting ancient Yi characters based on dual discriminator generative adversarial networks. Acta Automatica Sinica, 2022, 48(3): 853−864 doi: 10.16383/j.aas.c190752

一种双判别器GAN的古彝文字符修复方法

doi: 10.16383/j.aas.c190752
基金项目: 国家自然科学基金 (61603310), 国家社会科学基金(19BYY171), 重庆市自然科学基金(cstc2019jcyj-msxm2550), 模式识别国家重点实验室开放课题(201900010), 中央高校基本科研业务费(XDJK2018B020), 重庆市教育委员会科学技术研究计划青年项目(KJQN201801901, KJQN201801902)资助
详细信息
    作者简介:

    陈善雄:西南大学计算机与信息科学学院教授. 2013年获得重庆大学计算机科学与技术博士学位. 主要研究方向为数据挖掘, 模式识别, 压缩感知. E-mail: csxpml@163.com

    朱世宇:重庆工程学院计算机与物联网学院博士研究生. 2012年获得重庆大学电子与通信工程专业硕士学位. 主要研究方向为计算机视觉和自然语言处理. 本文通信作者.E-mail: herb_711@icloud.com

    熊海灵:西南大学商贸学院教授. 2007年获得西南大学土壤学博士学位. 主要研究方向为农业信息化与数字农业, 农业生态与环境过程计算机模拟.E-mail: xionghl@swu.edu.cn

    赵富佳:西南大学计算机与信息科学学院硕士研究生. 2017年获得盐城工学院计算机科学与技术学士学位. 主要研究方向为模式识别, 图像处理, 深度学习.E-mail: zfj_ai_python@163.com

    王定旺:西南大学计算机与信息科学学院软件学院硕士研究生. 2017年获得重庆理工大学计算机科学与工程学院学士学位. 主要研究方向为模式识别, 图像处理, 深度学习.E-mail: 18996632565@163.com

    刘云:贵州工程应用技术学院彝学研究院副研究员. 主要研究方向为彝文古籍整理, 彝文信息化. E-mail: ly_3692022@163.com

A Method of Inpainting Ancient Yi Characters Based on Dual Discriminator Generative Adversarial Networks

Funds: Supported by National Natural Science Foundation of China (61603310), National Social Science Foundation of China (19BYY171), Natural Science Foundation of Chongqing (cstc2019jcyj-msxm2550), Open Projects Program of National Laboratory of Pattern Recognition (201900010), Fundamental Research Funds for the Central Universities, of China (XDJK2018B020), and Youth Project of Science and Technology Research Program of Chongqing Education Commission(KJQN201801901, KJQN201801902)
More Information
    Author Bio:

    CHEN Shan-Xiong Professor at the College of Computer and Information Science, Southwest University. He received his Ph.D. degree in computer science from Chongqing University in 2013. His research interest covers data mining, pattern recognition, and compressed sensing

    ZHU Shi-Yu Ph.D. candidate at the College of Computer and Internet of Things, Chongqing Institute of Engineering College. He received his master degree in electronics and communication engineering from Chongqing University in 2012. His research interest covers machine vision and natural language processing (NLP). Corresponding author of this paper

    XIONG Hai-Ling Professor at the Business College, Southwest University. He received his Ph.D. degree in soil science from Southwest University in 2007. His research interest covers agricultural informatization and digital agriculture, computer simulation of agricultural ecology and environmental process

    ZHAO Fu-Jia Master student at the School of Computer and Information Science, Southwest University. He received his bachelor degree in computer science and technology from Yancheng Institute of Technology in 2017. His research interest covers pattern recognition, image processing, and deep learning

    WANG Ding-Wang Master student at the School of Computer and Information Science, Southwest University. He received his bachelor degree in computer science and engineering from Chongqing University of Technology in 2017. His research interest covers pattern recognition, image processing, and deep learning

    LIU Yun Associate researcher at the Resrarch Institute of Yi Nationality Studies, Guizhou University of Engineering Science. His research interest covers arrangement of Yi ancient books, and Yi language informatization

  • 摘要: 在中国, 彝文古籍文献日益流失而且损毁严重, 由于通晓古彝文的研究人员缺乏, 使得古籍恢复工作进展十分缓慢. 人工智能在图像文本领域的应用, 为古籍文献的自动修复提供可能. 本文设计了一种双判别器生成对抗网络(Generative adversarial networks with dual discriminator,D2GAN), 以还原古代彝族字符中的缺失部分. D2GAN是在深度卷积生成对抗网络的基础上, 增加一个古彝文筛选判别器. 通过三个阶段的训练来迭代地优化古彝文字符生成网络, 以获得古彝文字符的文字生成器. 根据筛选判别器的损失结果优化D2GAN模型, 并使用生成的字符恢复古彝文中丢失的笔画. 实验结果表明, 在字符残缺低于1/3的情况下, 本文提出的方法可使文字笔画的修复率达到77.3%, 有效地加快了古彝文字符修复工作的进程.
  • 稀土是由镧系元素, 钪和钇等17种元素组成, 且以共生矿形式存在.稀土元素素有"工业维生素"的美称, 已成为极其重要的战略资源.目前我国稀土萃取分离工艺技术已达到世界先进水平[1], 但稀土工业生产自动化水平普遍较低, 仍停留在"离线分析、手工调整、经验控制"状态, 导致企业生产效率低、资源消耗大、产品质量不稳定, 成为制约稀土工业发展的瓶颈.因此, 部分学者对稀土萃取过程建模与控制进行了大量研究, 并取得了一定成果.公锡泰等[2]在稀土萃取过程的测控敏感点安装在线分析仪, 并在"四出口"萃取生产线上开展自动反馈控制试验, 取得了较好的结果.但在线分析仪因频繁堵塞, 工作人员保养习惯等原因导致难以长时间使用. Giles等[3]提出了基于人工神经网络的稀土萃取控制模型.贾文君等[4]考虑稀土萃取过程的多段动态特性的基础上, 建立了一种具有状态滞后的双线性动态模型.该模型忽略了级间相互作用, 损失了部分动态特性, 难以实现萃取过程流量自动控制, 现场仍采用操作员手动调节模式.黄桂文[5]和Jia等[6]研究了利用稀土离子的特征颜色进行串级萃取工艺控制的可行性, 并总结归纳了萃取剂、料液、洗涤液流量调节的经验公式.

    Chai等[7]和杨辉等[8]提出了具有两层体系结构的基于案例推理技术的稀土萃取过程优化控制系统.该方法采用案例推理技术实现上层优化器的功能, 对底层回路值进行预设定, 并利用元素组分含量的离线化验值和预测值对预设定值进行修正, 底层流量回路采用PID控制器实现闭环控制.但是上述方法的优化层只采用了案例推理技术. Yang等[9]提出了稀土萃取过程组分含量分布控制方法, 首先根据每一级建立渐进阶梯式模型, 结合质量平衡模型和监测级组分含量变化趋势, 通过动态补偿萃取剂和洗涤剂流量, 使两端出口产品满足纯度要求, 该方法忽略了较多实际因素, 难以应用于生产实践.针对稀土萃取过程非线性和动态特性, 杨辉等[10]采用自适应神经模糊推理系统(ANFIS)对CePr/Nd萃取过程进行描述, 运用广义预测控制方法(GPC)实现各控制流量的优化控制.文献[11]提出一种稀土萃取过程多模型建模和组分含量预测控制方法, 基于组分含量预测值对CePr/Nd萃取过程中萃取剂和洗涤剂流量进行动态补偿, 以保证两端出口产品的纯度.但一旦工作环境发生变化, 给定参数的PID控制器无法自适应调整, 不能达到最优效果.很明显, 目前还没有较好的稀土萃取过程萃取量和洗涤量控制方法.

    针对CePr/Nd稀土萃取过程入矿条件各参数的重要性程度不一致特点, 本文首先提出基于特征属性加权的案例推理预设定模型以实现对萃取剂/洗涤剂的流量预设定, 然后针对镨/钕稀土萃取过程的稀土溶液颜色与组分含量密切相关的特点[12], 利用机器视觉技术获得稀土溶液颜色值, 建立基于LS-SVM的组分含量软测量模型, 在此基础上, 采用基于模糊的智能动态补偿模型补偿萃取剂和洗涤剂流量设定值.

    由于分馏萃取可同时得到二个高纯度、高收率的产品, 且单级萃取不能达到有效的分离效果, 因此稀土工业生产普遍采用串级分馏萃取方式(若干级萃取槽串联起来), 使含料水相与有机相多次接触, 并从混合稀土溶液中分离、富集、提取所需具有一定纯度和收率的单一稀土元素. CePr/Nd稀土萃取分离工艺如图 1所示.按照Ce、Pr和Nd分离元素的化学特性, 该萃取过程将Ce和Pr设为难萃组分$ B $, 将Nd设为易萃组分$ A $.在稀土元素Pr和Nd位置间切割, 料液进料为有机相进料方式.整个稀土萃取生产线包括$ {n} $级萃取段和$ m $级洗涤段共$ {n+m} $级萃取槽组成, 每级萃取槽上面溶液为有机相下面溶液为水相.该生产工艺流程第1级萃取槽加入萃取剂$ u_1 $, 第$ {n} $级萃取槽加入料液$ u_3 $; 第$ {n+m} $级萃取槽加入洗涤剂$ u_2 $.由于萃取槽体独特的结构设计和搅拌力的作用, 使得各级萃取槽中有机相从左向右顺向流动、水相从右向左逆向流动.稀土溶液在萃取剂和洗涤剂的作用下经过萃取段和洗涤段各级萃取槽的交换和纯化, 最后从萃取段第1级水相出口得到目标纯度为$ {Y_B} $的难萃$ B $产品, 洗涤段第$ {n+m} $级有机相出口得到目标纯度为$ {Y_A} $的易萃$ A $产品.

    图 1  CePr/Nd萃取生产过程工艺流程
    Fig. 1  CePr/Nd extraction process

    该萃取工艺包含料液、萃取剂和洗涤剂等进料, 水相、有机相等出料, 料液指混合稀土溶液, 萃取剂用于控制有机相出口易萃组分的纯度(即组分含量), 洗涤剂(酸液)用于控制水相出口难萃组分的纯度.

    CePr/Nd稀土萃取过程的控制目标是指产量、有机相出口易萃稀土元素$ A $纯度和水相出口难萃稀土元素$ B $纯度等指标满足生产要求.稀土萃取过程的入矿条件及扰动量包括萃取流程原料处理量、料液稀土元素配分、皂化度、萃取剂浓度、料液浓度和洗涤剂浓度等工艺参数.其控制变量包括稀土萃取过程添加的萃取剂、洗涤剂、料液等流量值(简称萃取量、洗涤量、料液量).当稀土萃取生产工艺产量指标不变时, 料液的流量设定值一般不作调整, 即料液流量与稀土处理量指标直接关联, 一般不作为工业稀土萃取生产过程的调整量.

    稀土萃取生产过程对关键工艺参量稀土元素组分含量缺乏有效的在线检测手段; 且萃取机理复杂、工艺流程长, 难以建立精确的萃取量(洗涤量)与稀土元素组分含量的关系数学模型, 并依据组分含量实现萃取量(洗涤量)优化控制.因此, 实际稀土萃取工业生产过程中工作人员一般根据萃取流程原料处理能力、料液稀土元素配分、皂化度、萃取剂浓度、料液浓度和酸液浓度等入矿条件预先设定萃取量和洗涤量, 然后当入矿条件或扰动量波动导致萃取溶液的元素组分含量变化时, 工作人员再根据CePr/Nd萃取溶液颜色或离线化验值凭借操作经验调整萃取量和洗涤量.例如, 当工况异常导致难萃组分$ B $纯度的化验值或易萃组分$ A $纯度的化验值低于目标值时, 操作人员需根据组分含量化验值调整萃取剂和洗涤剂的流量设定值.但因为萃取工艺流程长、离线化验流程(一班组等隔化验2次)等因素导致化验延时较长, 操作人员难以及时有效地调整萃取量和洗涤量, 导致CePr/Nd萃取工艺生产效率较差.当工况变化频繁时, 操作人员不能及时准确地判断工况和调整萃取量(洗涤量)设定值, 将会使控制系统性能变坏甚至瘫痪, 常常造成故障工况[13-14].

    根据稀土萃取工业生产操作实践经验可知, 操作人员首先根据入矿条件预先设定萃取剂和洗涤剂流量值, 在萃取生产过程中再根据稀土元素组分含量的变化对流量设定值动态微调补偿, 使萃取生产稳定运行.由于生产工艺指标稀土元素组分含量与底层控制回路的萃取量和洗涤量设定值密切相关, 它们间的静态、动态特性难以用精确模型描述.为此, 本文结合调节控制和过程运行优化思想, 采用两层结构的实时优化过程模型[15-16], 上层采用非线性专家经验模型优化基础控制回路的萃取量(洗涤量)设定值, 下层通过PID基础控制系统跟踪设定值, 尽可能使过程运行在经济优化状态.由于PID控制能够实时有效地调节萃取剂和洗涤剂流量在给定值上, 因此本文重点研究上层优化结构模型, 该模型首先模拟操作人员专家行为, 采用案例推理方法预设定萃取剂和洗涤剂流量值, 其次, 在萃取生产过程中根据稀土元素组分含量变化采用模糊推理技术对预设定值动态微调补偿, 从而实时优化萃取量和洗涤量设定值.

    综上所述, 本文提出一种案例推理静态设定和模糊推理动态补偿相结合的CePr/Nd萃取过程萃取量和洗涤量优化控制方法, 其总体方案结构如图 2所示, 主要包括萃取剂/洗涤剂流量预设定模型、稀土元素组分含量软测量模型及模糊推理智能补偿模型.借鉴人工操作专家行为, 预设定模型采用案例推理方法根据稀土生产入矿条件和生产工艺指标从历史操作案例库中按照一定的匹配原则选择最接近实时工况的流量作为流量预设定值; 然后根据稀土溶液组分含量的变化动态补偿流量预设定值, 包括基于稀土溶液颜色特征的组分含量软测量模型和模糊推理智能补偿模型; 组分含量软测量模型根据机器视觉系统采集到的稀土溶液图像颜色特征软测量稀土溶液的组分含量, 模糊推理智能补偿模型根据软测量模型计算的组分含量与目标组分含量的比较值对预设定的萃取量和洗涤量动态补偿; 最终得到一个满足实时工况要求的萃取剂和洗涤剂流量设定值, 作为底层控制回路的给定值, 由PID控制器实现对萃取剂和洗涤剂流量给定值的跟随控制.此外, 为使各个模型更加准确, 离线检测到的组分含量数据也要作为各个模型数据库更新、校正的参考数据源.

    图 2  萃取量和洗涤量优化控制方案
    Fig. 2  Optimal control scheme for the flow rate of the extractant and the detergent

    稀土萃取实际生产过程工作人员通常凭借操作经验根据萃取生产流程原料处理量、料液稀土元素配分、皂化度、萃取剂浓度、料液浓度和酸液浓度等入矿条件预先设定萃取剂和洗涤剂流量值.案例推理(Case based reasoning, CBR)是借鉴人类思维方式对问题进行处理的人工智能推理方法, 并从案例中获取专家知识用于解决新问题.基于此, 采用案例推理技术预设定稀土萃取过程萃取剂/洗涤剂的最初值.

    基于CBR技术的分类器主要根据已知类别中的历史案例来判断一个新案例属于哪个类.分类精度是评价案例推理技术的一个重要指标, 案例检索、案例重用、案例更新与删除都对分类精度有着重要的影响.由于每个入矿条件的重要程度不一致, 所以利用CBR检索时需要考虑给每个入矿条件属性分配不同的权重, 显然如何确定属性的权重显得尤为重要.此外, 由于实际稀土萃取生产每隔2分钟就要检索案例并产生新的案例, 如果没有合适的案例更新及退出机制, 则案例库中的案例会急剧膨胀, 导致案例检索的精度下降, 因此本文主要针对这两方面问题进行设计.

    1) 案例构造

    通过对稀土萃取生产过程工艺机理和过程工况数据的分析表明, 稀土料液配分$ {Y_A} $、料液浓度$ {C_F} $、萃取剂浓度$ {C_S} $、洗涤剂浓度$ {C_W} $、料液处理量$ {G} $、产品纯度指标$ {P_{B, 1}} $及$ {P_{A, n+m}} $等入料参数和工艺指标参数直接影响到CePr/Nd萃取过程萃取量和洗涤量的设定值.每条案例都由工况描述和案例解构成, 可以用式(1)表示:

    $$ \begin{equation} {\rm{Case}}(X_k, V_k) = {\rm{Case}}((x_{1, k}, \cdots, x_{7, k}), (v_{1, k}, v_{2, k})) \end{equation} $$ (1)

    其中, $ {X_k = (x_{1, k}, \cdots, x_{7, k})} $表示第$ {k} $条工况描述, $ {x_{1, k}, \cdots, x_{7, k}} $分别表示稀土料液配比$ {R_F} $、料液浓度$ {C_F} $、萃取剂浓度$ {C_S} $、洗涤剂浓度$ {C_W} $、料液处理量$ {G} $、产品纯度指标$ {P_{B, 1}} $及$ {P_{A, n+m}} $等; $ {V_k = (v_{1, k}, v_{2, k})} $表示第$ {k} $条工况描述所对应的解答, $ {v_{1, k}, v_{2, k}} $分别表示萃取量预设定值$ {V_S} $和洗涤量预设定值$ {V_W} $.

    2) 案例检索、匹配与重用

    CBR作为一种重要的问题求解范式, 基于"相似问题具有相似解"的认知假设, 通过检索存储于案例库中的相似案例来解决当前的新问题[17].案例检索策略可采用最近邻检索策略, 通过计算目标案例与源案例的特征向量之间的距离衡量二者的相似性, 距离越小说明两案例间的相似度越高. KNN (K-nearest neighbor)检索策略[18]是最近邻检索策略的扩展, 当$ {K = 1} $时就是最近邻策略; 当$ {K>1} $时, 需将历史案例库中的$ {p} $个案例与目标案例的相似度$ {S(X, X_1), \cdots, S(X, X_P)} $从大到小排列, 依次选出相似度最大的$ {K} $个近邻, 并从这$ {K} $个近邻的结果中取出最多的某类结果作为新问题的建议解答.该策略简单有效, 计算方便, 但这种策略对于噪声或者不相关的数据比较敏感, 其解决通常依赖于给案例属性分配不同的权重.因此在利用KNN计算案例之间的相似度之前, 须对工况描述属性的权重进行合理的分配[19].

    根据工况描述$ {X_k = (x_{1, k}, \cdots, x_{7, k})} $进行案例检索和匹配, 由入矿条件对稀土萃取过程的影响程度, 设定权重为$ {w_i\; (i = 1, \cdots, 7)} $.在检索过程中, 由入矿条件和工艺指标计算所有案例与该工况输入条件的相似度, 检索出满足设定阈值的所有案例.

    假设案例库中共有$ {n} $条案例, 其工况描述为$ {X_k = (x_{1, k}, \cdots, x_{7, k})} $, 对应案例解为$ {V_k = (v_{1, k}, v_{2, k})} $.设目标工况为$ {X = (x_{1}, \cdots, x_{7})} $, 其对应的案例解为$ {V = (v_{1}, v_{2})} $.定义单个工况输入条件$ {X_k} $和$ {X} $间的相似度函数为:

    $$ \begin{align} S(X, X_k) = &1-\left(\sum\limits_{i = 1}^m {{w_i}({x_i-x_{i, k})^2}}\right)^{\frac{1}{2}}, \\& k = 1, 2, \cdots, n \end{align} $$ (2)

    其中, $ {w_i} $是第$ {i} $个特征属性的权重, 满足$ {\sum\nolimits_{i = 1}^m{w_i = 1}} $.将案例库中与新工况输入条件的相似度$ {S(X, X_k)} $大于或等于阈值的所有历史案例都检测出来作为候选案例.

    假设系统检索出大于相似度阈值的$ {l\; (l<n)} $个候选案例, 其相似度值分别为$ {s_1, \cdots, s_l} $, 对应的案例解为$ {V_k = (v_{1, k}, v_{2, k})\; (k = 1, \cdots, l)} $, 则当前新工况的案例解按下式计算:

    $$ \begin{equation} \left\{ \begin{array}{l} v_1 = \dfrac{\sum\limits_{k = 1}^l{s_k}\times{v_{1, k}}}{\sum\limits_{k = 1}^l{s_k}}\\ v_2 = \dfrac{\sum\limits_{k = 1}^l{s_k}\times{v_{2, k}}}{\sum\limits_{k = 1}^l{s_k}} \end{array} \right. \end{equation} $$ (3)

    属性权重优化的方法较多[20-21], 由于萃取工艺生产过程的复杂多样性, 本文采用内省学习方法调整式(2)的属性权重, 根据两个案例相同属性的描述值, 需将案例属性分为匹配属性和不匹配属性, 给出其定义[22].

    定义1.  设历史案例库中的两案例$ C_I = (X_I;Y_I) $和$ {C_{II} = (X_{II};Y_{II})} $, 其中, $ X_I = (x_{1, I}, x_{2, I}, \cdots, $ $ x_{m, I}) $, $ X_{II} = (x_{1, II}, x_{2, II}, \cdots, $ $ x_{m, II}) $, 若满足$ 0\leq $ $ \mid\! x_{i, I}-x_{i, II}\!\mid\leq \xi_{\rm match} $, 则称案例$ {C_I} $与案例$ {C_{II}} $的第$ {i} $个属性为匹配属性, 否则, 称为不匹配属性.其中, $ {\xi_{\rm match}} $为判断属性是否匹配的阈值, 表示两个案例同一属性的接近程度.

    基于内省学习原理, 采用如下权重学习策略规则.

    规则1.  存在一条案例分类成功时, 匹配属性的权重增加, 同时不匹配属性的权重减少:当检索到的案例$ A $具有和目标案例$ C $相同的类别时(分类成功), 增大案例$ A $中与案例$ C $匹配的属性的权重, 同时减小案例$ A $中与案例$ C $不匹配的属性的权重.如此, $ A $与$ C $的相似度将增大, 保证在求解案例$ C $时, 相似案例$ A $更有可能被检索到.

    规则2.  存在一条案例分类失败时, 匹配属性的权重减少, 不匹配属性的权重增加:当检索到的案例$ B $具有和目标案例$ C $不同的类别时(分类失败), 减小案例$ B $中与案例$ C $匹配的属性的权重, 同时增大案例$ B $中与案例$ C $不匹配的属性的权重.如此, $ B $与$ C $的距离增大($ B $被"推离" $ C $), 保证在求解案例$ C $时, 案例$ B $不易被检索到.

    根据上述定义和规则, 考虑稀土入矿条件和参数属性对案例属性权重的影响, 对属性权重调整如下.

    权重增加:

    $$ \begin{equation} w_i(t+1) = w_i(t)+\frac{\triangle_i}{m}, \quad w_i(t+1)\geq0 \end{equation} $$ (4)

    权重减少:

    $$ \begin{equation} w_i(t+1) = w_i(t)-\frac{\triangle_i}{m}, \quad w_i(t+1)\geq0 \end{equation} $$ (5)

    其中, $ {w_i(t)} $是第$ {i} $个属性第$ {t} $次迭代的权重; $ {w_i(t+1)} $是第$ {i} $个属性第$ {t+1} $次迭代的权重; $ {\triangle_i/m} $决定权重的变化量, $ {m} $是案例属性的个数, 对于不同的分类问题, 如果属性个数不同, 即使$ {\triangle_i} $相同,权重的变化量也不相同.但是对于同一个案例库而言, $ {m} $是固定不变的, 因此权重的变化主要与$ {\triangle_i} $有关, 为实现属性权重有效调整, 且保证属性权重调整值不会超过属性权重之和, $ {\triangle_i} $在$ {(0, m)} $区间内取值, 并且需满足式(4)和式(5)更新后的权重不小于零.

    当一条案例所有属性的权重都被调整之后, 为保证所有属性的权重之和为1, 需进行归一化操作:

    $$ \begin{equation} w_i(t+1)' = \frac{w_i(t+1)}{\sum\limits_{i = 1}^m w_i(t+1)} \end{equation} $$ (6)

    其中, $ {w_i(t+1)'} $是第$ {i} $个属性第$ {t+1} $次迭代的最终权重.

    3) 案例修正、存储和删除

    经案例推理得到的案例结果与新工况不相匹配时, 则需应用其他领域相关知识对案例结果进行人工修正, 或更改参数和检索规则, 使得修正后的案例解(各回路流量设定值)与新工况相匹配.如果根据当前的工况描述没有被正确地分类成案例, 则说明当前案例库中没有该工况的案例, 则需要把当前的工况及案例加入案例库.如果根据当前的工况描述被正确地分类了, 但是其相似度太小, 则说明分类信息不完善, 则需要重新更新该类工况的案例, 并使得该类工况案例的相似度小于$ {\sigma[0, 1]}$并保存.

    很明显, 案例推理的案例库随生产工况变化不断增加, 不但增加了检索的难度, 而且降低检索的精确度.所以, 伴随着检索案例的增加, 同时应该删除一些已经不适合新工况的案例.

    假定每一个新的案例被保存时都被赋予一个初始遗忘因子, 根据内省学习原理, 其遗忘值的更新如下[23]:

    $$ \begin{equation} F_i(t+1) = F_i(t)+\beta\times{r_i}\times{F_i(0)} \end{equation} $$ (7)

    其中, $ {\beta} $为遗忘增强因子; $ {r_i} $为奖励函数, 当检索案例与新案例一致时为$ -1, $否则为1. $ {F_i(0)} $为案例的初始遗忘值, $ {F_i(t)} $为案例当前迭代的遗忘值, $ {F_i(t+1)} $为案例当前迭代后的遗忘值.很明显, 当检索案例每次都能成功地匹配当前工况时, 则遗忘值不断减少.当检索案例每次不能成功地匹配当前工况时, 则遗忘值不断增加, 当$ {F_i(t+1) > \varsigma} $时, 该案例需要删除.

    通常情况下, 采用基于案例推理的预设定模型可得到一个相对合理的萃取剂和洗涤剂流量预设定值, 但该模型没有考虑到实时外界干扰导致的工况变化.当稀土萃取生产过程遇扰动使工况变化时, 难以保证生产指标满足工艺要求.实际铈镨/钕萃取工业生产过程中, 工作人员周期性地根据溶液颜色或组分含量化验值等调整萃取剂和洗涤剂流量差值, 使得生产指标满足工艺要求.本质上, 领域专家或工作人员是根据稀土元素组分含量变化差和变化率调节药剂量流量, 因此, 基于铈镨/钕溶液颜色能够表征其组分含量, 本文采用基于稀土溶液颜色的在线补偿模型, 该模型主要包括基于稀土溶液颜色的组分含量软测量模型和基于模糊推理的智能补偿模型.

    稀土元素独特的电子层结构可使稀土离子具有丰富的发射光谱. CePr/Nd萃取生产工艺的易萃组分和难萃组分分别显现紫色和绿色颜色特征, 且颜色深浅与稀土离子溶度有关.由于难萃稀土组分含量和易萃稀土组分含量之和为1, 所以一般情况下只需要知道其中一个即可.将采集的稀土溶液颜色图像在HSI颜色空间下提取出H、S、I颜色特征分量.颜色H分量和S分量与组分含量的相关性分析如图 3所示, 图 3(a)为二维图, 图 3(b)为三维图.从图 3中可知, 溶液图像颜色与组分含量关系具有一定的规律, 且为非线性关系, 可表示如下[24]:

    $$ \begin{equation} y = f(C_h, C_s) \end{equation} $$ (8)
    图 3  颜色分量H、S特征值和组分含量相关性分析
    Fig. 3  Correlation between color feature and component content

    其中, $ {y} $表示组分含量, $ C_h $、$ C_s $分别表示HSI颜色空间的H和S分量.该非线性关系可采用非线性回归、神经网络等建模.由于最小二乘支持向量机(LS-SVM) 在求解二次规划问题时遇到的稀疏性及运算复杂问题时简化了采用二次规划求解优化问题的方式, 对小样本建模效果较好[25].因此, 式(8)采用LS-SVM建模.

    LS-SVM算法的最小化目标函数为[26]:

    $$ \begin{align} & \min\;J = \frac{1}{2}w^{ {\rm{T}}}w+\frac{1}{2}\gamma\sum\limits_{i = 1}^l{e_i^2} \\ & {\rm{s.t.}}\;y_i = w^ {{\rm{T}}}\cdot\varphi(x_i)+b+e_i, \;i = 1, \cdots, l \end{align} $$ (9)

    根据式(9), 构造其拉格朗日函数为:

    $$ \begin{equation} L = J-\sum\limits_{i = 1}^l \alpha_i[w^ {{\rm{T}}}\cdot\varphi(x_i)+b+e_i-y_i] \end{equation} $$ (10)

    其中, $ {\alpha_i} $为拉格朗日乘子, 根据KTT条件得:

    $$ \begin{equation} \begin{cases} \dfrac{\partial L}{\partial w} = 0 \to {w} = \sum\limits_{i = 1}^l \alpha_i\varphi(x_i)\\ \dfrac{\partial L}{\partial b} = 0 \to \sum\limits_{i = 1}^l \alpha_i = 0 \\ \dfrac{\partial L}{\partial e_i} = 0 \to \alpha_i = \gamma_{e_i}\\ \dfrac{\partial L}{\partial a_i} = 0 \to w^ {\rm{T}}\cdot\varphi(x_i)+b+e_i-y_i = 0\\ \end{cases} \end{equation} $$ (11)

    对于$ {i = 1, \cdots, l} $, 消去$ {w} $和$ {e} $, 得到如下线性方程:

    $$ \begin{equation} \left[ {\begin{array}{*{20}{c}} {0}&{{\bf{\mathord{\buildrel{\lower3pt\hbox{$ \rightharpoonup$}} \over 1} }}^ {\rm{T}}} \\ {\bf{\mathord{\buildrel{\lower3pt\hbox{$ \rightharpoonup$}} \over 1} }} &{ZZ^ {\rm{T}}+{\gamma^{-1}}I} \\ \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {{b}}\\ {{ \alpha}}\\ \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {0}\\ {y}\\ \end{array}} \right] \end{equation} $$ (12)

    其中, $ {\bf{\mathord{\buildrel{\lower3pt\hbox{$\scriptscriptstyle\rightharpoonup$}} \over 1} }} = [1\cdots1]^ {\rm{T}} $, $ {a = [\alpha_1, \cdots, \alpha_l]^ {\rm{T}}} $, $ {y = [y_1, \cdots, y_l]^ {\rm{T}}} $, $ {Z = [\varphi(x_1), \cdots, \varphi(x_l)]^ {\rm{T}}} $, 令 $ {K(x, x_i) = \varphi(x)^ {\rm{T}}\varphi(x_i)} $, 根据核函数方法, 引入$ {n\times n} $的高斯径向基函数为核函数:

    $$ \begin{equation} K(x_i, x_k) = {\rm{exp}}(-\frac{\|x_i-x_k\|}{2\sigma^2}) \end{equation} $$ (13)

    则采用支持向量机进行预测的输出为:

    $$ \begin{equation} y(x) = \sum\limits_{i = 1}^l a_iK(x, x_i)+b \end{equation} $$ (14)

    LS-SVM模型的参数$ {a, b} $由式(12)计算, 高斯核函数宽度$ {\sigma} $和误差惩罚系数$ {\gamma} $等由粒子群优化算法计算[27]:

    $$ \begin{align} v_{i, d}(t+1) = &wv_{i, d}(t)+c_1r_1(P_{i, d}-x_{i, d}(t))+ \\& c_2r_2(P_{g, d}-x_{i, d}(t)) \end{align} $$ (15)
    $$ \begin{equation} x_{i, d}(t+1) = x_{i, d}(t)+v_{i, d}(t+1) \end{equation} $$ (16)

    其中, $ {x_{i, d}} $, $ {v_{i, d}} $, $ {P_{i, d}} $, $ {P_{g, d}} $分别表示粒子当前位置、当前速度、当前最佳位置和粒子的最佳位置. $ {\beta} $ 称为内部系数, 为[0, 1]间常量; $ {c_1} $和$ {c_2} $为学习速率; $ {r_1} $和$ {r_2} $为[0, 1]随机系数. $ v_{id}\in[-v_{ \rm{max}}, v_{ \rm{max}}], {v_{ \rm{max}}}$ 为最大速率.

    稀土萃取生产过程中, 工作人员根据组分含量的变化量及变化率对当前的流量预设定值进行补偿, 使稀土元素组分含量满足生产指标.工作人员调节药剂量完全依赖于其经验知识, 且调节量是估计值.相比其他智能模型, 模糊推理是以模糊系统理论为基础、专家经验和领域知识为决策依据, 把定性模型转化成定量模型, 且不需要依据萃取机理建立关系模型.因此, 本节采用模糊推理建立萃取量和洗涤量补偿模型.该模型根据稀土元素组分含量软测量值, 补偿萃取量和洗涤量设定值.

    通过第3.1节得到的组分含量 ${Y_{B}(k-1)}$ , ${Y_{A}(k-1)} $与目标组分含量$ {P_{B, 1}} $, $ {P_{A, n+m}} $比较, 得到其偏差$ e_s = Y_A(k-1)-P_{A, n+m} $, $ e_w = Y_B(k-1)-P_{B, 1} $及其偏差的变化$ \triangle e_s = Y_A(k-1)-Y_A(k-2) $, $ \triangle e_w = Y_B(k-1)-Y_B(k-2) $作为模糊系统的输入量, 萃取剂和洗涤剂的补偿增量$ {\triangle v_s} $和$ {\triangle v_w} $作为模糊推理系统的输出量.

    基于模糊推理的智能补偿器如图 4所示, 主要包括模糊化/解模糊部分、模糊推理、知识库等; 知识库包括输入输出变量的比例(尺度)因子、隶属度函数和模糊语义控制规则.通常, 萃取剂流量模糊控制器的输入与输出之间的关系可表示为[28-29]:

    $$ \begin{equation} \triangle v_s(k) = f(e_s(k), \triangle e_s(k)|\psi) \end{equation} $$ (17)
    图 4  基于模糊的补偿器
    Fig. 4  Fuzzy based compensator

    其中, $ {f(\cdot|\psi)} $表示参数为$ {\psi} $ (如隶属度函数和模糊语义控制规则等)的模糊函数.假定萃取剂补偿量为$ {\triangle v_s(k)} $, 洗涤剂补偿量为$ {\triangle v_w(k)} $, 水相出口分数为$ {f_b'} $, 根据萃取机理模型[30], 可推导出下列公式:

    $$ \begin{equation} \triangle v_w(k) = \triangle v_s(k)+f_b' \end{equation} $$ (18)

    一般情况只需要模糊推理萃取剂或洗涤剂补偿量即可.

    式(17)所示的模糊补偿器需要确定模糊逻辑控制器的参数向量$ {\psi} $, 该参数向量包括模糊化方法、模糊推理方法、解模糊方法和知识库等.由于参数向量$ {\psi} $里面的参数相互作用影响输出控制量, 所以实际工业过程的模糊控制器要满足控制性能并非容易得到, 需要大量的时间和精力调整设计参数.基于模糊系统理论, 参数$ {\psi} $可写成[31]:

    $$ \begin{equation} \psi = [k_h, k_g, k_u;a_h, b_h, c_h;a_g, b_g, c_g;a_u, b_u, c_u] \end{equation} $$ (19)

    其中, $ {k} $表示模糊控制的比例(尺度)因子, $ {a} $、$ {b} $和$ {c} $表示三角形隶属度函数的参数集, 下标$ {h} $、$ {g} $和$ {u} $分别表示组分含量误差、组分含量变化差和萃取量.

    为实现对萃取量和洗涤量预设定值的补偿, 组分含量误差和组分含量变化差须模糊化成语义输入量, 经模糊推理得到的语义输出量须解模糊成实际补偿萃取量和洗涤量.假定组分含量误差的论域$ {e_s} $为$ [-e_{ \rm{smax}}, e_{ \rm{smax}}]$, 组分含量变化差$ {\triangle e_s} $的论域为 $ [-\triangle e_{ \rm{smax}}, \triangle e_{ \rm{smax}}]$, 补偿萃取量$ {\triangle v_s} $的论域为$ [-\triangle e_{ \rm{wmax}}, \triangle e_{ \rm{wmax}}]$. 语义变量 $ {E_s} $的论域为$ {-n, \cdots, 0, \cdots, n} $;语义变量$ {\triangle E_s} $的论域为$ {-m, \cdots, 0, \cdots, m-1, m} $; 语义变量$ {\triangle V_s} $ 的论域为$ {-l, \cdots, 0, \cdots, l-1, l} $ .因此, 比例因子表示为:

    $$ \begin{equation} k_g = \frac{n}{e_{ \rm{smax}}}, k_h = \frac{m}{\triangle e_{ \rm{smax}}}, k_u = \frac{\triangle v_{ \rm{smax}}}{l} \end{equation} $$ (20)

    $ {E_s} $和$ {\triangle E_s} $的论域定义七个模糊集为: "NB", "NM", "NS", "ZE", "PS", "PM", "PB". $ {\triangle V_s} $的论域定义七个模糊集为: "NL", "NB", "NM", "NS", "ZE ", "PS ", "PM ", "PB ", "PL".正如图 5(a)所示, 三角形隶属度函数主要有三个顶点参数$ {a} $、$ {b} $、$ {c} $确定.当两个相邻的隶属度函数确定后, 其交点$ {\lambda} $即确定.在一般情况下, 交点$ {\lambda} $ 应该位于两个隶属度函数中心位置, 能够对输入量快速响应及输入量小范围内变化精确响应.组分含量误差、组分含量变化差和萃取剂补偿量的隶属度函数如图 5(b)5(c)(d)所示.

    图 5  组分含量误差、组分含量变化差和萃取剂补偿量的隶属度函数
    Fig. 5  Membership functions of component content, change error of component content and compensator of the flow rate of the extractant

    根据专家经验和领域知识, 模糊规则表如表 1所示, 利用下面的语义控制规则构建模糊规则:

    表 1  语义规则表
    Table 1  Semantic rules table
    NB NM NS ZO PS PM PB
    \hline NB PL PL PB PM PM PS ZO
    NM PL PB PM PS PS ZO NS
    NS PB PM PS PS ZO NS NM
    ZO PM PM PS ZO NS NM NM
    PS PM PS ZO NS NS NM NB
    PM PS ZO NS NS NM NB NL
    PB ZO NS NM NM NB NL NL
    下载: 导出CSV 
    | 显示表格

    $ {R_i} $:如果$ {E_s} $为$ {\Gamma_i(E_s)} $和$ {\triangle E_s} $为$ {\Gamma_i(\triangle E_s)} $则$ {\triangle V_s} $为 $ {\Gamma_i(\triangle V_s)} $.

    其中, $ {R_i} $表示$ {{\rm if}\cdots {\rm then}\cdots} $规则, $ \Gamma_i(\cdot)\in $ $ \{{\rm NL, NB, NM, NS, ZO, PS, PM, PB, PL}\} $. $ {NL} $、$ {NB} $等代码为模糊简写符号. $ {E_s} $、$ {\triangle E_s} $ 和$ {\triangle V_s} $ 分别表示组分含量误差、组分含量变化差和萃取量.

    由于模糊推理的输出结果为一模糊集, 反映了模糊语义输出变量论域各元素隶属度大小的组合, 为了能够得到萃取剂补偿量的精确答案, 使用重心法解模糊:

    $$ \begin{equation} \triangle v_s = \frac{\sum\limits_{i = 1}^n(\Gamma_i(\triangle v_s)\cdot\mu_v(\triangle v_s))}{\sum\limits_{i = 1}^n\mu_v(\triangle v_s)} \end{equation} $$ (21)

    其中, $ {n} $指输出论域元素的总数, $ {\Gamma_i(\triangle v_s)} $表示输出论域的元素, 权系数$ {\mu_v(\triangle v_s)} $表示相应的隶属度函数.需要注意的是, $ {\triangle v_s} $需要乘以比例因子$ {k_u} $才能施加于预设定值.

    为了验证基于LS-SVM的组分含量软测量模型, 以某公司CePr/Nd萃取分离生产过程为对象, 在CePr/Nd萃取生产现场采集350份混合溶液样品, 287份作为训练样本, 63 份作为测试样本.溶液样本经机器视觉和化验两个基本预处理过程, 其提取出混合溶液图像的H、S颜色特征分量一阶矩作为软测量模型的输入变量, 化验后的CePr/Nd组分含量作为模型的输出变量.为了消除变量间由于数量级差异带来的影响, 将输入/输出样本组成的数据集进行归一化处理.在案例推理模型中, 案例阈值设为0.93, $ {\triangle_i} $设为2, 遗忘因子$ {\beta} $设为0.1, $ {\varsigma} $设为2.其案例推理权重属性通过学习策略计算分别为{0.12, 0.11, 0.21, 0.21, 0.11, 0.12, 0.12}. 在LS-SVM 模型中, 通过粒子群优化算法确定核宽度$ {\sigma} $ 为10.37, 惩罚系数(核参数) $ {\gamma} $为0.53, 287组样本数据集训练LS-SVM 模型, 利用63 组测试样本在已经训练好的模型上进行测试, 并与人工化验得到的63组样本所对应的组分含量化验值进行对比分析, 结果如图 6所示.结果表明该软测量模型基本满足要求, 为了说明该方法的有效性, 与RBF、SVM方法比较, 采用均方根误差(反映测量数据偏离真实值的程度)和均方差(反映一个数据集的离散程度)作为评价指标, 计算公式分别为: $ {{\rm RMSE} = \sqrt{\frac{1}{n}\sum\nolimits_{i = 1}^n(\hat{y_i}-y_i)^2}} $ 和 $ {{\rm MSE} = \sqrt{\frac{1}{n}\sum\nolimits_{i = 1}^n(y_i-\bar{y})^2}} $. 结果如表 2所示.本文采用的LS-SVM相比较于RBF和SVM方法效果更好.

    图 6  基于溶液颜色的组分含量软测量与化验值比较
    Fig. 6  Component content by soft sensing based on color feature of rare earth solution and lab test
    表 2  不同方法建模结果的性能比较
    Table 2  Performance comparison of different models
    Nd (%) CePr (%)
    $ {\rm Method} $ $ {\rm RMSE} $ $ {\rm MSE} $ $ {\rm RMSE} $ $ {\rm MSE} $
    $ {\rm RBF} $ 0.5829 0.623 0.6243 0.6109
    $ {\rm SVM} $ 0.6021 0.528 0.5820 0.652
    $ {\rm LSSVM} $ 0.5012 0.4891 0.5391 0.509
    下载: 导出CSV 
    | 显示表格

    模糊药剂量补偿模块主要确定特征参数$ {\psi} $的值, 即$ \psi = [k_h, k_g, k_u; a_h, b_h, c_h; a_g, b_g, c_g; $ $ a_u, b_u, c_u] $.根据工艺试验, Nd组分含量误差的论域为$ [-3, 3] $, 其预定值为0, 因此$ \triangle h $的论域为 $ [-3, 3] $, 则$ {k_h = 3/6 = 0.5} $. Nd组分含量变化率的论域为$ [-0.18, 0.18], $其预定值为0, 因此$ {\triangle g} $的论域为$ [-0.18, 0.18] $, 则$ {k_g = 0.18/6 = 0.03} $.萃取量范围为 [3ml/min, 5ml/min] 且其预定值为4ml/min, 则$ {\triangle u} $的论域[$ -\rm{ml}/min, 1ml/min]$, 则$ {k_u = 8/1 = 8} $. $ {\triangle h} $, $ {\triangle g} $和$ {\triangle u} $隶属度函数的特征参数如表 3所示.工业现场工艺实验可知水出口分数$ {f_b' = 0.9298} $.

    表 3  使用特征参数的隶属度函数
    Table 3  Membership functions of the fuzzy logic controller using characteristic parameters
    Nd组分含量误差 Nd组分含量变化率 萃取补偿量
    NO. $ {a_h} $ $ {b_h} $ $ {c_h} $ $ {a_g} $ $ {b_g} $ $ {c_g} $ $ {a_u} $ $ {b_u} $ $ {c_u} $
    1 -8 -6 -4 -8 -6 -4 -10 -8 -6
    2 -6 -4 -1.5 -6 -4 -1.5 -8 -6 -4
    3 -4 -1.5 0 -4 -1.5 0 -6 -4 -1.5
    4 -1.5 0 1.5 -1.5 0 1.5 -1.5 0 1.5
    5 0 1.5 4 0 1.5 4 -4 -1.5 0
    6 1.5 4 6 1.5 4 6 0 1.5 4
    7 4 6 8 4 6 8 1, 5 4 6
    8 $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ 4 6 8
    9 $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ $ {\times} $ 6 8 10
    下载: 导出CSV 
    | 显示表格

    采用人工方法、单独采用案例推理预设定模型与本文方法的控制效果相比较如何?为此, 结合图 1所示的铈镨/钕稀土萃取过程, 将图 2所提出的本文方法与人工设定方法、基于案例推理的单模型方法进行对比试验.三组萃取工艺生产流程来自于同一料液, 一组生产流程采用本文的智能设定方法; 另一组萃取生产流程采用人工设定方法; 第三组采用基于案例推理的单模型设定方法.三组萃取生产过程的料液相同, 萃取生产过程条件相同.每隔2小时记录萃取量和洗涤量、图像分析的组分含量及其对应的离线化验的组分含量, 各得到60组运行数据, 统计结果如图 7所示.由于图像采集安装在萃取段第20级, 在该级萃取槽中镨的组分含量应该在[26, 30]范围内, CePr的组分含量在[70, 74]范围内.

    图 7  人工方法、单模型方法与本文方法的组分含量
    Fig. 7  Component content obtained by manual method、CBR based method and the proposed method

    图 7方法1采用本文提出智能设定方法对其设定, 方法2表示人工方法设定, 方法3表示仅采用基于案例推理的单模型设定方法进行设定.图 8为三种方法的药剂量变化情况.图 7图 8反映了组分含量的变化及相应的萃取剂和洗涤剂流量变化情况.在61个样本点里, 当组分含量变化时, 其对应的萃取剂和洗涤剂流量也就相应地变化, 并反映到了下一个组分含量上面.由图 7可看出, 方法2的人工调节随意性较大, 调节时喜欢采用一步到位法, 导致组分含量波动较大.方法3的组分含量波动较大, 并且不满足要求的数据点较多, 其原因在于通过案例推理模型计算萃取量, 但没有利用反馈补偿机制.采用方法1的组分含量波动较小, 满足要求的数据点较多.通过计算, 单独采用预设定模型时精矿品位合格率为73.78%, 采用人工方法的合格率为75.4%, 而采用本文方法时的合格率为83.5%.这说明经过智能动态补偿模型的补偿作用, 所提方法克服了工况变化带来的扰动, 弥补了人工方法、基于案例推理的单模型设定方法的不足.

    图 8  人工方法、案例推理方法与本文方法的萃取量
    Fig. 8  Flow rate of the extractant and the detergent consumed by manual method、CBR based method and the proposed method

    在稀土萃取工业生产领域, 除了生产指标, 萃取量和洗涤量消耗也是萃取生产厂家最为关注的指标之一, 与工厂绩效密切相关.为评估本文方法的有效性, 比较人工设定、单模型方法及智能优化设定的生产指标和萃取量和洗涤量消耗量等.图 8显示本文方法的萃取量和洗涤量波动比另外两种方法要小, 说明药剂量控制相对比较平缓; 而且, 如表 4所示, 本文方法的萃取量和洗涤量消耗比另外两种方法要少.结果表明, 比起另外两种方法, 本文方法的三个模型为紧密联系的有机整体, 能够实现更好的萃取剂和洗涤剂流量控制, 并达到了降低药剂量消耗的目的.

    表 4  萃取量和洗涤量消耗统计表(升)
    Table 4  Sum of the extractant and detergent comsumend by three methods (L)
    人工方法 案例推理方法 本文方法
    洗涤量 82.67 82.29 81.54
    萃取量 28.63 28.26 27.54
    下载: 导出CSV 
    | 显示表格

    萃取生产过程具有机理复杂、非线性、时滞大等特点, 难以建立精确的萃取量和洗涤量优化控制模型, 当前稀土生产企业采用主观性强的人工设定方式, 易造成工艺指标波动较大.本文根据稀土溶液颜色与组分含量密切相关的特点, 提出了一种基于案例推理静态设定和机器视觉动态补偿的萃取量和洗涤量优化控制方法.首先根据案例推理确定萃取剂/洗涤剂预设定值, 然后根据软测量组分含量和目标组分含量的差值采用模糊推理方法智能补偿萃取量预设定值, 实现了萃取生产过程萃取溶液组分含量的优化控制, 试验表明该方法能够实现较好的生产指标及经济效益.


  • 收稿日期 2019-10-30 录用日期 2020-04-16 Manuscript received October 30, 2019; accepted April 16, 2020 国家自然科学基金 (61603310), 国家社会科学基金 (19BYY171), 重庆市自然科学基金 (cstc2019jcyj-msxm2550), 模式识别国家重点实验室开放课题 (201900010), 中央高校基本科研业务费(XDJK2018B020), 重庆市教育委员会科学技术研究计划青年项目(KJQN201801901, KJQN201801902) 资助 Supported by National Natural Science Foundation of China (61603310), National Social Science Foundation of China (19BYY171), Natural Science Foundation of Chongqing (cstc2019jcyj-msxm2550), Open Projects Program of National Laboratory of Pattern Recognition (201900010), Fundamental Research Funds for the Central Universities of China (XDJK2018B020), and Youth Project of Science and Technology Research Program of Chongqing Education Commission (KJQN201801901, KJQN201801902) 本文责任编委 金连文 Recommended by Associate Editor JIN Lian-Wen 1. 西南大学计算机与信息科学学院 重庆 400715 2. 重庆工程学院计算机与物联网学院 重庆 400056 3. 贵州工程应用技术学
  • 院彝学研究院 毕节 551700 1. College of Computer and Information Science, Southwest University, Chongqing 400715 2. College of Computer and Internet of Things, Chongqing Institute of Engineering College, Chongqing 400056 3. Institute of Yi Studies, Guizhou University of Engineering Science, Bijie 551700
  • 图  1  彝文残卷

    Fig.  1  The incomplete literature of the ancient Yi

    图  2  本文双判别器生成式对抗网络结构

    Fig.  2  Generative adversarial networks with double discriminator in the paper

    图  3  文献[25]中D2GN结构

    Fig.  3  The structure of D2GN in [25]

    图  4  生成式对抗网络的结构图

    Fig.  4  The strutrue of generative adversarial networks

    图  5  古彝文字符判别器模型详细结构

    Fig.  5  Detailed structure of the ancient Yi character discriminator model

    图  6  古彝文字符生成器模型详细结构

    Fig.  6  Detailed structure of the ancient Yi character generator model

    图  7  古彝文字符筛选判别器模型

    Fig.  7  Selecting discriminator model for ancient Yi character

    图  8  原始样本

    Fig.  8  The original sample

    图  9  待修复古彝文

    Fig.  9  Ancient Yi character need to be restored

    图  10  通过生成器模型输出图像$G(\boldsymbol z)$

    Fig.  10  Output images $G(\boldsymbol z)$ from the generator model

    图  11  原始样本

    Fig.  11  The original sample

    图  12  古彝文硬笔(上)和软笔(下)

    Fig.  12  Ancient Yi hard pen (upper) and soft pen (down)

    图  13  待修复古彝文

    Fig.  13  Ancient Yi character need to be restored

    图  14  古彝文手写数据集样例

    Fig.  14  The handwritten sample of ancient Yi

    图  15  学习率0.2, 0.02, 0.002的损失值变化曲线

    Fig.  15  The loss variation of the learning rate involving 0.2, 0.02 and 0.002

    图  16  学习率为0.0002, 0.001, 0.002的损失值变化曲线

    Fig.  16  The loss variation of the learning rate involving 0.0002, 0.001 and 0.002

    图  17  不同训练次数下生成器生成图像

    Fig.  17  The generator generates image under different training times

    图  18  10000次训练后损失值变化曲线

    Fig.  18  Loss curve after 10000 training

    图  19  生成器输出图像

    Fig.  19  Output images by generator

    图  20  筛选判别器的训练过程中损失值变化曲线

    Fig.  20  The loss curve in process of the training of the selecting discriminator

    图  21  训练得到$\boldsymbol z',$ 然后输入$\boldsymbol z'$到生成器得到的图像

    Fig.  21  After trainning, $\boldsymbol z'$ is generated, and then input $\boldsymbol z'$ to the generator to get the image

    图  22  修复后的图像

    Fig.  22  The restored image

    图  23  部分古彝文修复结果

    Fig.  23  The repair effect of some ancient Yi

    图  24  多形状残缺修复结果

    Fig.  24  The repair effect of ancient Yi characterof multiple shape occlusion

    图  25  彝文古籍文献中残缺字符修复效果

    Fig.  25  The repair effect of incomplete characters in ancient Yi literature

    图  26  古彝文残缺字符修复失败效果

    Fig.  26  The failed repair effect of ancient Yi incomplete characters

    表  1  判别器模型参数表

    Table  1  Parameter table of the discriminator model

    层信息卷积核
    个数
    卷积核
    大小
    步长特征图
    大小
    参数个数
    C1层 (卷积层)1285×5232×323328
    C2层 (卷积层)2565×5216×16819456
    C3层 (卷积层)5125×528×83277312
    C4层 (卷积层)10245×524×413108224
    OUTPUT层 (输出层)1×116385
    下载: 导出CSV

    表  2  生成器模型参数表

    Table  2  Parameter table of the generator model

    层信息卷积核
    个数
    卷积核
    大小
    步长特征图
    大小
    参数个数
    C1层 (卷积层)1654784
    C2层 (卷积层)5125×528×813107712
    C3层 (卷积层)2565×5216×163277056
    C4层 (卷积层)1285×5232×32819328
    OUTPUT层 (输出层)15×5264×643201
    下载: 导出CSV

    表  3  古彝文字符修复比例

    Table  3  Restoration proportion of ancient Yi characters

    个数占比 (%)时间消耗 (ms)
    完全修复5235252
    部分修复2462561
    未完成修复2312356
    合计1000100
    下载: 导出CSV
  • [1] 王正贤. 彝文金石述略. 贵州民族研究, 2002, 20(2): 156-163 doi: 10.3969/j.issn.1002-6959.2002.02.028

    Wang Zhen-Xian. Brief introduction to Yiwen and gold stone language. Guizhou Ethnic Studies, 2002, 20(2): 156-163 doi: 10.3969/j.issn.1002-6959.2002.02.028
    [2] 吴勰, 禄玉萍, 王明贵. 论贵州古彝文编码字符集构建. 中文信息学报, 2014, 28(4): 153-158 doi: 10.3969/j.issn.1003-0077.2014.04.021

    Wu Xie, Lu Yu-Ping, Wang Ming-Gui. On the construction of ancient Yi coded character sets in Guizhou. Journal of Chinese Information, 2014, 28(4): 153-158 doi: 10.3969/j.issn.1003-0077.2014.04.021
    [3] 史金波. 当代少数民族文字历史文献整理与研究刍议. 青海民族研究, 2018, 29(3): 108-115

    Shi Jin-Bo. On the collation and research of historical literature of contemporary minority characters. Qinghai Nationalities Research, 2018, 29(3): 108-115
    [4] 滇川黔桂彝文协作组, 滇川黔桂彝文字典. 云南民族出版社, 2001, 1(1): 278−286

    Yunnan Sichuan Guizhou Guangxi Yi Language Cooperation Group. Yunnan Sichuan Guizhou Guangxi Yi Dictionary. Yunnan Nationalities Publishing House, 2001, 1(1): 278−286
    [5] 陈善雄, 王小龙, 韩旭等. 一种基于深度学习的古彝文识别方法, 浙江大学学报(理学版), 2019, 46(3): 261-269

    Chen Shan-Xiong, Wang Xiao-Long, Han Xu. A recognition method of ancient Yi based on deep learning. Journal of Zhejiang University (SCIENCE EDITION), 2019, 46(3): 261-269
    [6] Bouillon M, Ingold R, Liwicki M. Grayification: A meaningful grayscale conversion to improve handwritten historical documents analysis. Pattern Recognition Letters, 2018, 121(15): 46-51
    [7] Cao XY, Zhou F, Lin X, et al. Hyperspectral image classification with markov random fields and a convolutional neural network. IEEE Transactions on Image Processing, 2018, 27(5): 2354-2367 doi: 10.1109/TIP.2018.2799324
    [8] Danilla C, Persello C, Tolpekin V, Bergado J R. Classification of multitemporal sar images using convolutional neural networks and Markov random fields. In: Proceedings of the 2017 Geoscience and Remote Sensing Symposium (IGARSS), IEEE International Symposium. Fort Worth, TX, USA: IEEE, 2017. 2231−2234
    [9] Huang X, Belongie S. Arbitrary style transfer in real-time with adaptive instance normalization. In: Proceedings of the 2017 IEEE International Conference on Computer Vision (ICCV). Venice, Italy: IEEE, 2017(1): 1510−1519
    [10] Cireşan D C, Meier U, Schmidhuber J. Transfer learning for Latin and Chinese characters with deep neural networks. In: Proceedings of the 2012 International Joint Conference on Neural Networks. Brisbane, QLD, Australia: IEEE, 2012. 1−6
    [11] Goodfellow L, Pouget-Abadie J, Mirza M, Xu B, Warde-Farley D, Ozair S, Courville A, Bengio Y. Generative adversarial nets. In: Proceedings of the 27th International Conference on Neural Information Processing Systems. Kuching, Malaysia: MIT Press, 2014, 2672−2680
    [12] Chen S M, Wang W J, Xia B H, You X G, Peng Q M, Cao Z H, et al. CDE-GAN: Cooperative dual evolution-based generative adversarial network, IEEE Transactions on Evolutionary Computation, 2021, 25(5): 986−1000
    [13] Radford A, Metz L, Chintala S. Unsupervised representation learning with deep convolutional generative adversarial networks. In: Proceedings of the 4th International Conference on Learning Representations (ICLR ). San Juan, Puerto Rico: OpenReview.net, 2016.
    [14] 赵树阳, 李建武. 基于生成对抗网络的低秩图像生成方法. 自动化学报, 2018, 44(5): 829-839

    Zhao Shu-Yang, LI Jian-Wu. Generative adversarial network for generating low-rank images. ACTA AUTOMATICA SINICA, 2018, 44(5): 829-839
    [15] Wei P, Li H, Hu P. Inverse discriminative networks for handwritten signature verification. In: Proceedings of the 2019 IEEE Conference  on  Computer  Vision  and  Pattern  Recognition (CVRP). Long Beach, CA, USA: IEEE, 2019. 5764−5772
    [16] Li X W, Shen H Q, Yu M, Wei X, Han J, Zhu J L, et al. Blind image inpainting using pyramid GAN on thyroid ultrasound images. In: Proceedings of the 2019 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). San Diego, CA, USA: IEEE, 2019. 678−683
    [17] Chang  J,  Gu  Y  J,  Zhang  Y,  Wang  Y  F.  Chinese  handwriting imitation  with  hierarchical  generative  adversarial  network.  In: Proceedings  of  the  2018  British  Machine  Vision  Conference (BMVC). Newcastle, UK: BMVA, 2018.
    [18] Chang B, Zhang Q, Pan S Y, Meng L L. Generating handwritten Chinese characters using cyclegan. In: Proceedings of the 2018 IEEE Winter Conference on Applications of Computer Vision (WACV). Lake Tahoe, NV, USA: IEEE, 2018. 199−207
    [19] Yeh R A, Chen C, Lim T Y, Schwing A G, Hasegawa-Johnson M, Do M N. Semantic image inpainting with deep generative models. In: Proceedings of the 2007 IEEE Conference on Computer Vision and Pattern Recognition (CVPR). Honolulu, HI, USA: IEEE, 2017. 6882−6890
    [20] 付晓, 沈远彤, 李宏伟, 程晓梅. 基于半监督编码生成对抗网络的图像分类模型. 自动化学报, 2020, 46(3): 531-539

    Fu Xiao, Shen Yuan-Tong, Li Hong-Wei, Cheng Xiao-Mei. A Semi-supervised encoder generative adversarial networks model for image classification. ACTA AUTOMATICA SINICA, 2020, 46(3): 531-539
    [21] Gao Y P, Gao L, Li X Y. A generative adversarial network based deep learning method for low-quality defect image reconstruction and recognition. IEEE Transactions on Industrial Informatics, 2021, 17(5): 3231−3240
    [22] You C Y, Cong W X, Vannier M W, Saha P K, Hoffman E A, Wang G, et al. CT super-resolution GAN constrained by the identical, residual, and cycle learning ensemble (GAN-CIRCLE). IEEE Transactions on Medical Imaging, 2020, 39(1): 188−203
    [23] Liu M, Ding Y K, Xia M, Liu X, Ding E R, Zuo W M, Weng S L. Stgan:  A  unified  selective  transfer  network  for  arbitrary  image attribute editing. In: Proceedings of the 2019 IEEE Conference  on  Computer  Vision  and  Pattern  Recognition  (CVRP). Long Beach, CA, USA: IEEE, 2019. 3668−3677
    [24] Zhong Z L, Li J, Clausi D A, Wong A. Generative adversarial networks and conditional random fields for hyperspectral image classification. IEEE Transactions on Cybernetics, 2020, 50(7): 3318−3329
    [25] Nguyen T D, Le T, Vu H T, Phung D. Dual discriminator generative adversarial nets.  In: Proceedings of the 31st Internation Conference on Neural Information Processing Systems (NIPS). Long Beach, CA, USA: Curran Associates, 2017. 2671−2681
    [26] Memisevic R. Gradient-based learning of the 2011 higher-order image features. In: Proceedings of the 2011 International Conference on Computer Vision. Barcelona, Spain: IEEE, 2011, 1591−1598
    [27] Zhang S, Wang M, Xiong J, Liu S, Chen P Y. Improved linear convergence of training CNNs with generalizability guarantees: A one-hidden-layer case, IEEE Transactions on Neural Networks and Learning Systems, 2021, 32(6): 2622−2635
    [28] Shi Y, Davaslioglu K, Sagduyu Y E. Generative adversarial network in the air: deep adversarial learning for wireless signal spoofing. IEEE Transactions on Cognitive Communications and Networking, 2021, 7(1): 294−303
    [29] Chang J, Scherer S. Learning representations of emotional speech with deep convolutional generative adversarial networks. In: Proceedings of the 2017 IEEE  International Conference on Acoustic, Speech and Signal Processing. New Orleans, LA USA: IEEE, 2017. 2746−2750
    [30] Meng Y Y,  Kong  D Q, Zhu Z F, Zhao Y. From night to day: GANs based low quality image enhancement. Neural Processing Letters, 2019, 50(1): 799−814
    [31] 马锦卫, 彝文起源及其发展考论 [博士学位论文], 西南大学, 中国, 2010.

    Ma Jin-Wei. On the Origin and Development of Yi Language [Ph.D. dissertation]. Southwest University, China, 2010.
    [32] Arcos-García, álvaro, álvarez-García, Juan A, Soria-Morillo L M. Deep neural network for traffic sign recognition systems: an analysis of spatial transformers and stochastic optimisation methods. Neural Networks, 2018, 99: 158−165
    [33] Shin H C, et al. Deep convolutional neural networks for computer-aided detection: cnn architectures, dataset characteristics and transfer learning. IEEE Transactions on Medical Imagings, 2016, 35(5): 1285-1297 doi: 10.1109/TMI.2016.2528162
    [34] Chen S X, Han X, Gao W Z, Liu X X, Mo B F. A classification method of oracle materials based on local convolutional neural network framework. IEEE Computer Graphics and Applications, 2020, 40(3): 32−44
  • 期刊类型引用(3)

    1. 王子一,李光亚,张志艺,李旭卿,简丽. 基于D2GAN和颜色迁移的侯马盟书虚拟修复. 国外电子测量技术. 2023(09): 186-192 . 百度学术
    2. 杨能普,周苗,王文昆,张娟,孙永奎. 基于R-P图像注意融合网络的列车轮对踏面缺陷识别. 铁道科学与工程学报. 2023(12): 4811-4822 . 百度学术
    3. 刘栋斌,王慧琴,王可,王展,甄刚. 基于内容风格迁移的残缺稀疏文字图像盲修复方法. 激光与光电子学进展. 2022(24): 106-117 . 百度学术

    其他类型引用(11)

  • 加载中
  • 图(26) / 表(3)
    计量
    • 文章访问数:  987
    • HTML全文浏览量:  472
    • PDF下载量:  179
    • 被引次数: 14
    出版历程
    • 收稿日期:  2019-10-30
    • 录用日期:  2020-04-16
    • 网络出版日期:  2022-02-17
    • 刊出日期:  2022-03-25

    目录

    /

    返回文章
    返回