A stochastic programming approach for EOL electric vehicle batteries recovery network design under uncertain conditions

With the development of the electric vehicle industry, the number of power batteries has increased dramatically. Establishing a recycling EOL (end-of-life) battery network for secondary use is an effective way to solve resource shortage and environmental pollution. However, existing networks are challenging due to the high uncertainty of EOL batteries, e.g., quantity and quality, resulting in a low recycling rate of the recovery network. To fill this gap, this paper proposes a stochastic programming approach for recovery network design under uncertain conditions of EOL batteries. Firstly, a multi-objective model for battery recovery network is established, considering carbon emissions and economic benefits. Secondly, a stochastic programming approach is proposed to clarify the model. Subsequently, the genetic algorithm is employed to solve the proposed model. Finally, a recovery network case of Region T is given to verify the credibility and superiority of the proposed method. The results demonstrate that the proposed model reduces carbon emissions by 20 metric tons and increases overall economic benefits by 10 million yuan in Region T compared to the deterministic model. Furthermore, the two portions affecting the optimization results are also discussed to provide a reference for reducing carbon emissions and improving economic efficiency in recycling networks.


A stochastic programming approach for EOL electric vehicle batteries recovery network design under uncertain conditions
Wei Yan 1,2 , Xiao Wang 1* , Ying Liu 3 , Xu-mei Zhang 1 , Zhi-gang Jiang 2 & Lin Huang 1 With the development of the electric vehicle industry, the number of power batteries has increased dramatically.Establishing a recycling EOL (end-of-life) battery network for secondary use is an effective way to solve resource shortage and environmental pollution.However, existing networks are challenging due to the high uncertainty of EOL batteries, e.g., quantity and quality, resulting in a low recycling rate of the recovery network.To fill this gap, this paper proposes a stochastic programming approach for recovery network design under uncertain conditions of EOL batteries.Firstly, a multiobjective model for battery recovery network is established, considering carbon emissions and economic benefits.Secondly, a stochastic programming approach is proposed to clarify the model.Subsequently, the genetic algorithm is employed to solve the proposed model.Finally, a recovery network case of Region T is given to verify the credibility and superiority of the proposed method.The results demonstrate that the proposed model reduces carbon emissions by 20 metric tons and increases overall economic benefits by 10 million yuan in Region T compared to the deterministic model.Furthermore, the two portions affecting the optimization results are also discussed to provide a reference for reducing carbon emissions and improving economic efficiency in recycling networks.The quantity of EOL batteries flowed between battery processing center j and reuse center m Q jn

List of symbols
The quantity of EOL batteries flowed between battery processing center j and material recycling center n Q ij The quantity of EOL batteries flowed between battery collection center i and battery processing center j p 11 Unit economic benefits of energy storage over the life cycle p n Unit economic benefits of material recycling center n g Unit economic benefits of government subsidy for recycling of EVB d ij Distance between battery collection center i and battery processing center j d jm Distance between battery processing center j and reuse center m d jn Distance between battery processing center j and material recycling center n f j Fixed cost of building battery processing center j f m Fixed cost of building reuse center m f n Fixed cost of building material recycling center n c Unit transportation cost between centers c j The unit processing cost of battery processing center j www.nature.com/scientificreports/c m The unit processing cost of reuse center m c n The unit processing cost of material recycling center n p The unit price of battery collection center purchasing EOL batteries a Unit CO 2 emissions per distance generated during EOL battery transport b CO2 emissions from building a center e j Unit CO 2 emissions from EOL batteries processing at battery processing center j e m Unit CO 2 emissions from EOL batteries processing at reuse center m e n Unit CO 2 emissions from EOL batteries processing at material recycling center n β Electric emission factor for energy storage scenario h Carbon tax

MaxJ
The upper limit of processing capacity for battery processing center j MaxM The upper limit of processing capacity for reuse center m

MaxN
The upper limit of processing capacity for material recycling center n Decision variables y j ∈ {0, 1} If a battery processing center is established at the location j ∈ J is 1, otherwise 0 y n ∈ {0, 1} If a material recycling center is established at the location n ∈ N is 1, otherwise 0 y m ∈ {0, 1} If a reuse center is established at the location m ∈ M is 1, otherwise 0

