伸缩保差法求解矿井风网电费优化问题

胡章地，倪小山，李卫，胡文军，徐亮，邱丹丹，刘豪

1.湖北虎渡河科技有限公司，湖北荆州 434300；

2.湖北兴发化工集团股份有限公司，湖北宜昌 443700；

3.宜昌市夷陵区安全生产监督管理局，湖北宜昌 443100；

4.武汉工程大学资源与土木工程学院，湖北武汉 430074

伸缩保差法求解矿井风网电费优化问题

胡章地1，倪小山2，李卫3，胡文军4*，徐亮2，邱丹丹4，刘豪4

1.湖北虎渡河科技有限公司，湖北荆州 434300；

2.湖北兴发化工集团股份有限公司，湖北宜昌 443700；

3.宜昌市夷陵区安全生产监督管理局，湖北宜昌 443100；

4.武汉工程大学资源与土木工程学院，湖北武汉 430074

选用伸缩保差法作为算法，并开发便于操作的软件，求解矿井风网通风电费优化问题.以部分自然分风、部分按需分风通风系统的主扇及辅扇年总电费最小为优化目标，把非固定风量分支的风量及分支增阻、增压调节值作为变量，进行优化计算的初始点由风网解算结果进行调整而得.将风网解算软件和风网优化计算程序进行集成开发出C#应用软件，软件考虑了节点号不连续、自然风压作用等实际问题.软件通过一简例予以调试验算，计算结果表明风网电费优化效果显著.

矿井通风系统；优化；伸缩保差法；软件

随着矿山生产不断发展和变化，我国一些矿山的通风系统一方面存在主要运输行人巷道和采掘工作面通风效果差、污风循环使用现象严重等问题，另一方面又存在风量分配不合理、漏风现象严重、井下辅扇及通风构筑物的设置缺乏科学计算而不合理导致通风电能的浪费［1-5］.因此，对矿井风流进行优化控制已势在必行.本文旨在选用一种较为有效的非线性规划求解方法，并开发便于操作的软件，对矿井风流控制进行优化，在满足用风需要、满足风流基本规律、不妨碍生产运输、技术上可行的条件下，使通风电费最少［6-7］.

1 风网电费优化数学模型

以部分自然分风、部分按需分风通风系统的主扇及辅扇年总电费最小为优化目标，满足节点风量平衡定律、回路风压平衡定律、同一分支（增阻、增压）调节唯一性、调节量上限约束，将非固定风量分支的风量及分支阻力调节值作为未知量，得到式（1）的非线性规划问题［8］：

2 风网电费优化数模的求解方法

大多数约束非线性最优化方法基本围绕下列三个基本思路之一：借助反复的线性逼近把线性方法扩展到非线性规划问题中，采用罚函数把约束非线性问题变换为一系列无约束问题，采用伸缩保差法以便同时容纳可行和不可行的自变量矢量.考虑到解非线性规划问题的相对有效性、计算逻辑的简单性及实际的适用性，本文采用伸缩保差法（又称可变容差法）.

伸缩保差法用于求解带有等式约束和不等式约束（约束皆可为线性或非线性函数）的最优化问题［9］.伸缩保差法的原理是将多约束求优问题转变为单约束求优问题，即：

满足于Φ(k)-T(X)≥0.式中：Φ(k)为公差准则函数，表示在进行第k步探索时给定的可变允许公差，Φ(k)是一个下降序列，随迭代次数的增加逐渐趋近到0；T（X）为约束破坏估计量，表示变量不满足约束的程度.

伸缩保差法的运算过程：给定初始点X（0）和初始单纯形边长l.从初始点X（0）出发，按照无约束的单纯形加速法，对f（X）进行下降迭代.将问题（2）的约束区域适当伸缩，并在近似可行的基础上，借用非线性单纯形法求解非线性规划问题.

3 风网电费优化软件

3.1 风网电费优化软件结构

本文采用全球使用最多的.NET用户界面控件套包DXperience，开发出C#应用软件，软件分3大模块：工程基础模块，用于工程管理与系统帮助；风网解算模块，用于参数设置、风网解算、结果展示；风网优化模块，用于优化参数、风网优化、优化结果.

项目主工程为C#工程，项目核心算法为FORTRAN工程，WANGFOR为风网解算算法，VENTOPT为风网优化算法［10］，项目采用C#与FORTRAN外部调用的方式进行集成，VENTOPT.exe为风网优化程序，系统将直接对FORTRAN程序进行输入输出控制.

3.2 风网电费优化软件几个实际问题的处理

3.2.1 软件支持的数据格式由于实际风网分支数较多（有的可达几百条），数据量庞大，为方便技术人员进行数据管理，设计本软件支持.xml和.xlsx两种数据输入输出格式，利用excel强大的编辑功能方便地实现数据的批量处理.

