Next Article in Journal
Lake Level Evolution of the Largest Freshwater Lake on the Mediterranean Islands through Drought Analysis and Machine Learning
Next Article in Special Issue
A Heuristic for the Two-Echelon Multi-Period Multi-Product Location–Inventory Problem with Partial Facility Closing and Reopening
Previous Article in Journal
Exploring the Relationship between Data Analytics Capability and Competitive Advantage: The Mediating Roles of Supply Chain Resilience and Organization Flexibility
Previous Article in Special Issue
Spatial Structure and Development Patterns of Urban Traffic Flow Network in Less Developed Areas: A Sustainable Development Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Bi-Level Programming Model for the Integrated Problem of Low Carbon Supplier Selection and Transportation

1
School of Business Administration, Shandong Technology and Business University, Yantai 264005, China
2
School of Electrical Engineering, Shandong University, Jinan 250061, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(16), 10446; https://doi.org/10.3390/su141610446
Submission received: 15 July 2022 / Revised: 16 August 2022 / Accepted: 18 August 2022 / Published: 22 August 2022

Abstract

:
In this paper, we investigate an integrated problem of low-carbon supplier selection and transportation. The supplier selection decision depends on the location and energy consumption level of batching plants at the manufacturing stage. Meanwhile, ready-mixed concrete is allocated and delivered to construction sites by concrete mixer trucks at the transportation stage. A bi-level programming model for the integrated problem is established. The bi-level optimization problem is transformed into a single-level problem by KKT (Karush–Kuhn–Tucker) optimality conditions. In order to validate the proposed model, a case study is conducted based on real-world problems. Experimental results show that the proposed method efficiently solves the integrated problem and the model can not only reduce carbon emissions but also optimize transportation time.

1. Introduction

