最小二乘法在造紙堿回收的黑液濃度軟測(cè)量技術(shù)中的應(yīng)用

1、基于造紙堿回收的黑液濃度軟測(cè)量技術(shù)摘要摘要：研究了造紙工業(yè)堿回收蒸發(fā)工段黑液濃度的在線測(cè)量，用數(shù)據(jù)擬合的方法建立了數(shù)學(xué)模型，提出了黑液濃度控制的一種有效方案。關(guān)鍵詞關(guān)鍵詞：黑液濃度，數(shù)據(jù)擬合，串級(jí)控制，前饋反饋控制黑液濃度是紙機(jī)堿回收過程的一個(gè)重要參數(shù)，必須進(jìn)行在線測(cè)量與控制，而國(guó)內(nèi)目前還沒有成熟可靠的黑液濃度在線測(cè)量?jī)x，國(guó)外的又特別昂貴，因此必須研究黑液濃度的在線自動(dòng)檢測(cè)方法。在此我們采用溫度傳感器和壓力傳感器，對(duì)黑液溫度和壓力兩個(gè)參

2、數(shù)進(jìn)行了監(jiān)測(cè)，在此基礎(chǔ)上通過采用數(shù)據(jù)擬合技術(shù)間接求出了黑液濃度，并建立了黑液濃度的數(shù)學(xué)模型。并且提出了黑液濃度的一種有效控制方案。1黑液濃度在線測(cè)量黑液濃度在線測(cè)量1.1數(shù)據(jù)分析數(shù)據(jù)分析為了獲得濃度Q（波美度）、溫度和壓力（MPa）的關(guān)系，我們先對(duì)所采集到的大量數(shù)據(jù)進(jìn)行簡(jiǎn)單分析，我們首先在壓力一定的情況下分析濃度與溫度的關(guān)系：在壓力（MPa）分別為0.0530.0650.071時(shí)取若干組數(shù)據(jù)并繪制成散點(diǎn)圖，從圖中可以看出在壓力一定的情況

3、下濃度與溫度的關(guān)系為非線性關(guān)系，Q隨T變化的趨勢(shì)近似二次函數(shù)關(guān)系，所以我們用二次函數(shù)對(duì)其進(jìn)行擬合：設(shè)Q(PiT)=q(T)=C2T2C1TC0，使用最小二乘法則需要求能使誤差平方和達(dá)到最小值的。即：。由三元函數(shù)極值)(210cccr210ccc????402210))(()(miniiiqTqcccr理論，現(xiàn)在令：，得到三元一次正規(guī)方程組(1)，解方程組：000210?????????crcrcr（1）?????????????????

4、????????????????????????????????????????????????????40240402104044034024034024040240)()()(1iiiiiiiiiiiiiiiiiiiiiiiiTqTTqTTqCCCTTTTTTTT把表1中的原始數(shù)據(jù)代入（1）式得到在不同壓力下的擬合結(jié)果分別為：MPaP053.0?0002.00572.00006.02???TTQMPaP065.0?0003.0064

5、6.0001.02???TTQMPaP071.0?0004.00095..00019.02???TTQ表1壓力一定溫度與濃度的原始數(shù)據(jù)表P0（MPa）T（c?）q（波美度）P1（MPa）T（c?）q（波美度）P2（MPa）T（c?）q（波美度）0.05393.510.50.06596.915.30.07193.617.50.05398.711.30.065100.316.30.07196.618.30.053100.311.70.065

6、102.316.70.07199.819.60.053101.611.90.065104.617.30.071105.821.70.053103.412.20.065107.318.10.071107.222.9數(shù)據(jù)擬合曲線和散點(diǎn)圖如圖1所示小曲面投影區(qū)域內(nèi)，即T0≤T≤T4，P0≤P≤P2，我們就可以利用Q（P，T）求出濃度Q。所以當(dāng)接收到現(xiàn)場(chǎng)傳送來的數(shù)據(jù)P，T時(shí)，我們先判斷P，T所在的范圍，即Pi1≤P≤Pi1，Ti1≤T≤Ti1，

7、然后再在這一范圍內(nèi)利用以上所闡述的思想求出Q（P，T），從而求出P，T所對(duì)應(yīng)的Q。下面就具體說明一下我們這一建模思想。利用數(shù)據(jù)擬合，我們先求出當(dāng)P一定時(shí)Q與T的函數(shù)關(guān)系Q（Pi，T）。設(shè)Q（Pi，T）=a2T2a1Ta0為了方便起見，我們以前面的數(shù)據(jù)為例，闡述二元函數(shù)的建模過程，在此我們直接引用。MPaP053.00?0002.00572.00006.02???TTQMPaP065.01?0003.00646.0001.02???TTQ

8、MPaP071.02?0004.00095..00019.02???TTQ對(duì)于任一個(gè)T從以上求出的三個(gè)表達(dá)式Q（P0，T），Q(P1T)，Q(P2T)可以求出三個(gè)不同的濃度Q而這些濃度值又是同一溫度下不同壓力所對(duì)應(yīng)的濃度，利用這些數(shù)據(jù)，通過最小二乘法擬合確定Q(PT)=b2P2b1Pb0，在求能使誤差平方和達(dá)到最小值的（)(210bbbr）的過程中我們得到正規(guī)方程組如下:210bbb（3）????????????????????????

9、?????????????????????????????????????????????20220202102042032022032022020220)()()(1iiiiiiiiiiiiiiiiiiiiiiiiTPQPTPQPTPQbbbPPPPPPPP（4）??????????????????????????????????????????????????????????????????????20220201204203202

10、2032022020220210)()()(1iiiiiiiiiiiiiiiiiiiiiiiiTPQPTPQPTPQPPPPPPPPbbb從（4）式可知b2b1b0是T的函數(shù)，解方程組可以求出b2b1b0關(guān)于T的解析函數(shù)式從而可以求出黑液濃度關(guān)于溫度和壓力的二元函數(shù)Q（P，T）。2控制方案控制方案在工程中，我們還采用了前饋反饋控制（如圖4所示）。在蒸發(fā)工段，主要干擾是稀黑液的流量，當(dāng)干擾強(qiáng)烈時(shí)，就不能采用上述的反饋調(diào)節(jié)系統(tǒng)，否則就會(huì)引起

11、黑液濃度發(fā)生大幅度波動(dòng)。由于稀黑液流量可以預(yù)先測(cè)量出來，所以此時(shí)可以采用前饋控制。根據(jù)對(duì)象的物質(zhì)平衡條件，計(jì)算出適應(yīng)該擾動(dòng)的調(diào)節(jié)量，設(shè)計(jì)一個(gè)合適的前饋調(diào)節(jié)器即可進(jìn)行控制。無論何時(shí)，只要干擾出現(xiàn)，就立即進(jìn)行校正使得在它影響黑液濃度之前就被抵消掉。由于根據(jù)稀黑液流量設(shè)計(jì)的前饋控制只能克服這一種干擾，當(dāng)其它干擾出現(xiàn)時(shí)，它無能為力，所以我們?cè)谇梆伒幕A(chǔ)上還加上了反饋，構(gòu)成前饋反饋控制系統(tǒng)。這樣既發(fā)揮了前饋控制作用及時(shí)的優(yōu)點(diǎn)，又保持了反饋控制能克