呂陳秋

(中國昆侖工程公司，北京100037)

隨著全球化纖生產(chǎn)進一步向中國轉(zhuǎn)移，我國已經(jīng)成為世界最大的化纖生產(chǎn)國，我國化纖產(chǎn)量占據(jù)全球總量的60%以上，在“十一五”期間更是得到了空前發(fā)展。2011年我國化纖產(chǎn)量達到33 900.7 kt，占國內(nèi)紡織纖維加工總量的70%以上［1］?；w產(chǎn)業(yè)發(fā)展迅速，科學準確地預測化纖產(chǎn)量顯得尤為迫切，且準確程度將直接影響到投資、產(chǎn)業(yè)布局和運行的合理性，不僅對化纖行業(yè)的健康發(fā)展，甚至對整個國民經(jīng)濟的發(fā)展有著重要的意義［2－4］。作者通過對化纖原料消費量的分析，找出化纖產(chǎn)量與各種原料消費量之間的內(nèi)在聯(lián)系，從而預測出化纖產(chǎn)量在未來時期的變化趨勢。

1 模型建立

1.1 方法

采用定量預測方法中的多元線性回歸法建立化纖產(chǎn)量預測模型。多元回歸分析含有2個或2個以上的自變量，由多個自變量的最優(yōu)組合共同來預測或估計因變量，比只用一個自變量進行預測更有效。因此多元線性回歸法比一元線性回歸法的實用意義更大［5－6］。

多元線性回歸的數(shù)學模型一般表示為［5－8］:

式中:ε為隨機誤差;β0為常數(shù)項;βi為偏回歸系數(shù)，i=1，2，3，…，n 。

自變量y的變化由兩部分引起，一是n個因變量x的變化;二是其他隨機因素。模型建立后，對給定的預測點 x1，x2，x3…xn，便可求出對應的預測點的預測值和預測置信區(qū)間［6］。

1.2 模型變量的選取

影響化纖產(chǎn)量的因素有很多，有原料供應、進出口關稅、國家匯率、原油價格、政策與法規(guī)等，選取原料消耗量作為影響因素，因為消耗多少確定產(chǎn)出多少，而且消耗量已經(jīng)包含了原油價格、進出口關稅、原料供應等影響因素，所以以原料消耗量作為模型的因變量是合適的。2001—2011年我國化纖產(chǎn)量分別為 8 431.8，9 921.0，11 811.5，14 245.4，16 647.9，20 731.8，24 137.8，24 050.0，27 478.7，30 897.0，33 900.7 kt，主要化纖的產(chǎn)量如表1所示。

表1 2001—2011年我國主要化纖產(chǎn)量Tab.1 China chemical fiber outptut over 2001 －2011 kt

從表1可以看出，化纖各類產(chǎn)品中滌綸占的比例最大，其次是粘膠、錦綸和腈綸，而丙綸、維綸和氨綸所占比例最少，其相加值不足4%。因此忽略丙綸、維綸和氨綸對化纖產(chǎn)量的影響，只考慮滌綸、錦綸、腈綸和粘膠對化纖產(chǎn)量的影響。對苯二甲酸(PTA)、己內(nèi)酰胺(CPL)、己二酸(AA)、丙烯腈(AN)、漿粕(CP)作為化纖行業(yè)的重要原料，主要用于生產(chǎn)滌綸、錦綸、腈綸和粘膠，因此以PTA，CPL，AA，AN，CP 的消耗量作為自變量，建立五元線性回歸方程進行模擬。

1.3 數(shù)據(jù)處理

原始數(shù)據(jù)來源于中國統(tǒng)計年鑒2011［1］和中國化纖經(jīng)濟信息網(wǎng)(見表2)，采用統(tǒng)計分析軟件PASW/SPSS Statistics 18對原始數(shù)據(jù)進行處理［9］(見表3)。

表2 2001—2011年化纖原料消費量Tab.2 Raw material consumption for chemical fiber over 2001－2011 kt

表3 回歸系數(shù)Tab.3 Regression coefficient

2 結(jié)果與討論

2.1 回歸方程

根據(jù)表3所示的回歸方程的偏回歸系數(shù)，可以得到完整的多元回歸方程:

式中:y為化纖產(chǎn)量;x1為PTA消費量;x2為CPL消費量;x3為AN消費量;x4為AA消費量;x5為CP消費量。

2.2 相關性檢驗

R為多元相關系數(shù)，R2為0.998，調(diào)整R2為0.996，標準估計誤差為568.836。由此進行相關性即擬合優(yōu)度檢驗［10］。參考調(diào)整R2數(shù)值，其值接近1，表明擬合度較高，說明化纖產(chǎn)量的99.6%由原料 PTA，CPL，AA，AN，CP 的消耗量決定。

2.3 顯著性檢驗

由表4可以看出，F(xiàn)檢驗統(tǒng)計量的觀測值較大，為461.896，大于F值概率為0，結(jié)果表明最終的回歸方程應包含CP，AN，AA，CPL，PTA消費5個變量，顯著性很好。