3.2.2 风网电费优化数模的几点说明调节方式考虑增阻（风窗）调节、增压（辅扇）调节，可调分支为可进行增阻调节或增压调节的分支（二者居其一）［11-12］.

漏风风量视为固定风量，漏风风量分支为不可调节分支，其风阻可参考实际情况选定.

对自然风压的处理：把各分支的自然风压作为其恒压头代入计算.

3.2.3 风网节点号不连续问题的处理实际风网比较复杂，在进行节点编号（从1开始）和修改时，往往不能保证节点连续编号.本软件采取了措施，使得在节点不连续编号时不致产生计算和输入数据混乱.

3.2.4 伸缩保差法初始点的选取根据前述风网解算提供的结果，对其中的Cg（可调性标识符）、Q（z非固定风量）、S、p等参数进行适当修改，并调整l、Ep（控制精度）等参数，作为伸缩保差法初始点，进行反复试算.

3.2.5 伸缩保差法自变量的数值变换由于本文数学模型中未知变量的量级相差悬殊，即分支j的增阻值Sj与分支j的增压值pj和分支j的风量Qj相差很大，对初始多面体大小的选择造成困难，故根据各参数上限值对其进行数值变换（除以相应的上限值），使其均在0至1之间变化［13-15］.

4 简例优化计算过程

简例如图1所示.经计算，风网分支数为7，风网节点数为5，风网最大节点号为5，固定风量分支数为2，主扇分支数为1，独立网孔数为3，可调分支数为2.变量数为9，等式约束数为8，全部约束数为16.网孔1的分支：5，-4，-2，3；网孔2的分支：6，-2，3；网孔3的分支：7，1，2，4.可调分支：（5），（6）；主扇分支：（7）.

图1 风网简例Fig.1One case of ventilation network

本例构成问题（3）的非线性规划问题：

上述数模共有未知量：Qz1、Qz2、Qz3、Qz4、Qz7、S5、p5、S6、p6.数模约束共有16个（其中等式约束8个）.本例采用伸缩保差法进行优化计算，取l=0.010 0，Ep=0.003，=0.098 1 N·s2/m8，= 289.40 Pa，C=0.600元/千瓦时，tv=8 760.00 h.优化初始点以风网解算结果为基础进行修改而得.

优化计算总次数为106次，收敛精度为1.72×10-3，目标函数值为1.02×105.简例风网解算和优化结果见表1，其中，NO.为分支号，JA为始节点号，JB为末节点号，PAHP为风网中风机功率.

表1 简例风网解算和优化结果对照表Tab.1Comparison of results of network calculation and network optimization

从表1看出，简例风网经过优化后主、辅扇功耗之和（19.38 kW），小于一般风网解算结果（24.45 kW）.

5 结语

1）本文选择伸缩保差法求解矿井风网电费优化非线性规划问题，其优化计算的初始点由矿井风网解算结果进行调整而得；考虑到自变量参数的量级悬殊可能导致迭代失败，根据自变量的上限值及下限值，对自变量进行数值变换，使其均在0～1之间变化.

2）本文开发的C#应用软件从界面和程序框架2个方面对软件后期的维护升级预留接口，将矿井风网解算软件和风网优化计算程序进行集成；软件支持.xml和.xlsx两种数据输入输出格式，利用excel强大的编辑功能可以方便实现数据的批量处理；软件考虑了节点编号不连续、自然风压作用等实际问题.

3）本文开发的软件以一简例予以调试、验算通过，其优化计算结果表明，简例风网经过优化后主、辅扇功耗之和小于一般风网解算结果，优化效果较为明显.

4）下一步的研究工作方向：对某磷矿通风系统进行实际工程应用，得出该矿通风系统风流优化控制方案.

［1］陈国平，阮德超，胡文军，等.大冶铁矿尖林山车间井下通风系统改造［J］.武汉工程大学学报，2013，35（6）：80-86.

CHEN G P，RUAN D C，HU W J，et al.Reconstruction of Jianlinshan mine ventilation system in Daye iron mine［J］.Journal of Wuhan Institute of Technology，2013，35（6）：80-86.

［2］文永胜.金属矿山地下开采通风系统优化的重要性［J］.金属世界，2009（2）：63-66.

WENYS.Theimportanceofventilationsystem optimization for underground mining at metal mines［J］. Metal World，2009（2）：63-66.

［3］赵千里，刘剑.用矿井通风仿真系统（MVSS）确定通风系统优化改造方案［J］.中国安全科学学报，2002，12（2）：14-18.

ZHAO Q L，LIU J.Optimizing the reform protocol for ventilation system by mine ventilation simulation system［J］.China Safety Science Journal，2002，12（2）：14-18.

［4］杨文杰，贺昶友.昭通市铅锌矿井下通风系统改造设计风量计算［J］.云南冶金，2009，38（3）：8-12.