Fuzzy variable Qi
The fuzzy quantity of EOL batteries recycled from battery collection center i μ The fuzzy percentage of EOL batteries flowed to reuse center m ω The fuzzy percentage of EOL batteries flowed to material recycling center n In response to the global energy crisis and the imperative to strengthen environmental protection, the National Development and Reform Commission (NDRC) formulated the "New Energy Automobile Production Admittance Management Rule" to promote the marketization of new energy vehicles 1 .As a result, new electric vehicles (EVs), serving as an alternative to traditional vehicles, are experiencing rapid development 2 .According to a report released by the International Energy Agency (IEA), global sales of EVs increase by 55% to over 10 million in 2022.And the share of EVs in total vehicle sales rises to 14%, making its highest historical proportion 3 .In such a case, the number of batteries, serving as the core components of EVs, is on the rise 4 .As outlined in China's 14th Five-Year Plan, notable enhancements have been made in the technical indicators of energy batteries, resulting in an extension of battery service life from 8-10 to 20 years 5 .Consequently, a substantial quantity of batteries is anticipated to reach their end-of-life (EOL) in the coming decade.It is predicted that the number of EOL electric vehicle batteries is expected to expand to around 7.08 million tons, and the EV batteries recycling market may even exceed 100 billion yuan by 2030 6 .Furthermore, EV batteries contain metals such as cobalt, nickel, manganese, and lithium, possessing significant recoverable value.If these metals are hydrolyzed, they will cause acidification of the ground and produce toxins and other harmful substances 7,8 .Failure to promptly dispose of EOL batteries leads to both resource wastage and environmental harm.Therefore, an effective recovery network of EOL batteries is the key enabler for the sustainable development of the EV industry.However, there are some unavoidable uncertainties in the battery recycling process.The battery recycling process involves several steps and statistics indicate that irregularities in the operation and management of these steps can cause uncertainties in the quantity and quality of batteries recycled.For instance, there is some risk of leakage or damage during logistics and transportation as well as a lack of standardization in the battery recycling and disposal processes.Simultaneously, there are problems of insufficient regulation and enforcement in some areas, such as illegal dumping or the re-entry of untreated batteries into the market.Additionally, the inability of customers to replace EV batteries on time can result in a diversity of battery capacity at the EOL stage.The uncertainties in the recycling process may affect the location and number of network nodes.A wellfunctioning recovery network is a crucial component of a used battery recycling strategy, as it reduces costs, minimizes environmental impact, maximizes profits, and enhances resource efficiency.Therefore, it is crucial to take into account the uncertainty involved in the process of recycling batteries at the end of their life cycle, and then optimize the network accordingly.
The objective of this paper is to address this gap through the design of a recycling network for EOL batteries.We examine the quality variations among EOL batteries and incorporate the secondary energy storage echelon utilization of batteries, along with the recycling of precious metal materials.Additionally, under the policies of "carbon neutrality" and "peak carbon dioxide emissions", and given the requirements of the global low-carbon strategy, we also fully consider the economic and environmental benefits of the entire logistics network.To be closer to the actual situation, we have also considered the uncertainty related to the quantity and quality of batteries in the recycling process.
Thus, this paper proposes a stochastic programming approach for recovery network design under uncertain conditions of EOL batteries.The study can effectively deal with the quantity and quality uncertainty of EOL batteries, thereby enhancing recycling efficiency.The primary contributions of this paper include: (1).This paper presents a multi-objective optimization model for EOL battery networks, taking into account uncertainties, which is not only a theoretical attempt but more importantly fits the real-world scenarios.
It not only facilitates resource conservation and environmental protection but also promotes the adoption and application of sustainable energy.
(2).This paper employs a blend of stochastic programming methods and triangular fuzzy numbers to address the uncertainties in both quantity and quality during the battery recycling process.Intelligent algorithms are employed to tackle this intricate problem.It not only expands the scope of the solution, but also improves the feasibility of solving the problem and its practical applicability.(3).The model is validated through a real case and the impact of uncertainties on the results is analyzed.In addition, the paper uses the actual distance of the nodes in the case identified in the AutoNavi to estimate the actual transportation costs.
The rest of this paper is organized as follows.Section "Review" reviews the relevant literature.Section "Problem statement" describes the battery recycling network and related assumptions and parameters.Section "Methodology" presents the constructed stochastic programming fuzzy model.And in this section, triangular fuzzy numbers are introduced to make the model clear and easy to solve with the help of a fuzzy chance-constrained programming model.Section "Case study" validates the model through a real-life example.The research is ended in Sect.6 with a conclusion and recommendations based on the analysis of the results.

