黑龍江省農作物生物質能能值分析與物流網絡路徑優化

任永泰，孫慧偉，郗通通，孫阿夢

（1.東北農業大學 理學院，黑龍江 哈爾濱 150030；2.東北農業大學 工程學院，黑龍江 哈爾濱 150030）

黑龍江省農作物生物質能能值分析與物流網絡路徑優化

任永泰1，孫慧偉2，郗通通2，孫阿夢2

（1.東北農業大學 理學院，黑龍江 哈爾濱 150030；2.東北農業大學 工程學院，黑龍江 哈爾濱 150030）

面對傳統能源儲量減少、能源需求量越來越多、環境污染日益加劇的問題，開發利用新能源已成為全世界共同關注的焦點。通過引入生命周期法對傳統能值分析法加以改進，對生物質能系統做能值分析，以實際定量結果反映出植物生長過程中每種能值的貢獻，為生產勞作提供建議。對建三江區域農作物秸稈的生物質原料物流網絡進行規劃設計，以重心法結合層次分析法優化原料收儲點。并以典型月份的原料運輸為例，基于ArcGIS優化得到最優車輛路徑和原料運輸方案，結果顯示優化后的方案減少了車輛投入、降低了運輸成本。

能值分析；生命周期分析；物流網絡優化；黑龍江；農作物

1 引言

能源植物是生物質能的重要來源之一[1]，傳統化石能源儲量日益減少，環境污染問題越來越嚴重，開發生物質能已成為各國政府關注的焦點，植物生物質能不僅可以帶來顯而易見的經濟效益[2]，還具有緩解能源危機、減少溫室氣體排放量等環境效益。國外對生物質能的研究起步較早，1996年Mitchell和Bridgwater[3]通過建立生物質能供應鏈系統，第一次將生物質原料和能源作為系統的輸入和輸出。2001年Debarata Das,T.NegatVeziroglu[4]提出將植物產生的氫蒸汽重組作為能源的來源。2007年Lynch[5]提出發展“可再生能源+生物質能源”的路線。2008年Rubo Leng[6]將生命周期應用于評價環境影響、能量消耗以及能量效率。相比國外生物質能領域的發展研究，我國對生物質能源的研究與開發應用相對較晚，不過目前我國生物質能源已步入商業化階段。我國北方農村“四位一體”能源生態模式已經取得了顯著的效益[7]，實現了產氣、積肥同步，種植、養殖并舉，能流、物流良性循環的生態目標。張永、陳曉嬌等[8]運用系統動力學原理對生物質能供應鏈的運行進行模擬，通過構建模型，完善農戶、政府與企業的三方協調機制，優化了生物質能供應鏈，降低了成本，提高了供應鏈的運行效率。張宗蘭、劉輝利等人[9]基于我國生物質能的發展現狀，對生物質能的利用與開發路線與技術提出建議，并積極倡導開展國際合作與交流，引進國外先進技術。劉志雄，何曉嵐[10]將技術水平和政府支持等方面與國外相比較，認為我國應加強政策支持，引進先進技術，鼓勵自我研發，積極發展低碳技術[11]。

已有研究在對植物進行能值分析時往往忽略了植物本身的生長周期，為了準確計算植物生物質能的能量投入，本文引入生長周期對傳統能值分析加以改進，定量的給出植物生長過程中每種能值的貢獻，并對原料運輸物流網絡進行設計優化，為生產勞作提出更合理的意見。

2 基本理論與模型構建

2.1 能值分析法及改進

能值分析可以將系統內各生態流都轉換為統一的單位[12]，其提供了衡量和比較各種類別能值投入的統一標準，開辟了定量研究的新方法。能值分析基本指標的計算公式如下：

（1）能值產出率：

（2）環境負荷率：

（3）能值投資率：

（4）可持續發展指數：

其中，R—可再生能值；N—不可再生能值；F—人類經濟社會反饋投入的能值；FN—不可更新資源能值；FR—可更新資源能值；Y—產出的能值。

能值分析的與眾不同之處在于考慮了自然環境資源對生態系統的能值投入。各種環境資源的能量計算公式如下：

其中，S—太陽光能，A—土地面積，I—年均太陽輻射量。

其中,RC是雨水化學能，P是年均降水量，D是雨水密度，ΔG是吉布斯自由能。

