This research includes the following technical limitations: infrastructure restrictions on the maximum permitted train speed; restrictions on the maximum number of trains on railway lines; limitations on the available number of rolling stock for train service; restrictions for implementation of the rehabilitation of the sections of the railway lines, according to the operational program transport. From the software point of view, there are no limitations. It can work with hundreds of alternatives and criteria.
4.1. Determination of the Initial Decision Matrix. Performing the SIMUS Procedure
Table 1 presents the initial decision matrix (in the case of the Bulgarian railway network) for the analysis and selection of a suitable policy for railway operators. It is formed with alternatives in columns and criteria in rows.
Alternative A1 presents the current situation and includes fast and accelerated fast passenger trains. The accelerated fast trains have mandatory seat reservations and serve major cities and transport nodes. The trains are composed of wagons. The average technical speed of movement of passenger trains is one of the lowest in Europe. At a design speed of 120–130 km/h, the movement of trains is achieved at 75–80 km/h, and in certain areas, it is limited to 40–60 km/h in order to ensure traffic safety [
44]. A significant part of the rolling stock does not meet European standards regarding comfort and quality. The existing structure of trains by type could be improved through the implementation of new types of trains.
Alternative A2 offers three categories of trains—fast trains, accelerated fast trains and direct express trains. The direct express trains are intercity trains that have also mandatory seat reservations and serve major transport and administrative centres. They have a lower number of stops in comparison to accelerated fast trains. Express trains and accelerated fast trains in alternatives A2 and A3 are composed of novel EMUs. As Alternatives A2 and A3 include fewer stops the value of criterion W6 decreases. The number of trains for alternatives A2 and A3 increases due to the replacement of old rolling stock with new ones. The new electric multiple unit trains (EMU) are characterized by lower electricity consumption compared to trains composed of wagons and locomotives due to their lower mass. The inclusion in the scheme of transport of a new generation of trains leads to a reduction of electricity costs and, accordingly, of carbon dioxide emissions for the production of electricity by the power plants. The increase in electricity consumption for alternative A3 is due to the movement of trains with increased maximum speed 100–120 km in the direction Sofia–Plovdiv–Burgas (this is part of the core TEN-T network). The introduction of novel EMUs leads to a reduction of direct operating costs by 3.44% for A2, and upon completion of the rehabilitation of the railway infrastructure by 1.25% for A3 (
Table 1—criterion S1).
Figure 2 shows the scheme of the itineraries of the passenger trains in the Bulgarian railway network.
The values of criteria W5, O1, O3, T1–T4 are determined using scale 0, 1. The values of criteria S2, S3, W1, W7, O2 and O3 are determined using a scale of 0, 1, 2, 3.
The next step in the model is to form the normalized matrix. The normalization could be performed in different ways. This study uses the Sum of All Values method.
Table 2 shows the normalized matrix, the type of actions for each criterion, the type of the operator for the restrictive conditions, the limits called “Right Hand Side” (RHS). The values of the RHS are obtained from the left normalized values and the type of action. In the case of minimum, RHS is equal to the minimum value of the row; in the case of maximum, the RHS value is equal to the maximum value of the row in the normalized matrix. The type of operator depends on the type of objective function. In the case of maximum, the operator is “≤”; in the case of minimum, the operator is “≥”.
The linear optimization models are performed by using the data in
Table 2. For example, the first optimization linear model is formed for the first objective Z1(criterion S1) as follows:
where:
represents the score of each alternative,
.
The restrictive conditions for the optimization model are formed by using the others rows of the ERM matrix. For example, for criterion S2, the restrictive condition is:
The restrictive conditions are formed successively using all other rows in the Normalized Sum Matrix. The final restrictive condition for the first optimization model is performed by criterion T4 (objective Z17) based on the data in the last row in the Normalized Sum Matrix, as follows:
Similar optimization linear models are performed for all other criteria.
Table 3 shows the results of the linear optimization of each criterion. Each row represents the values of the scores of the alternatives according to the optimization models by applying each criterion as an objective function. The Efficient Results Matrix (ERM) is formed. The values of the objective functions for each linear model (LHS) are given in the last column of the table.
4.2. Preliminary Analysis
Table 4 shows a preliminary analysis and makes a comparison between the desired values for each objective for SWOT criteria, which are identified as RHS, and the optimum values of the objective functions obtained through SIMUS (LHS). The optimal values for the criteria must not be confused with the lack of optimal values for the strategies. When there is equality of the LHS and RHS, it means that the objective is satisfied 100%.
Comparing the LHS and RHS columns, it can be seen that a series of criteria coincide with the same series of objectives, regarding the mathematical symbol. For instance, S3, calling for maximization, and with an LHS value of 0.50, must be less than or equal (≤) to RHS with also a value of 0.50. Other criteria, such as W2 show, also for maximization, that LHS is less or equal than objective Z5. Another criterion such as W4, aimed for minimization, has an LHS value of 0.32, whilst it is lower than the RHS values by 0.33.
However, it is worth noting that there is also a series of criteria that shows no coincidence between LHS and RHS. While a coincidence indicates that the objective is 100% satisfied, the opposite, a lack of coincidence, shows that it is not satisfied, either in deficit or in excess. The last case appears to be contradictory, but it is not.
Assume that a criterion such as investment requires maximization, and whose LHS is greater than RHS which is the maximum amount of funding available. Consequently, this would indicate that the value of the objective surpasses the maximum amount of money available, and thus, making the objective not feasible. This is the case of W2 (frequency of trains), whose value is 0.42, and is thus larger than the objective which is 0.34. The meaning is that at present, there are more frequencies than those recommended by the software when this criterion is compared with others, with which it interacts in both ways. For criterion W6 (average number of stops), it is the opposite, since it aims for a minimum number of 0.33 and the computation shows 0.32. Remember that all these values are normalized and then correspond to integer values in the initial decision matrix. Consequently, the minimum desired number of stops determined by studies and surveys is not satisfied. It can be seen how this procedure allows for the identification of facts that perhaps were not considered at first sight, and that are also a consequence of the interaction of the criteria. For instance, it shows that there could be a link between higher frequency (W2) and operating costs (W4). Perhaps an analysis by stakeholders may show the necessity to decrease frequency, to save in operating costs. It can be seen how this procedure may pinpoint aspects that need to be examined.
Table 4 is prepared on the ERM, where each objective combines with its option, Max or Min, and LHS and RHS values. Out of the 17 objectives, there are 6 significant ones. They are objectives Z5, Z7, Z9, Z14, Z15, Z16 and are in bold.
Objective Z5 refers to the frequency of trains, which naturally must be maximum. Observe that the computed value or LHS (0.42) is greater than the maximum value established (0.34) in the original table. This means that the actual frequency is higher than expected, which is obviously something positive for passengers, but possibly not so good for railway finances, and possibly related to a strain on operating costs.
Objective Z7, operating costs, which naturally must be minimized, shows that the railway is spending a little less (0.32) than the desired (0.33). This could be related to the diminishing number of passengers or less funding available, or too high frequencies as objective 5 shows. It is interesting to have a look at what statistics say. The reference [
45], shows a steady decline in the number of passengers, 2422 million passengers.km in 2006, to 1438 in 2017, with a slight recovery in 2018 and 2019 (1524 million passengers.km). That is, these statistics may confirm the lower expenditure shown in operating costs in
Table 1 due to a lower number of passengers transported. As a matter of fact, and from the financial point of view, [
46] shows a set of five financial indicators as IRR and NPV for investment in seven railway projects; four of them show negative returns, while IRR for capital, is the only one that shows a maximum value of 0.03. This lack of investment incentive may lead to not enough funding to have a reasonable operating cost level.
For objective Z9 (Average number of stops), interpreting the numbers from the table we find that apparently the minimum number of stops is not reached. This can be interpreted as positive or negative. Positive—as it means lower operating costs and increasing speed. Negative—as it may produce a reduction in the number of passengers.
For objective Z14 (Delayed purchase of rolling stock), it shows that RHS = 0; it then suggests that no long-term plans are in study to renovate the equipment, however, the computed value of 0.50 suggests that this must be considered as very important and in a high value. This is confirmed by the study in [
47] which says, about Bulgaria, that more than 75% of all transport funding should be allocated for: 1—Public urban transport system; 2—Integrated regional and suburban public transport systems 3—Railways (infrastructure and passenger rolling stock), 4—Intermodal infrastructure for shifting freight from road to rail, 5—Bicycle lanes and paths, and 6—Traffic management systems. We believe that for the development of the sustainable transport system in the country and for reversing the current dangerous trends of the predominant use of road transport there is a need for a new approach in transport investment.
For objective Z15 (Increasing Road quality), Bulgaria has about 20,000 km of paved routes, which seems very high in a country of 7,300,000 inhabitants with only one major city—Sofia, with about 1,300,000 people, and several low medium-size cities of about 350,000 inhabitants each. The territory of Bulgaria is not very large; it is about 111,000 km2 and ranks in 16th place among 51 European countries. This fact means that its density is 266 people/km2, which is low, and decreasing. Thus, it would be worth examining the occupancy rate of the wagons which is linked with criterion W1 in the transportation market, and W4, operating cost. This objective is ambitious, and speaks very well of the road infrastructure whose value of 0.50 is higher than the RHS value of 0 (or no plans). Nevertheless, railways have the advantage of the larger quantity of passengers and freight they can transport in comparison with busses, from the point of view of weight (and thus road conservation), passengers’ comfort, fuel consumption and contamination.
Objective Z16 (Increasing speed of road transportation), with a value greater than RHS, suggests the advantage of the road system.
Objective Z17 (Decline in traffic), shows no plans since RHS = 0, and then that whatever plan comes is OK.
The results show that the areas that must be improved are related to service (speed, number of stops, travel time, directness, costs, etc.) as well as those areas where it is necessary to work together with the country road development program, due to the competition between the two systems in those zones where the road and the rails run in parallel. The case is interesting since Bulgaria has some highways with a maximum speed of 140 km/h, while new rolling stock can reach 160 km/h.
4.3. Determination of the Weights of Criteria
The SIMIUS method does not use the weights of criteria for ranking the alternatives. This method is based on Linear Programming (LP) and the use of the Simplex algorithm. Linear Programming does not use any type of weights, and for that reason, they are not needed in the SIMUS method. LP works following an iterative process until it finds the optimal solution (if it exists). LP uses a vector procedure to establish a ranking between criteria, and thus, for each iteration, every criterion gets a new significance. In this way, the criteria significance is dependent on the set of alternatives to evaluate. It is a similar concept to using Shannon entropy, to evaluate criteria weights. There is no subjectivity here, everything is objective, and then, it does not matter who or how many decision-makers intervene. The result is always the same, and based on the initial data.
In this study, the Efficient Results Matrix (ERM) is used to determine the weights of criteria with the aim to assess the impact of the SWOT criteria. For this purpose, the maximum of each row is determined.
Table 5 presents the weights of the criteria for the SWOT group determined by the proposed approach.
The global weights are calculated by dividing these values by the sum of the maximum of the rows. Then, the weights of the SWOT groups are calculated. The results from the different SWOT clusters clearly show that weaknesses and threats are prevalent in the Bulgarian rail network. Both SWOT factors show sharp quantitative differences with their counterparts, i.e., strength and opportunities, both with the same importance or weight of 0.180, as shown in
Table 5. However, this is also significant since it shows that both favorable conditions exist in the network. Consequently, there is room for improvement. If the criteria are considered, it can be seen that the weights in
Table 5 are all practically the same, about 0.060, and therefore, no discrimination by importance is possible. However, observe that it is shown that although clusters Weakness and Threats are dominant, their criteria have the same weights as those of Opportunities and Strength.
4.4. Ranking of the Alternatives
The SIMUS algorithm produces an efficient matrix, where each of its rows is a Pareto Efficient solution, and thus, optimal. From this matrix, which is a mapping of the original matrix, SIMUS obtains two results. The first, following the Weighted sum procedure, and the second, based on Outranking. Both results coincide in their rankings, and then they mutually check their results.
The results obtained in the Efficient Results Matrix are applied to start the procedure of ranking the alternatives. The normalization sum method was applied to the Efficient Results Matrix. The results of the normalization are given in the first part of
Table 6. The second part presents the steps in obtaining the ranking. First, the sum of the column is determined. Then, the participation factors, which indicate the number of satisfactions of each alternative by each objective, are determined. The normalization of the participation factor is carried out by dividing the number of criteria. The final score of the alternatives is calculated by multiplying the sum of the columns by the normalized participation factor.
As can be seen, the best strategy is A3 (some reconstruction of railway infrastructures and new rolling stock on some lines), with the highest score of 3.76, followed by A2 (new rolling stock on some lines), with a score of 2.71. Note that the status-quo strategy has a very low score of 0.43, and then, should be discarded, which is in line with stakeholders’ opinion.
The SIMUS Method uses the outranking approach to validate the results obtained using ERM ranking.
Table 7 presents the results of the outranking approach where a new matrix is formed, called Project Dominance Matrix (PDM), that determines the best alternative. The number of columns and the rows in PDM is equal to the number of alternatives. The data of the ERM matrix are used. Starting from the highest value in the first row the difference between values in the same row of normalized ERM is calculated. The procedure is repeated with all the values. The net dominance is calculated as the difference between row sum and column sum. The alternatives are ranked according to the maximal value of the net dominance. It can be seen that the ranking formed using both procedures is the same.