The present invention relates to a field of power system engineering cost management, in particular, to a post-evaluation and risk management and control method of power transmission engineering cost.
Power transmission engineering is an important part of power network construction, and its amount of investment is usually large. There are many uncertainty factors in the power transmission engineering implementation process, such as design changes, schedule delays, equipment price change, etc. These uncertainty factors may cause fluctuation of power transmission engineering cost, resulting in cost and investment risks. Thus, it is necessary to consider the influence of uncertainty factors when evaluating the cost of the power transmission, analyze its possible risks, and study appropriate risk control measures, so as to predict and determine the costs in various stages of project reasonably, and thus control the total cost within a reasonable limit.
The total cost of power transmission engineering consists of a number of independent sub-item costs of unit project. Since it is complicated and difficult to analyze the effect of uncertainty factors on the total cost, and relatively simple to analyze the effect of sub-item cost, the evaluation and risk control on total cost can be converted to evaluation and risk control on sub-item costs. At present, usually the weight analysis method is used to evaluate the effect of different factors on power transmission engineering cost. Commonly-used methods include AHP (Analytic Hierarchy Process, AHP), principal component analysis and entropy weight method. These methods have the following drawbacks when used to evaluate the effect of each sub-item cost on total cost: usually AHP method requires expert knowledge and experiences, which is quite subjective; principal component analysis requires dimension reduction on various sub-item cost or factors, and the final principal component can not fully reflect the independent weight of single cost or single factor; entropy weight method is easily interfered by abnormal fluctuations, that is, when there are unreasonable data with great deviation in the samples, the results by weight calculation have great errors. Moreover, the three methods are unable to measure the risk interval of sub-item cost fluctuation.
So far, a variety of methods for risk measurement have been proposed, including Tail Conditional Expectation (Tail Conditional Expectation, TCE), VaR (Value at Risk, VaR) and Conditional Value at Risk (Conditional Value at Risk, CVaR), etc; TCE is similar to CVaR, and under a continuous distribution function, TCE is equivalent to CVaR. VaR is commonly used in the financial field, which means the maximum potential loss of a portfolio of investment products in a specific time period in the future under normal market conditions and a given confidence level. The definition of CVaR is the conditional mean that the loss exceeds VaR under a confidence level, which reflects the mean level of excess losses. Generally it is believed that tail loss measurement of VAR has non-sufficiency, while CVaR can reflect the potential risk of investment portfolio appropriately. Since the natural and market environment for power transmission engineering in the same region are the same, and some experiences have accumulated from historical engineering, when risk measurement of VaR and CVaR is used to actual power transmission engineering, its tail risk is usually caused by special reasons (such as major policy mistakes), which can be usually ignored; but if CVaR is used, which considers this tail risk, the actual risk is overstated.
The object of the invention is to overcome the above problems and provide a VAR-based post-evaluation and risk management and control method of power transmission engineering cost.
To achieve the object, the invention employs the following technical solutions:
A VAR-based post-evaluation and risk management and control method of power transmission engineering cost is provided. Through introducing the VAR theory into the evaluation and risk measurement of costs, a VAR-based post-evaluation and risk management and control method of power transmission engineering cost is proposed.
A VAR-based post-evaluation method of power transmission engineering cost, including:
(1) breaking down the power transmission engineering into a plurality of unit projects, including earthwork, foundation engineering, tower engineering, overhead line engineering and accessory engineering according to the contents of the power transmission engineering projects, namely five sub-item costs. Each unit project has a strong independence, to convert the evaluation of total cost into evaluation of sub-item costs, so as to have a detailed evaluation on the engineering cost, and identify the root cause of sub-item cost volatility based on the evaluation results.
(2) Calculating the weight of sub-item cost of each unit project based on VAR theory.
A weight evaluation model of unit project sub-item cost is established based on VaR theory.
The sub-item cost of the power transmission engineering is divided by unit project, when X represents unit project vector, vector YεRm (Rm represents m-dimensional real number space) represents uncertainty factor, each sub-item cost can be expressed as f(X,Y). The probability distribution function of f(X,Y) is calculated as follows:
Assuming that the JPDF of Y is p(Y), the probability of f(X,Y) that does not exceed a given critical value α (α represents a given sub-item cost level) for a definite X is:
where, X=[X1, X2, . . . , Xn], n is the number of unit project; φ(X,α) is the distribution function of sub-item cost.
It is usually difficult to directly obtain the JPDF of random variables affecting the sub-item cost, and thus it is unable to get the distribution function of sub-item cost through the above formula (1). However, since the changes in sub-item cost caused by random factors have features of volatility clustering, it is generally simulated by normal distribution, namely, obtaining the distribution curve of φ(X,α) by direct simulation based on historical data.
The corresponding VaR value of each sub-item cost f(X,Y) is represented by αβ(X) when the confidence that f(X,Y) does not exceed the critical value α is β, then αβ(X) can be represented by the following formula (2):
αβ(X)=min{αεR:φ(X,α)≧β} (2)
u(X) represents the mean of sub-item costs in φ(X,α) distribution curve; under a given confidence, the greater the deviation of VaR of sub-item cost from the mean, the greater the risk of sub-item cost fluctuations, thus, the overall volatility weight of sub-item cost can be described by the mean of deviation of VaR under a given confidence. Assuming that Xi represents an unit project, the mean of corresponding sub-item costs is represented by u(Xi), the corresponding VaR value under the confidence β is represented by αβ(Xi), assuming that θβ(Xi) represents weight, it can be represented by the ratio in the formula (3):
θβ(Xi)=αβ(Xi)/u(Xi) (3)
where, Xi represents a unit project; i=1, 2, . . . , n, n is the number of unit project;
θβ(Xi)
θβ(Xi) can be used to qualitatively assess the risk of fluctuation of each unit project sub-item cost, thus it is required to focus on the unit project with greater weight in the management and control process.
(3) Establishing a post-evaluation method for cost fluctuation according to the weight of sub-item cost of each unit project.
Using the weight result under confidence β as a reference, the post-evaluation of the level of risk of sub-item cost for a specific power transmission engineering is performed, assuming that c(Xi) represents a sub-item cost corresponding to a unit project Xi of a particular project; θ(Xi) represents the degree of deviation from the population mean, i.e. deviation coefficient, which is measured according to the ratio in the formula (4).
θ(Xi)=c(Xi)/u(Xi) (4)
Assuming that σ(Xi) represents the risk assessment score of Sub-item cost c(Xi), its value is measured by difference of overall volatility weight θβ(Xi) between θ(Xi) and sub-item cost under confidence β, namely:
σ(Xi)=θ(Xi)−θβ(Xi) (5)
The smaller σ(Xi) is, the smaller the risk of fluctuations of sub-item cost c(Xi) for the particular project.
A VAR-based risk management and control method of power transmission engineering, including:
(1) The sub-item cost of each unit project is different, and its influence on total cost is different. Here, an index of contribution degree is proposed to measure the effect of sub-item cost of each unit project on total cost, and determine the unit project that should focus on management and control according to the sequence of contribution degree.
The contribution degree consists of two parts: (1) overall volatility level of all sub-item costs, i.e. the weight as stated above; (2) the proportion of each sub-item cost among total cost. Firstly, if only analysis of the ratio, for a sub-item cost with high ratio but less fluctuation, the sub-item cost can be approximated as a constant, and its influence on overall cost volatility is less; secondly, if only analysis of the volatility level, for a sub-item cost with great volatility but less ratio, its influence on overall cost volatility is less. By comprehensively considering the weight and ratio, the contribution degree index can be calculated by the following formula:
k(Xi)=kp(Xi)×θβ(Xi) (6)
Where, kp(Xi) is the proportion of sub-item cost c(Xi) among total cost; i=1, 2, . . . , n, n is the number of unit projects.
(2) Screening the main sub-item cost according to the sequence of contribution degree, establishing an optimized stochastic linear programming model of sub-item cost.
Unit projects with great contribution can be identified according to the sequence of contribution degree, so as to focus on their risk management and control.
Many factors will influence the cost, including voltage level, the number of loops, wind speed, icing and other technical factors; besides, many social factors such as project management level, equipment development level, inflation rate, etc. The program design is based on technical factors; after formation of a specific design program, the designs may be changed due to the design, construction and quality problems during implementing power transmission projects. The cost of design change risk can be expressed as the product of general cost level multiplied by the rate of design changes according to the formula (7):
c
1=γ1×y0 (7)
Where, γ1 is the rate of design change risk, y0 is the general cost level; γ1 is random but its value should be within a range of [γ1−, γ1+], if exceeding the range, the engineering is feasible; γ1− and γ1+ are determined by historical data or experiences. The design change cost under a given risk rate can be used to guide the formation and comparison of design change program.
Duration management risks exist during the implementation of power transmission engineering. The project implementation cost mainly consists of two parts: direct cost and indirect cost. The direct cost refers to the sum of direct costs of all processes for the power transmission engineering plan, which include the costs of raw materials, machinery and equipment and labors for the processes. The indirect cost mainly includes the costs for management, supervision, inspection and coordination during the engineering project implementation, which is correlated to the project duration. The longer the duration, the higher the indirect cost; for an actual project, an approximate indirect cost within a unit time can be given. By this way, the cost for the duration management risk can be expressed by the formula (8):
c
2
=T×i
c (8)
where, T represents time exceeding the estimated duration, ic is the rate of indirect cost within unit time, which can be estimated by experts or calculated from the indirect cost in the budget statement by the duration. Assuming that γ2 represents the rate of duration management risk, γ2=T/T0, T0 is the estimated duration, then the formula (8) can be expressed as:
c
2=γ2×T0×ic (9)
The rate of duration management risk is random and varied, and its range of variation is limited by the longest duration. The range of values can be expressed as [γ2−, γ2+]. We can facilitate the scientific management of duration by calculating the duration management cost within a risk rate.
There are uncertainty factors of equipment price changes and labor cost changes in the span of construction duration. The two costs can be expressed by the formula (10):
c
3
=c
31×γ31+c32×γ32 (10)
Where, c31 represents the estimated cost of equipment, γ31 represents the rate of equipment price risk; c32 represents the estimated cost of labor, γ32 represents the rate of labor cost risk. γ31 and γ32 vary with the fluctuation of supply-demand relationship of equipments and labor markets. It can be believed that their value intervals can be predicted, represented by [γ31−, γ31+] and [γ32−, γ32+] respectively. The calculation of the equipments and labor costs within a given rate of risks can facilitate the procurement and employment negotiations, so as to control the cost within a reasonable risk interval.
Based on above work, a stochastic linear programming model can be established to control the cost within a reasonable range of risk:
y
β−≦(1+γ1)×y0+γ2×T0×ic+c31×γ31+c32×γ32≦yβ+ (11)
γ1−≦γ1≦γ1+ (12)
γ2−≦γ2≦γ2+ (13)
γ31−≦γ31≦γ31+ (14)
γ32−≦γ32≦γ32+ (15)
Where, β is confidence, yβ− and yβ+ represent the left boundary value and right boundary value of VaR at the confidence β; y0 represents a general cost level under the design scheme; γ1, γ2, γ31 and γ32 represent the rate of design change risk, rate of duration management risk, rate of equipment price risk and rate of labor cost risk respectively; [γi−, γi+] represents the interval of γi, i.e. the constraints of the optimization model.
Assuming that y=(1+γ)×y0+γ2×T0× ic+c31×γ31+c32× γ32, the target of the above optimization problem is to control y within the interval [yβ−, yβ+], the solution of optimization model is the risk intervals of γ1, γ2, γ31 and γ32, and the above problem is converted to obtain the combination of optimal solutions of γ1, γ2, γ31 and γ32 when y is the right boundary yβ+; here, the combination of optimal solution refers to the maximum of combination (γ1, γ2, γ31, γ32), that is, as long as the risk is controlled within the range of optimal combination, it can ensure that y is within [yβ−, yβ+], namely:
max(γ1,γ2,γ31,γ32) (16)
(1+γ1)×y0+γ2×T0×ic+c31×γ31+c32×γ32=yβ+ (17)
γ1−≦γ1≦γ1+ (18)
γ2−≦γ2≦γ2+ (19)
γ31−≦γ31≦γ31+ (20)
γ32−≦γ32≦γ32+ (21)
{circle around (3)} Random variables in the above optimization model are simulated using Monte Carlo method, to obtain a group of optimal solutions based on samples in each group, and for the combination of 50 groups of optimal solutions within the solution domain, their mean is used as the optimal solution.
The VAR-based post-evaluation method and risk management and control method of power transmission engineering cost provided in the invention can implement effective evaluation on the cost fluctuation risks; and through screening the key unit project by the contribution degree index, it can effectively reduce the risk control dimension, enhance risk management and control efficiency; and finally through establishing a stochastic programming model, it can achieve control over the cost within a risk interval.
The present invention is further described in combination with drawings and specific embodiments.
For the post-evaluation and risk management and control method of power transmission engineering cost, an engineering example in the same region is given to illustrate the effectiveness of the method.
Forty groups of historical power transmission engineering data in a region are selected for analysis. Normal distribution is employed to simulate all sub-item cost data. Based on VaR theory, the weights of unit projects are calculated at confidence levels 0.05, 0.1 and 0.15. In
As shown from
The standard confidence of the risk is set at 0.1 and a post-evaluation on the sub-item costs of two specific projects is performed, as shown in
The contribution degree index of each sub-item can be calculated in combination with the weights and ratio of sub-item costs at the standard confidence, as shown in
Taking tower engineering as an example, we perform calculation using the aforementioned optimization model. Due to diversity of common devices in the market, the overall price level of devices changes little. Assuming that γ31=0 and given γ32 ε[0, 0.2], we can calculate the certain linear optimization model to get the results as shown in
The specific embodiments described above, as preferred embodiments, are used to explain rather than limit the invention. Any modification, equivalent replacement and improvement made within the spirit and claims of the invention shall fall within the scope of protection of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
201610262139.2 | Apr 2016 | CN | national |