表4 離散分析數(shù)據(jù)Tab.4 Discrete analysis data

2.4 適合性檢驗

模型的適合性檢驗主要是殘差分析［9］，見表5。通過殘差分析，可以確定數(shù)據(jù)的可靠性、周期性或其他干擾。從圖1可看出，各觀測點都分布在對角線上，據(jù)此可初步判斷殘差符合正態(tài)分布。

表5 殘差統(tǒng)計量Tab.5 Residual statistics

圖1 回歸標準化殘差的標準P-PFig.1 Standard P-P plot of standard regression residual

從圖2可看出，各散點分布在標準化殘差為0，并以此為中心的橫帶中，證明模型是適合的。

圖2 散點圖Fig.2 Scatter diagram

3 預測

3.1 原料消費量預測

從表6可看出，PTA的消費主要用于生產(chǎn)聚酯，而聚酯產(chǎn)品中占比例最多的還是滌綸，瓶片和膜片所占比例較少?！笆晃濉逼陂g，非纖用聚酯比例雖然增加不大，但產(chǎn)量增加很大，尤其是瓶級聚酯。據(jù)統(tǒng)計，2011年，我國瓶級聚酯有接近5 000 kt的產(chǎn)能，但國內(nèi)需求卻不足3 000 kt，產(chǎn)能嚴重過剩。隨著內(nèi)需的擴大、中西部的開發(fā)和啤酒瓶市場的開拓，瓶級聚酯在“十二五”期間仍會保持增長，但增長速度應有所緩慢，所占聚酯比例也應有所下降。膜級聚酯將是“十二五”期間聚酯產(chǎn)品中的重要產(chǎn)品，將會有突飛猛進的增長趨勢，預計2015年在聚酯產(chǎn)品中的比例將突破5%?！笆濉逼陂g，纖維用聚酯將保持持續(xù)增長趨勢，雖然常規(guī)滌綸目前已處于飽和，但功能性、差別化滌綸目前還不能滿足市場需求，高阻燃、抗熔滴、高導濕、抗靜電、抗紫外線、抗菌防臭、防輻射等單一功能和復合功能的新型滌綸是滌綸行業(yè)的發(fā)展方向［2，10－11］。

表6 “十一五”期間我國聚酯產(chǎn)品產(chǎn)量Tab.6 Output of China polyester products during the eleventh five-year period

根據(jù)以上分析，PTA消費量大致保持每年穩(wěn)步增長趨勢［12－13］，平均每年增長 1 735 kt(見圖3)，根據(jù)曲線模擬可以預測得到2015年PTA消費量為 29 021.0 kt。

圖3 PTA消費量增長曲線Fig.3 PTA consumption growth curve

近年來，我國各錦綸企業(yè)進行了產(chǎn)業(yè)調(diào)整，逐漸形成規(guī)?；硗?，在市場需求的強烈刺激下，錦綸聚合產(chǎn)能進入快速增長期［14－18］，使主要原料CPL和AA的消費量急劇上升，特別是“十二五”指出重點發(fā)展CPL的原料建設［2］，根據(jù)Origin擬合曲線(見圖4)，預計到2015年CPL和AA的消費量分別能達到 1 626.6 kt和 1 411.5 kt。

圖4 CPL與AA和AN消費量增長曲線Fig.4 CPL，AA and AN consumption growth curves

國內(nèi)腈綸需求旺盛，而生產(chǎn)能力遠低于需求，一直靠進口彌補，因此從2001開始，腈綸產(chǎn)量高速增長，然而這種跳躍式擴張僅維持到2005年。從2006年開始，國際原油價格一路飆升，使腈綸原料AN的價格一路上漲，另外，常規(guī)腈綸已處于飽和狀態(tài)［19－22］，國內(nèi)企業(yè)規(guī)模偏小，均使腈綸行業(yè)受到多方面壓力，使得AN的消耗量增長緩慢，利用Origin進行曲線擬合得2015年AN消費量為 1 722.2 kt。

雖然粘膠制造行業(yè)屬于高能耗高污染行業(yè)［23，24］，但粘膠纖維的需求仍不斷上漲，而且“十二五”指出要加快發(fā)展人纖木漿的原料建設［2，25］，因此漿粕的消費量在未來幾年仍然會保持持續(xù)增長趨勢，如圖5所示，預計2015年CP的消耗量將達到 4 108.6 kt。

圖5 CP消費量增長曲線Fig.5 CP consumption growth curves

3.2 化纖產(chǎn)量點預測及區(qū)間預測

根據(jù)預測的2015年各原料的消耗量，再通過以上擬合的五元回歸模型(式2)，即可預測得到2015年化纖產(chǎn)量為41 121.4 kt。

由于現(xiàn)實情況的變化和各種環(huán)境因素的影響，預測的實際值總會與預測值產(chǎn)生或大或小的偏移，因此在點預測的基礎上，通過PASW Statistics 18進行分析，得出區(qū)間預測值，如表7所示。