YANG W J，HE C Y.Air quantity calculation for the reformation design of underground ventilation system in Zhaotong lead＆zinc mine［J］.Yunnan Metallurgy，2009，38（3）：8-12.

［5］姜勇国.矿井通风管理及其信息管理系统［J］.煤炭技术，2008，27（5）：169-170.

JIANG Y G.Mine ventilation management and informa-tion management system［J］.Coal Technology，2008，27（5）：169-170.

［6］叶巍巍，谢贤平，李镇.金属矿山井下通风系统网络结构优化分析［J］.矿冶，2014，23（2）：40-48.

YE W W，XIE X P，LI Z.Analysis of optimization of underground ventilation system of metal mine［J］. Mining&Metallurgy，2014，23（2）：40-48.

［7］梅林芳，吴冷峻.大冶铁矿井下通风系统改造设计［J］.矿业快报，2007，23（4）：63-64.

MEI L F，WU L J.Design for reconstruction of underground ventilation system in Daye iron mine［J］.Express Information of Mining Industry，2007，23（4）：63-64.

［8］胡文军.一种风量调节优化数模及其解法［J］.暖通空调，1997（1）：10-12.

HU W J.Mathematical model for optimal air volume regulating in mine ventilation［J］.Heating Ventilating &Air Conditioning，1997（1）：10-12.

［9］罗雄麟，杨斌，杨洁琦，等.全局校正的可变容差法及其在化工过程在线优化中的应用［J］.化工进展，2010，29（10）：1832-1839.

LUO X L，YANG B，YANG J Q，et al.Global-guided flexible tolerance method and its application to on-line optimization of chemical process［J］.Chemical Industry and Engineering Progress，2010，29（10）：1832-1839.

［10］万耀青，梁庚荣，陈志强.最优化计算方法常用程序汇编［M］.北京：工人出版社，1983：220-257.

［11］王德明.矿井通风与安全［M］.徐州：中国矿业大学出版社，2007：199-218.

［12］支学艺，张红婴.矿井通风与防尘［M］.北京：化学工业出版社，2013：81-103.

［13］王丽铮，席龙飞.可变容差法在船舶最优化设计中的应用［J］.武汉造船，1996（3）：1-4. WANG L Z，XI L F.Application of flexible tolerance method in optimal ship design［J］.Wuhan Shipbuilding，1996（3）：1-4.

［14］刘敏华，郑洁，严彩梅.基于可变容差法的电力系统最优潮流计算［J］.电力自动化设备，2003，23（11）：34-36.

LIU M H，ZHENG J，YAN C M.Optimal power flow computationinelectricpowersystemsbasedon flexible tolerance method［J］.Electric Power Automation Equipment，2003，23（11）：34-36.

［15］王从陆.复杂矿井通风网络解算及参数可调度研究［D］.长沙：中南大学，2003.

本文编辑：苗变

Optimization of Electricity Cost of Mine Ventilation Network by Flexible Tolerance Method

HU Zhangdi1，NI Xiaoshan2，LI Wei3，HU Wenjun4*，XU Liang2，QIU Dandan4，LIU Hao4
1.Hubei Huduhe Technology Co.，LTD，Jingzhou 434300，China；
2.Hubei Xingfa Chemical Group Co.，LTD，Yichang 443700，China；
3.Administration of Work Safety in Yiling District of Yichang，Yichang 443100，China；
4.School of Resource and Civil Engineering，Wuhan Institute of Technology，Wuhan 430074，China

s：We chose the flexible tolerance method as solution method，and developed a software to optimize the electricity cost of mine network ventilation.We took the annual minimum electricity cost of main and auxiliary fans in ventilation system，where natural and on-demand air splitting both exist，as the objective.The air quantity of unfixed-quantity branches，regulated value of increased resistance and air pressure of branches were selected as variables，and the initial values were obtained through adjusting calculation results of ventilation network in mine.A software based on C#was developed by integrating calculation software of ventilation network with ventilation optimization program.This software also takes discontinuous node numbers and natural air pressure into consideration.The software was debugged successfully by a case and the calculation result proves the obvious optimization effect of electricity consumption.

mine ventilation system；optimization；flexible tolerance method；software

TD725

A

10.3969/j.issn.1674-2869.2017.02.014

1674-2869（2017）02-0182-04

2016-09-14

胡章地，硕士.E-mail：hzd62514931@163.com

*通讯作者：胡文军，硕士，教授.E-mail：huwenjunwuhan@126.com

胡章地，倪小山，李卫，等.伸缩保差法求解矿井风网电费优化问题［J］.武汉工程大学学报，2017，39（2）：182-185.

HU Z D，NI X S，LI W，et al.Optimization of electricticy cost of mine ventilation network by flexible tolerance method［J］.Journal of Wuhan Institute of Technology，2017，39（2）：182-185.