With the rapid development of urbanization and industrialization, the construction of infrastructure has a severe impact on air pollution [1]. According to an analysis of the International Energy Agency (IEA), greenhouse gas emissions in buildings and the transport sector together represented in 2019 over half of global emissions, with transport comprising 27%, followed by buildings (25%). The increasing requirements of construction lead to a large demand for building materials. As the most widely used man-made building material, concrete consumption is approximately 1 ton for every human being each year, which produces tremendous carbon dioxide emissions [2].
Concrete CO2 emissions result from energy consumption not only in the production stage but also in the transport phase to the construction site [3]. Generally, concrete is made up of cement, water, and aggregates. Concrete batching and mixing is commonly carried out at concrete batching plants, located in different positions around the city. The mixing equipment is powered by electricity, which is the main source of emissions. Therefore, the energy consumption level of mixing equipment is the most crucial factor of emissions generated during the concrete production stage. After the production process, the concrete is transported to construction sites by a concrete mixer truck. The diesel fuel consumption of concrete transport is also a contributor to CO2 emissions [1]. In order to reduce concrete CO2 emissions, it is necessary to optimize the integrated problem of concrete production and transport.
There are two significant decisions: supplier selection and transportation planning, which are interdependent. The selection of concrete batching plants is a strategic decision that will have a long-term impact on the performance. This decision to choose the energy consumption level and facility location can influence not only the environment cost but also transportation times.
In the past few decades, facility location problems have attracted more and more attention. Farahani et al. conducted a comprehensive review on covering problems in facility locations [4]. Extending the classical facility location problems, Ortiz-Astorquiza et al. presented a survey of formulations, algorithms, and applications for multi-level facility location problems [5]. Puerto et al., provided an updated literature review of extensive facility location problems in networks [6]. A lot of research branches have been intensively investigated. Emirhüseyinoğlu and Ekici considered the facility location decision problem with supplier selection [7]. Arıoğlu et al. studied supplier selection and evaluation for the production of concrete in a large construction project [8]. Another research focus is the transportation or dispatching problem. Ready-mixed concrete is supplied to construction sites from batching plants along planned roads. Narayanan et al. introduced the delivery planning problem of ready-mixed concrete, which is formulated as a mixed-integer programming problem. Lagrangian relaxation is used to solve the concrete dispatching problem [9]. Maghrebi et al., addressed the large-scale ready-mixed concrete delivery problem and developed a column generation-based algorithm [10].
In the above-mentioned studies, facility location and transportation problems were considered separately. Additionally, it is of great importance to investigate the integrated problem of facility location and transportation [11]. There have been a few related works proposed from the supply chain perspective. Tezenji et al. presented an integrated model for supplier location-selection and order allocation with distance-based transportation cost [12]. Das et al., explored an integration problem between the facility location problem and transportation with the goal of determining optimum places for facilities and minimizing transportation cost [13]. Carlo et al., proposed an extended transportation location problem where the number of new facilities is a decision variable [14]. It is noted that the studies focused on the economic cost.
Recently, researchers have attempted to establish the integrated model from a low-carbon perspective. Leng et al. proposed a regional low-carbon location-routing problem with a bi-objective model. The total logistics cost and vehicle waiting time were both considered as objectives [15]. Hemmati and Pasandideh investigated a mixed-integer nonlinear program for the suppler location, supplier selection, and order allocation problem with green constraints. A bi-objective model was established by taking the total cost as the first objective and CO2 emissions as the second objective [16]. Kang et al. designed a novel bi-objective model for the supplier selection and inventory allocation planning problem with carbon trading [17]. Leng et al., developed a low-carbon location-routing problem for cold chain logistics. A novel bi-objective model was established with the goal of economic costs and environmental benefits [18].
However, the facility location problems are strategic decisions from the management point of view. Meanwhile, transportation planning issues are tactical decisions. Bi-level programming is applicable to solve these issues [19,20]. Wang et al., explored a bi-level programming approach to the location-routing problem with cargo splitting under low-carbon policies [21].
From the perspective of model building, the supplier selection decisions for the integrated problem are related to the number and location of plants, warehouses, or other facilities without determining the level of production energy consumption. Accordingly, we investigated the low-carbon integrated problem of concrete supplier selection and transportation planning in terms of selecting the energy consumption level of supplier equipment, determining the location of the concrete supplier, and planning the transportation of ready-mixed concrete. As mentioned above, a bi-level programming model [22] was established to deal with the strategic and tactical decision problems.
The main contributions of this paper are threefold. First, we propose an integrated problem of supplier selection and transportation from a low-carbon perspective. The CO2 emissions are generated not only in the production stage but also in the transportation stage. Second, a bi-level programming model is formulated to optimize both strategic and tactical decisions. The upper level is to minimize total CO2 emissions according to environmental criteria, whereas the lower-level model is designed to minimize the transportation time from the perspective of economy. Third, the bi-level optimization problem is transformed into a single-level problem through KKT conditions and dual theory.
The remainder of the paper is structured as follows. In Section 2, we describe the problem and present the mathematical model. In Section 3, we provide the solution method. In Section 4, a case study is conducted. In Section 5, we provide some concluding remarks.

2. Problem Description and Mathematical Model

2.1. Problem Description

The problem is described as follows. An integrated problem of low-carbon concrete supplier selection and transportation (IPLCCSST) is considered. At the strategic decision level, the construction enterprise needs to select several concrete batching plants for serving construction sites, which are located in different positions around the city. Generally, the concrete prices of different plants in a certain area are the same. The supplier selection decision depends on the location and energy consumption level of batching plants. At the tactical decision level, the batching plants are allocated to construction sites to satisfy their demands. After production, the allocated ready-mixed concrete can be delivered to corresponding construction sites by a concrete mixer truck. According to the concrete setting time, it is not allowed to deliver all concrete at once.
The following assumptions are considered:
  • The locations of the concrete batching plants and construction sites are known.
  • The capacities of the batching plants are known. The demand of each construction site is known, which should be supplied by one batching plant.
  • There are certain concrete mixer trucks initially available at each batching plant.
  • The distance and travelling time between plant and construction site is known. Due to the different speed limit, distance and time are not in direct proportion.
  • The truck unloading time is ignored.
  • We define volume as a unit term to measure the amount of concrete, which equals the capacity of one concrete mixer truck.
  • The energy consumption level of the batching plant is known.
  • The carbon dioxide emissions generated at the production stage are related to the energy consumption level and capacity.
  • The carbon dioxide emissions in the transportation phase are associated with distance.