Review Reverse logistics network design
The significance of reverse logistics extends beyond cost reduction; it encompasses enhancing resource utilization and environmental protection.In most of the literature studies, it aims to portray real life constraints in a more refined way by modeling the targets.For example, Kilic et al. 9 developed a cost-conscious reverse logistics network work for waste electrical and electronic equipment (WEEE) considering the constraints of the recycling rate of WEEE; Alumur et al. 10 proposed a profit maximizing modeling framework for reverse logistics network design and finally applied it to a realistic case of washing machines and tumble dryers in Germany; Sasikumar 11 developed a mixed-integer nonlinear programming model for Indian tires that considered maximizing the profit of the reverse logistics network; Momenitabar and Mattson et al. 12 proposed a multi-objective mixed-integer linear programming model and used it in a logistics network in the tire industry, which successfully reduced the total cost of ownership and promoted the sustainability of the supply chain; Al-Salem et al. 13 used an optimization model for minimizing the total cost of the reverse logistics network and treated the nonlinear objective as linear by segmentation method, and the results showed that a rational reverse logistics network can achieve cost savings.
While numerous studies exist on reverse logistics network designs for various objects, the number of studies focusing directly on EOL electric vehicle batteries is limited.In the early days, Yoda et al. 14 argued that the combustion of fossil fuels can have adverse effects on the environment.Consequently, the recycling of efficient and long-lasting batteries holds significant advantages for both environmental preservation and energy sustainability.And Beer et al. 15 explored the potential opportunities of using EOL electric vehicle batteries for secondary applications and concluded that centralized storage and reuse of EOL batteries could create significant economic value.Therefore, effective recycling of EOL batteries can bring economic as well as environmental benefits.More and more scholars have begun to focus on the topic of building a reverse logistics network for power batteries.Among them, Kannan et al. 16 established a multi-cycle and multi-product battery reverse logistics network model; Tadaros et al. 17 developed a mixed-integer programming model considering the total cost of the network and applied it to the recycling of used lithium-ion batteries in Sweden.However, from a circular economy and environmental point of view, the degree of environmental impact needs to be taken into account in the optimization objectives as well.Later, through a comparative exploration of greenhouse gas emissions from the cradle-togate process of ordinary internal combustion engines and EVs 18 and an estimation of greenhouse gas emissions from EV batteries 19 , it was found that carbon emissions generated during the life cycle of EVs or EV batteries are not negligible.Therefore, many articles indicate the extent of the environmental impact by calculating the carbon emissions generated during the recycling process.In exploring ways to reduce environmental impact, many scholars have found that effective recycling and remanufacturing strategies will reduce carbon emissions to some extent: Xiong et al. 20 found that recycling and remanufacturing EOL batteries could lower greenhouse gas emissions by 6.62%; Hao et al. 21learned from their study that recycling can reduce greenhouse gas emissions by about 10%.Therefore, Chen et al. 22 argued that carbon dioxide emissions should be prioritized in a closedloop supply chain, and thus proposed a model with a dual objective of total cost and carbon emissions; Wang et al. 6 constructed a recovery network considering remanufacturing strategies for batteries and built a model that considered both total costs and carbon emissions; Wang et al. 23 think that reverse logistics should not only consider the economic benefits but also focus on environmental protection.To this end, they have developed a reverse logistics network optimization model for waste batteries, with the aim of reducing logistics costs and CO2 emissions; M. Momenitabar et al. 24 designed a sustainable bioethanol supply chain network considering costs, environmental impacts, and employment opportunities, and finally used machine learning methods and meta-heuristic algorithms for solving the problem.Although carbon emissions have been a significant concern, many papers generally consider only the carbon emissions and costs during transportation and processing in the battery recycling network optimization process, and few of them also consider the echelon utilization of EOL batteries after recycling.Considering that batteries also generate carbon emissions as well as profitability issues in the reuse process, this paper considers the costs, profits, and carbon emissions in the battery recycling process from consumers to the battery echelon utilization process for secondary energy storage in building the model.

Uncertainty in reverse logistics network optimization
Compared to models of forward logistics, there is a higher level of uncertainty in the recycling process of the product 25 .In other words, uncertainty is an inherent characteristic of reverse logistics (RL).Uncertainty could be generally found in market demand, the capacity of the facility to process the products, the return quantity of www.nature.com/scientificreports/ the recycled products, the quality variability of the return flow, and the associated cost parameters 26 .In 2021, Wang et al. 27 investigated a reverse logistics network dedicated to household hazardous waste, considered a mixed integer deterministic and stochastic model, and set the objective function in the model to be an expectation value; Amin et al. 28 designed a closed-loop supply chain network for remanufacturing tires and solved the problem of uncertainty in the demand and return of the tires by using a decision tree approach; Momenitabar and Mattson et al. 12 proposed a fuzzy multi-objective hybrid linear programming model for the tire industry and solved the model using fuzzy robust optimization method as well as goal programming method.An extensive reading of the literature reveals that there are still relatively few articles on the impact of uncertainty in recycling on reverse logistics siting models in the EVB industry.Bao et al. 29 considered the total cost of a reverse logistics network for retired power batteries and developed a fuzzy programming model considering the uncertainty in the percentage of end-of-life, so as to determine the location of the network facilities; Jafari and Abharian 30 proposed a mathematical model for solving the closed-loop supply chain site selection problem in the automotive battery industry under uncertainty in recycling rate, facility capacity, and cost, and minimized the cost and impact on the environment by introducing trapezoidal fuzzy numbers.
For reverse logistics, uncertainty on the input side can complicate the structure of the network and create unmeasured cost issues.Certainly, there is always uncertainty in the quantity and quality of the batteries to be recycled.This paper explores recovery networks dealing with uncertainties in both the quality and quantity of recycled EOL batteries.And yet, fuzzy programming is highly subjective, requires iterative adjustments, and has many limitations in practice.And the results of robust optimization are too pessimistic and too rigid conservative.Meanwhile, Liang 31 states that triangular membership functions are prone to deviations from the central value and are relatively flexible.Hence, this paper chooses to convert the constructed model into a fuzzy chance-constrained programming model based on stochastic programming.Accordingly, the article uses the triangular membership functions to represent the fuzzy numbers in the constraints and objectives, enabling it to be combined with a fuzzy chance-constrained programming model to transform the uncertainty problem into an equivalent deterministic programming model.