其中，Rg為雨水重力勢能，E為平均海拔高度，g為重力加速度。

其中，W為風能，r為空氣密度，c為風阻系數，vg為地轉風速度，G為風速梯度。其中，vg=v·10/6，v為年均風速。

其中，N為表土凈損失，E為風蝕速率，O為土壤有機質含量,EO為土壤有機質的能量。

為了準確計算植物生物質能的能量投入，應考慮植物的生長周期。生命周期分析是分析產品“從誕生到消亡”全過程的一種監督管理方法，也是衡量制品在它的一生中對環境造成的負荷大小的一種評價方法。不同的植物有不同的生命生長周期，同時，要考慮太陽光反射率和雨水蒸發率，以便更準確地計算能量投入。此時各種環境資源投入植物中的能量計算公式應為：

其中，α為太陽光反射率。

其中，β為雨水蒸發率。

其中各種能源的轉化率均可知。

2.2 重心法的收儲點選址模型

重心法就是將物流系統中需求點和資源點看成是一個物體系統分布在平面上，各個點的供需量作為物體的重量，則此物流系統的最佳設施點就是物體系統的重心。

使用重心法構建模型如下：

設有n個供應點，坐標表示為(xi,yi)(i=1,…，2n)待定物流設施的地址坐標(x0,y0)，總運輸費用為C：

化簡得到最佳位置坐標：

2.3 層次分析法

層次分析法是指將一個復雜的多目標決策問題作為一個系統，把單個目標分解為多個目標或準則，進而分解為多指標（或準則、約束）的若干層次，通過定性指標量化算出層次單排序（權數）和總排序，以作為目標（多指標）、多方案優化決策的系統方法。層次分析法是將決策問題按總目標、各層子目標、評價準則直至具體的備選方案的順序分解為不同的層次結構，然后用求解判斷矩陣特征向量的辦法，求得每一層次的各元素對上一層次某元素的優先權重，最后再用加權和的方法求得最終權重，最終權重最大者即為最優方案。

備選收儲點的選取，主要由原料、運輸及相關費用這三種因素決定。由此將備選收儲點分為三個層次結構：目標層、標準層和決策方案層，并引入相對重要性的標度以便比較各備選收儲點的重要性。

運用規范列平均法，在收儲點影響因素兩兩比較矩陣的基礎上求各因素的權重。

2.4 物流網絡路徑優化模型構建

模型基本參數為：C—總成本；C1—運輸成本；C2—庫存成本；C3—運營成本；q—運輸量；K—庫存總量；k'—階段庫存新增量；x—表示工廠的階段消耗量。

生物質原料物流網絡模型即：

（1）運輸成本模型為：

約束為：

（2）階段運輸成本模型為：

約束為：

（3）庫存成本模型為：

3 案例分析

3.1 能值分析

根據建三江農墾統計年鑒的數據，運用改進后的能值分析計算公式，可以得到玉米秸稈的各項能值指標，見表1。

表1 玉米秸稈的能值指標

運用同樣的方法對大豆和水稻秸稈生產系統進行能值分析。

表2顯示，建三江農墾管理局轄下農場的幾種主要農作物中，相比玉米和水稻，大豆的能值轉換率較高，也就是說每單位的大豆蘊含較多的能量。我們可以通過提高產量和減少種植過程中的能值投入兩種方法來降低其能值轉換率。通過對主要農作物能值指標的分析可以看出：三種農作物的ELR值適中，在未來的生產種植過程中，要適當控制能值投資率EIR，以達到既不依賴環境資源又不會影響經濟發展的目的。水稻的可持續發展系數高于其他兩種農作物，黑龍江作為全國糧食主產地，水稻種植系統的可持續發展性良好，適合生物質能產業大力發展。

表2 農作物秸稈的能值指標

生物質資源供應具有季節性，只在每年農作物收獲期有豐富的生物質資源，這就導致高額的收購和運輸費用。且農作物種植區域呈分散分布，擴大了收購范圍，使企業投入較大的人力、財力和物力。依據能值分析結果規劃不同種類農作物的種植面積，得到每個農場的農作物秸稈產量，為生物質原料收儲點的選擇與優化提供依據。

3.2 收儲點選址