As described above, the upper-level decision-maker develops the supplier selection strategy to optimize the environmental criteria by minimizing carbon dioxide emissions caused in the production and transportation stage from the perspective of society, whereas the low-level decision-maker formulates the transportation schedule to improve delivery efficiency by minimizing concrete delivery time.

2.2. Evaluation of Carbon Dioxide Emissions

The carbon dioxide emissions are generated in both the concrete production and transportation stages. The evaluation of CO2 emissions is presented as follows:
(1)
Carbon Dioxide Emissions from Concrete Production
In the concrete production stage, the energy consumption is based on the energy consumption level and the capacity of the production batching plants. Hence, the amount of energy consumption in the production stage can be calculated by the following equation [23]:
EC = i E L i × z i
In this equation, EC represents the energy consumption of the concrete batching plant, whose units are standard coal equivalent corresponding to the energy produced by the combustion of 1 kg of coal (in kgce). E L i denotes the energy consumption level of concrete batching plant i. z i is amount of concrete produced at batching plant i.
Hence, the CO2 emissions in the production stage can be calculated by the following equation [24]:
E CO 2 p = E C × E F p
where EC is the total energy consumption of the concrete batching plant, and EFp is the CO2 emission factor in the production stage.
(2)
Carbon Dioxide Emissions from Concrete Transportation
In the transportation phase, the CO2 emissions are generated when the ready-mixed concrete is transported to construction sites by a concrete mixer truck. Due to the special type of vehicle, the following calculation equation is adopted [24]:
E CO 2 t = d i j × a c × E F t
where dij is the distance between concrete batching plant i and construction site j, ac is the average fuel consumption of the concrete mixer truck, and EFt is the CO2 emission factor in the transportation stage.

2.3. The Bi-Level Programming Model for IPLCCSST

To develop the model, notations are defined as follows:
  • Sets:
  • I: set of candidate concrete batching plants
  • J: set of construction sites
  • Parameters:
  • dij: distance between concrete batching plant i and construction site j
  • p: number of selected concrete batching plants
  • dj: demand of construction site j
  • si: capacity of concrete batching plant i
  • c: load capacity of the concrete mixer truck
  • tij: time of a shipment of concrete transported from batching plant i to construction site j by concrete mixer truck
  • Decision variables:
  • xi: 1 if candidate batching plant i is selected to supply concrete
  • 0 otherwise
  • yij: amount of shipments of concrete transported from batching plant i to construction site j by concrete mixer truck
  • zi: amount of shipments of concrete supplied from batching plant i
With the notations above, we can formulate a bi-level programming model for the integrated problem of low-carbon concrete supplier selection and transportation (BIIPLCCSST) as follows:
The upper-level model (UBIIPLCCSST):
Objective function:
min F 1 = i E L i × z i × c × E F p + i j y i j × c × d i j × a c × E F t
subject to:
i x i p
i z i = j d j
z i s i
x i { 0 , 1 }
z i 0 , z i
Equation (4) is the objective function of the upper model, which minimizes the total CO2 emissions. The first term indicates the CO2 emissions generated in the concrete production stage. The second term implies the CO2 emissions in the transportation phase. Constraint (5) limits the number of selected concrete batching plants to p. Constraint (6) ensures that the sum of transportation volumes from all batching plants is equal to the total demand of all constructions site. Constraint (7) enforces that concrete batching plant capacities are not violated. Constraint (8) is a binary restriction. Constraint (9) defines the domains of the decision variables.
The lower-level model (LBIIPLCCSST):
Objective function:
min F 2 = i j t i j y i j
subject to:
i y i j = d j , j = 1 , 2 , n
j y i j = z i , i = 1 , 2 , m
y i j M 1 x i , i = 1 , 2 , m , j = 1 , 2 , n
y i j 0 , y i j
Equation (10) is the objective function of lower-level model, which minimizes the total transportation time. Constraint (11) ensures that for each construction site, the total transportation volumes from all batching plants are equal to its demand, whereas constraint (12) guarantee that the total transportation volumes from each batching plant are equal to the supply of the plant. Constraint (13) imposes that construction sites can be served only from selected batching plants where M1 represents a large integer number. Constraint (14) is the variable domain constraints.

