高 賀,柴士超
(東北林業(yè)大學(xué) 土木工程學(xué)院,黑龍江 哈爾濱 150040)
?
基于模糊數(shù)學(xué)理論的工程成本-工期研究
高 賀,柴士超
(東北林業(yè)大學(xué) 土木工程學(xué)院,黑龍江 哈爾濱 150040)
為了提高工程在不確定環(huán)境下對工程成本的控制,通過對工程成本分析,并依據(jù)工程成本與工程規(guī)模、工期的關(guān)系,建立以工程成本為目標(biāo),人工費(fèi)、材料費(fèi)、機(jī)械費(fèi)為不確定因素影響條件的優(yōu)化模型,該模型通過模糊數(shù)學(xué)理論和方法,結(jié)合四邊形模糊數(shù)的解析方法,將廣義Lagrange乘子法與反證法相結(jié)合,得出模型最優(yōu)解的求解過程,以期為工程成本控制提供一定的借鑒。
工程造價(jià);不確定性;梯形模糊數(shù);廣義拉格朗日乘子法
隨著經(jīng)濟(jì)的快速發(fā)展和客觀環(huán)境的日益復(fù)雜,我國施工企業(yè)必須在項(xiàng)目管理中正確處理成本、工期、質(zhì)量三者的關(guān)系,運(yùn)用技術(shù)、經(jīng)濟(jì)、管理等措施對“三大方向”實(shí)行有效管制,從而獲得較好的經(jīng)濟(jì)效益和社會效益。其中成本管理無疑成為施工企業(yè)最主要的任務(wù)之一。暴克琦[1]通過工程規(guī)模、工期和成本的關(guān)系建立相應(yīng)的數(shù)學(xué)模型,最終得出工期的最優(yōu)解。黃松顯等[2]在不確定因素和人為因素對成本控制的影響與對策中,通過理論闡述了不確定因素下的成本管理措施。雄鷹等[3]采用蟻群算法對成本進(jìn)行優(yōu)化。KENZO等[4]針對工程項(xiàng)目中的不確定性問題,建立了GERT網(wǎng)絡(luò)計(jì)劃模型。LUIS等[5]運(yùn)用蒙特卡洛仿真方法估計(jì)工程成本與工期的可能風(fēng)險(xiǎn)。上述學(xué)者對工程成本的優(yōu)化主要集中在工程項(xiàng)目成本管理的概念、內(nèi)容、特點(diǎn)和項(xiàng)目成本結(jié)構(gòu)描述,或是對工程項(xiàng)目成本管理中成本管理一些方法的應(yīng)用介紹。
筆者首先明確工程成本構(gòu)成,總結(jié)工程成本內(nèi)在規(guī)律;其次分析成本形成影響因素的定性研究;最后通過模糊數(shù)學(xué)理論建立以項(xiàng)目成本為目標(biāo)的數(shù)學(xué)模型,并應(yīng)用廣義Lagrange乘子法求最優(yōu)解,得出工期為最優(yōu)解時(shí)的成本最優(yōu)解。該模型能夠?qū)こ坛杀具M(jìn)行預(yù)測,有利于及時(shí)發(fā)現(xiàn)成本管理中的問題和薄弱環(huán)節(jié),合理控制工程進(jìn)度,繼而控制工程成本。
工程成本是指建造一個(gè)工程項(xiàng)目所需花費(fèi)的全部費(fèi)用,即從工程項(xiàng)目確定建設(shè)意向直至建成、竣工驗(yàn)收的整個(gè)建設(shè)期間所支出的總費(fèi)用。這是保證工程正常施工的資金,是建設(shè)項(xiàng)目投資的重要組成部分。工程成本主要包括兩部分:直接費(fèi)和間接費(fèi)。直接費(fèi)包括人工費(fèi)、材料費(fèi)、機(jī)械使用費(fèi)。間接費(fèi)包括其他直接費(fèi)和間接費(fèi)。其中:人工費(fèi)=用工日×日平均工資;材料費(fèi)=材料耗用量×材料單價(jià);機(jī)械使用費(fèi)=機(jī)械臺班數(shù)×臺班單價(jià);其他直接費(fèi)=工程直接費(fèi)×其他直接費(fèi)費(fèi)率;間接費(fèi)=規(guī)費(fèi)+企業(yè)管理費(fèi),規(guī)費(fèi)=計(jì)算基數(shù)×規(guī)費(fèi)費(fèi)率(以直接費(fèi)為計(jì)算基數(shù)),企業(yè)管理費(fèi)=計(jì)算基數(shù)×企業(yè)管理費(fèi)費(fèi)率(以直接費(fèi)為計(jì)算基數(shù))。根據(jù)工程成本的主要成分,建立以成本為目標(biāo)的基本函數(shù)關(guān)系式:
C=GWX+MQX+TZX+(GW+MQ+TZ)HX+
(GW+MQ+TZ)NX+(GW+MQ+TZ)FX
(1)
式中:X為工程的建筑面積;M為每平方米耗用的工程原料;Q為材料的現(xiàn)實(shí)單元價(jià);G為工程每平方米的工作日;W為日平均工資;T為機(jī)械的每平方米需用臺數(shù);Z為機(jī)械的臺班成本;H為該工程的其他直接費(fèi)的取費(fèi)率;N為該工程規(guī)費(fèi)費(fèi)率的取費(fèi)率;F為該工程的企業(yè)管理費(fèi)費(fèi)率。
羅桂平[6]通過對工程項(xiàng)目成本構(gòu)成與影響因素進(jìn)行研究,依據(jù)工程項(xiàng)目成本組成中各部分成本占總成本比例為:①材料費(fèi)是工程項(xiàng)目成本的主要組成部分,占總成本的50%以上,最高達(dá)到65%。②人工費(fèi)所占比例都在12%以上。③機(jī)械成本所占比例在7%~17%之間。勞動(dòng)力成本、材料成本和機(jī)械成本占總成本的85%~90%,而其他直接費(fèi)和間接費(fèi)占總成本的5%左右,最高不超過7%。
2.1 工程成本、工期、規(guī)模分析
工程建設(shè)項(xiàng)目的目標(biāo)主要包括進(jìn)度目標(biāo)、成本目標(biāo)、質(zhì)量目標(biāo),這三大目標(biāo)既相互對立又相互統(tǒng)一。BADU等[7]在論述進(jìn)度目標(biāo)、成本目標(biāo)和質(zhì)量目標(biāo)的關(guān)系時(shí)指出:工程項(xiàng)目施工速度加快,會影響到工程完成的質(zhì)量,并創(chuàng)建了線性模型來研究進(jìn)度、成本和質(zhì)量之間的平衡關(guān)系。EI-RAYES等[8]論述了進(jìn)度和成本的綜合優(yōu)化問題,并綜合了進(jìn)度和成本的關(guān)系,提出了可行的優(yōu)化方法。HSIEH[9]在論述進(jìn)度-成本-質(zhì)量的優(yōu)化模型中,提出了采用離散的工程項(xiàng)目模型進(jìn)行求解。上述學(xué)者根據(jù)進(jìn)度、成本、質(zhì)量三者之間的關(guān)系,找出了其平衡點(diǎn)。
圖1所示為工程成本與工期關(guān)系圖。對于一個(gè)特定的工程項(xiàng)目,工程總成本會隨著工期的延長先下降后上升,即工期-工程總成本函數(shù)為U形函數(shù)。由圖1可知:直接成本的費(fèi)用是一條曲線,在一般情況下,直接成本會隨著工期的縮短而增加,隨著工期的延長而減少。但是,如果工期無限制地延長,也可能會引起直接成本的增加。間接成本的費(fèi)用近似一條與工期成正比的直線,會隨著工期的縮短而減少,隨著工期的延長而增加。
由于工程項(xiàng)目成本構(gòu)成復(fù)雜,受到宏觀經(jīng)濟(jì)環(huán)境、工程自身特征、施工組織管理、天氣、施工環(huán)境和意外事故等眾多不確定因素的影響,因此為了簡化該理論模型,參考文獻(xiàn)[6]工程項(xiàng)目成本構(gòu)成與影響因素研究中引入的變量,將工程進(jìn)度和工程規(guī)模作為控制變量。工程成本與施工工期的關(guān)系如式(2)所示。
圖1 工程成本與工期關(guān)系圖
C(X,D)=α+α1X/D+α2D
(2)
式中:C為工程總成本;X為工程規(guī)模;D為工程工期;α、α1、α2為常數(shù)且大于0。
2.2 道路工程成本主要影響因素不確定分析
建筑工程項(xiàng)目領(lǐng)域廣,建設(shè)周期長,耗費(fèi)大量人力、物力和財(cái)力,影響因素眾多且情況復(fù)雜,致使工程造價(jià)高,即使工程成本變化幅度較小,也會導(dǎo)致工程造價(jià)出現(xiàn)比較大的偏差。
工程成本的不確定性主要表現(xiàn)在3個(gè)方面:①項(xiàng)目本身存在著不確定性;②項(xiàng)目規(guī)模及其消耗和占用的資源數(shù)量存在著不確定性;③項(xiàng)目所消耗和占用的資源價(jià)格存在不確定性。
影響工程成本不確定性的因素主要有:
(1)工期的不確定性。工期的不確定性直接導(dǎo)致人工費(fèi)與機(jī)械使用費(fèi)的不確定性。根據(jù)項(xiàng)目進(jìn)度和成本的關(guān)系,當(dāng)質(zhì)量和安全在一定的條件下,存在最佳的項(xiàng)目進(jìn)度,項(xiàng)目成本是最低的。
(2)材料價(jià)格的不確定性。價(jià)格水平只對個(gè)別的間接成本費(fèi)用有一定的影響,對大多數(shù)施工項(xiàng)目成本,特別是原材料成本并沒有顯著的影響,這主要是由于項(xiàng)目部在整個(gè)施工期內(nèi),對材料價(jià)格的變動(dòng)有準(zhǔn)確的判斷并且制定了相應(yīng)的措施,使得整體材料費(fèi)用受物價(jià)波動(dòng)的影響不大。
(3) 材料消耗量的不確定性,即生產(chǎn)率的不確定性,工程量清單計(jì)價(jià)中材料的凈用量一定,變化的是材料的損耗量。
(4)實(shí)際工程量的不確定性,例如不可預(yù)見工程量、冬雨季施工的施工附加量,以及作業(yè)空間增加的工程量,若采用工程量清單計(jì)價(jià)模式,投標(biāo)報(bào)價(jià)時(shí)清單工程量己給出,此時(shí),工程量的不確定性體現(xiàn)在綜合單價(jià)變化上。
筆者在影響工程成本不確定因素的分析中,假設(shè)質(zhì)量和安全為一定條件,通過模糊數(shù)學(xué)理論,對工程成本主要組成中的人工費(fèi)、材料費(fèi)、機(jī)械費(fèi)進(jìn)行模糊化處理,結(jié)合四邊形模糊數(shù)解析法,得出工程進(jìn)度的最優(yōu)解,進(jìn)而得出工程成本的最優(yōu)方案。
2.3 模糊數(shù)學(xué)理論的應(yīng)用
從20世紀(jì)90年代起,許多學(xué)者開展了模糊數(shù)學(xué)的應(yīng)用研究,其中MORTEZA[10]在研究短缺生產(chǎn)-庫存模型時(shí),假定全部生產(chǎn)因素為模糊數(shù),使用一套模糊數(shù)的運(yùn)算法則創(chuàng)建無欠缺生產(chǎn)-庫存模型,并用廣義Lagrange乘子法求最優(yōu)解,大大簡化了建模與求解的進(jìn)程。在建筑領(lǐng)域,ABALFAZL等開展掙值法在模糊理論的研究,將梯形模糊數(shù)轉(zhuǎn)化為掙值法,與項(xiàng)目成本和進(jìn)度相結(jié)合,在模糊數(shù)學(xué)理論基礎(chǔ)上預(yù)測項(xiàng)目的竣工時(shí)段。
由于在工程施工中對不確定的因素進(jìn)行定量處理非常復(fù)雜,因此在難以定量處理的方面可用語言變量表示,如目標(biāo)權(quán)系數(shù)可表示為“很低”、“低”、“中”、“高”、“很高”等語言變量的形式。因此按照廣義模糊數(shù)的線性取值圖將這些語言變量值的從屬函數(shù)轉(zhuǎn)變?yōu)椴灰?guī)則四邊形模糊數(shù)。通過對模糊值的大小估計(jì)重心法,則A=(α1,α2,α3,α4)的實(shí)際大小為:
(3)
圖2 廣義模糊數(shù)線性取值圖
其中,0 2.4 工程成本模糊四邊形 根據(jù)式(1)及四邊形模糊數(shù)的運(yùn)算法則,以工程模糊成本為目標(biāo)的模糊四邊形為: C=(Wbg+Qag+Zcg)(1+H+N+F)X(5) 根據(jù)式(3)重心法估計(jì)模糊數(shù)大?。?/p> (6) 其中,Lg=(Wbg+Qag+Zcg)×(1+H+N+F)X,g=1,2,3,4。 2.5 模型的建立 基于工程造價(jià)的定義,建立以成本為目標(biāo)的優(yōu)化模型,根據(jù)變量之間的線性關(guān)系、項(xiàng)目的規(guī)模,結(jié)合梯形模糊數(shù)的理論和方法,采用廣義拉格朗日乘子法和反證法,得到了模型的最優(yōu)解。但對工程造價(jià)與工程規(guī)模之間關(guān)系的分析,表明兩者的線性關(guān)系,并且呈遞增趨勢。因此筆者采用工程成本與工程進(jìn)度兩者之間的關(guān)系建立相應(yīng)的模型。聯(lián)立式(2)和式(5),建立工程成本模糊因素-工期-工程總成本模型: (7) Ig=(Wbg+Qag+Zcg)(1+H+N+F) (8) g=1,2,3,4 在工程工期為模糊狀況下,建立最終工程模糊成本數(shù)學(xué)模型: (9) 利用最優(yōu)乘子法的廣義拉格朗日解的約束極值問題,即: 由于有4個(gè)變量,直接求出解的具體形式非常困難,故在求解過程中,不妨假設(shè)所求問題總有解。求解過程如下: (10) 根據(jù)廣義拉格朗日乘子法的最優(yōu)解求解方法,可得:D2-D1=0,D3-D2=0,D4-D3=0,D1>0,即當(dāng)0 為了找到L(D1,D2,D3,D4,λ1,λ2,λ3)的最小值,采取求偏導(dǎo)的方法,對L(D1,D2,D3,D4,λ1,λ2,λ3)分別對D1,D2,D3,D4,λ1,λ2,λ3求偏導(dǎo)并令其導(dǎo)數(shù)等于零,即可得該模型的最優(yōu)解。 由式(9)可知: (11) 式中:α,α1,α2均為常數(shù)且大于0;B為工程工期;A為工程成本主要組成因素及工程規(guī)模。 總之,以數(shù)學(xué)建模作為基本指導(dǎo)思想和方法論,從多元素集成成本管理系統(tǒng)的角度,充分考慮工程項(xiàng)目成本的主要成分,通過對項(xiàng)目成本與進(jìn)度的關(guān)系進(jìn)行適當(dāng)簡化,結(jié)合廣義拉格朗日乘子法和反證法,得出該模型的最優(yōu)解。 哈爾濱市某市政工程維修一條道路。該項(xiàng)目合同總價(jià)為105萬元,工期約為50天,并根據(jù)該項(xiàng)目的施工特點(diǎn)編制施工組織設(shè)計(jì),該工程的其他直接費(fèi)的取費(fèi)率為0.15,規(guī)費(fèi)費(fèi)率為0.041,企業(yè)管理費(fèi)費(fèi)率為0.059。 為了簡化計(jì)算該項(xiàng)目的最佳工期,根據(jù)哈爾濱市2005年1月1日到2013年12月31日期間開工的工程項(xiàng)目成本核算數(shù)據(jù),經(jīng)過篩選后,剔除不完整數(shù)據(jù),如尚在施工未結(jié)算的數(shù)據(jù)及拆遷類數(shù)據(jù)等。由于代表宏觀經(jīng)濟(jì)因素的PPI同樣會影響施工總成本及工程項(xiàng)目的材料費(fèi),因此用PPI將施工總成本與工程規(guī)模分別標(biāo)準(zhǔn)化,以剔除PPI的干擾,作為施工總成本的替代變量。利用SPSS軟件對數(shù)據(jù)進(jìn)行多元線性回歸檢驗(yàn),得出α1=28.932;α2=2 104 277.216;通過Matlab對式(2)進(jìn)行擬合,得出α=5.321×109。 將項(xiàng)目成本的主要組成費(fèi)用換算為均值。①人工費(fèi)。810個(gè)工作日,其日平均工資為120元/日。②材料費(fèi)。材料的耗用量為110噸,其材料平均價(jià)格為5 000元/噸。③機(jī)械費(fèi)。臺班數(shù)為150個(gè)工作日,其日平均臺班工資為1 300元/日。為了簡化計(jì)算該項(xiàng)目的最佳工期,根據(jù)式(4),使用一個(gè)通用的語言規(guī)則轉(zhuǎn)移數(shù)據(jù),梯形模糊數(shù)為:“比X更大或更小”=(0.9X,0.95X,1.05X,1.1X)?!按蠹sX”=(0.95X,X,X,1.05X)。 根據(jù)工程項(xiàng)目的特性對該項(xiàng)目的人工費(fèi)、材料費(fèi)、機(jī)械費(fèi)、間接費(fèi)進(jìn)行評判得出:人工費(fèi)的不確定性為“大約X”=(9.234,9.72,9.72,10.206);材料費(fèi)的不確定性為“比X更大或者更小”=(49.50,52.25,57.75,60.50);機(jī)械費(fèi)的不確定性為“比X更大或者更小”=(17.550,18.525,20.475,21.450)。根據(jù)式(8)得出A=37.689;由式(11)可知該工程的最佳工期為44.832天。 在實(shí)際工程施工環(huán)境中,施工工期、材料價(jià)格、材料消耗量及實(shí)際工程量的不確定性構(gòu)成了工程成本的不確定性,從而導(dǎo)致工程成本的預(yù)算與工程實(shí)際成本會出現(xiàn)較大的偏差。通過基于模糊數(shù)學(xué)理論對工程成本分析,根據(jù)道路工程成本、工程規(guī)模和工程工期的關(guān)系,建立了以工程成本為目標(biāo)的理論數(shù)學(xué)模型,并將該理論數(shù)學(xué)模型與工程成本基本函數(shù)關(guān)系相結(jié)合,將其因素模糊化,得出最終工程成本模型。在模糊化的工程成本模型基礎(chǔ)上,依據(jù)廣義模糊數(shù)的運(yùn)算與廣義Lagrange乘子法,得出在已知的工程規(guī)模下工程項(xiàng)目工期的最優(yōu)解。 [1] 暴克琦.工程規(guī)模、工期與施工總成本:基于某大型國有施工企業(yè)工程項(xiàng)目成本數(shù)據(jù)的實(shí)證研究[D].成都:西南財(cái)經(jīng)大學(xué),2012. [2] 黃松顯,黃軍周.不確定因素和人為因素對成本控制的影響與對策[J].中州煤炭,2010(10):141-142. [3] 雄鷹,匡亞萍.基于蟻群算法的施工項(xiàng)目工期-成本優(yōu)化[J].系統(tǒng)工程理論與實(shí)踐,2007(3):105-111. [4]KENZOK,NOBUYUKIN.EfficientMonteCarlosimulationmethodofGERT-typenetworkforprojectmanagement[J].Computers&IndustrialEngineering,2002,42(2-4):521-531. [5]LUISFA,DAVIDBA,ANGELIQUESH,etal.Riskplanningandmanagementforthepanamacanalexpansionprogram[J].JournalofConstructionEngineeringandManagement,2011,137(10):111-113. [6] 羅桂平.工程項(xiàng)目成本構(gòu)成與影響因素研究:基于工程項(xiàng)目成本的數(shù)據(jù)挖掘[D].成都:西南財(cái)經(jīng)大學(xué),2011. [7]BADUAJG,SURESHN.Projectmanagementwithtime,costandqualityconsiderations[J].EuropeanJournalofOperationalResearch,1996,88(2):320-327. [8] EI-RAYES K, KANDIL A. Time-cost-quality trade-off analysis for highway construction[J]. Journal of Construction Engineering and Management,2005, 131 (4): 477-478. [9] HSIEH C H. Optimzation of fuzzy production inventory models[J].Information Sciences,2002,146(1):29-40. [10] MORTEZA B.Designing a control mechanism using earned value analysis:an application to production environment[J].The International Journal of Advanced Manufacturing Technology,2010 (49): 419-429. GAO He:Assoc. Prof.; School of Civil Engineering, Northeast Forestry University, Harbin 150040, China. [編輯:王志全] Project Cost and Time Limit Based on Fuzzy Mathematics Theory GAOHe,CHAIShichao In order to improve the control on engineering cost of project in uncertain environment, the project cost was analyzed; and an optimization model was established based on the relationship of project cost and project scale and construction period while engineering cost was taken as the goal. The uncertainty conditions included labor fee, materials fee and machinery fee. The fuzzy mathematics theory and method, combining with trapezoid fuzzy number analysis method were used in the model. The generalized Lagrange multiplier method was combined with reduction to absurdity to obtain the solving process of optimal solution. It provides cost control reference for construction enterprise. engineering cost; uncertainty; trapezoidal fuzzy number; generalized Lagrange multiplier method 2015-05-14. 高賀(1974-),男,黑龍江龍江人,東北林業(yè)大學(xué)土木工程學(xué)院副教授. 黑龍江省交通廳科研基金資助項(xiàng)目(MJ20110034). 2095-3852(2015)06-0786-04 A TU723 10.3963/j.issn.2095-3852.2015.06.0273 模型的求解
4 算例驗(yàn)證
5 結(jié)論