根據各農場的地理位置坐標和每個農場向生產工廠提供原料的數量，由重心法模型計算得到四個備選收儲點A、B、C、D，四個收儲點的地理位置分別在前進農場管理區、富興東部、勝利農場三隊以及七星農場第八管理區附近。在考慮收儲點建設的固定費用、經濟管理變動費用以及實際情況是否允許的情況下，應用層次分析法得到各收儲點間的兩兩比較矩陣，并計算各個因素的實際貢獻權重及評價指標。評價指標越高，越適于建立收儲點，為生產工廠供應原料。評價結果見表3。

表3 各備選收儲點各個因素的實際貢獻權重及評價指標

由表3知各收儲點評價指標排名為A、C、D、B，對比各個收儲點的地理位置、交通情況可知：

（1）建三江農墾管理局的收儲點遍布在生產工廠的周圍，有利于及時供應生物質原料。

（2）由于鴨綠河農場和前哨農場的原料供應量相對較少，B收儲點到生產工廠的路程較長，所以B的收儲點廠房建設投入資金可以適當減少。

（3）將收儲點C作為戰略儲備收儲點，以保證生產工廠在遇到緊急情況時可以迅速做出應對。

3.3 制定典型月份運輸方案

工廠從7月開始正式投入使用，為了達到年產量10萬t的目標，每日生產需要的原料量為980t，自然月以30d進行計算。各個收儲點初始庫存量均為10 000t，7月份在收貨期內，收購量為12 300t/月，各收儲點分別收購4 800t/月，則7月的約束條件為：

月末各收儲點的原料存儲量表示為：

各收儲點的庫存成本和運營成本是已知的，分別為45元/t、4.08元/t，則c'=45+4.08=49.08元/t也是常量，則為常量。為達到成本最小化，就要使最小。表示為函數即：minf=12.15×x1+

將約束條件代入運籌學軟件中進行計算，得出7月份最少運輸費用為473 156元。同時可以得到每個收儲點的原料送出量以及月底的原料庫存量，見表4。

表4 7月底各收儲點的原料送出量及原料庫存量

3.4 全年運輸方案制定

以7月份為初始月份制定全年運輸方案。根據各階段運輸成本模型以及每個收儲點每個月份的原料收儲范圍，運用線性規劃模型可求得全年原料的運輸方案以及每個月的最優運輸成本，結果見表5。

表5 全年原料運輸安排

3.5 車輛路徑優化

運輸路徑不同決定了各個收儲點車輛工作時間的不同，為了使運營成本最低，且盡量減少貨車裝載時的等待時間，現對貨車的運輸時間做詳細合理的規劃安排。

收儲點A：對運送原料的12輛貨車進行排號，且每周依次輪換。工作時間貨車依次裝載，每輛車裝載時間為10min，全部裝載完成需要用時120min，貨車從收儲點A到工廠往返需用時2h，可知在第一輛出發的貨車回來時可以全部完成裝載。同理，其他收儲點也是一樣的情況。假設各收儲點貨車統一時間開始運輸工作，用甘特圖表述各貨車到達工廠的時間，便于對工廠卸載機器進行合理安排。由圖1可以看出，同一時間上班，貨車運送原料到工廠的時候，在一定的時段會出現3輛貨車同時到達的情況，因此造成擁堵，且需要3臺卸載車同時工作才可以保證物流暢通。

圖1 貨車到達工廠的時間甘特圖

各收儲點貨車工作時間相同，為了更快更合理的運送原料到工廠，避免不同收儲點貨車同時到達造成卸載擁擠，可以采取不同時間上班制。通過網絡優化，得出：收儲點A的貨車工作開始時間較收儲點B的貨車開始工作時間推遲120min，B、C是同一時間開始工作；B收儲點的送貨車往返一次后休息110min；收儲點C完成兩次往返之后休息40min。工作時間重新安排規劃后，各收儲點的運貨車到達工廠后，卸載原料所用時間用甘特圖表示，如圖2所示。

圖2 各收儲點貨車在工廠的卸載時間

由圖2可以看出，采用不同工作時間制度，可以一定程度的減少等待，而且只需要配備2輛卸載車即可完成卸載工作，節約了設備購置的費用。

4 結論