3. Solution Method

In order to solve the lower-level model, we relaxed the lower-level constraint (14) to the following:
y i j 0
The problem (10)–(13) and (15) is denoted as relaxation-LBIIPLCCSST, which is a relaxation of LBIIPLCCSST. Let A1, A2 be the coefficient matrix of constraints (10)–(12), respectively. Let b1, b2 denote the right-hand side of constraints (10)–(12), respectively. Then, we present some results, which will be used in subsequent sections.
Lemma 1.
The polyhedron {constraints (11)–(13) and (15)} is integral.
Proof. 
Recall that A1 is totally unimodular, and A2 is I. We can obtain that [ A 1 A 2 ] is total unimodular [25]. As mentioned above, b1 and b2 are integral vectors. Due to Hoffman and Kruskal’s theorem, it is easy to yield the result. □
Proposition 1.
The relaxation-LBIIPLCCSST and its duality equation have integral optimum solutions.
Proof. 
Note that c is an integral vector. Direct from the fundamental properties of total unimodularity and integer linear programming [26], we can obtain the property of LBIIPLCCSST. □
Thereby, we provide a straightforward property for LBIIPLCCSST. The proof is omitted.
Corollary 1.
The LBIIPLCCSST can be optimally solved by the relaxation-LBIIPLCCSST.
By Corollary 1, the lower-level model can be optimally solved by its relaxation problem. We obtain the dual problem of the relaxation-LBIIPLCCSST:
max i z i u i + j d j v j + i j M 1 x i r i j
subject to:
u i + v j + r i j t i j
r i j 0
u i , f r e e
v j , f r e e
As the LBIIPLCCSST problem is linear, the bi-level optimization problem can be transformed into a single-level problem using the Karush–Kuhn–Tucker (KKT) optimality conditions. The equivalent single-level problem is formulated as follows.
min i E L i * z i * E F p + i j y i j * d i j * a c * E F t
subject to:
i x i p
i z i = j d j
z i s i
i y i j = d j
j y i j = z i
y i j M 1 x i
u i + v j + r i j t i j
r i j ( y i j M 1 x i ) = 0
x i { 0 , 1 }
z i 0 , z i
y i j 0
r i j 0
u i , f r e e
v j , f r e e
Equation (21) is the objective function of the equivalent single-level model. Constraints (22)–(24), (30) and (31) are the constraints of the upper-level problem. Constraints (25)–(27) and (32) are the primal constraints of the relaxed lower-level problem, whereas constraints (28) and (33)–(35) are the dual constraints of the relaxed lower-level problem. Constraint (29) is the complementary slackness conditions, which guarantee the optimality of the relaxed lower-level problem.
However, the complementary slackness conditions are non-linear. In order to linearize the complementary conditions, the big-M reformulation is used as follows [27]: Add a binary variable δ i j to the problem and append two additional constraints instead of the complementary conditions. We constructed the following big-M constraints:
r i j M 2 δ i j
y i j M 1 x i M 2 ( 1 δ i j )
where M 2 is a large positive constant. Then the resulting linearized single-level formulation is obtained as follows:
min i E L i * z i * E F p + i j y i j * d i j * a c * E F t
subject to:
i x i p
i z i = j d j
z i s i
i y i j = d j
j y i j = z i
y i j M 1 x i
u i + v j + r i j t i j
r i j M 2 δ i j
y i j M 1 x i M 2 ( 1 δ i j )
x i { 0 , 1 }
δ i j { 0 , 1 }
z i 0 , z i
y i j 0
r i j 0
u i , f r e e
v j , f r e e

