曹三海
摘要 以引黃灌區(qū)管道輸水灌溉工程為研究對象,灌溉工程費(fèi)用的計(jì)算包括工程投資費(fèi)用和工程年運(yùn)行費(fèi)用,采用分?jǐn)傁禂?shù)法計(jì)算灌區(qū)效益。以工程投資函數(shù)、工程年運(yùn)行費(fèi)用函數(shù)、灌區(qū)效益函數(shù)為基礎(chǔ),以灌溉面積為參數(shù),灌區(qū)年凈現(xiàn)值最大為目標(biāo)建立引黃灌區(qū)規(guī)模優(yōu)化模型,對引黃灌區(qū)水資源優(yōu)化配置具有一定的指導(dǎo)意義。
關(guān)鍵詞 引黃灌區(qū);優(yōu)化模型;年凈現(xiàn)值
中圖分類號 S277 文獻(xiàn)標(biāo)識碼 A 文章編號 1007-5739(2015)10-0196-04
Study on Optimization of Pipeline Irrigation Area of the Yellow River Irrigation Project
CAO San-hai
(Gansu Dayu Conserving Water Group Co.,Ltd.,Tianjin 301712)
Abstract The paper is focused on pipeline irrigation area of Yellow River irrigation project. Irrigation engineering cost included the annual operation cost engineering cost and investment projects,the irrigation benefit was calculated by sharing coefficient method.It was based on project investment function,the project annual operation cost function and irrigation benefit function. Yellow River irrigation area scale optimization model with the parameters,the maximum annual net present value as the goal. Irrigation area water resources optimal allocation had certain guiding significance.
Key words Yellow River irrigation areas;optimization model;annual net present value
隨著黃河流域農(nóng)業(yè)生產(chǎn)不斷發(fā)展,農(nóng)業(yè)用水供需矛盾十分突出,提高用水效率、減少水資源浪費(fèi)是保證農(nóng)業(yè)可持續(xù)發(fā)展的根本途徑。目前,黃河流域灌溉平均用水的效率均低于0.5[1]。例如,山東打漁張灌區(qū)灌溉水利用效率僅為0.25;河南人民勝利渠灌溉水利用率為0.4;陜西涇惠灌區(qū)灌溉水的利用效率為0.505;河南新鄉(xiāng)平均灌溉水的利用系數(shù)0.41[2]。
為了實(shí)現(xiàn)水、土、肥的合理利用,緩解灌溉用水的供需矛盾,提高水的利用效率,促進(jìn)農(nóng)業(yè)可持續(xù)發(fā)展,保障國家糧食安全,研究引黃灌區(qū)管道灌溉規(guī)模優(yōu)化模型對提高引黃灌區(qū)水的利用效率、老舊灌區(qū)改造、建設(shè)節(jié)水農(nóng)業(yè)、實(shí)現(xiàn)農(nóng)業(yè)生產(chǎn)現(xiàn)代化具有一定的指導(dǎo)意義。
1 工程費(fèi)用的計(jì)算
工程費(fèi)用主要包括2個方面:①工程投資,即灌溉系統(tǒng)在建設(shè)中所產(chǎn)生的所有費(fèi)用。②系統(tǒng)年運(yùn)行費(fèi)用,即灌溉系統(tǒng)在正常使用過程中所產(chǎn)生的費(fèi)用。
1.1 工程投資費(fèi)用計(jì)算
工程投資費(fèi)用由兩大部分組成,即:水源工程投資(泵站投資)、管網(wǎng)工程投資。
(1)水源工程投資函數(shù)。水源工程投資主要由四部分組成:水泵、電機(jī)、土建、變配電系統(tǒng)。新建泵站,應(yīng)當(dāng)根據(jù)工程建設(shè)的目的,結(jié)合灌區(qū)的具體情況,正確地處理好近期和長期、局部和整體的關(guān)系,因地制宜、因時制宜的進(jìn)行統(tǒng)籌規(guī)劃。泵站投資可用如下公式表示:
K1=Kp·Wp=■(1)
式中:K1—泵站總投資(元);Kp—泵站投資參數(shù)(元/kW);Wp—泵站裝機(jī)容量(kW);η裝—水泵電機(jī)效率;γ—水容重(N/m3);Q—水泵設(shè)計(jì)流量(m3/h);H—水泵揚(yáng)程(m)。
(2)管網(wǎng)投資函數(shù)。管網(wǎng)是管道輸水灌溉工程的重要組成部分,一般可分為輸水管道投資、配水管道投資、附件三大部分。輸配水管道組成了管網(wǎng)的主體,管網(wǎng)中的附件對主干管、干管、支干管等起連接作用,是管網(wǎng)的必要組成部分。管徑和單價的關(guān)系可用如下公式表示:
M=a·Db (2)
式中:M—管道單價(元/m);D—管道內(nèi)徑(mm);a、b—回歸參數(shù),可根據(jù)管道單價與管徑,利用最小二乘法求的。由式(2)可知管網(wǎng)中某一管徑的投資為:
Kj=Mj·Lj=aj·D■■·Lj (3)
式中,Kj—系統(tǒng)中第j種管道投資(元);Mj—系統(tǒng)中第j種管道單價(元/m);Lj—系統(tǒng)中第j種管道長度(m);Dj—系統(tǒng)中第j種管道內(nèi)徑(mm);aj、bj—系統(tǒng)中第j種管道回歸參數(shù),可根據(jù)管道單價與管徑,利用最小二乘法求的。
由式(3)可知管網(wǎng)總投資為:
K2=■Kj·(1+χ)=■(Mj·Lj)·(1+χ)=■(aj·D■■·Lj)·(1+χ) (4)
式中,K2—管網(wǎng)總投資(元);n—管網(wǎng)中所用管道種類數(shù);χ—管道附件投資占管道投資比例。
(3)工程總投資。整個系統(tǒng)投資為各部分投資之和。對于灌溉管網(wǎng)來說,是由水源投資、管網(wǎng)投資兩部分組成,所以總投資K為K1、K2兩部分投資之和。計(jì)算公式如下:
K=K1+K2=■+■(aj·D■■·Lj)·(1+χ)(5)
式中,K—工程總投資(元)。endprint
1.2 工程年運(yùn)行費(fèi)用計(jì)算
管道灌溉系統(tǒng)的年運(yùn)行費(fèi)包括燃料動力費(fèi)、維護(hù)費(fèi)、管理費(fèi)、管理人員工資、其他五部分,其中灌區(qū)維護(hù)費(fèi)、管理費(fèi)、管理人員的工資會隨著系統(tǒng)總投資的增加而增加,可用占系統(tǒng)總投資的比例來進(jìn)行計(jì)算。系統(tǒng)年運(yùn)行費(fèi)用可以用下式表示:
C=C1+C2+C3+C4+C5 (6)
式中:C—管道灌溉系統(tǒng)年運(yùn)行費(fèi)用(元);C1—燃料動力費(fèi)(元);C2—維護(hù)費(fèi)(元);C3—管理費(fèi)(元);C4—管理人員工資(元);C5—其他(元)。
(1)燃料動力費(fèi)的計(jì)算。根據(jù)本文的研究對象,燃料動力費(fèi)一般為提水電費(fèi),即泵站向管網(wǎng)輸水并且滿足灌溉所需壓力所花的動力費(fèi)用。水泵的提水電費(fèi)可以按照以下公式計(jì)算[3]:
C1=■ (7)
式中,E—電費(fèi)單價(元/kW·h)。
水泵的設(shè)計(jì)流量可采用如下公式計(jì)算:
Q=■(8)
式中,A—灌溉面積(hm2);m凈,綜—凈綜合灌水定額(m3/hm2);T—灌溉周期(d);t—每天灌溉時間(h);η水—灌溉水利用系數(shù)。
將式(8)帶入式(7)中得:
C1=■(9)
水泵揚(yáng)程為水泵的提水高度、管網(wǎng)水頭損失、出水口所需工作水頭之和,一般取水泵揚(yáng)程可以按照以下公式計(jì)算:
H=H0+hw+Δ (10)
式中,H0—提水高度(m);hw—管網(wǎng)水頭損失(m);Δ—管道分水口(給水栓)工作水頭(m)。
管網(wǎng)總水頭損失包括各管段沿程水頭損失和局部水頭損失??偹^損失可表示為:
hw=Σhf+Σhj (11)
式中:Σhf—管網(wǎng)沿程水頭損失(m);Σhj—管網(wǎng)局部水頭損失(m)。對于管道灌溉工程來說,沿程水頭損失可以用以下公式計(jì)算:
hf =■(12)
式中,f—摩擦系數(shù),隨水流的雷諾數(shù)Re而變化;L—管道長度(m);Q—管道流量(m3/h);D—管道內(nèi)徑(mm);m—流量指數(shù),與摩阻損失有關(guān);b—管徑指數(shù),該指數(shù)與摩阻損失有關(guān)。
管道灌溉系統(tǒng)管線較長,水頭損失以沿程為主,在計(jì)算中可用占沿程水頭損失的比例計(jì)算局部水頭損失,一般按沿程水頭損失的10%計(jì)算。
hw=1.1Σhf=1.1■■(13)
式中,fj—第j種管道摩擦系數(shù);Lj—第j種管道長度(m);Qj—第j種管道流量(m3/h);Dj—第j種管道內(nèi)徑(mm);mj—第j種管道流量指數(shù);bj—第j種管道管徑指數(shù)。
將式(13)代入式(10)的水泵揚(yáng)程為:
H=H0+1.1■■+Δ (14)
將式(3-13)代入式(3-8)的提水電費(fèi)為:
C1=■·H0+1.1■■+Δ (15)
(2)其他管理費(fèi)用計(jì)算。維護(hù)費(fèi)、管理費(fèi)、管理人員工資、其他的總和可用λK表示,其中K為系統(tǒng)總投資(元);λ為年運(yùn)行費(fèi)系數(shù)(不包括提水電費(fèi))。
(3)工程年運(yùn)行費(fèi)用計(jì)算。計(jì)算公式如下:
C=C1+λK (16)
將式(5)、(14)代入式(15):
C=■·H0+1.1■■+Δ
+λ■+[■(aj·D■■·Lj)]·(1+χ)(17)
(4)工程年費(fèi)用函數(shù)。根據(jù)工程經(jīng)濟(jì)學(xué)理論,系統(tǒng)年費(fèi)用等于投資本利攤還與年運(yùn)行費(fèi)之和,即示為:
Mc=K0+C=K■+λ+C1 (18)
式中,Mc—系統(tǒng)年費(fèi)用(元);K0—系統(tǒng)投資年本利攤還值(元);ic—年利率(%);r—工程使用年限(年)。
2 工程效益的計(jì)算
灌溉工程的效益主要反映在有無灌溉或者在現(xiàn)有的灌溉耕地經(jīng)過改造后農(nóng)作物產(chǎn)品質(zhì)量、產(chǎn)量的增加所帶來產(chǎn)值的增加。產(chǎn)量是各種農(nóng)業(yè)技術(shù)綜合作用的結(jié)果,灌溉是農(nóng)業(yè)增產(chǎn)的手段之一。常用計(jì)算灌區(qū)效益的方法有3種[4]:即分?jǐn)傁禂?shù)法、扣除農(nóng)業(yè)生產(chǎn)費(fèi)用法、以灌溉保證率為參數(shù)推求多年平均增產(chǎn)效益。
本文計(jì)算中采用分?jǐn)傁禂?shù)法,可采用以下公式計(jì)算:
B=A·ε(Y-Y0)·P (19)
式中,B—灌區(qū)水利工程措施分?jǐn)偟亩嗄昶骄旯喔刃б妫ㄔ籄—作物種植面積(hm2);ε—灌溉效益分?jǐn)傁禂?shù);Y—采取灌溉措施后,作物單位面積的多年平均產(chǎn)量(kg/hm2);Y0—未采取灌溉措施前,作物單位面積的多年平均產(chǎn)量(kg/hm2);P—農(nóng)作物產(chǎn)品價格(元/hm2)。
在灌區(qū)內(nèi),作物的種類是比較復(fù)雜的而不是單一的某一種農(nóng)作物,因此灌區(qū)的效益是各種作物效益之和??捎孟率接?jì)算灌區(qū)的效益:
B=■A·αk·εk(Yk-Y0k)·Pk (20)
式中,S—作物種類數(shù);αk—第k種作物種植系數(shù);k—第k種作物灌溉效益分?jǐn)傁禂?shù);Yk—第k種作物采取灌溉措施后,作物單位面積的多年平均產(chǎn)量(kg/hm2);Y0k—第k種作物未采取灌溉措施前,作物單位面積的多年平均產(chǎn)量(kg/hm2);Pk—第k種作物農(nóng)作物產(chǎn)品價格(元/kg)。
3 建立優(yōu)化模型
3.1 目標(biāo)函數(shù)
在確定使用年限的情況下,以灌區(qū)年凈現(xiàn)值最大為優(yōu)化目標(biāo),目標(biāo)函數(shù)為:
NAVmax=max{B-Mc}=max{■A·αk·εk(Yk-Y0k)·Pk
-K[■+λ]-C1}(21)
將式(5)、(12)代入式(18),令F=■+λ得:
NAVmax=max{■A·αk·εk(Yk-Y0k)·Pk-F[■+[■(aj·D■■·Lj)]·(1+χ)]-■·(H0+1.1■■+Δ)}(22)endprint
目標(biāo)函數(shù)(21)分別是由工程效益函數(shù)(20)、系統(tǒng)投資函數(shù)(5)、工程年運(yùn)行費(fèi)用函數(shù)(18)組成。
3.2 約束條件
目標(biāo)函數(shù)(22)受以下條件約束:
(1)灌溉面積。灌溉面積應(yīng)介于最小灌溉面積和最大灌溉面積之間。最小灌溉面積即為灌區(qū)設(shè)計(jì)時所規(guī)劃的最小面積,最大灌溉面積是指受到灌溉區(qū)域地物地貌及氣候等自然條件限制,適合農(nóng)業(yè)灌溉的最大耕地面積。
Mc=0.111 9×{1.4×106A(3.2+5×106■+1.12×107■+4.4×104V32.385A0.385)+[107.4×(■)1.862+171.1×(■)1.862+0.12A×(■)1.862]}+150 839.1×(3.2+5×106■+1.12×107■+4.4×104V32.385A0.385) (元)
(5)工程效益計(jì)算。灌區(qū)小麥、夏玉米和春玉米的灌溉效益分?jǐn)傁禂?shù)εk分別取0.5,0.4,0.40,小麥、春玉米、夏玉米的產(chǎn)量Yk分別取10 500、15 000、15 000 kg/hm2,Y0k分別取6 000、11 070、11 070 kg/hm2,糧食價格Pk分別取2.3、3.2、3.2元/kg。由此,工程灌溉效益為:
B=A×0.75×0.5×(10 500-6 000)×2.3+A×0.25×0.4×(15 000-11 070)×3.2+A×0.75×0.4×(15 000-11 070)×3.2+A×0.75×0.4×(15 000-11 070)×3.2=8 911.65A
(6)目標(biāo)函數(shù)。由以上費(fèi)用、效益函數(shù),本例優(yōu)化模型的目標(biāo)函數(shù)為:
HAV=8 911.65A-0.111 9×{(1.4×106A3.2+5×106■+1.12×107■+4.4×104V32.385A0.385)+[107.4×(■)1.862+171.1×(■)1.862+0.12A×(■)1.862]}-150 839.1×(3.2+5×106■+1.12×107■+4.4×104V32.385A0.385)
現(xiàn)代農(nóng)業(yè)科技2015年10期