應用改進后的能值分析定量化得出建三江墾區水稻、玉米、大豆的能值指標，由此得出大豆的能值轉換率最高，通過提高產量和減少種植過程的能值投入來降低其能值轉換率，提高經濟效益。水稻在三種作物中的可持續發展系數最高，其種植系統可持續發展性良好，有利于大力發展生物質能產業。玉米的可持續性發展指數較低。根據研究結果，為生產勞作提供合理的意見，降低成本投入，提高經濟效益。并對生物質企業的原料供應網絡進行設計和優化。計算得到四個最佳收儲點，并選定原料供應量最大、交通運輸便利、地理位置優勢明顯和各方面條件成熟的C（即勝利農場）收儲點作為戰略儲備收儲點。并優化車輛運輸路徑，減少了車輛的投入數量和等待時間，制定了全年的運輸方案，節約了企業的資金投入。

[1]吳創之,周肇秋,陰秀麗,等.我國生物質能源發展現狀與思考[J].農業機械學報,2009,(1):91-99.

[2]丁建勛.我國能源強度的最優規模實證研究[J].天津財經大學學報,2007,(2):50-54.

[3]Mitchell C P,Bridgewater A V,Mackie K L.Bioenergy systems[A]. In:Bridgewater A V,Boocock D G B,eds.Developments in Thermochemical Biomass Conversion[C].London:Chapman& Hall, 1996.

[4]郭昌倫.乙醇汽油汽車燃料消耗量的計算[J].機械工業標準化與質量,2011,(457):27-28.

[5]Lynch J.Copenhagen Biofuel Assessment Conference[J].Bioenergy Policy Analysis,2007,(6):45.

[6]胡志遠,譚丕強,樓狄明,等.不同原料制備生物柴油生命周期能耗和排放評價[J].農業工程學報,2006,(11).

[7]Laurent Lardon,Arnaud Helias,Breno Sialve,et al.Life-Cycle Assessment of Biodiesel production from Microalgae[J].Environ Sci Technol,2009,43:6 475-6 481.

[8]張永,陳曉嬌,景月明.生物質能供應鏈的系統動力學建模研究[J].物流技術,2009,(10):114-116.

[9]張永,陳曉嬌,景月明.生物質能供應鏈協調機制研究[J].物流技術,2009,(3):119-122.

[10]張宗蘭,劉輝利,朱義年.我國生物質能利用現狀與展望[J].中外能源,2009,(4):27-32.

[11]劉志熊,何曉嵐.低碳經濟背景下我國生物質能發展分析及比較[J].生態經濟,2012,(1):117-121.

[12]石惠春,孫宏山,甘世民.西北干旱區能值分析理論及應用[J].安徽農業科學,2010,(15):8 125-8 127,8 199.

Emergy Analysis and Logistics Network Routing Optimization of Agricultural Biomass in Heilongjiang

Ren Yongtai1,Sun Huiwei2,Xi Tongtong2,Sun Ameng2
(1.School of Science,Northeast University of Forestry,Harbin 150030; 2.School of Engineering,Northeast University of Forestry,Harbin 150030,China)

In this paper,we introduced the life circle method to improve the traditional emergy analysis process in its application with the biomass energy to quantitatively reflect every emergy contribution in the growth circle of the plant.Then we designed and programmed the biomass raw material logistics network for the Three-river Area,optimized the location of the raw material purchasing outlets using the gravity-AHP method and at the end,with the raw material transportation in a typical month as the example,obtained the optimal vehicle routing and raw material transportation plan based on ArcGIS.

emergy analysis;life circle analysis;logistics network optimization;Heilongjiang;agricultural product

TK6

A

1005-152X(2016)12-0055-05

10.3969/j.issn.1005-152X.2016.12.014

2016-10-28

黑龍江省教育廳科學技術研究項目“黑龍江省農作物生物質能物流能流分析與循環利用功能評價”（12541038）

任永泰（1973-），通訊作者，男，黑龍江安達人，東北農業大學理學院教授，碩士生導師，研究方向：水資源優化與利用、農業系統工程理論與方法、應用統計學等；孫慧偉，男，黑龍江鶴崗人，東北農業大學工業工程專業碩士研究生，研究方向：工業系統分析方法與優化技術；郗通通，男，山東濟寧人，東北農業大學工業工程專業碩士研究生，研究方向：物流工程；孫阿夢，女，山東煙臺人，東北農業大學工業工程專業碩士研究生，研究方向：物流工程。