A brief summary
From the literature, scholars have increasingly turned their attention to the economic and environmental benefits of EOL battery recovery, as demonstrated in Table 1.Some of the recovery papers about EOL batteries have been categorized in Table 1 based on research target, reverse network, objective function, whether echelon utilization is considered, whether a specific scenario about echelon utilization is considered for reused batteries and uncertainty modeling approach.However, it also faces several challenges.
A fundamental challenge lies in the high degree of uncertainty regarding the quantity and quality of EOL batteries.Existing networks struggle to harness their full potential due to this uncertainty, which affects factors such as location decisions and distribution planning.Furthermore, many EOL batteries can undergo secondary energy storage after recycling, creating additional economic benefits.
In summary, while numerous researchers have started to focus on EOL battery recycling and the uncertainties within reverse logistics networks, literature that comprehensively addresses uncertainties in the EOL battery recycling process, taking into account factors such as economic benefits, costs, carbon emissions, and the potential for secondary battery usage after recycling, remains limited.Therefore, this paper, recognizing the potential for echelon utilization of EOL batteries, establishes a multi-objective model for designing economically and environmentally efficient recycling networks.It proposes a stochastic programming approach and employs triangular fuzzy numbers to address uncertainty within the problem.

Description of the recovery network
The EOL batteries recovery network includes a battery collection center, battery processing center, reuse center, and material recycling center, as shown in Fig. 1.
In order to recycle EOL batteries more efficiently, the proposed model would purchase the batteries voluntarily recovered by consumers.The battery collection center is the place where the EOL batteries are traded and is responsible for sending them to the battery processing center.And the battery processing center needs to test and screen the quality of the batteries shipped to determine whether they can continue to be used.Furthermore, the EOL batteries eligible for continued use can be reassembled in new batteries, thereby being utilized in other scenarios.Otherwise, the EOL batteries will be sent to a material recycling center for deep dismantling and recycling of heavy metals.The proposed model will be able to determine the location of every type of center constructed in this network and the quantity of EOL batteries which be flowed from one center to another.

Main modeling assumptions
The proposed model has ten constraints.To simplify the model, the following assumptions are proposed: • The flow of the EOL batteries only occurs between different centers, and it does not flow between centers of the same kind.• The price of EOL batteries in the battery collection center is the same.Each center's construction cost is fixed.
• The locations of centers are uncertain, but the number of alternative centers and the capacity constraints of centers are known in advance.www.nature.com/scientificreports/

Notation
In this paper, we introduce the notation, encompassing sets, parameters, decision variables, and fuzzy variables used to design the recovery network.In order to achieve the desired results, determining the accurate values of the notations is crucial.Therefore, within the limits of the mathematical model, the values of the following main notations will be determined in List of symbols section.

Methodology
In this section, the paper constructs a stochastic programming-based recovery network model for EOL batteries.First, the stochastic programming approach will be introduced, followed by a description of the model constructed.Then, the existing model will be converted into a deterministic model.Finally, the model-solving method will be identified.

Stochastic programming
Stochastic programming is a frequently used method of optimization under uncertainty.As it can be seen in the above literature, robust optimization and stochastic programming are two common optimization methods under uncertainty.Due to the intractability and higher costs of robust optimization 33 and the assumption that the probability distribution of uncertain parameters is known in stochastic programming as proposed by Dantzig (1955) 34 , the more tractable stochastic programming method specifically designed to deal with modeling problems involving uncertainty was chosen for this paper 35 .In this paper, a fuzzy chance-constrained programming model in stochastic programming is introduced.

Model construction
The paper expresses the uncertainty of the quantity of EOL batteries recycled in terms of its ambiguity, using the ambiguity of the probability of reuse and material recycling to express the uncertainty of the quality of EOL batteries recycled.