表7 2015年化纖產(chǎn)量預測及置信區(qū)間Tab.7 Prediction and confidence limit of chemical fiber output in 2015

4 結(jié)論

a.通過PASW/SPSS Statistics對化纖產(chǎn)量進行多元回歸分析，擬合的化纖產(chǎn)量和各原料消費量的多元模型相關性、顯著性、適合性均明顯。

b.根據(jù)多元模型預測，2015年化纖產(chǎn)量為41 121.4 kt，接近“十二五”規(guī)定產(chǎn)量目標41 000 kt，其95%置信度的預測值和個別值上限分別為46 907.2 kt和47 089.2 kt，與“十二五”規(guī)定產(chǎn)能目標46 000 kt接近。

［1］中華人民共和國國家統(tǒng)計局.中國統(tǒng)計年鑒［EB/OL］.［2013－01 － 01］.http://www.stats.gov.cn/tjsj/ndsj/2012/indexch.htm.

［2］工業(yè)和信息化部.化纖行業(yè)“十二五”發(fā)展規(guī)劃［EB/OL］.［2012－01－19］.http://www.gov.cn/gzdt/2012－01/19/content_2049031.htm.

［3］金離塵.我國GDP增長率走高將促進化纖工業(yè)持續(xù)增長［J］.合成纖維工業(yè)，2004，27(2):31－34.

［4］張神勇.中國化纖業(yè)發(fā)展:需求拉動與資源約束［J］.東華大學學報:自然科學版，2010，36(3):332－337.

［5］王惠文，孟潔.多元線性回歸的預測建模方法［J］.北京航空航天大學學報，2007，33(4):500－504.

［6］彭鵬，彭家紅.基于多元線性回歸模型的電力負荷預測研究［J］.中國安全生產(chǎn)科學技術，2011，7(9):158－161.

［7］宋云雪，張科星，史永勝.基于多元線性回歸的發(fā)動機性能參數(shù)預測［J］.航空動力學報，2009，24(2):427－431.

［8］王光輝，范程，田文中.焦炭熱態(tài)性質(zhì)預測模型的研究［J］.燃料與化工，2009，40(1):1－3.

［9］李志輝，羅平.PASW/SPSS Statistics中文版統(tǒng)計分析教程［M］.北京:電子工業(yè)出版社，2010:300－319.

［10］陳向玲，宋睿，趙萬金，等.國內(nèi)外滌綸產(chǎn)品與技術發(fā)展趨勢［J］.合成纖維工業(yè)，2011，34(2):47－49.

［11］姜兆輝，白瑛，金劍，等.差別化聚酯纖維技術研究進展［J］.合成纖維工業(yè)，2011，34(1):52－46.

［12］汪英枝.我國對苯二甲酸產(chǎn)業(yè)鏈市場現(xiàn)狀與展望［J］.化學工業(yè)，2011，29(10):19－23.

［13］楊愛民.近期 PTA行情的走勢［J］.聚酯工業(yè)，2011，24(1):7－9.

［14］錢伯章.中國錦綸發(fā)展現(xiàn)狀［J］.合成纖維，2011，40(7):55.

［15］侯湘湘.己內(nèi)酰胺的供需現(xiàn)狀及發(fā)展趨勢［J］.合成纖維工業(yè)，2011，34(4):44－47.

［16］錢伯章.己內(nèi)酰胺的市場分析［J］.合成纖維工業(yè)，2009，32(2):23－23.

［17］宋超，文夢君，余毅.聚酰胺纖維生產(chǎn)現(xiàn)狀及發(fā)展展望［J］.合成纖維工業(yè)，2012，35(1):49－53.

［18］金離塵.我國聚酰胺纖維產(chǎn)業(yè)鏈的發(fā)展［J］.合成纖維工業(yè)，2009，32(1):36－39.

［19］鄒志量，夏延致，記全，等.中國腈綸工業(yè)的現(xiàn)狀和發(fā)展方向的探討［J］.合成纖維，2009，38(6):1－5.

［20］盧東，李青山，楊秀珍.中國腈綸行業(yè)的發(fā)展現(xiàn)狀與展望［J］.合成纖維，2006，35(6):28－31.

［21］錢伯章.我國腈綸生產(chǎn)現(xiàn)狀與發(fā)展［J］.合成纖維工業(yè)，2009，32(5):44 －46.

［22］金離塵.我國產(chǎn)業(yè)用腈綸的現(xiàn)狀與發(fā)展［J］.合成纖維工業(yè)，2009，32(4):37 －40.

［23］趙利娜，楊波.粘膠與染整廢水綜合處理工程與實例［J］.水處理技術，2010，36(10):119－121.

［24］彭志，周永超.粘膠短纖行業(yè)發(fā)展分析［J］.經(jīng)濟研究導刊2012(6):191－195.

［25］蘆長椿.“十二五”期間國內(nèi)化纖工業(yè)與可持續(xù)發(fā)展［J］.合成纖維，2012，41(1):1－8.