4. Case Study

In this section, we present an application of the bi-level programming model. Firstly, we conducted a case study of a subway station construction project based on real-world problems. Then, we verified the effectiveness of the proposed solution method. The experimental results were obtained by the method. Finally, we analyzed the results, and some managerial implications are provided.

4.1. Description of Subway Station Construction Project

In order to illustrate and validate the proposed bi-level programming model, we considered a case study based on real-world problems in Yantai city, Shandong Province, China. The subway stations are planned to be built. Concrete batching plants need to be selected and concrete is allocated and transported to subway stations under construction. Taking the High-Tech Industrial Development Zone as an example, there are six concrete batching plants and 14 planned stations. The locations of the candidate concrete batching plants and construction sites are presented in Figure 1.
To conduct the test instance of the problem, detailed information is given as follows. There are three energy consumption levels of the concrete batching plants: 0.3 kgce/m3, 0.7 kgce/m3, and 1.1 kgce/m3. The number of selected concrete batching plants p = 2 L/km. A concrete mixer truck of 8 cubic meters is used to transport concrete from batching plants to stations under construction. A concrete mixer truck’s speed is 40 km/h with an average fuel consumption of 0.37 L/km. The demand of a station under construction is 500 shipments. The CO2 emission factor in the production stage is EFp = 2.6604 kg CO2/kgce, whereas the CO2 emission factor in the transportation stage is EFt = 3.1212 kg CO2/kg. The energy consumption level of the concrete batching plants is EL = [1.1, 0.7, 0.7, 0.3, 0.7, 0.3]. The distances from the candidate concrete batching plants to the construction sites are shown in Table 1.

4.2. Experimental Results

The case study was solved by the proposed method, which was conducted with Matlab and Gurobi on an Intel Core Quad PC with 3.3 GHz CPU and 8.00 GB RAM. The results were obtained within seconds. The optimal solution of the case is stated as follows:
x 3 = x 4 = 1 ; x 1 = x 2 = x 5 = x 6 = 0 ;   z 3 = z 4 = 3500 ; z 1 = z 2 = z 5 = z 6 = 0 . It is noted that Plant 3 and Plant 4 were selected for producing concrete with a supply amount of 3500 shipments. The detailed transportation plan is listed in Table 2.
The optimal objective value of the upper-level model was 301,348.76 kg, which is the minimum CO2 emissions for this case. The optimal objective value of the lower-level model was 613.8750 h, which implies the minimum total transportation time.
In order to analyze the impact of the number of selected concrete batching plants and concrete demand, a series of computational experiments were conducted. Consider that the number of selected concrete batching plants was 1, 2, 3, 4, 5, 6, respectively and the demands were 100, 200, 300, 400, 500, 600, 700, 800, 900, and 1000 shipments respectively. The results are shown in Table 3 and Table 4.

4.3. Discussion

Based on the results, the following can be obtained:
(1)
As seen in this model, the upper-level decision-maker (leader) should select concrete batching plants to minimize total CO2 emissions, taking into account the lower-level decision-maker’s plan. The lower-level decision-maker (follower) reacts to the leader’s action, and then carries out the transportation planning depending on the leader’s decision. Therefore, it is suitable to solve the low-carbon integrated problem of supplier selection and transportation by bi-level programming.
(2)
From the optimal solution, it can be seen that Plant 3 and Plant 4 were selected for producing concrete. Unlike the separated model, the proposed integrated model obtained a global optimal solution. Generally speaking, the batching plant with the lowest energy consumption level should be selected in the production model by optimizing the carbon emissions, and the closest demand sites should be allocated to the facility in the transportation model by minimizing transportation time. In this model, the solution was achieved by jointly optimizing both production and transportation.
(3)
According to the experimental results, we saw that the optimal objective value of the upper-level model was 301,348.76, whereas the optimal objective value of the lower-level model was 613.8750. By implication, that is not only a game between leader and follower but also a global optimum of the problem.
(4)
Table 3 and Figure 2 show that with increases in the number of selected batching plants, the total CO2 emissions and transportation time decreased at the first stage. When the number of selected batching plants was increased to four, the total CO2 emissions and transportation time did not decrease. This implies that it is feasible to adjust CO2 emissions and transportation time by increasing the number of selected batching plants at first. Up to the boundary, it will not be improved.
(5)
As shown in Table 4, it is notable that the demands of the subway stations under construction affected the objective values directly. The total CO2 emissions and transportation time had a linear correlation with the demands.