Objective function
The model encompasses two objectives.According to the above assumptions, a multi-objective model is developed to optimize site selection in the entire recovery network to minimize carbon emissions and maximize the economic benefits of the whole network.
where Eqs.(1) and (2) represent the objective function of the model.As shown below, the economic benefits of the entire network comprise six parts, i.e., benefits from batteries used for echelon utilization and material recycling F1; government subsidy F2; fixed cost F3; inter-center transportation cost F4; processing cost F5; and purchase cost F6.In Eq. ( 2), CO2 emissions are summed from inter-center transportation, the center facility establishment, processing, and echelon utilization scenarios, which are represented by P1, P2, P3, and P4.
Economic benefits.The economic benefit from reuse and material recycling in Eq. ( 3) refers to the fact that after recycled EOL batteries, they are usually sent to reuse centers or material recycling centers for processing in order to realize their residual value.In general, the echelon utilization scenario for recycled EOL batteries is battery energy storage system (BESS), based on Ref. [36][37][38][39][40] .This paper considered the benefit generated by the energy storage scenario over the life cycle, not limited to the first year.Therefore, EOL batteries that can continue to be used will be sent to the energy storage system for energy storage to create income.EOL batteries that cannot be reused in the network will be sent to a material recycling center for precious metal recycling.
Government subsidy means that as the government attaches importance to the recycling of EOL batteries, the state will introduce relevant policies to achieve compensation for the recycling of unit batteries, in Eq. ( 4).
The fixed construction cost refers to the fixed cost needed to build the processing center, reuse center, and material recycling center of the EVB recovery network, in Eq. (5).
Transport cost refers to the cost incurred in transporting the recycled EOL batteries between the various centers of the network, in Eq. ( 6).
(1) www.nature.com/scientificreports/ The processing cost refers to the cost incurred by each center in processing recycled EOL batteries.The processing cost of each center is calculated by Eq. (7).
The purchase cost refers to the recycling cost of the EOL batteries from the consumption area in the network, in Eq. ( 8).It means that the recycling of EOL batteries can be promoted by providing subsidies to consumers.
In this paper, the first objective is to maximize the total economic benefits of the network of EOL batteries for EVs, and the following formula is established.
Carbon emissions.Environmental costs are generally calculated by monitoring greenhouse gases in the environment, including CO2 generated during the transportation and processing of batteries related to the network.Furthermore, this paper takes into account the carbon emissions produced by the reuse of batteries.
P1 refers to the carbon emissions from transportation in the network.
P2 refers to the carbon emissions from the process of center construction.
P3 refers to the carbon emissions from each center in the processing of the battery.
P4 refers to the carbon emissions from the Echelon utilization scenario.
The second objective is to minimize the total carbon emissions of the network of EOL batteries for EVs.

