1     引言

    工業(yè)實現信息化發(fā)展是工業(yè)現代化的內在要求和必然趨勢。對于水泥工業(yè)而言，信息化技術是水泥工業(yè)邁向先進制造業(yè)的基石，同時也為水泥生產工藝節(jié)能提供了技術支撐[1]。

    水泥工業(yè)從原料礦山開采、生料制備、熟料煅燒、水泥粉磨各生產過程都廣泛結合現代科技成果采用了信息化技術，目前對生產時消耗的能源情況也逐步開展了信息化管理工作。

    在一個實際工業(yè)系統中，當一些變量影響目標時，決策變量（參數、生產條件）與目標間的函數關系必定存在一個以上的峰（谷），從而可以尋求合理的決策變量來達到優(yōu)化目標的過程稱為優(yōu)化[2]。水泥生產各過程中信息化建設離不開軟件系統的支撐。軟件系統的基礎在于選取合適的策略和數學模型進行描述，并且有高效的軟件運算工具。

2    LINGO軟件簡介

    LINGO(Linear Interactive and General Optimizer)軟件，其中文意思是“交互式的線性和通用優(yōu)化求解器”，是由美國LINDO系統公司研制開發(fā)的，可以求解線性規(guī)劃、二次規(guī)劃、非線性規(guī)劃和組合優(yōu)化（圖論）等問題。最大的特色在于它允許優(yōu)化模型中的決策變量為整數（即整數規(guī)劃），而且執(zhí)行速度快。它是一種優(yōu)秀的建模語言，包括許多常用的函數可供使用者調用，并提供了與其他數據文件的接口，易于方便地輸入、輸出、求解和分析大規(guī)模最優(yōu)化問題，并能與ACCESS、EXCLE進行數據交換[3]。

    LINGO模型以MODEL開始，以END結束。中間為語句，，分為四大部分(SECTION)：

   （1）集合部分(SETS)。這部分以“SETS：”開始，以“ENDSETS結束。這部分的作用在于定義必要的變量，便于后面編程進行大規(guī)模計算。

   （2）目標與約束。這部分定義了目標函數、約束條件等，一般要用到LINGO的內部函數。

   （3）數據部分（DATA）。這部分以“DATA：”開始，以“ENDDATA”結束。其作用在于對集合的屬性（數組）輸入必要的數值，主要是方便數據的輸入。

   （4）初始化部分（INIT）。這部分以“INIT：”開始，以“END INIT”結束。其作用在于對集合的屬性（數組）定義初值。對于非線性規(guī)劃求解時，通常得到的是局部最優(yōu)解，而局部最優(yōu)解受輸人的初值影響。通�？筛淖兂跚閬淼玫讲幌虻慕�，從而發(fā)現更好的解[4]。 
 
3    具體實例

3.1 最優(yōu)的水泥生產計劃

    合理安排水泥生產計劃，使得利潤最大。設某廠425#水泥年產量為x1萬噸，525#水泥的年產量為x2萬噸，其費用、售價、生產能力等條件見文獻[5]。

    設L為年度利潤總和，則目標函數為：max L=(100-65.40)x1+(160-75)x2-1052

    約束條件為：

    熟料能力：0.552x1+0.965x2≤38

    水泥能力：131.06x1+138.89x2≤7884

    水渣的可供量：1.028x1+0.656x2≤57

    可比產品工業(yè)總產值：55x1+70x2≥3400

    在LINGO Model里編寫程序見圖1所示。

  
圖1 實例1LINGO程序

    LINGO優(yōu)化后的結果見圖2所示。

    即在該年度內生產425#水泥43.02萬噸，525#水泥14.77萬噸，將獲得最佳經濟效益1692萬元。

 圖2 實例1LINGO程序運行結果