5. Conclusions

In order to achieve sustainable development, a low-carbon integrated problem of supplier selection and transportation planning is presented. From a low-carbon perspective, the evaluation of carbon dioxide emissions from the concrete production and concrete transportation stage are provided. From an integrated perspective, a bi-level programming model is formulated. A test instance based on real-world problems in Yantai city is conducted. According to experimental results, the proposed method achieved the optimal solution by optimizing both production and transportation. Sensitivity analysis on the number of selected concrete batching plants and demands is presented.
There are more complex constraints in real world, i.e., traffic congestion and uncertainty. For future research directions, it would be interesting to extend the integrated model with traffic congestion or uncertainty. Moreover, another future research direction is to develop more suitable algorithms for complex problems.

Author Contributions

Conceptualization, H.Z.; methodology, H.Z. and C.L.; software, H.Z. and C.L.; validation, H.Z., C.L. and Y.S.; formal analysis, H.Z. and C.L.; investigation, H.Z. and Y.S.; resources, H.Z. and Y.S.; data curation, H.Z. and C.L.; writing—original draft preparation, H.Z.; writing—review and editing, H.Z., C.L. and Y.S.; visualization, H.Z. and C.L.; supervision, H.Z.; project administration, H.Z.; funding acquisition, H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Humanities and Social Sciences Foundation of the Ministry of Education in China (grant number 17YJC630238).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the editor and the reviewers for their constructive suggestions and kind help.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Behera, M.; Bhattacharyya, S.K.; Minocha, A.K. Recycled aggregate from C&D waste & its use in concrete—A breakthrough towards sustainability in construction sector: A review. Constr. Build. Mater. 2014, 68, 501–516. [Google Scholar]
  2. Nazari, A. Handbook of Low Carbon Concrete; Elsevier: Oxford, UK, 2017; pp. 389–412. [Google Scholar]
  3. Goswein, V.; Gonalves, A.B.; Silvestre, J.D.; Freire, F.; Habert, G.; Kurda, R. Transportation matters—Does it? GIS-based comparative environmental assessment of concrete mixes with cement, fly ash, natural and recycled aggregates. Resour. Conserv. Recycl. 2018, 137, 1–10. [Google Scholar] [CrossRef]
  4. Farahani, R.Z.; Asgari, N.; Heidari, N.; Hosseininia, M.; Goh, M. Covering problems in facility location: A review. Comput. Ind. Eng. 2012, 62, 368–407. [Google Scholar] [CrossRef]
  5. Ortiz-Astorquiza, C.; Contrerasn, I.; Laporte, G. Multi-level facility location problems. Eur. J. Oper. Res. 2018, 267, 791–805. [Google Scholar] [CrossRef] [Green Version]
  6. Puerto, J.; Ricca, F.; Scozzari, A. Extensive facility location problems on networks: An updated review. TOP 2018, 26, 187–226. [Google Scholar] [CrossRef]
  7. Emirhüseyinoğlu, G.; Ekici, A. Dynamic facility location with supplier selection under quantity discount. Comput. Ind. Eng. 2019, 134, 64–74. [Google Scholar] [CrossRef]
  8. Arıoğlu, M.; Sarkis, J.; Dhavale, D.D. Selection of suppliers using Bayesian estimators: A case of concrete ring suppliers to Eurasia Tunnel of Turkey. Int. J. Prod. Res. 2020, 5, 1–12. [Google Scholar] [CrossRef]
  9. Narayanan, P.K.; Rey, D.; Maghrebi, M.; Waller, S.T. Using Lagrangian Relaxation to Solve Ready Mixed Concrete Dispatching Problems. Transp. Res. Rec. 2015, 2498, 84–90. [Google Scholar] [CrossRef]
  10. Maghrebi, M.; Periaraj, V.; Waller, S.T.; Sammut, C. Column Generation-Based Approach for Solving Large-Scale Ready Mixed Concrete Delivery Dispatching Problems. Comput. Aided Civ. Infrastruct. Eng. 2016, 31, 145–159. [Google Scholar] [CrossRef]
  11. Yan, S.; Wang, W.; Chang, G.; Lin, H.-C. Effective ready mixed concrete supply adjustments with inoperative mixers under stochastic travel times: Transportation Letters. Transp. Lett. 2016, 8, 286–300. [Google Scholar] [CrossRef]
  12. Tezenji, F.R.; Mohammadi, M.; Pasandideh, S.; Koupaei, M.N. An integrated model for supplier location-selection and order allocation under capacity constraints in an uncertain environment. Sci. Iran. Trans. E: Ind. Eng. 2016, 23, 3009–3025. [Google Scholar]
  13. Das, S.K.; Roy, S.K.; Weber, G.W. Heuristic approaches for solid transportation-p-facility location problem. Cent. Eur. J. Oper. Res. 2020, 28, 939–961. [Google Scholar] [CrossRef]
  14. Carlo, H.J.; David, V.; Salvat, G. Transportation-Location Problem with Unknown Number of Facilities. Comput. Ind. Eng. 2017, 112, 212–220. [Google Scholar] [CrossRef]
  15. Leng, L.; Zhao, Y.; Wang, Z.; Zhang, J.; Wang, W.; Zhang, C. A Novel Hyper-Heuristic for the Biobjective Regional Low-Carbon Location-Routing Problem with Multiple Constraints. Sustainability 2019, 11, 1596. [Google Scholar] [CrossRef] [Green Version]
  16. Hemmati, M.; Pasandideh, S. A bi-objective supplier location, supplier selection and order allocation problem with green constraints: Scenario-based approach. J. Ambient Intell. Humaniz. Comput. 2021, 12, 8205–8228. [Google Scholar] [CrossRef]
  17. Kang, K.; Pu, W.; Ma, Y.; Wang, X. Bi-objective inventory allocation planning problem with supplier selection and carbon trading under uncertainty. PLoS ONE 2018, 13, e0206282. [Google Scholar]
  18. Leng, L.; Zhang, C.; Zhao, Y.; Wang, W.; Zhang, J.; Li, G. Biobjective Low-carbon location-routing problem for cold chain logistics: Formulation and heuristic approaches. J. Clean. Prod. 2020, 273, 122801. [Google Scholar] [CrossRef]
  19. Jacobsen, S.K.; Madsen, O.B.G. A comparative study of heuristics for a two-level routing-location problem. Eur. J. Oper. Res. 1980, 5, 378–387. [Google Scholar] [CrossRef]
  20. Zhou, L.; Lin, Y.; Wang, X.; Zhou, F. Model and algorithm for bilevel multisized terminal location-routing problem for the last mile delivery. Int. Trans. Oper. Res. 2017, 26, 131–156. [Google Scholar] [CrossRef]
  21. Wang, C.; Peng, Z.; Xu, X. A Bi-Level Programming Approach to the Location-Routing Problem with Cargo Splitting under Low-Carbon Policies. Mathematics 2021, 9, 2325. [Google Scholar] [CrossRef]
  22. Vallejo, J.F.; Rodríguez, E.G.; Almaguer, F.J.; Gonzalez-Ramirez, R.G. A bi-level optimization model for aid distribution after the occurrence of a disaster. J. Clean. Prod. 2015, 105, 134–145. [Google Scholar] [CrossRef]
  23. The Norm of Energy Consumption per Unit Product of Ready-Mixed Concrete, GB36888-2018. Available online: http://std.samr.gov.cn/gb/search/gbDetailed?id=7AC3075CD354B4B1E05397BE0A0AC4DE (accessed on 19 November 2018).
  24. Paustian, K.; Ravindranath, N.H.; Amstel, A.V. 2006 IPCC Guidelines for National Greenhouse Gas Inventories; IPCC: Geneva, Switzerland, 2006. [Google Scholar]
  25. Nemhauser, G.; Wolsey, L. Integer and Combinatorial Optimization; Wiley: Hoboken, NJ, USA, 1988. [Google Scholar]
  26. Schrijver, A. Theory of Linear and Integer Programming; Wiley: Hoboken, NJ, USA, 1998. [Google Scholar]
  27. Fortuny-Amat, J.; Mccarl, B. A Representation and Economic Interpretation of a Two-Level Programming Problem. J. Oper. Res. Soc. 1981, 32, 783–792. [Google Scholar] [CrossRef]