Constraints
To solve the model, the following constraints are considered.
Demand constraints.The demand for batteries at collection centers is limited.Constraint (15) indicates the demand restriction constraints of the collection centers.
( Capacity constraint.The processing capacity of the battery processing center, the remanufacturing capacity of the reuse center, and the recycling capacity of the material recycling center are limited.Constraints (17) represent restrictions on the battery processing center, reuse center, and material recycling center.
Decision variables constraints.Constraints (18) are related to the corresponding decision variables.
Other constraints.Constraint (19) indicates that the sum of the proportion of batteries shipped to the reuse center and material recycling center is 1.

Model transformation
Survey statistics show that the carbon tax is an essential tool to promote remanufacturing and reduce carbon emissions 41 .According to the carbon tax policy, the processes associated with carbon emissions will involve the cost of the carbon tax 42 .As a result, the article combines the multi-objective optimization problem with the carbon tax, treating carbon emissions as environmental costs and converting the model into an easily solvable single objective function.As the above model contains related parameters of recycling quantity and recycling quality Qi , μ, ω , it cannot be solved directly.The paper applies stochastic programming theory and triangular fuzzy numbers to convert the model into a deterministic one.The triangular fuzzy numbers used in this paper is one of the most commonly used fuzzy numbers 43 , which can more accurately express the fuzziness and uncertainty of complex data.The fuzzy chance-constrained programming mentioned above is a kind of uncertain mathematical planning based on possibility theory and fuzzy set theory, which means that when there are fuzzy parameters in the decision environment, the decision result may not satisfy the constraint limit under certain circumstances.In order to make the objective function and the constraints can meet the constraints, it is necessary to set a certain confidence level, that is, solving for the optimal solution of the model under certain probability conditions.
Typically, the fuzzy chance-constrained programming model with fuzzy parameters is used for modeling as follows: where x is the decision variable; ξ is the fuzzy covariate; f is the objective function; g is the constraint function; ∂ is the confidence level of the objective function; β is the confidence level of the ith constraint; pos[•] denotes the probability of the event.
The key to solving fuzzy chance-constrained optimal problems is to deal with chance constraints.One method is fuzzy simulation, but the solution process is time-consuming and imprecise.Another approach is to convert the chance constraint into a clear equivalence class and then solve the model.Given the uncertainty of the objective situation and the vagueness of human thinking, fuzzy set theory has become the main method for decision making.In particular, triangular fuzzy numbers or fuzzy languages have been widely used in fuzzy control and decision making 44,45 .For example, Wang et al. 46 proposed a comprehensive evaluation method for group support systems using triangular fuzzy number measurement indicators.Thus, in order to simplify the process of solving the model, the paper introduces triangular fuzzy numbers to clarify the fuzzy chance constraints.
Assume that the triangular fuzzy number is ∂ = (l i , m i , r i ) , and the membership function is µ(x) , then: . For any given confidence level ∂(0 www.nature.com/scientificreports/if and only if the conditions for (1 − ∂)l i + ∂m i ≤ z is satisfied; for any given confidence level β(0 ≤ β ≤ 1) , pos( ∂ = z) ≥ β is valid if and only if the conditions for Where m i is called the central value of ∂ , l i and r i are the left and right boundaries of ∂ .When ∂ = l i = m i = r i , ∂ is an exact value.The membership function µ 47 represents the likelihood measure of the fuzzy factor and can be defined as , the image of membership function is shown in Fig. 2.
According to the above model, the transformed fuzzy chance-constrained programming model is as follows: St.
The paper uses triangular fuzzy numbers to represent the fuzzy recycling quantity Qi and fuzzy recycling quality µ and ω in the model, and uses the confidence levels to limit the magnitude of risk due to uncertainty, making the above model and constraints clearer.Deriving Eqs. ( 20) ~ (23) gives Eqs. ( 24) ~ ( 27). ( 20)

Model solution
The above model falls within the category of network optimization problems.In addressing such problems, conventional solving algorithms have traditionally been favored in the literature.However, these conventional methods are associated with significant computational costs, susceptibility to local optima, and low solving efficiency, as noted in reference 48 .Therefore, many scholars use intelligent algorithms, such as genetic algorithms (GA) and taboo search.Compared with traditional algorithms, intelligent algorithms have better stability.In the 1970s, the genetic algorithm was first mentioned as a technical term 49 .GA is able to find the global optimal solution rather than just the local optimal solution by searching multiple solution spaces in the problem space.Furthermore, GA uses probabilistic search methods to automatically obtain the optimal search space, which is very suitable for network optimization problems 50 including routing problems, resource allocation problems, and scheduling problems, etc. 51 , and plays a great role in solving network optimization problems.Given the robust theoretical foundation and extensive success in both theoretical and practical applications 52 , along with GA's simplicity, ease of use, and high efficiency, this paper opts for GA as the solution for the model.The calculation steps are shown in Fig. 3.
Step 1: Import of basic data and encode.Firstly, the parameters in the model and the base data of the network are imported to complete the parameter initialization.The coding method is determined according to the characteristics of the variables, and the corresponding initial populations are randomly generated by the coding operation.
Step 2: Population initialization.According to the encoding method selected in Step 1, N random initial string structures are generated to form the population, each string structure is an individual of the population, and N is the size of the generated population.
Step 3: Calculate the fitness value.In a genetic algorithm, the value of the fitness function is the criterion for distinguishing the goodness of individuals in a population.The fitness function is developed based on the objective function, and the population is decoded using certain rules to calculate individual fitness values.(24)   [− i∈I j∈J www.nature.com/scientificreports/ Step 4: Genetic operation.Genetic operators perform the recombination and mutation of genetic algorithms, including selection, crossover, mutation, etc.The genetic material drives the optimization of the population in these three ways, moving on to Step 2, finding the superior offspring.
Step 5: Determine termination conditions.One of the reasons why a genetic algorithm is so computationally efficient is that it can set conditions to control the algorithm process, thus reducing useless computations.When the algorithm satisfies the target termination and iteration termination conditions, the operation is terminated and the optimal result is output; if not, the above steps are continued.

Case study Background
The paper takes the Region T in China as an example.As a major province in the manufacture and consumption of new energy vehicles in China, new EVs in Region T are growing rapidly and contribute significantly to the national total 53 .In response to policies and regulations supporting battery recycling, and considering the uncertainty of the quantity and quality of EOL batteries in the recycling process, the paper proposes an EOL batteries recovery network.The output of Lithium iron phosphate (LFP batteries) batteries in China reached 125.4 GWh in 2021, accounting for 57.1% of the total power batteries output in China.And in 2021, the National Energy Administration (NEA) issued the "Twenty-five Key Requirements for Preventing Electricity Production Accidents", which states that ternary lithium-ion batteries (NMC batteries) will no longer be used for energy storage.The market share of LFP batteries increased further in 2022, and its production accounted for 61% of the total production of power batteries.Therefore, this paper focuses on the recycling of LFP batteries in Region T 54 .As shown in Fig. 4, ten battery collection centers, five battery processing centers, five reuse centers, and three material recycling centers are involved in the battery recycling industry chain in Region T. In particular, the scale of the map in Fig. 4 is 1:7000,000, where 1 cm on the map corresponds to 70 km in the actual physical space.

Date description
According to the survey statistics, it is assumed that the ratio of the batteries flowing to the reuse center and material recycling center is generally 0.88 and 0.12 based on the reference 55 .From other literature, it is known that fuzzy affiliation generally maps 10% of the number of fuzzy centers 56 .Thus, the article makes that the fuzzy affiliation maps within 10% of the number of fuzzy centers.Based on relevant information and actual enterprise application data research results, Table 2 shows the number of EOL batteries of each region, expressed in kWh 8, 57 .The values are presented in the original units (kg) and for clarity, based on the calculations from the Chinese LCA Database also in kWh using the conversion 1 kWh = 8.643 kg.Other uncertain parameters involved in the network are shown in Table 3.The unit-relevant parameters are shown in Table 4.Meanwhile, according to the "Notice on the key work related to the management of enterprise greenhouse gas emission report in 2022" issued by the Ministry of Ecology and Environment of the People's Republic of China (MEEPRC), the latest grid emission factor is 0.5810tCO2/ MWh 58 .As stated by the "2021 China Carbon Price survey report", the average carbon price in the national carbon market in 2022 is expected to be 49 yuan/ton 59 .www.nature.com/scientificreports/ The processing center connects the collection center and reuse center, as well as the material recycling center, so its location directly affects the network.Region T can be roughly divided into 10 regions.Based on the number of logistics parks or new energy companies in each region, the corresponding region is selected as an alternative location for the battery processing centers, reuse centers, and material recycling centers.And the location coordinates of alternative centers are calculated from AutoNavi.
The related emission data and capacity limits are shown in Table 5.The paper assumes that each alternative center covers 200 square meters.Based on the relevant data and actual enterprise application data finding, the confidence level of uncertain parameters is set, in which suppose ∂ 1 = ∂ 2 = 0.9,β 1 = β 2 = β 3 = 0.8.

Results
The paper utilized MATLAB2022Ra software to perform this algorithmic example with a standalone run time of 120 s.A laptop with Intel Core i5 processor and 16 GB RAM was used as the hardware device.The case study was solved using the Improved Hybrid Collaborative Genetic Algorithm, in which the population size was assumed to be 200, the maximum number of iterations to be 600, the maximum crossover probability to be 0.8; and the minimum crossover probability to be 0.2.When the constraint is satisfied, it can be seen from Tables 6 and 7 that the logistics network is most environmentally friendly and economically beneficial when J1, J2, and J4 are selected for the battery processing center, M1, M2, and M3 for the reuse center, and N1 and N3 for the material recycling centers.The optimal total economic benefits value of the model is calculated to be 0.229845 trillion yuan, of which the optimal environmental cost is 15,999 yuan and the economic cost is 0.457549 trillion yuan.As shown in Fig. 5, nearly 46% of total carbon emissions arise from the batteries processing process, followed by the transportation process.And the batteries processing process generates the most economic costs and found it accounted for 75% of total economic costs, followed by the transportation process.

Discussion
Comparison of model results.Generally speaking, considering uncertainty in a model can have a substantial impact on the design of a network.After satisfying the constraints, the paper compares uncertain and deterministic conditions.By varying the amount of battery recycling, three scenarios are set up. Figure 6 shows the

Parameters Value
Processing cost in the processing center 17,000 yuan per metric ton 8,18 Processing cost in reuse center 1,28,000 yuan per metric ton 8,18 Processing cost in a material recycling center 42,000 yuan per metric ton 8 Fixed construction costs in the processing center 0.8 million yuan 6 Fixed construction costs in the reuse center 1 million yuan www.nature.com/scientificreports/comparison of the three scenarios in terms of economic benefits and environmental costs under uncertain and deterministic conditions.The deterministic model has about 20,000 kg more carbon emissions than the uncertainty model.The total benefits of the deterministic model are on average 10 million less than those of the uncertain model.It can be observed from the figures that the uncertain condition gives significantly better results than the deterministic condition for all three scenarios, both in terms of economic benefits and environmental costs.Moreover, the difference in economic benefits is distinct between the scenarios.This indicates that with an increase in the quantity of recycling, the economic benefits will be enhanced, which also encourages companies to focus on battery recycling.The results suggest that uncertainty needs to be taken into account in the model to obtain greater benefits and reduce the impact on the environment.Solutions under uncertainty not only provide greater flexibility for the network but also have a significantly smaller and more meaningful impact on the environment.Region T only represents a fraction of the country.Based on the distribution of registered enterprises for EV battery recycling in China, it is evident that the East and Central regions of  68 .Based on the findings from the case study in Region T, if the entire nation or even the global community focuses on EOL battery recycling and promotes similar models, it could yield more substantial potential environmental and economic benefits.We recommend the formulation of a nationwide policy for the re-cycling and utilization of EOL batteries, in collaboration with local governments, industries, and various sectors of society, to ensure feasibility and sustainability.This initiative would contribute to a perspective on global sustainable development.

Sensitivity analysis
It is clear from the research that many factors influence the total economic benefits of the network, such as the time of recycling, the location of recycling, the quantity of batteries recycled, the quality of the batteries, the unit purchase cost, and relevant policies.This paper considers the uncertainty of the quantity of batteries recycled and the quality of the batteries.Consequently, this section intends to analyze the impact on the model results and the value of the objective function by varying the values of the parameters with uncertainty.
Firstly, in order to analyze the changes in the model results, the amount of battery recycled was changed and five scenarios were set up.The results show that the change in recycling quantity has little effect on site selection results.And from Fig. 7, it can be intuitively seen that the difference between the unit economic benefits and unit environmental costs under the five scenarios is also not significant.But due to the scale effect, as the proportion of EOL electric vehicle batteries recycled increases, the total profit of the reverse logistics network tends to increase, so the number of nodal facilities required will also increase.Therefore, the decision maker needs to decide on the number or scale of alternative points for recycling nodal facilities based on the quantity of EOL electric vehicle batteries to be recycled.
Furthermore, the paper also changes the battery recycling rate, similar to changing the recycling quality, which is the probability of moving the battery from the battery processing center to the reuse center.Since the recycling rate in this case is already large enough, this paper analyzes the impact on the network by reducing the remanufacturing recycling rate.As can be seen from Fig. 8  the total benefit is reduced by 60.205 million yuan and the carbon emission is increased by 116,250 kg.As the recycling rate decreases, the total economic benefits and unit economic benefits decrease significantly, while the environmental costs increase significantly.It is reasonable to assume that an increase in recycling rates should contribute to an increase in total benefits and a decrease in environmental costs.The conclusion suggests that the model results would be better if more EOL batteries were reused.This also proves the advantages of echelon utilization, being more environmentally friendly and sustainable.The lower the quality of the product, the lower the network benefits.Likewise, if the quality of the product is improved, network benefits increase, and then more reuse centers will be needed.Therefore, by comparing the results of the above scenarios, it is found that the changes in recycling rate have a greater impact on the network than the changes in recycling quantity.

Conclusions
This study introduces an EOL electric vehicle battery recovery network design approach that accounts for uncertainty.In this study, we developed a battery recycling model that incorporates environmental costs and benefits within energy storage scenarios, taking into account the uncertain quality and quantity of EOL batteries.To enhance model clarity, we introduced a stochastic programming approach, combining uncertainty programming and a fuzzy mathematical method.The GA was applied for solving this model.These models and approaches offer a rational decision-making framework for new energy enterprises involved in battery recycling.The results from the Region T case study demonstrate the effectiveness of the proposed method in enhancing network benefits and lowering network costs.In contrast to the deterministic model, the model that accounts for uncertainty is expected to reduce carbon emissions by 20 metric tons and boost overall benefits by approximately 10 million yuan.Furthermore, the paper concludes that battery processing has the greatest influence on environmental and eco-economic costs.Hence, the primary focus of network optimization is now on developing and implementing cost-effective and environmentally friendly recycling technologies and transportation methods.Moreover, the data above reveals that a 10% reduction in the recycling rate significantly impacts the model's outcomes, leading to a decrease of 60.205 million yuan in total revenue and an increase of 116.25 tons in carbon emissions.This highlights the strong correlation between the recycling rate, economic benefits, and environmental impacts.Increasing the recycling of used batteries not only enhances the likelihood of reuse but also boosts economic benefits, with a reduced negative impact on the environment.Consequently, this study promotes increased active involvement of companies in the recycling of hazardous used batteries.While alterations in recycling volume have a minor influence on the network, they do impact the overall network scale.
To deepen our understanding of EOL battery recycling, further research on uncertainty is necessary.While the uncertainty in the quantity and quality of EOL batteries significantly affects the results of network optimization, it's essential not to overlook other factors, such as collection time and prices, for a more precise establishment of recycling network models.Dealing with uncertainty is crucial to address existing uncertain factors, and methods for handling uncertainty, such as scenario analysis and Monte Carlo simulations, can be applied.Furthermore, there are various optimization algorithms available for model solving.Moreover, future research can compare the impact of different processing and transportation methods on recycling since they have a profound influence on EOL battery recycling.

Figure 2 .
Figure 2. The image of membership function.

FFigure 7 .
Figure 7.Comparison of results in five scenarios.

Figure 8 .
Figure 8.Comparison of results in three scenarios.

Table 1 .
Summary of previous studies.

Table 2 .
Quantity of recycled at each collection center of Region T (unit: kWh).

Table 6 .
Flow Distribution from the battery collection centers to the processing centers.

Table 7 .
Flow Distribution from the battery processing centers to the next centers.Percentage of Carbon emissions and Economic costs by category.China show a greater emphasis on the recycling of EOL batteries.According to the registered enterprise data in Region T which belongs to the East regions, it comprises only 10% of the national total 67 .Additionally, the data from the case study indicates that the quantity of recycling in Region T is approximately 2% of the national total Vol.:(0123456789) Scientific Reports | (2024) 14:876 | https://doi.org/10.1038/s41598-024-51169-6www.nature.com/scientificreports/ , when the recycling rate is reduced by 10%, Comparison of economic benefits and environmental in two conditions.