3.2   水泥生產運營計劃

    某水泥廠年生產能力為110萬噸，其主要產品有425#水泥、525#水泥、75℃油井水泥。其各水泥產品的價格、成本、單位利潤見文獻[6]中表1。在水泥生產中，生料磨、窯系統、水泥磨是三大生產環(huán)節(jié)，其各自的臺時產量與運轉率見文獻[6]中表2。同時原材料的年供應能力及單位產品的能源消耗也有其約束。

    設425#水泥、525#水泥、75℃油井水泥的年產量分別為x1、x2、x3噸；425#水泥合同外產量為x4噸，525#水泥合同外產量為x5噸。

    設Z為年度產品的目標利潤總額，則目標函數為：max Z=65(x1-x4)+84(x2-x5)+125x3+100x4+99x5

    約束條件為：

    電：88.52x1+106.75x2+103.42x3≤107100000

    煤：230.6xl+234.2x2+237.8x3≤250195320

    石灰石：729xl+719x2+754x3≤8800000000

    礦渣：44xl+43x2+45x3≤52500000

    原料磨：0.0252x1+0.0253x2+0.0256x3≤26671.75

    熟料窯：0.0408x1+0.0410x2+0.0414x3≤43834.33

    水泥磨：0.057x1+0.057x2+0.057x3≤57120.06

    425#水泥產量：xl≥430000

    525#水泥產量：x2≥400000

    75℃油井水泥產量：100000≤x3≤170000

    425#水泥合同外產量：x4=x1-430000

    525#水泥合同外產量：x5=x2-400000

    在LINGO Model里編寫程序見圖3所示。

圖3 實例2LINGO程序

    LINGO優(yōu)化后的結果見圖4所示。

    即在該年度內生產425#水泥432106.3噸，525#水泥400000噸，75℃油井水泥170000噸，425#水泥合同外產量2106.3噸，將獲得最佳經濟效益8301.063萬元。（注：文獻[4]計算結果有誤。）

圖4 實例1LINGO程序運行結果

3.3  水泥生料配料過程優(yōu)化控制系統

    水泥行業(yè)一直在嘗試運用優(yōu)化技術來對水泥生料配料進行優(yōu)化，起初應用最小二乘法（包括廣義最小二乘法），現在發(fā)展到應用線性和非線性規(guī)劃的措施來進行解決。眾所周知，要求得非線性規(guī)劃的精確解是件非常困難的事情，基本上是利用各種措施轉化成線性規(guī)劃來求得近似解。

    而LINGO的強大功能能較好的處理非線性規(guī)劃問題，也能對此類優(yōu)化問題進行解決。文獻[7]通過對影響水泥配料的數學參數進行分析，認為可采用非線性規(guī)劃數值方法來進行優(yōu)化計算，其實質是參數優(yōu)化問題。參數優(yōu)化就是求出一組參數x={x1,x2,…,xn},以滿足在某種意義上的最優(yōu)。

    通過尋求一組配方，使得三個實際率值對理想率值KH0、SM0、AM0，各自對應的差值最下，建立目標函數如下：

min F(x)=k1∣KH-KH0∣2+k2|SM-SM0|2+k3|AM-AM0|2
   

    約束條件：

    總平衡關系條件：

    熟料中有害含量限制條件：g1(x)=W5-MgOu≤0

    變量非負條件：g1+i(x)=-xi≤0(i=1,2,3,4)

    針對模型中的多目標函數以及非線性約束特點，利用LINGO編制算法程序并通過應用程序接口(API)嵌入到ABB France2000的DCS系統中。

    工業(yè)實驗表明該系統能有效地對配料過程進行優(yōu)化控制，是解決非線性多目標優(yōu)化問題的一種可行方法。

4    結束語

    LINGO具有良好的線性和非線性優(yōu)化性能，并且具有編程方法簡便、收斂性好、可靠性高等優(yōu)點，在工業(yè)生產領域優(yōu)化控制方面將有廣闊前景。目前水泥行業(yè)應用此軟件的還很少，因此有必要加強LINGGO軟件在水泥行業(yè)的推廣應用。

參考文獻

[1]汪瀾，開發(fā)利用信息化技術促進水泥工業(yè)節(jié)能減排[J].中國水泥，2009，（6）：33-35。

[2]何小榮，化工過程優(yōu)化[M].清華大學出版社，2003：1-14。

[3]薛毅，數學建�；A[M].科學出版社，2010。

[4]肖華勇，基于MATLAB和LINGO的數學實驗[M].西北工業(yè)大學出版社，2009。

[5]孫仲英，應用線性規(guī)劃制定最優(yōu)的水泥生產計劃[J].水泥，1988，（6）：19-20。

[6]陸曉恒，淺談運籌學在水泥生產技術中的應用[J].銅陵學院學報，2003，（1）：96-97。

[7]張志剛等，水泥生料配料過程優(yōu)化控制系統[J].計算技術與自動化，2007，26（1）：30-32。

[8]任競爭等，LINGO及其在化工過程優(yōu)化中的應用[J].計算機與應用化學，2010，27（7）：
975-978。

[9]杜宇，優(yōu)化數學模型及LINGO軟件求解[J].黔西南民族師范高等�？茖W校學報，2008，（3）：111-114。