Figure 1. Locations of the candidate concrete-batching plants and construction sites.
Figure 1. Locations of the candidate concrete-batching plants and construction sites.
Sustainability 14 10446 g001
Figure 2. Impact of the number of selected concrete batching plants.
Figure 2. Impact of the number of selected concrete batching plants.
Sustainability 14 10446 g002
Table 1. Distances from candidate concrete batching plants to stations (km).
Table 1. Distances from candidate concrete batching plants to stations (km).
Station 1Station 2Station 3Station 4Station 5Station 6Station 7Station 8Station 9Station 10Station 11Station 12Station 13Station 14
Plant 15.554.84.77.37.19.33.40.881.12.23.856
Plant 25.65.55.24.97.57.19.13.71.21.52.53.856
Plant 35.95.773.86.66.17.83.71.50.311.334.25.2
Plant 48.28.46.35.64.93.53.97.44.93.82.92.92.23.2
Plant 55.65.85.35.27.36.88.641.51.82.53.54.75.7
Plant 69.58.27.46.664.53.37.85.44.23.33.32.42.5
Table 2. Detailed transportation plan of the case study.
Table 2. Detailed transportation plan of the case study.
Station 1Station 2Station 3Station 4Station 5Station 6Station 7Station 8Station 9Station 10Station 11Station 12Station 13Station 14
Plant 35005000500000500500500500000
Plant 40050005005005000000500500500
Table 3. Results of different numbers of selected concrete batching plants.
Table 3. Results of different numbers of selected concrete batching plants.
p = 1p = 2p = 3p = 4p = 5p = 6
F1(t)359.27301.35295.34290.82290.82290.82
F2(h)851.25613.875597.625573.875573.875573.875
Table 4. Results of different demands.
Table 4. Results of different demands.
1002003004005006007008009001000
F1(t)60.27120.54180.81241.08301.35361.62421.89482.16542.43602.68
F2(h)122.775245.55368.325491.1613.875736.65859.425982.21105.001227.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhu, H.; Liu, C.; Song, Y. A Bi-Level Programming Model for the Integrated Problem of Low Carbon Supplier Selection and Transportation. Sustainability 2022, 14, 10446. https://doi.org/10.3390/su141610446

AMA Style

Zhu H, Liu C, Song Y. A Bi-Level Programming Model for the Integrated Problem of Low Carbon Supplier Selection and Transportation. Sustainability. 2022; 14(16):10446. https://doi.org/10.3390/su141610446

Chicago/Turabian Style

Zhu, Hongli, Congcong Liu, and Yongming Song. 2022. "A Bi-Level Programming Model for the Integrated Problem of Low Carbon Supplier Selection and Transportation" Sustainability 14, no. 16: 10446. https://doi.org/10.3390/su141610446

APA Style

Zhu, H., Liu, C., & Song, Y. (2022). A Bi-Level Programming Model for the Integrated Problem of Low Carbon Supplier Selection and Transportation. Sustainability, 14(16), 10446. https://doi.org/10.3390/su141610446

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop