<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.3 20210610//EN" "JATS-journalpublishing1-3-mathml3.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:ali="http://www.niso.org/schemas/ali/1.0/" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" dtd-version="1.3" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Comput. Sci.</journal-id>
<journal-title-group>
<journal-title>Frontiers in Computer Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Comput. Sci.</abbrev-journal-title>
</journal-title-group>
<issn pub-type="epub">2624-9898</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fcomp.2025.1729185</article-id>
<article-version article-version-type="Version of Record" vocab="NISO-RP-8-2008"/>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Original Research</subject>
</subj-group>
</article-categories>
<title-group>
<article-title>A hierarchical optimization model for off-peak battery swapping scheduling of electric trucks in open-pit mines</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Mao</surname> <given-names>Chaoli</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Conceptualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/conceptualization/">Conceptualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Data curation" vocab-term-identifier="https://credit.niso.org/contributor-roles/data-curation/">Data curation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Formal analysis" vocab-term-identifier="https://credit.niso.org/contributor-roles/formal-analysis/">Formal analysis</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Methodology" vocab-term-identifier="https://credit.niso.org/contributor-roles/methodology/">Methodology</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Software" vocab-term-identifier="https://credit.niso.org/contributor-roles/software/">Software</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Validation" vocab-term-identifier="https://credit.niso.org/contributor-roles/validation/">Validation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Visualization" vocab-term-identifier="https://credit.niso.org/contributor-roles/visualization/">Visualization</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; original draft" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-original-draft/">Writing &#x2013; original draft</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
<uri xlink:href="https://loop.frontiersin.org/people/3249460"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Tan</surname> <given-names>Yonghong</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<xref ref-type="corresp" rid="c001"><sup>&#x0002A;</sup></xref>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Funding acquisition" vocab-term-identifier="https://credit.niso.org/contributor-roles/funding-acquisition/">Funding acquisition</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Supervision" vocab-term-identifier="https://credit.niso.org/contributor-roles/supervision/">Supervision</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Xie</surname> <given-names>Shuangbo</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhou</surname> <given-names>XueBin</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Zeng</surname> <given-names>Xianren</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Linhui</given-names></name>
<xref ref-type="aff" rid="aff1"/>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Investigation" vocab-term-identifier="https://credit.niso.org/contributor-roles/investigation/">Investigation</role>
<role vocab="credit" vocab-identifier="https://credit.niso.org/" vocab-term="Writing &#x2013; review &amp; editing" vocab-term-identifier="https://credit.niso.org/contributor-roles/writing-review-editing/">Writing &#x2013; review &#x00026; editing</role>
</contrib>
</contrib-group>
<aff id="aff1"><institution>Hunan Provincial Higher Education Institutions Key Laboratory of Small and Micro Intelligent Agricultural Machinery Equipment and Application, Hunan Engineering Research Center for Smart Agriculture (Fruits and Vegetables) Information Perception and Early Warning, School of Intelligent Manufacturing, Hunan Institute of Intelligent Manufacturing and Big Data Modern Industry, Hunan University of Science and Engineering</institution>, <city>Yongzhou, Hunan</city>, <country country="cn">China</country></aff>
<author-notes>
<corresp id="c001"><label>&#x0002A;</label>Correspondence: Yonghong Tan, <email xlink:href="mailto:tyh2977@huse.edu.cn">tyh2977@huse.edu.cn</email></corresp>
</author-notes>
<pub-date publication-format="electronic" date-type="pub" iso-8601-date="2025-12-10">
<day>10</day>
<month>12</month>
<year>2025</year>
</pub-date>
<pub-date publication-format="electronic" date-type="collection">
<year>2025</year>
</pub-date>
<volume>7</volume>
<elocation-id>1729185</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>10</month>
<year>2025</year>
</date>
<date date-type="rev-recd">
<day>14</day>
<month>11</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>11</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2025 Mao, Tan, Xie, Zhou, Zeng and Wang.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Mao, Tan, Xie, Zhou, Zeng and Wang</copyright-holder>
<license>
<ali:license_ref start_date="2025-12-10">https://creativecommons.org/licenses/by/4.0/</ali:license_ref>
<license-p>This is an open-access article distributed under the terms of the <ext-link ext-link-type="uri" xlink:href="https://creativecommons.org/licenses/by/4.0/">Creative Commons Attribution License (CC BY)</ext-link>. The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</license-p>
</license>
</permissions>
<abstract>
<p>This study addresses the queuing inefficiencies caused by synchronized battery-swapping demands for electric trucks in open-pit mines. Through Discrete Event Simulation (DES), we identified systemic bottlenecks stemming from this synchronization. To mitigate this, we propose a hierarchical off-peak battery-swapping scheduling framework comprising an inner-layer Mixed-Integer Linear Programming (MILP) and an outer-layer Bayesian Optimization (BO) mechanism. Validated through three large-scale case studies, the model achieved 65% and 80% reductions in queuing times for single and dual loading platform scenarios, respectively, with 5.2%&#x02013;5.7% improvements in transport throughput. Notably, expanding battery-swapping stations to four achieved equivalent efficiency gains (667 trips) as the optimization strategy (665 trips), highlighting the cost-effectiveness of intelligent scheduling over infrastructure scaling. Furthermore, in the third case study, by increasing loading platforms to alleviate constraints from upstream processes, the optimized model boosts transportation trips by up to 10%, demonstrating its capability to eliminate battery-swapping bottlenecks and fully unlock the potential of energy replenishment workflows.</p></abstract>
<kwd-group>
<kwd>hierarchical optimization model</kwd>
<kwd>off-peak battery swapping scheduling</kwd>
<kwd>open-pit mines</kwd>
<kwd>electric trucks</kwd>
<kwd>discrete event simulation</kwd>
</kwd-group>
<funding-group>
  <funding-statement>The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Natural Science Foundation of Hunan Province, China (Grant Nos. 2024JJ6227, 2024JJ7189, and 2024JJ7197).</funding-statement>
</funding-group>
<counts>
<fig-count count="14"/>
<table-count count="3"/>
<equation-count count="15"/>
<ref-count count="23"/>
<page-count count="16"/>
<word-count count="9541"/>
</counts>
<custom-meta-group>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Networks and Communications</meta-value>
</custom-meta>
</custom-meta-group>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<label>1</label>
<title>Introduction</title>
<p>Traditional fuel-powered mining trucks that burn diesel or gasoline emit a large amount of pollutants such as carbon dioxide, nitrogen oxides, and particulate matter into atmosphere (<xref ref-type="bibr" rid="B1">Bai et al., 2024</xref>), which exacerbates the ongoing concerns about global climate change. In contrast, electric mining trucks are almost zero emissions as their eletricity can be from renewable energy sources (<xref ref-type="bibr" rid="B21">Zhang et al., 2024</xref>). Besides, they also have an advantage in terms of haulage costs (<xref ref-type="bibr" rid="B10">Lindgren et al., 2022</xref>). Hence, transport electrification in open-pit mines is becoming a trend (<xref ref-type="bibr" rid="B12">Teng et al., 2024</xref>).</p>
<p>Despite this trend, the electrification rate of open-pit mining trucks remains the lowest across all vehicle segments (<xref ref-type="bibr" rid="B5">GlobalData, 2025</xref>). This is primarily attributed to the fact that, in the context of open-pit mines, electric-powered mining trucks necessitate batteries of substantially larger capacity, specifically in the range of 700&#x02013;800 KWh. Even when employing fast charging technology, it requires more than one hour for such a large battery to reach fully charged status. Such extend charging time significantly hinder the production efficiency of open-pit mines. Given the number of mining trucks, excavators, and the capacity of crushing plant, open-pit mines aim to maximize the number of transport trips per unit of time. If electric-powered mining trucks spend excessive time on charging, their avaiable time for transportation is inevitably reduced, leading to a decrease in the number of daily transport trips.</p>
<p>To address the efficiency decline caused by prolonged charging times, battery-swapping electric-powered mining trucks can be deployed as an alternative (<xref ref-type="bibr" rid="B13">Vallera et al., 2021</xref>). In battery-swapping mode, robotic arms automatically replace depleted batteries with fully charged ones in an average of 6-8 minutes, which is comparable to refueling times. Thus, this mode is particularly well-suited for open-pit mining operations (<xref ref-type="bibr" rid="B23">Zhu et al., 2023</xref>), where stringent demands for rapid energy replenishment are critical. The Chinese open-pit mining sector currently has thousands of battery-swap-capable electric haul trucks in active service.</p>
<p>However, battery-swapping electric mining trucks differ significantly from conventional diesel haulers in terms of operational constraints such as range limitations, the need for battery-swapping infrastructure, and fleet management logistics. Traditional dispatch strategies designed for fuel-based fleets are thus ill-equipped for battery-swapping operations, potentially resulting in suboptimal productivity and elevated haulage costs. Therefore, an optimized dispatch model tailored to these electric trucks is essential (<xref ref-type="bibr" rid="B6">He et al., 2017</xref>; <xref ref-type="bibr" rid="B18">Xu et al., 2019</xref>).</p>
<p>Recognizing this challenge, researchers have proposed various multi-objective scheduling models to address diverse operational demands. <xref ref-type="bibr" rid="B20">Zhan et al. (2022)</xref> conducted a comprehensive review of optimal charging scheduling strategies for battery-swapping stations. The main optimization objectives involved are minimizing operational costs (e.g., charging expenses and battery degradation), maximizing grid stability (peak shaving, renewable energy integration), ensuring user satisfaction (reduced waiting time, service availability) and managing uncertainties in battery swapping demand. <xref ref-type="bibr" rid="B4">Deng et al. (2023)</xref> provided a cost-effective framework for the design of battery-swapping station tailored to logistic fleets of electric trucks, integrating battery degradation dynamics and practical charging constraints. The proposed management strategy enhances battery lifespan and reduces lifecycle costs, supporting sustainable electrification of freight transport. <xref ref-type="bibr" rid="B11">Sun et al. (2024)</xref> proposed an effective battery-swapping dispatch framework for a self-sustained highway energy system, focusing on balancing battery supply-demand dynamics between battery-charging and batteryswapping stations, while minimizing total transportation costs. To begin with, a deep-learning-based spatiotemporal traffic flow network was designed to accurately forecast the demands, thereby providing a crucial foundation for the subsequent model decisions. Subsequently, a battery-swapping dispatch model was formulated, as an extension of vehicle routing problem. It comprehensively considers multiple cost factors, including the routing cost, which is determined by the distance traveled and a fixed cost per kilometer; the cost associated with the degree of satisfaction, which is contingent upon the accuracy of demand prediction; and the cost of using battery trucks. Concerning the short-term prediction of battery-swapping demands for passenger electric vehicle, <xref ref-type="bibr" rid="B16">Wang et al. (2023)</xref> also proposed a series of deep learning modles. A real-world dataset containing 2,529 battery-swapping events collected from 36 battery-swapping stations in Beijing was utilized to train these modles. <xref ref-type="bibr" rid="B15">Wang et al. (2024)</xref> pointed out that high investment costs and suboptimal operational strategies hinder the profitability of battery-swapping stations for heavy-duty electric trucks. To address these challagens, they proposed a novel scheduling method and a two-layered optimization framework to enhance the performance of battery-swapping stations, while integrating photovoltaic systems for cost-effectiveness. <xref ref-type="bibr" rid="B19">Yang et al. (2024)</xref> presented an online scheduling framework to dispatch batteries between battery-charging stations and battery-swapping stations efficiently, combining partial delivery logistics and real-time optimization to enhance operational flexibility and reduce costs. The proposed solving algorithm balances computational efficiency and solution quality, making it suitable for practical electric vehicle infrastructure management. Besides, the solving algorithm is embedded in a rolling-horizon framework with the introduction of dummy copies, which enables the algorithm to tackle future uncertainties in the prediction of battery demands. Regarding the same type of problems within this battery-swapping-charging system, <xref ref-type="bibr" rid="B22">Zhang and Wang (2016)</xref> proposed a two-direction model which is solved by particle swarm optimization(PSO) method to reduce transportation cost. <xref ref-type="bibr" rid="B7">Huang et al. (2021)</xref> established a nonlinear Mixed-Integer Programming model considering dynamic electricity prices to determine the optimal charging schedules for batteries in battery-charging stations and transportation schedules for the dispatch of batteries. <xref ref-type="bibr" rid="B2">Ban et al. (2021)</xref> focused on battery-swapping-charging system composed of multiple nanogrids and battery-swapping stations. The system aims to enhance energy supply cleanliness and promote transportation electrification by aggregating distributed renewable energy sources. The authors established a joint optimal scheduling model based on mixed-integer linear programming, which takes into account battery charging/discharging in the battery-swapping stations, battery swapping among individual units, and the vehicle routing problem of battery transporters. To simplify the model, several pre-processing technologies and assumptions are employed, such as setting dummy copies and discretizing the planning interval. Besides, <xref ref-type="bibr" rid="B9">Jordehi et al. (2021)</xref> investigated the placement of battery-swapping stations in microgrids with various renewable energy sources. The authors formulated a nonlinear Mixed-Integer Programing model considering multiple constraints. Case studies show that location of battery-swapping stations significantly impacts operation cost of microgrids. From the perspective of the battery-swapping scenario of electric trucks in open-pit mines, the work by <xref ref-type="bibr" rid="B17">Xiao et al. (2024)</xref> is the most related one. The authors proposed a multi-objective scheduling optimization model, aiming to minimize total haulage cost and total waiting time during a single shift(about 8 hours). The model takes into account unique features like battery-swapping alerts, station selection, and the impact of ambient temperature on battery capacity.</p>
<p>While prior studies have established foundational frameworks for developing battery dispatch strategies in open-pit mine electric haulage systems, our comprehensive literature review reveals a critical gap in real-time scheduling optimization for battery-swapping operations. The most related research (<xref ref-type="bibr" rid="B17">Xiao et al., 2024</xref>) appears predominantly focused on static scheduling models, leaving dynamic operational scheduling with fluctuating energy demands substantially underexplored. This gap becomes particularly acute in scenarios characterized by high-density truck deployment coupled with insufficient swapping infrastructure, where limited battery inventory at individual stations creates operational bottlenecks. Such resource-constrained conditions inevitably lead to queuing inefficiencies during peak demand periods, adversely impacting overall productivity through extended equipment downtime. In response to these operational challenges, this study proposes a hierarchical scheduling optimization model with queuing time minimization as its primary objective. The core architecture comprises of an inner-layer of Mixed-Integer Linear Programming (MILP) and an out-layer of Bayesian optimization (BO). This two-layer design resolves multi-truck coordination complexity in battery-constrained scenarios while keeping computations feasible via reduced parameter space. To formulate the inner MILP, the planning interval is discretized (<xref ref-type="bibr" rid="B2">Ban et al., 2021</xref>). In a sort of way, this pre-processing technology can also enables decision-making under uncertainly.</p>
<p>The remainder of this paper is structured as follows. Section 2 deconstructs the operational workflow of battery-swapping electric trucks in open-pit mines, followed by a Discrete Event Simulation (DES) to quantify bottlenecks under varying truck densities and station capacities. Section 3 formulates the hierarchical optimization framework. Section 4 validates the framework through several case studies, benchmarking its efficacy against conventional scheduling strategies via several metrics (such as total trips and waiting time at the battery-swapping stations). Section 5 concludes with theoretical and practical implications.</p></sec>
<sec id="s2">
<label>2</label>
<title>The workflow in open-pit mines</title>
<sec>
<label>2.1</label>
<title>System description</title>
<p>The transportation system with battery-swapping trucks in open-pit mines comprises multiple operational nodes, including loading platforms, unloading zones, and battery-swapping stations. Electric haul trucks cycle continuously between these nodes to transport materials while periodically requiring battery exchanges at designated swapping stations, as depicted in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig position="float" id="F1">
<label>Figure 1</label>
<caption><p>The workflow in open-pit mines.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0001.tif">
<alt-text content-type="machine-generated">Diagram of a mine transportation system showing two loading platforms, a battery-swapping station area, and an unloading zone. Loading platform 1 includes shovel 1 and 2 with trucks entering and exiting. Loading platform 2 mirrors this arrangement. The battery-swapping area features two stations for trucks. The unloading zone has a crusher with trucks entering and exiting. Arrows depict the travel path for loading and unloading.</alt-text>
</graphic>
</fig>
<p>In line with the workflow depicted in <xref ref-type="fig" rid="F1">Figure 1</xref>, a simple scheduling method for mining truck battery swapping functions as follows. Once unloading is completed, the system assesses whether the truck&#x00027;s remaining battery capacity can support one more haul cycle. If it can, the truck is sent to continue its hauling operations. If not, it is routed to the battery-swapping station with a shorter queue for battery replacement. At the station, trucks in the queue are served strictly on a first-come-first-served (FCFS) basis, without considering any operational priority levels. However, this straightforward scheduling approach may give rise to two issues. Firstly, when the battery depletion moments of a large number of mining trucks approach simultaneously, these trucks will flock to the battery-swapping stations for battery replacement. Due to the limited capacity of the battery-swapping stations and the finite number of fully charged batteries, long queues will be formed, which will affect the operational duration of the mining trucks and cause traffic congestion around the battery-swapping stations. Secondly, FCFS basis is suboptimal and as a result, the total queuing time for battery replacement of the entire fleet may be excessively long.</p>
</sec>
<sec>
<label>2.2</label>
<title>Discrete event simulation of the system</title>
<p>Discrete Event Simulation (DES) is a computational modeling approach that mimics the dynamic behavior of complex systems by simulating discrete events (e.g., task arrival, service completion) occurring at specific time points (<xref ref-type="bibr" rid="B3">Banks et al., 2013</xref>). Unlike continuous simulation, which models systems with smooth state changes (e.g., fluid flow), DES focuses on abrupt state transitions triggered by events. For example, in mining truck scheduling, DES tracks when trucks arrive at battery-swapping stations, how long they queue, and when they resume hauling, considering random battery depletion times and station capacities (<xref ref-type="bibr" rid="B8">Huayanca et al., 2023</xref>). Such a simulation can help identify bottlenecks and improve resource allocation.</p>
<p>In this paper, we conduct a discrete event simulation (DES) to model the workflow in open-pit mines. The parameter configurations for the simulation are summarized in <xref ref-type="table" rid="T1">Table 1</xref>. These settings are based on an actual open-pit coal mine and the work by <xref ref-type="bibr" rid="B8">Huayanca et al. (2023)</xref>. The simulation framework encompasses three interconnected processes: loading, unloading, and battery-swapping. Initially, all trucks&#x02014;equipped with a preset state of charge (SoC)&#x02014;depart from the unloading zone and travel to designated loading platforms. Upon arrival, trucks queue at the loading platform until they are loaded with material. Once loaded, they return to the unloading zone for material discharge. After completing the unloading process, the system evaluates whether the truck&#x00027;s remaining SoC is sufficient for another hauling cycle. If the SoC falls below the required threshold, the truck is redirected to a battery-swapping station to replace its depleted battery with a fully charged one. Given the inherent uncertainties in travel time, spotting time, battery-swapping duration, and other operational variables, the simulation&#x02014;modeling an 18-hour production period&#x02014;was executed 100 times to account for stochastic variations. The results yielded an average of 628 total trips completed across all iterations, while the average waiting time for battery-swapping stood at 86.6 minutes. <xref ref-type="fig" rid="F2">Figure 2</xref> illustrates the relationship between total trips and battery-swapping station wait times through a scatter plot. As shown, a distinct negative correlation emerges between these two metrics, suggesting that prolonged battery-swapping delays directly reduce the number of achievable haulage cycles. This inverse relationship positions battery-swapping operations as a potential critical bottleneck in the production workflow.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p>The parameter settings for DES.</p></caption>
<table frame="box" rules="all">
<thead>
<tr>
<th valign="top" align="left"><bold>Parameter</bold></th>
<th valign="top" align="left"><bold>Description</bold></th>
<th valign="top" align="left"><bold>Value/range/expression</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Simulation duration</td>
<td valign="top" align="left">Total time span of the simulation</td>
<td valign="top" align="left">18 h</td>
</tr>
<tr>
<td valign="top" align="left">Truck fleet size</td>
<td valign="top" align="left">Number of mining trucks in the simulation</td>
<td valign="top" align="left">36</td>
</tr>
<tr>
<td valign="top" align="left">Battery depletion rate for unloaded travel</td>
<td valign="top" align="left">State of charge consumption per minute, normal distribution</td>
<td valign="top" align="left">Mean: 0.25%/min, standard error: 3e-4%/min</td>
</tr>
<tr>
<td valign="top" align="left">Battery depletion rate for loaded travel</td>
<td valign="top" align="left">State of Charge consumption per minute, normal distribution</td>
<td valign="top" align="left">Mean: 0.33%/min, standard error: 5e-4%/min</td>
</tr>
<tr>
<td valign="top" align="left">Swap station capacity</td>
<td valign="top" align="left">Number of parallel service bays per station</td>
<td valign="top" align="left">1</td>
</tr>
<tr>
<td valign="top" align="left">Number of stations</td>
<td valign="top" align="left">Number of battery-swapping stations</td>
<td valign="top" align="left">2</td>
</tr>
<tr>
<td valign="top" align="left">Service time</td>
<td valign="top" align="left">The time required to replace a battery and prepare to service next truck, normal distribution</td>
<td valign="top" align="left">Mean: 10 min, standard error: 1 min</td>
</tr>
<tr>
<td valign="top" align="left">Full battery number</td>
<td valign="top" align="left">Spare fully charged batteries per station</td>
<td valign="top" align="left">5</td>
</tr>
<tr>
<td valign="top" align="left">Service rule of battery swapping</td>
<td valign="top" align="left">Queueing discipline at battery-swapping stations</td>
<td valign="top" align="left">FCFS</td>
</tr>
<tr>
<td valign="top" align="left">Service rule of loading</td>
<td valign="top" align="left">Queueing discipline at shovel</td>
<td valign="top" align="left">FCFS</td>
</tr>
<tr>
<td valign="top" align="left">Service rule of unloading</td>
<td valign="top" align="left">Queueing discipline at crusher</td>
<td valign="top" align="left">FCFS</td>
</tr>
<tr>
<td valign="top" align="left">Unloading travel time from unloading zone to loading platform 1</td>
<td valign="top" align="left">The time required to travel from unloading zone to loading platform 1, normal distribution</td>
<td valign="top" align="left">Mean: 12 min, standard error: 0.6 min</td>
</tr>
<tr>
<td valign="top" align="left">Unloading travel time from unloading zone to loading platform 2</td>
<td valign="top" align="left">The time required to travel from unloading zone to loading platform 2, normal distribution</td>
<td valign="top" align="left">Mean: 15 min, standard error: 0.9 min</td>
</tr>
<tr>
<td valign="top" align="left">loading travel time from loading platform 1 to unloading zone</td>
<td valign="top" align="left">The time required to travel from loading platform 1 to unloading zone, normal distribution</td>
<td valign="top" align="left">Mean: 22 min, standard error: 1.2 min</td>
</tr>
<tr>
<td valign="top" align="left">loading travel time from loading platform 2 to unloading zone</td>
<td valign="top" align="left">The time required to travel from loading platform 2 to unloading zone, normal distribution</td>
<td valign="top" align="left">Mean: 25 min, standard error: 1.5 min</td>
</tr>
<tr>
<td valign="top" align="left">unloading time</td>
<td valign="top" align="left">The time required for one truck to be loaded, normal distribution</td>
<td valign="top" align="left">Mean: 1 min, standard error: 0.02 min</td>
</tr>
<tr>
<td valign="top" align="left">loading time</td>
<td valign="top" align="left">The time required for one truck to be unloaded, normal distribution</td>
<td valign="top" align="left">Mean: 5 min, standard error: 0.5 min</td>
</tr>
<tr>
<td valign="top" align="left">Threshold of State of Charge</td>
<td valign="top" align="left">After unloading, if state of charge is under this value, go to exchange battery</td>
<td valign="top" align="left">30%</td>
</tr>
<tr>
<td valign="top" align="left">Spotting time at shovel</td>
<td valign="top" align="left">The time required to spot around shovel, beta distribution</td>
<td valign="top" align="left">0.73&#x0002B;0.84*BETA(1.37, 1.46) min</td>
</tr>
<tr>
<td valign="top" align="left">Spotting time at unloading spot</td>
<td valign="top" align="left">The time required to spot arround unloading spot, normal distribution</td>
<td valign="top" align="left">Mean: 1.55 min, standard error: 0.06 min</td>
</tr>
<tr>
<td valign="top" align="left">Spotting time at battery-swapping bay</td>
<td valign="top" align="left">The time required to spot arround battery-swapping bay, normal distribution</td>
<td valign="top" align="left">Mean: 1.0 min, standard error: 0.05 min</td>
</tr>
<tr>
<td valign="top" align="left">charging rate</td>
<td valign="top" align="left">The charging rate for depleted battery in battery-swapping station</td>
<td valign="top" align="left">0.95%/min</td>
</tr></tbody>
</table>
</table-wrap>
<fig position="float" id="F2">
<label>Figure 2</label>
<caption><p>Total trips vs average waiting time for battery-swapping.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0002.tif">
<alt-text content-type="machine-generated">Scatter plot showing the relationship between average waiting time for battery swapping in minutes and total trips. Blue dots represent data points. A red trend line indicates a negative correlation, with total trips decreasing as waiting time increases.</alt-text>
</graphic>
</fig>
<p>To investigate the root cause of prolonged battery-swapping delays, <xref ref-type="fig" rid="F3">Figure 3</xref> visualizes the state of charge (SoC) time series for all trucks during one representative simulation run. The data reveals a critical synchronization issue: during the first battery-swapping cycle, all trucks converge on the swapping stations approximately 350 min after production initiation. This synchronization creates a demand surge that overwhelmes the stations limited capacity, forcing some trucks in this cohort to endure wait times exceeding two hours for battery replacement. Notably, subsequent swapping cycles exhibite significantly reduced delays, as evidenced by the dispersed timing of later SoC troughs in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig position="float" id="F3">
<label>Figure 3</label>
<caption><p>SoC time series for all trucks during one simulation.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0003.tif">
<alt-text content-type="machine-generated">Line graph showing the state of charge (SoC) over time in minutes. Multiple colored lines represent different simulations, demonstrating fluctuating SoC levels, mostly between 0.2 and 1.0.</alt-text>
</graphic>
</fig>
<p>This stark contrast between the first and subsequent cycles strongly suggests that synchronized battery-swapping demand&#x02013;rather than inherent station inefficiency&#x02014;is the primary driver of bottlenecks. The clustering of trucks at the stations during peak intervals directly correlates with the observed productivity loss (i.e., reduced total trips). These findings imply that decentralizing battery-swapping activity through demand-leveling strategies, such as staggered battery replacement schedules or predictive SoC-based dispatching, could mitigate congestion and unlock throughput gains. Consequently, designing an off-peak battery-swapping protocol emerges as a critical priority for optimizing production efficiency in open-pit mining operations.</p>
<p>To further validate these findings, we designed a proactive staggered battery-swapping strategy to mitigate synchronization-induced bottlenecks. The strategy involves dynamically selecting Z trucks with state of charge (SoC) below Y% at X-minute intervals for prioritized battery replacement, contingent on the availability of fully charged batteries at the swapping stations. To prevent overloading station capacity, the parameter Z caps the maximum number of trucks dispatched per interval&#x02013;even when surplus charged batteries are available, only Z trucks are proactively rerouted. To identify the optimal parameter combination (X, Y, Z), we performed a grid search over the hyperparameter space: <italic>X</italic>:[10, 15, 20, 25] &#x000D7; <italic>Y</italic>:[40, 50, 60, 70] &#x000D7; <italic>Z</italic>:[2, 3, 4, 5]. For each combination, the simulation was replicated 100 times, with results averaged to compute total trips completed and average battery-swapping wait time. The optimal combination (<italic>X</italic> &#x0003D; 25, <italic>Y</italic> &#x0003D; 70, <italic>Z</italic> &#x0003D; 3) yielded 665 total trips&#x02014;a 5.7% improvement over the baseline (628 trips)&#x02014;while reducing average battery-swapping wait time to 32.5 min (from 86.5 min). <xref ref-type="fig" rid="F4">Figure 4</xref> illustrates the SoC time series under this configuration, demonstrating the near-elimination of prolonged delays during the initial battery-swapping cycle. Meanwhile, <xref ref-type="fig" rid="F5">Figure 5</xref> plots total trips against average wait time across all parameter combinations, reaffirming the negative correlation observed earlier. This consistent trend underscores battery-swapping efficiency as a leverage point for systemic productivity gains, with suboptimal scheduling directly constraining haulage throughput.</p>
<fig position="float" id="F4">
<label>Figure 4</label>
<caption><p>SoC time series for all trucks corresponding to the best parameter combination.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0004.tif">
<alt-text content-type="machine-generated">Line graph showing the state of charge (SoC) of a battery over simulation time in minutes. The graph features multiple overlapping, colorful lines indicating various scenarios or simulations. SoC ranges from 1.0 to below 0.2 across times up to 1200 minutes.</alt-text>
</graphic>
</fig>
<fig position="float" id="F5">
<label>Figure 5</label>
<caption><p>Total trips vs. average waiting time in staggered battery replacement schedules.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0005.tif">
<alt-text content-type="machine-generated">Scatter plot showing the relationship between average waiting time for battery-swapping in minutes (x-axis) and total trips (y-axis). Blue dots represent data points, and a red line indicates a negative correlation.</alt-text>
</graphic>
</fig>
<p>Another way to assess whether battery-swapping efficiency constitutes the dominant bottleneck is by analyzing the relationship between total haulage trips and the number of battery-swapping stations. As demonstrated in <xref ref-type="fig" rid="F6">Figure 6</xref>, increasing the number of stations initially drives a near-linear improvement in total trips, followed by a plateau (even though the waiting time decreases sharply) as the system transitions to being constrained by other workflow components&#x02013;such as truck fleet size, shovel loading capacity, or crusher throughput. Notably, expanding the number of stations to 4 achieves a throughput (667 trips) parity with the optimized staggered battery-swapping strategy proposed earlier (665 trips).</p>
<fig position="float" id="F6">
<label>Figure 6</label>
<caption><p>The effect of increasing station number on total trips and waiting time at battery-swaping stations.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0006.tif">
<alt-text content-type="machine-generated">Bar and line graph showing total trips and average waiting times across six stations. Blue bars represent total trips, peaking at Station 2. Green lines with points denote average waiting times, decreasing after Station 1.</alt-text>
</graphic>
</fig>
<p>While both approaches alleviate battery-swapping congestion, their cost implications diverge significantly. Deploying additional stations necessitates substantial capital investment in infrastructure and battery inventory. In contrast, the staggered scheduling strategy&#x02013;requiring no hardware expansion&#x02013;achieves comparable efficiency gains through operational optimization. This contrast underscores that smarter dispatching protocols, rather than infrastructure scaling, represent a cost-effective lever for mitigating bottlenecks in open-pit mining operations.</p>
<p>Building on the operational insights derived from the simulation experiments, the subsequent section proposes an optimization-driven framework for battery-swapping scheduling. This framework systematically reduces synchronization induced delays by integrating real-time state of charge (SoC) monitoring, station capacity constraints, and dynamic truck dispatching, thereby establishing a generalizable methodology to alleviate bottlenecks in open-pit mining operations.</p></sec>
</sec>
<sec id="s3">
<label>3</label>
<title>Optimization-based heriarchical model</title>
<p>In this section, we propose a hierarchical optimization framework for off-peak battery-swapping scheduling, with the primary objective of minimizing queuing delays at battery-swapping stations. The framework integrates a Mixed-Integer Linear Programming (MILP) model and a Bayesian Optimization layer to address both operational and systemic bottlenecks. At its core, the MILP model optimizes short-term resource allocation within stations by resolving battery availability, truck sequencing, and service time constraints, thereby maximizing the inherent efficiency of individual stations. Simultaneously, the Bayesian Optimization layer governs the outer loop of the hierarchy, dynamically determining optimal battery-swapping timings for each truck to avoid synchronization-induced congestion&#x02013;specifically targeting the prolonged delays observed during the first battery-swapping cycle (Section 2). By coordinating these two layers, the framework ensures that real-time scheduling decisions align with long-term system-wide efficiency goals: the MILP model fine-tunes station-level operations, while Bayesian Optimization proactively redistributes battery-swapping demand across temporal and spatial dimensions. The following subsection details the formulation of the MILP model first, which serves as the computational foundation for this hierarchical approach.</p>
<sec>
<label>3.1</label>
<title>Mixed-Integer Linear Programming model</title>
<p>Given the predicted battery-swapping demand moments during the first battery-swapping cycle (Section 2), the scheduling of electric mining trucks across multiple stations is formulated as a Mixed-Integer Linear Programming (MILP) model. The formulation begins with a temporal discretization step, where continuous time is partitioned into fixed intervals. This approach converts the intractable continuous-time optimization problem into a discrete, computationally tractable counterpart. Crucially, temporal discretization enhances robustness against demand prediction uncertainties by enabling probabilistic adjustments to resource allocations within each interval&#x02014;a mechanism that mitigates timing discrepancies between forecasted and actual demand. For illustration, <xref ref-type="fig" rid="F7">Figure 7</xref> demonstrates a 5-minute time-slice configuration. This discretization eliminates the need for exact demand timing predictions (e.g., 16:15). Instead, when a demand surge is projected for the third time slice, any occurrence within the corresponding 16:15 &#x02212; 16:20 window remains operationally viable. Such flexibility ensures resilience to minor temporal prediction errors while maintaining scheduling feasibility.</p>
<fig position="float" id="F7">
<label>Figure 7</label>
<caption><p>An example of time discretization.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0007.tif">
<alt-text content-type="machine-generated">Timeline showing an index of time slices labeled from \( t = 0 \) to \( t = 5 \), corresponding to moments in 5-minute intervals starting at 16:00 and ending at 16:30. It also displays earlier moments in descending order from \( T - 7 \) at 08:25 to \( T - 1 \) at 09:00.</alt-text>
</graphic>
</fig>
<p>Within this discretized framework, the decision variable in the MILP model is defined as a binary parameter:</p>
<disp-formula id="EQ1"><mml:math id="M1"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mn>1</mml:mn><mml:mtext class="textrm" mathvariant="normal">&#x02003;if truck i is scheduled for battery swapping at station</mml:mtext></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mtext class="textrm" mathvariant="normal">&#x02003;&#x000A0;&#x000A0;j during time slice t</mml:mtext><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mn>0</mml:mn><mml:mtext class="textrm" mathvariant="normal">&#x02003;otherwise</mml:mtext><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(1)</label></disp-formula>
<p>Here, <italic>x</italic><sub><italic>i, j, t</italic></sub> &#x0003D; 1 indicates that the dispatching system allocates truck <italic>i</italic> to the <italic>j</italic>&#x02212;<italic>th</italic> station for battery swapping at the <italic>t</italic>&#x02212;<italic>th</italic> time slice, while <italic>x</italic><sub><italic>i, j, t</italic></sub> &#x0003D; 0 denotes no such assignment. For a fleet of 36 electric mining trucks with projected battery-swapping demands over a 4-hour horizon, the MILP model&#x00027;s total number of decision variables is calculated as 36(trucks) &#x000D7; 48(time slices) &#x000D7; 2(stations) &#x0003D; 3, 456, where the temporal resolution is discretized into 5-min time slices (yielding <inline-formula><mml:math id="M2"><mml:mfrac><mml:mrow><mml:mn>4</mml:mn><mml:mo>&#x000D7;</mml:mo><mml:mn>60</mml:mn></mml:mrow><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:mfrac><mml:mo>=</mml:mo><mml:mn>48</mml:mn></mml:math></inline-formula> intervals) and two battery-swapping stations are available. This parameterization ensures the model systematically encodes all potential scheduling permutations through the binary variable <italic>x</italic><sub><italic>i, j, t</italic></sub>&#x02014;defined in <xref ref-type="disp-formula" rid="EQ1">Equation 1</xref>&#x02014;while maintaining computational tractability under the discretized temporal and spatial constraints outlined in Section 2.</p>
<p>The model&#x00027;s objective function explicitly minimizes the total queuing time, addressing the synchronization-induced bottlenecks identified in Section 2, while enforcing operational constraints such as station capacity limits, battery availability, and truck SoC thresholds. The objective function is defined as:</p>
<disp-formula id="EQ2"><mml:math id="M3"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>I</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>A</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(2)</label></disp-formula>
<p>where <italic>I</italic> is the total number of trucks, <italic>S</italic><sub><italic>i</italic></sub> denotes the start time of battery-swapping for truck <italic>i</italic> and <italic>A</italic><sub><italic>i</italic></sub> represents the arrival time of truck <italic>i</italic>, which corresponds to the predicted battery-swapping demand time (an input to the MILP model). The start time <italic>S</italic><sub><italic>i</italic></sub> is derived from the binary decision variable <italic>x</italic><sub><italic>i, j, t</italic></sub>&#x02014;defined in <xref ref-type="disp-formula" rid="EQ1">Equation 1</xref>, which assigns truck <italic>i</italic> to station <italic>j</italic> at time slice <italic>t</italic>, and is calculated as:</p>
<disp-formula id="EQ3"><mml:math id="M4"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>J</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(3)</label></disp-formula>
<p>where <italic>J</italic> is the total number of battery-swapping stations, <italic>T</italic> is the total number of discrete time slices, and <italic>t</italic> indexes the time slices (e.g., <italic>t</italic> &#x0003D; 0, 1, &#x02026;, <italic>T</italic>&#x02212;1). This formulation ensures that each truck is assigned to exactly one station within a single time slice, translating the continuous-time scheduling problem into a tractable discrete optimization framework.</p>
<p>To ensure the MILP model aligns with operational realities, four critical constraints are incorporated. First, <xref ref-type="disp-formula" rid="EQ4">Equation 4</xref> enforces that each truck undergoes exactly one battery-swapping event within the planning horizon, reflecting the focus on the first battery-swapping cycle:</p>
<disp-formula id="EQ4"><mml:math id="M5"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>j</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>J</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>=</mml:mo><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(4)</label></disp-formula>
<p>where, <italic>x</italic><sub><italic>i, j, t</italic></sub> is defined in <xref ref-type="disp-formula" rid="EQ1">Equation 1</xref>. Second, <xref ref-type="disp-formula" rid="EQ5">Equation 5</xref> ensures causality by requiring battery-swapping to occur after arrival:</p>
<disp-formula id="EQ5"><mml:math id="M6"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mi>S</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>&#x02265;</mml:mo><mml:msub><mml:mrow><mml:mi>A</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02200;</mml:mo><mml:mi>i</mml:mi><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:mi>I</mml:mi></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(5)</label></disp-formula>
<p>where <italic>S</italic><sub><italic>i</italic></sub> and <italic>A</italic><sub><italic>i</italic></sub> are defined in <xref ref-type="disp-formula" rid="EQ3">Equations 3</xref> and <xref ref-type="disp-formula" rid="EQ2">2</xref>, respectively. Third, station throughput limitations are modeled via <xref ref-type="disp-formula" rid="EQ6">Equation 6</xref>, which restricts each station <italic>j</italic> to servicing at most one truck within any sliding time window of length <italic>T</italic><sub><italic>e</italic></sub> (battery-swapping and preparation time):</p>
<disp-formula id="EQ6"><mml:math id="M7"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>I</mml:mi></mml:mrow></mml:munderover></mml:mstyle><mml:mstyle displaystyle="true"><mml:munderover accentunder="false" accent="false"><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:munderover></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x02264;</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x02200;</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:mi>J</mml:mi></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x02200;</mml:mo><mml:mi>t</mml:mi><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>0</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:mi>T</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>e</mml:mi></mml:mrow></mml:msub><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mo>}</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(6)</label></disp-formula>
<p>Fourth, battery availability constraints (<xref ref-type="disp-formula" rid="EQ7">Equation 7</xref>) approximate the finite battery inventory and charging cycle without explicitly tracking individual battery states&#x02014;a simplification to maintain model tractability. For station <italic>j</italic> with capacity <italic>N</italic><sub><italic>j</italic></sub> (maximum stored batteries), the cumulative swaps within the charging time window <italic>T</italic><sub><italic>c</italic></sub> (time to recharge a depleted battery) must not exceed <italic>N</italic><sub><italic>j</italic></sub>:</p>
<disp-formula id="EQ7"><mml:math id="M8"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mtable columnalign="left"><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>I</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mn>0</mml:mn></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x02264;</mml:mo><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02200;</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:mi>J</mml:mi></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mtext class="textrm" mathvariant="normal">if</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x0003C;</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>I</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mstyle displaystyle="true"><mml:msubsup><mml:mrow><mml:mo>&#x02211;</mml:mo></mml:mrow><mml:mrow><mml:mi>t</mml:mi><mml:mo>-</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub><mml:mo>&#x0002B;</mml:mo><mml:mn>1</mml:mn></mml:mrow><mml:mrow><mml:mi>t</mml:mi></mml:mrow></mml:msubsup></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>x</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>j</mml:mi><mml:mo>,</mml:mo><mml:mi>t</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x02264;</mml:mo><mml:msub><mml:mrow><mml:mi>N</mml:mi></mml:mrow><mml:mrow><mml:mi>j</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02200;</mml:mo><mml:mi>j</mml:mi><mml:mo>&#x02208;</mml:mo><mml:mrow><mml:mo>{</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:mi>J</mml:mi></mml:mrow><mml:mo>}</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd><mml:mtd><mml:mtext class="textrm" mathvariant="normal">if</mml:mtext><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>t</mml:mi><mml:mo>&#x02265;</mml:mo><mml:msub><mml:mrow><mml:mi>T</mml:mi></mml:mrow><mml:mrow><mml:mi>c</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:mrow></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(7)</label></disp-formula>
<p>This formulation implicitly enforces that each battery used at time <italic>t</italic> becomes available again at <italic>t</italic>&#x0002B;<italic>T</italic><sub><italic>c</italic></sub>, aligning with the charging cycle while avoiding combinatorial complexity from per-battery tracking.</p>
<p>Batteries are charged in batches at a fixed rate (0.95%/minute), so the charging time window <italic>T</italic><sub><italic>c</italic></sub> (time to fully recharge a depleted battery) can be precomputed as a deterministic parameter (e.g., 105 minutes for a 100% SoC recovery, consistent with the 700&#x02013;800 kWh battery capacity). Spare batteries are maintained as a shared pool (5 per station) rather than assigned to specific trucks. Thus, the cumulative number of swaps within <italic>T</italic><sub><italic>c</italic></sub> directly reflects the maximum number of batteries that can be cycled (used &#x02192; charged &#x02192; reused), making the constraint (<xref ref-type="disp-formula" rid="EQ7">Equation 7</xref>) a reasonable proxy for actual inventory limits. Furthermore, the constraint (<xref ref-type="disp-formula" rid="EQ7">Equation 7</xref>) could be reformulated as a dynamic one. For a battery-swapping station with multiple charging batteries, their remaining charging times can be estimated and sorted in ascending order. As operational time progresses (e.g., across the MILP&#x00027;s discrete time slices; Section 3.1), each time the time exceeds a remaining charging time in the sorted list, the station&#x00027;s available battery count increases by 1&#x02013;aligning with the study&#x00027;s operational context.</p>
<p>The proposed MILP model (<xref ref-type="disp-formula" rid="EQ1">Equations 1</xref>&#x02013;<xref ref-type="disp-formula" rid="EQ7">7</xref>) is solved using the open-source solver SCIP (<xref ref-type="bibr" rid="B14">Vigerske and Gleixner, 2016</xref>), which outputs the battery replacement schedule&#x02013;specifically, the optimal battery-swapping station assignment and service order for each mining truck. To generate the model&#x00027;s input (i.e., the predicted battery-swapping demand times during the first cycle), we propose a state-based prediction method tailored to the operational logic of mining trucks. Unlike passenger electric vehicles, where battery-swapping demand can arise continuously, mining trucks can only request battery replacement after completing unloading tasks. This discrete demand pattern arises because trucks must first finish unloading and travel from the unloading zone to a battery-swapping station. The state-based method estimates available battery-swapping moments based on a truck&#x00027;s current operational state:</p>
<list list-type="bullet">
<list-item><p>Unloading or queuing for unloading: The available time is calculated as <italic>T</italic><sub><italic>re</italic></sub>&#x0002B;<italic>T</italic><sub><italic>UL</italic>&#x02192;<italic>S</italic></sub> (where <italic>T</italic><sub><italic>re</italic></sub> is the remaining unloading time and <italic>T</italic><sub><italic>UL</italic>&#x02192;<italic>S</italic></sub> is the travel time to a station) or <italic>n</italic>&#x000B7;<italic>T</italic><sub><italic>ul</italic></sub>&#x0002B;<italic>T</italic><sub><italic>UL</italic>&#x02192;<italic>S</italic></sub> (where <italic>n</italic> is the number of trucks queuing ahead and <italic>T</italic><sub><italic>ul</italic></sub> is the average unloading time).</p></list-item>
<list-item><p>Loading, traveling to loading platform, or returning to unloading zone: The truck&#x00027;s state is first reduced to the equivalent unloading-area state (e.g., estimating remaining time until unloading completion), after which the above logic applies.</p></list-item>
</list>
<p>This approach ensures that battery-swapping demands align with the discrete operational milestones inherent to mining workflows, thereby enabling accurate input generation for the MILP model. The state-based method generates outputs such as those illustrated in <xref ref-type="table" rid="T2">Table 2</xref>, where mandatory battery-swapping times correspond to instances when the battery is depleted (i.e., the remaining state of charge, SoC, is insufficient to sustain another operational cycle), and <italic>n</italic><sub>1</sub>, <italic>n</italic><sub>2</sub>, &#x02026;, <italic>n</italic><sub><italic>I</italic></sub> denote the number of available battery-swapping times for trucks 1, 2, &#x02026;, <italic>I</italic>, respectively. The final input to the MILP model is an aggregated set of predicted demand times, structured as [<italic>t</italic><sub>1, 2</sub>, <italic>t</italic><sub>2, 1</sub>, &#x022EF;&#x02009;, <italic>t</italic><sub><italic>i, l</italic></sub>, &#x022EF;&#x02009;, <italic>t</italic><sub><italic>I, k</italic></sub>] from <xref ref-type="table" rid="T2">Table 2</xref>. Here, <italic>t</italic><sub><italic>i, l</italic></sub> represents the <italic>l</italic>&#x02212;<italic>th</italic> battery-swapping time scheduled for truck <italic>i</italic>, and <italic>t</italic><sub><italic>I, k</italic></sub> denotes the <italic>k</italic>&#x02212;<italic>th</italic> time allocated to truck <italic>I</italic>, where <italic>I</italic> is the total number of trucks in the fleet.</p>
<table-wrap position="float" id="T2">
<label>Table 2</label>
<caption><p>The battery-swapping demands for all trucks.</p></caption>
<table frame="box" rules="all">
<thead>
<tr>
<th valign="top" align="left"><bold>Truck ID</bold></th>
<th valign="top" align="center" colspan="3"><bold>Available moments</bold></th>
<th valign="top" align="center"><bold>Mandatory moments</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">truck 1</td>
<td valign="top" align="center"><italic>t</italic><sub>1, 1</sub></td>
<td valign="top" align="center"><italic>t</italic><sub>1, 2</sub></td>
<td valign="top" align="center">&#x022EF;&#x02009;</td>
<td valign="top" align="center"><italic>t</italic><sub>1,<sub><italic>n</italic></sub><sub>1</sub></sub></td>
</tr>
<tr>
<td valign="top" align="left">truck 2</td>
<td valign="top" align="center"><italic>t</italic><sub>2, 1</sub></td>
<td valign="top" align="center"><italic>t</italic><sub>2, 2</sub></td>
<td valign="top" align="center">&#x022EF;&#x02009;</td>
<td valign="top" align="center"><italic>t</italic><sub>2,<sub><italic>n</italic></sub><sub>2</sub></sub></td>
</tr>
<tr>
<td valign="top" align="left">&#x022EE;</td>
<td valign="top" align="center">&#x022EE;</td>
<td valign="top" align="center">&#x022EE;</td>
<td valign="top" align="center">&#x022F1;</td>
<td valign="top" align="center">&#x022EE;</td>
</tr>
<tr>
<td valign="top" align="left">truck <italic>I</italic></td>
<td valign="top" align="center"><italic>t</italic><sub><italic>I</italic>, 1</sub></td>
<td valign="top" align="center"><italic>t</italic><sub><italic>I</italic>, 2</sub></td>
<td valign="top" align="center">&#x022EF;&#x02009;</td>
<td valign="top" align="center"><italic>t</italic><sub><italic>I</italic>,<sub><italic>n</italic></sub><sub><italic>I</italic></sub></sub></td>
</tr></tbody>
</table>
</table-wrap>
<p>First, we evaluate the potential of the MILP model to mitigate queuing delays for mining trucks at battery-swapping stations. To this end, input parameters for the MILP model are generated using a normal distribution with a mean of 210 (i.e., 3.5 hours after the current time). By systematically varying the standard deviation (&#x003C3;), we simulate distinct temporal distributions of battery-swapping demand times across a fleet of 36 mining trucks. Smaller &#x003C3; values produce more clustered demand profiles, emulating scenarios with synchronized battery depletion&#x02013;a critical bottleneck in mining operations. <xref ref-type="fig" rid="F8">Figure 8</xref> compares the average waiting times achieved by MILP-optimized scheduling against the baseline First-Come-First-Served (FCFS) strategy. The results reveal a maximum reduction of 15 minutes in queuing time at &#x003C3; &#x0003D; 35, whereas the minimum reduction of 12 minutes occurs under the most concentrated demand scenario (&#x003C3; &#x0003D; 0.5). This demonstrates that merely enhancing intrinsic station throughput (e.g., faster battery swapping) cannot resolve synchronization-induced congestion, underscoring the necessity of off-peak scheduling to address systemic inefficiencies.</p>
<fig position="float" id="F8">
<label>Figure 8</label>
<caption><p>The improvement by MILP with FCFS as baseline.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0008.tif">
<alt-text content-type="machine-generated">Bar and line chart comparing average waiting times across different standard deviations. Blue and green bars represent FCFS and MILP respectively. Red and orange lines indicate trends in waiting times.</alt-text>
</graphic>
</fig>
</sec>
<sec>
<label>3.2</label>
<title>Bayesian optimization</title>
<p><xref ref-type="fig" rid="F8">Figure 8</xref> reveals that greater dispersion in battery-swapping demand distributions correlates with shorter average waiting times, regardless of scheduling strategy (MILP or FCFS). This observation suggests an intuitive off-peak strategy: by systematically exploring all feasible battery-swapping time combinations derived from <xref ref-type="table" rid="T2">Table 2</xref>, one could theoretically identify the optimal schedule with minimal queuing time. However, this approach is computationally intractable due to combinatorial explosion. To illustrate, consider a fleet of 40 mining trucks, each with 4 candidate battery-swapping times before battery depletion. The total number of combinations is 4<sup>40</sup>&#x02248;10<sup>24</sup>. Assuming an MILP solver requires 1 second per combination, exhaustive evaluation would demand 10<sup>24</sup> seconds (&#x02248;10<sup>16</sup> years)&#x02013;a timescale exceeding practical feasibility.</p>
<p>Bayesian optimization is a data-driven method for efficiently optimizing black-box functions (e.g., complex systems with unknown analytical forms). Its core mechanism involves constructing a probabilistic surrogate model&#x02014;typically a Gaussian process&#x02014;to approximate the objective function. By iteratively updating this model with prior observations, Bayesian optimization computes a posterior distribution that balances exploration (sampling uncertain regions) and exploitation (refining known promising regions), thereby guiding the search toward globally optimal solutions. To address the combinatorial explosion inherent in scheduling problems, Bayesian optimization circumvents exhaustive enumeration of all combinations. Instead, it leverages the surrogate model&#x00027;s probability estimates to prioritize candidate solutions with higher likelihoods of optimality. This is achieved through an acquisition function (e.g., Expected Improvement), which intelligently selects the next evaluation point based on the posterior distribution. By focusing computational resources on high-potential regions, Bayesian optimization dramatically reduces the effective search space, mitigating the computational burden caused by combinatorial explosion (e.g., reducing 10<sup>24</sup> combinations to a tractable subset).</p>
<p>For the off-peak battery-swapping scheduling problem under study, the objective function approximated via Bayesian optimization is formally defined as:</p>
<disp-formula id="EQ8"><mml:math id="M9"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mo>:</mml:mo><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow><mml:mo>&#x02192;</mml:mo><mml:mi>&#x0211D;</mml:mi></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(8)</label></disp-formula>
<p>where <inline-formula><mml:math id="M10"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula> denotes the feasible solution space of battery-swapping time combinations, and &#x0211D; maps to the real-valued output of either the MILP model (<xref ref-type="disp-formula" rid="EQ1">Equations 1</xref>&#x02013;<xref ref-type="disp-formula" rid="EQ7">7</xref>) or the FCFS queuing time calculation. One sample in <inline-formula><mml:math id="M11"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula> can be expressed as <inline-formula><mml:math id="M12"><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle><mml:mo>=</mml:mo><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>1</mml:mn><mml:mo>,</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn><mml:mo>,</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mo>,</mml:mo><mml:mi>l</mml:mi></mml:mrow></mml:msub><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msub><mml:mrow><mml:mi>t</mml:mi></mml:mrow><mml:mrow><mml:mi>I</mml:mi><mml:mo>,</mml:mo><mml:mi>k</mml:mi></mml:mrow></mml:msub></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:math></inline-formula>, with <italic>t</italic><sub><italic>i, l</italic></sub> representing the <italic>l</italic>&#x02212;<italic>th</italic> battery-swapping time scheduled for truck <italic>i</italic> as defined in <xref ref-type="table" rid="T2">Table 2</xref>. By treating <italic>f</italic> as a computationally expensive black-box function, Bayesian optimization constructs a probabilistic surrogate model (e.g., Gaussian process) trained on prior evaluations to guide the search toward high-potential regions of <inline-formula><mml:math id="M13"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula>, thereby avoiding exhaustive enumeration of all combinations (e.g., 4<sup>40</sup>) and balancing exploration-exploitation trade-offs through acquisition functions like Expected Improvement (EI) and Upper Confidence Bound (UCB).</p>
<p>At the core of the Gaussian process (GP) lies the assumption that the function values of any finite set of points follow a joint Gaussian distribution, formally expressed as:</p>
<disp-formula id="EQ9"><mml:math id="M14"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0007E;</mml:mo><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>GP</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msub><mml:mrow><mml:mi>&#x003BC;</mml:mi></mml:mrow><mml:mrow><mml:mi>f</mml:mi></mml:mrow></mml:msub><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>,</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(9)</label></disp-formula>
<p>where &#x003BC;<sub><italic>f</italic></sub>(<bold>c</bold>) is the mean function (set to 0 due to the absence of prior knowledge, ensuring reliance on the kernel for localized feature capture) and <italic>k</italic>(<bold>c</bold>, <bold>c</bold>&#x02032;) is the covariance kernel quantifying similarity between combinations <bold>c</bold> and <bold>c</bold>&#x02032;. We employ the Mat&#x000E9;rn 5/2 kernel, defined as:</p>
<disp-formula id="EQ10"><mml:math id="M15"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mn>1</mml:mn><mml:mo>&#x0002B;</mml:mo><mml:msqrt><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msqrt><mml:mi>r</mml:mi><mml:mo>&#x0002B;</mml:mo><mml:mfrac><mml:mrow><mml:mn>5</mml:mn><mml:msup><mml:mrow><mml:mi>r</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mrow><mml:mn>3</mml:mn></mml:mrow></mml:mfrac><mml:msup><mml:mrow><mml:mi>e</mml:mi></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:msqrt><mml:mrow><mml:mn>5</mml:mn></mml:mrow></mml:msqrt><mml:mi>r</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(10)</label></disp-formula>
<p>where <inline-formula><mml:math id="M16"><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mo>&#x02225;</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle><mml:mo>-</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>&#x02032;</mml:mi></mml:mrow></mml:msup><mml:msub><mml:mrow><mml:mo>&#x02225;</mml:mo></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msub></mml:math></inline-formula> is the L2-norm distance between <bold>c</bold> and <bold>c</bold>&#x02032;. This kernel&#x00027;s twice-differentiable properties balance smoothness and flexibility, critical for modeling queuing time trends and localized demand variations in battery-swapping scheduling. Another advantage of this kernel function is that there is no need for hyperparameter optimization.</p>
<p>With the Gaussian process framework (<xref ref-type="disp-formula" rid="EQ9">Equations 9</xref>&#x02013;<xref ref-type="disp-formula" rid="EQ10">10</xref>), the queuing time for a mining truck fleet under a new battery-swapping schedule can be probabilistically predicted using observed historical data. Let <inline-formula><mml:math id="M17"><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mo>&#x02282;</mml:mo><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula> denote a set of <italic>m</italic> observed combinations of battery-swapping times, and <bold>f</bold> &#x0003D; [<italic>f</italic>(<bold>c</bold><sup>1</sup>), <italic>f</italic>(<bold>c</bold><sup>2</sup>), &#x02026;, <italic>f</italic>(<bold>c</bold><sup><italic>m</italic></sup>)]<italic>T</italic> represent their corresponding queuing times (computed via MILP or FCFS). For a candidate combination <bold>c</bold><sup>&#x0002A;</sup>, Bayesian inference yields the posterior distribution of its predicted queuing time::</p>
<disp-formula id="EQ11"><mml:math id="M18"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>P</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>|</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>f</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>N</mml:mi></mml:mstyle></mml:mrow><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mi>&#x003BC;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(11)</label></disp-formula>
<p>where the posterior mean &#x003BC;(<bold>c</bold><sup>&#x0002A;</sup>) and variance &#x003C3;(<bold>c</bold><sup>&#x0002A;</sup>) are derived as:</p>
<disp-formula id="EQ12"><mml:math id="M19"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:mi>&#x003BC;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>k</mml:mtext></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>K</mml:mtext></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003C3;</mml:mi></mml:mrow><mml:mrow><mml:mi>o</mml:mi><mml:mi>b</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mstyle mathvariant="bold"><mml:mtext>E</mml:mtext></mml:mstyle></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mstyle mathvariant="bold"><mml:mtext>f</mml:mtext></mml:mstyle></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>-</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>k</mml:mtext></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>K</mml:mtext></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:msubsup><mml:mrow><mml:mi>&#x003C3;</mml:mi></mml:mrow><mml:mrow><mml:mi>o</mml:mi><mml:mi>b</mml:mi></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msubsup><mml:mstyle mathvariant="bold"><mml:mtext>E</mml:mtext></mml:mstyle></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mrow><mml:mrow><mml:mo>-</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mstyle mathvariant="bold"><mml:mtext>k</mml:mtext></mml:mstyle><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(12)</label></disp-formula>
<p>Here:</p>
<list list-type="bullet">
<list-item><p><bold>k</bold>(<bold>c</bold><sup>&#x0002A;</sup>, <bold>C</bold>) &#x0003D; [<italic>k</italic>(<bold>c</bold><sup>&#x0002A;</sup>, <bold>c</bold><sup>1</sup>), &#x02026;, <italic>k</italic>(<bold>c</bold><sup>&#x0002A;</sup>, <bold>c</bold><sup><italic>m</italic></sup>)] is a 1 &#x000D7; <italic>m</italic> row vector of kernel similarities between <bold>c</bold><sup>&#x0002A;</sup> and observed combinations;</p></list-item>
<list-item><p><bold>K</bold>(<bold>C</bold>, <bold>C</bold>) is the <italic>m</italic>&#x000D7;<italic>m</italic> kernel matrix with entries <inline-formula><mml:math id="M21"><mml:msub><mml:mrow><mml:mi>k</mml:mi></mml:mrow><mml:mrow><mml:mi>p</mml:mi><mml:mo>,</mml:mo><mml:mi>q</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mi>k</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>p</mml:mi></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>q</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:math></inline-formula>;</p></list-item>
<list-item><p>&#x003C3;<sub>ob</sub> is the observation noise standard deviation, set to 0 in this work to reflect deterministic queuing time evaluations;</p></list-item>
<list-item><p><bold>E</bold> is a <italic>m</italic>&#x000D7;<italic>m</italic> identity matrix;</p></list-item>
<list-item><p><bold>k</bold>(<bold>C</bold>, <bold>c</bold><sup>&#x0002A;</sup>) &#x0003D; <bold>k</bold>(<bold>c</bold><sup>&#x0002A;</sup>, <bold>C</bold>)<italic>T</italic> is a column vector.</p></list-item>
</list>
<p>Once the Gaussian process surrogate model is established, the Upper Confidence Bound (UCB) acquisition function guides the selection of the next candidate solution to evaluate. The UCB is defined as:</p>
<disp-formula id="EQ13"><mml:math id="M22"><mml:mtable class="eqnarray" columnalign="left"><mml:mtr><mml:mtd><mml:mi>U</mml:mi><mml:mi>C</mml:mi><mml:mi>B</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>=</mml:mo><mml:mi>&#x003BC;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>&#x0002B;</mml:mo><mml:mi>&#x003BA;</mml:mi><mml:mo>&#x000B7;</mml:mo><mml:mi>&#x003C3;</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mo>*</mml:mo></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(13)</label></disp-formula>
<p>where &#x003BC;(<bold>c</bold><sup>&#x0002A;</sup>) and &#x003C3;(<bold>c</bold><sup>&#x0002A;</sup>) are the posterior mean and standard deviation from <xref ref-type="disp-formula" rid="EQ12">Equation 12</xref>, and &#x003BA; &#x0003D; 2.5 is a hyperparameter balancing exploration (prioritizing uncertain regions) and exploitation (refining known optima). At each Bayesian optimization iteration:</p>
<list list-type="bullet">
<list-item><p>A batch of candidate combinations {<bold>c</bold><sup>&#x0002A;</sup>} is sampled from <inline-formula><mml:math id="M23"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula>.</p></list-item>
<list-item><p>Their UCB values are computed via <xref ref-type="disp-formula" rid="EQ15">Equation 14</xref>.</p></list-item>
<list-item><p>The combination <bold>c</bold><sup>iter &#x0003D; 1</sup> with the highest UCB is selected, evaluated (via MILP or FCFS), and added to the observation set:</p></list-item>
</list>
<disp-formula id="EQ14"><mml:math id="M24"><mml:mtable columnalign='center'><mml:mtr><mml:mtd><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>C</mml:mi></mml:mstyle><mml:mrow><mml:mtext>iter</mml:mtext><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy='false'>[</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mn>1</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mn>2</mml:mn></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mi>m</mml:mi></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mrow><mml:mtext>iter</mml:mtext><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo stretchy='false'>]</mml:mo><mml:mo>,</mml:mo><mml:mtext>&#x02003;</mml:mtext><mml:msub><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>f</mml:mi></mml:mstyle><mml:mrow><mml:mtext>iter</mml:mtext><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mo stretchy='false'>[</mml:mo><mml:mi>f</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mn>1</mml:mn></mml:msup><mml:mo stretchy='false'>)</mml:mo><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mn>2</mml:mn></mml:msup><mml:mo stretchy='false'>)</mml:mo><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mi>m</mml:mi></mml:msup><mml:mo stretchy='false'>)</mml:mo><mml:mo>,</mml:mo><mml:mtext>&#x02009;</mml:mtext><mml:mi>f</mml:mi><mml:mo stretchy='false'>(</mml:mo><mml:msup><mml:mstyle mathvariant='bold' mathsize='normal'><mml:mi>c</mml:mi></mml:mstyle><mml:mrow><mml:mtext>iter</mml:mtext><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo stretchy='false'>)</mml:mo><mml:msup><mml:mo stretchy='false'>]</mml:mo><mml:mi>T</mml:mi></mml:msup><mml:mo>.</mml:mo></mml:mtd></mml:mtr></mml:mtable></mml:math></disp-formula>
<p>After N iterations, the augmented datasets become:</p>
<disp-formula id="EQ15"><mml:math id="M25"><mml:mtable class="eqnarray" columnalign="center"><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>C</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo>]</mml:mo></mml:mrow></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:msub><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>f</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:msub><mml:mo>=</mml:mo><mml:mrow><mml:mo>[</mml:mo></mml:mrow><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x02026;</mml:mo><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>m</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>1</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mn>2</mml:mn></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mo>,</mml:mo><mml:mo>&#x022EF;</mml:mo><mml:mspace width="0.3em" class="thinspace"/><mml:mo>,</mml:mo></mml:mtd></mml:mtr><mml:mtr><mml:mtd><mml:mi>f</mml:mi><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mrow><mml:msup><mml:mrow><mml:mstyle mathvariant="bold"><mml:mtext>c</mml:mtext></mml:mstyle></mml:mrow><mml:mrow><mml:mi>i</mml:mi><mml:mi>t</mml:mi><mml:mi>e</mml:mi><mml:mi>r</mml:mi><mml:mo>=</mml:mo><mml:mi>N</mml:mi></mml:mrow></mml:msup></mml:mrow><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:msup><mml:mrow><mml:mo>]</mml:mo></mml:mrow><mml:mrow><mml:mi>T</mml:mi></mml:mrow></mml:msup></mml:mtd></mml:mtr></mml:mtable></mml:math><label>(14)</label></disp-formula>
<p>The optimal battery-swapping schedule <bold>c</bold><sup>best</sup> corresponds to the combination with the minimum queuing time in <bold>f</bold><sub>iter &#x0003D; <italic>N</italic></sub>. This iterative process efficiently navigates the combinatorial space <inline-formula><mml:math id="M28"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula>, avoiding exhaustive evaluation of all 4<sup>40</sup> combinations while balancing data-driven exploration and exploitation.</p>
<p>With both the inner MILP model (<xref ref-type="disp-formula" rid="EQ1">Equations 1</xref>&#x02013;<xref ref-type="disp-formula" rid="EQ7">7</xref>) and the outer Bayesian optimization (<xref ref-type="disp-formula" rid="EQ9">Equations 9</xref>&#x02013;<xref ref-type="disp-formula" rid="EQ15">14</xref>) rigorously defined, the proposed hierarchical optimization framework is systematically established (detailed in <xref ref-type="other" rid="algorithm_1">Algorithm 1</xref>), as illustrated in <xref ref-type="fig" rid="F9">Figure 9</xref>. To adapt Bayesian optimization&#x02014;originally designed for continuous decision variables&#x02014;to the discrete combinatorial space <inline-formula><mml:math id="M29"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula>, a binary search-based discretization module (detailed in <xref ref-type="other" rid="algorithm_2">Algorithm 2</xref>) is integrated. This module maps continuous samples generated by the Gaussian process surrogate model to the nearest valid discrete battery-swapping times (<xref ref-type="table" rid="T2">Table 2</xref>), ensuring feasibility while preserving the probabilistic exploration-exploitation balance. The hierarchical integration of MILP (exact queuing time calculation) and Bayesian optimization (global search guidance) enables near-optimal scheduling with drastically reduced computational effort.</p>
<statement content-type="algorithm" id="algorithm_1">
<label>Algorithm 1</label>
<p>Bayesian optimization for off-peak battery swapping scheduling.
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0015.tif"/>
</p>
</statement>
<fig position="float" id="F9">
<label>Figure 9</label>
<caption><p>The Optimization-based heriarchical framework, <italic>its</italic> denotes the total iteration steps.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0009.tif">
<alt-text content-type="machine-generated">Flowchart illustrating a process that starts with a &#x0201C;state-based prediction method.&#x0201D; It includes steps like &#x0201C;available battery-swapping times,&#x0201D; &#x0201C;initial sampling,&#x0201D; &#x0201C;map to discrete space,&#x0201D; and &#x0201C;fitting Gaussian process.&#x0201D; The diagram incorporates Bayesian optimization, loops involving &#x0201C;MILP or FCFS,&#x0201D; and decisions with arrows indicating progression. The process concludes at &#x0201C;end.&#x0201D; </alt-text>
</graphic>
</fig>
<statement content-type="algorithm" id="algorithm_2">
<label>Algorithm 2</label>
<p>Binary search-based discretization module.
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0016.tif"/>
</p>
</statement>
</sec>
</sec>
<sec id="s4">
<label>4</label>
<title>Case studies</title>
<sec>
<label>4.1</label>
<title>Case study 1: single loading platform and unloading zone</title>
<p>The first scenario simulates battery-swapping operations in an open-pit mine with one loading platform (equipped with two excavators) and one unloading zone (serving a single crusher). Key operational parameters include:</p>
<list list-type="bullet">
<list-item><p>Haul Road Configuration: A 4-km bidirectional route connecting the loading and unloading zones.</p></list-item>
<list-item><p>State of Charge (SoC) consumption:
<list list-type="simple">
<list-item><p>&#x02013; Unloading travel (empty haul): 5% SoC consumed over 10 min.</p></list-item>
<list-item><p>&#x02013; Loading travel (loaded haul): 10% SoC consumed over 15 min.</p></list-item></list></p></list-item>
<list-item><p>Processing Times:
<list list-type="simple">
<list-item><p>&#x02013; Unloading: 3 min (average) at the crusher.</p></list-item>
<list-item><p>&#x02013; Loading: 5 min (average) per excavator.</p></list-item>
</list></p></list-item>
<list-item><p>Initial fleet state:
<list list-type="simple">
<list-item><p>&#x02013; Each excavator is actively loading one haul truck, with 3 trucks queued per excavator (total 6 trucks waiting to load).</p></list-item>
<list-item><p>&#x02013; 28 trucks are in transit (empty-hauling) between zones.</p></list-item>
</list></p></list-item>
</list>
<p>These settings are based on an actual open-pit coal mine. Using these inputs, the state-based prediction module (Section 3) generates a discrete battery-swapping schedule (<xref ref-type="table" rid="T2">Table 2</xref>) by simulating truck trajectories, energy consumption, and queuing dynamics. This table provides all feasible battery-swapping time windows for the fleet, serving as input to the hierarchical optimization framework (<xref ref-type="fig" rid="F9">Figure 9</xref>) for minimizing queuing delays.</p>
<p>A systematic sensitivity analysis evaluates the impact of initial sampling points (init_points) and Bayesian optimization iterations (n_iter) on solution quality. As shown in <xref ref-type="table" rid="T3">Table 3</xref>, increasing init_points enhances solution optimality initially but plateaus beyond a threshold (e.g., init_points=30), suggesting sufficient exploration of the combinatorial space <inline-formula><mml:math id="M48"><mml:mrow><mml:mstyle mathvariant="script"><mml:mi>C</mml:mi></mml:mstyle></mml:mrow></mml:math></inline-formula>. In contrast, increasingn_iter beyond 5 yields diminishing returns, indicating rapid convergence of the Gaussian process surrogate model. Based on these findings, subsequent experiments adoptinit_points=30 and n_iter=5, balancing computational efficiency with solution quality. This configuration ensures robust exploration-exploitation trade-offs while aligning with the hierarchical framework&#x00027;s goal of mitigating combinatorial explosion (Section 3).</p>
<table-wrap position="float" id="T3">
<label>Table 3</label>
<caption><p>The effect of the number of initial sampling points (init_points) and iteration steps (n_iter) on solution in case 1.</p></caption>
<table frame="box" rules="all">
<thead>
<tr>
<th valign="top" align="left"><bold>init_points</bold></th>
<th valign="top" align="center"><bold>n_iter</bold></th>
<th valign="top" align="center"><bold>average queuing time (minutes)</bold></th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">20</td>
<td valign="top" align="center">5</td>
<td valign="top" align="center">45.3</td>
</tr>
<tr>
<td valign="top" align="left">30</td>
<td valign="top" align="center">5</td>
<td valign="top" align="center">41.3</td>
</tr>
<tr>
<td valign="top" align="left">50</td>
<td valign="top" align="center">5</td>
<td valign="top" align="center">41.3</td>
</tr>
<tr>
<td valign="top" align="left">100</td>
<td valign="top" align="center">5</td>
<td valign="top" align="center">41.3</td>
</tr>
<tr>
<td valign="top" align="left">100</td>
<td valign="top" align="center">10</td>
<td valign="top" align="center">41.3</td>
</tr>
<tr>
<td valign="top" align="left">100</td>
<td valign="top" align="center">20</td>
<td valign="top" align="center">41.3</td>
</tr>
<tr>
<td valign="top" align="left">100</td>
<td valign="top" align="center">40</td>
<td valign="top" align="center">41.3</td>
</tr></tbody>
</table>
</table-wrap>
<p><xref ref-type="fig" rid="F10">Figure 10</xref> illustrates the scheduling outcomes generated by the Hierarchical Framework (Section 3). The results demonstrate that the framework proactively schedules more than 50% of trucks for battery swaps prior to SoC depletion, effectively dispersing demand peaks and reducing total queuing time. The optimized schedule achieves an average queuing time of 41.3 min at the battery-swapping station, representing a 65% reduction compared to the MILP baseline (116.4 min) and a 68% improvement over the FCFS strategy (128.4 min) when both use mandatory battery-depletion moments (<xref ref-type="table" rid="T2">Table 2</xref>) as input. This performance underscores the framework&#x00027;s ability to balance proactive scheduling (exploiting demand dispersion) with computational efficiency (avoiding the 4<sup>40</sup> combinatorial search space), validating its superiority over conventional rule-based or exhaustive optimization approaches.</p>
<fig position="float" id="F10">
<label>Figure 10</label>
<caption><p>Case 1: the scheduling results generated by the hierarchical framework, the red markers indicate the battery-swapping moments, the orange markers denote the arrival times of trucks at the station, and the blue markers represent the battery depletion moments (when the battery&#x00027;s State of Charge reaches critical levels).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0010.tif">
<alt-text content-type="machine-generated">Graph displaying a grid with truck IDs on the vertical axis and time slice indices on the horizontal axis. The background is green with scattered red squares indicating data points, which are more concentrated between time slices 35 and 45.</alt-text>
</graphic>
</fig>
</sec>
<sec>
<label>4.2</label>
<title>Case study 2: dual loading platforms and unloading zone</title>
<p>The second scenario models battery-swapping operations in an open-pit mine with two loading platforms (each equipped with two excavators) and one unloading zone (serving a single crusher). Key operational parameters include:</p>
<list list-type="bullet">
<list-item><p>Haul road configuration:
<list list-type="simple">
<list-item><p>&#x02013; A 3-km bidirectional route connects the first loading platform to the unloading zone.</p></list-item>
<list-item><p>&#x02013; A 5-km bidirectional route links the second loading platform to the unloading zone.</p></list-item>
</list></p></list-item>
<list-item><p>State of Charge (SoC) consumption:
<list list-type="simple">
<list-item><p>&#x02013; Unloading Travel (Empty Haul) to the first loading platform: 3.75% SoC consumed over 7.5 min.</p></list-item>
<list-item><p>&#x02013; Unloading Travel (Empty Haul) to the second loading platform: 6.25% SoC consumed over 12.5 min.</p></list-item>
<list-item><p>&#x02013; Loading Travel (Loaded Haul) from the first loading platform: 7.5% SoC consumed over 11.25 min.</p></list-item>
<list-item><p>&#x02013; Loading Travel (Loaded Haul) from the second loading platform: 12.5% SoC consumed over 18.75 min.</p></list-item>
</list></p></list-item>
<list-item><p>Processing times:
<list list-type="simple">
<list-item><p>&#x02013; Unloading: 1 min (average) at the crusher.</p></list-item>
<list-item><p>&#x02013; Loading: 5 min (average) per excavator.</p></list-item>
</list></p></list-item>
<list-item><p>Initial Fleet State:
<list list-type="simple">
<list-item><p>&#x02013; Each excavator is actively loading one haul truck, with 3 trucks queued per excavator (total 12 trucks waiting to load).</p></list-item>
<list-item><p>&#x02013; 10 trucks are in transit (empty-hauling) between the unloading zone and the first loading platform.</p></list-item>
<list-item><p>&#x02013; 10 trucks are in transit (empty-hauling) between the unloading zone and the second loading platform.</p></list-item>
</list></p></list-item>
<list-item><p>Production scheduling:
<list list-type="simple">
<list-item><p>&#x02013; Trucks alternate loading platforms for consecutive trips: if a trip loads at the first platform, the subsequent trip must load at the second platform to balance resource utilization.</p></list-item>
</list></p></list-item>
</list>
<p>Using these inputs, the state-based prediction module (Section 3) simulates truck trajectories, energy consumption, and queuing dynamics to generate a discrete battery-swapping schedule (<xref ref-type="table" rid="T2">Table 2</xref>). This schedule enumerates all feasible battery-swapping time windows, which serve as input to the hierarchical optimization framework (<xref ref-type="fig" rid="F9">Figure 9</xref>) for minimizing queuing delays through proactive demand dispersion.</p>
<p><xref ref-type="fig" rid="F11">Figure 11</xref> illustrates the scheduling outcomes generated by the Hierarchical Framework (Section 3). Compared to Case Study 1 (single loading platform), the framework proactively schedules 75% of trucks (versus 50% in Case 1) for battery swaps prior to SoC depletion, achieving enhanced demand dispersion across the dual-platform configuration. The optimized schedule attains an average queuing time of 26.4 min at the battery-swapping station, demonstrating a 78% reduction compared to the MILP baseline (122.5 min) and an 80% improvement over the FCFS strategy (132.1 min), both of which rely on mandatory battery-depletion moments (<xref ref-type="table" rid="T2">Table 2</xref>).</p>
<fig position="float" id="F11">
<label>Figure 11</label>
<caption><p>Case 2: the scheduling results generated by the hierarchical framework, the red markers indicate the battery-swapping moments, the orange markers denote the arrival times of trucks at the station, and the blue markers represent the battery depletion moments (when the battery&#x00027;s State of Charge reaches critical levels).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0011.tif">
<alt-text content-type="machine-generated">Chart displaying truck activity over time slices. The y-axis represents truck IDs, and the x-axis represents time slices. Green indicates inactivity and red indicates activity, with clusters of red appearing at various points, notably around time slice 55.</alt-text>
</graphic>
</fig>
</sec>
<sec>
<label>4.3</label>
<title>Case study 3: dual loading platforms with DES-driven inputs</title>
<p>The third scenario extends the battery-swapping operations to a Discrete Event Simulation (DES)-integrated environment, maintaining the dual loading platforms (each with two excavators) and a single unloading zone (serving one crusher). Key enhancements include:</p>
<list list-type="bullet">
<list-item><p>Input Generation: Available battery-swapping moments are derived from DES instead of the state-based prediction module, enabling realistic production scheduling with stochastic variations (e.g., delays, resource contention).</p></list-item>
<list-item><p>Scalability Validation: The hierarchical framework&#x00027;s scheduling results are tested within DES to verify robustness under dynamic, large-scale operational conditions.</p></list-item>
</list>
<p>Key operational parameters include:</p>
<list list-type="bullet">
<list-item><p>State of Charge (SoC) consumption:
<list list-type="simple">
<list-item><p>&#x02013; Empty Haul to Platform 1: 3.0% SoC consumed over 12 minutes (3 km route).</p></list-item>
<list-item><p>&#x02013; Empty Haul to Platform 2: 3.75% SoC consumed over 15 minutes (5 km route).</p></list-item>
<list-item><p>&#x02013; Loaded Haul from Platform 1: 7.33% SoC consumed over 22 minutes.</p></list-item>
<list-item><p>&#x02013; Loaded Haul from Platform 2: 8.33% SoC consumed over 25 minutes.</p></list-item>
</list></p></list-item>
<list-item><p>Processing times:
<list list-type="simple">
<list-item><p>&#x02013; Unloading: 3 min (average) at the crusher.</p></list-item>
<list-item><p>&#x02013; Loading: 5 min (average) per excavator.</p></list-item>
</list></p></list-item>
<list-item><p>Initial fleet state:
<list list-type="simple">
<list-item><p>&#x02013; 36 trucks are staged at the unloading zone, simulating a high-demand scenario with synchronized fleet initialization.</p></list-item>
</list></p></list-item>
</list>
<p>Instead of the state-based prediction module (Section 3), DES generates the discrete battery-swapping schedule (<xref ref-type="table" rid="T2">Table 2</xref>), enumerating feasible time windows while accounting for real-world variability (e.g., traffic, equipment downtime). This DES-derived schedule serves as input to the hierarchical optimization framework (<xref ref-type="fig" rid="F9">Figure 9</xref>), which minimizes queuing delays through proactive demand dispersion.</p>
<p><xref ref-type="fig" rid="F12">Figure 12</xref> illustrates the scheduling outcomes generated by the Hierarchical Framework (Section 3) for Case Study 3. The framework converges to a high-quality solution, reducing the total waiting time to 30 minutes by proactively scheduling 100% of trucks for battery swaps prior to SoC depletion during the first operational cycle, demonstrating a 99% reduction compared to the MILP baseline (79.6 minutes) and an 99% improvement over the FCFS strategy (96.0 minutes). It is worth noting that DES yields an average waiting time of 81.3 minutes and 599 total trips completed. It can be inferred that the solution results based on the above hierarchical framework will significantly increase transport trip throughput by the DES.</p>
<fig position="float" id="F12">
<label>Figure 12</label>
<caption><p>Case 3: the scheduling results generated by the hierarchical framework, the red markers indicate the battery-swapping moments, the orange markers denote the arrival times of trucks at the station, and the blue markers represent the battery depletion moments (when the battery&#x00027;s State of Charge reaches critical levels).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0012.tif">
<alt-text content-type="machine-generated">Chart showing a grid with truck IDs on the vertical axis and time slice indices on the horizontal axis. Red squares appear sporadically over a green background, indicating specific data points.</alt-text>
</graphic>
</fig>
<p>The scheduling solution derived by the Hierarchical Framework (<xref ref-type="fig" rid="F12">Figure 12</xref>, Section 3) is coupled into the Discrete Event Simulation (DES) environment to validate its real-world applicability. Prior to integration, the solution undergoes an operational translation: The framework&#x00027;s time-based battery-swapping assignments (e.g., &#x0201C;swap at 10:15 AM&#x0201D;) are converted into trip-triggered thresholds (e.g., &#x0201C;swap after completing 3 trips&#x0201D;), aligning the proactive scheduling logic with DES&#x00027;s discrete, event-driven architecture. This transformation ensures seamless integration of the hierarchical schedule into DES, enabling rigorous testing of its robustness under dynamic operational variability. The scheduling solution derived by the Hierarchical Framework achieves 630 total transport trips&#x02014;a 5.2% throughput increase over the DES baseline (599 trips)&#x02014;while reducing the average battery-swapping wait time to 4.5 min, a 94.5% reduction from the DES benchmark of 81.3 minutes. <xref ref-type="fig" rid="F13">Figure 13</xref> illustrates the SoC time series with and without the framework&#x00027;s solution, demonstrating near-elimination of prolonged queuing delays during the initial battery-swapping cycle through proactive demand dispersion (<xref ref-type="fig" rid="F12">Figure 12</xref>). The nonlinear scalability of total transport trips&#x02014;despite the 94.5% reduction in battery-swapping queuing time&#x02014;highlights that the bottlenecks inherent to open-pit mining operations have changed from battery-swapping to other subsystems (e.g., Truck fleet throughput, Excavator loading rates and Crusher processing capacity). To validate the hypothesis that mechanical subsystems&#x02014;not battery-swapping efficiency&#x02014;govern ultimate productivity, we tested the hierarchical framework&#x00027;s scalability by incrementally expanding the number of loading platforms and unloading zones. <xref ref-type="fig" rid="F14">Figure 14</xref> (left) shows that increasing loading platforms from 2 to 3 elevates total transport trips by 10.6% (from 688 to 761 trips), significantly exceeding the 5.2% gain observed in the dual-platform configuration. This nonlinear scalability confirms that excavator loading capacity (fixed at 5 min per truck) becomes the dominant bottleneck as platforms scale, overriding energy management optimizations. Conversely, <xref ref-type="fig" rid="F14">Figure 14</xref> (right) reveals marginal diminishing returns when expanding unloading zones: increasing zones from 1 to 2 raises trip gains by only 1.1% (10.6% &#x02192; 11.7%, 11.7% is obtained by the increment from 695 to 776 trips), with no further improvement when adding a third zone. This plateau&#x02014;attributable to the crusher&#x00027;s fixed processing capacity (3 min per unload)&#x02014;demonstrates that unloading zone availability ceases to constrain productivity beyond a critical threshold. While these results clarify the hierarchical framework&#x00027;s role in resolving energy-time bottlenecks (e.g., battery queuing), holistic optimization of mechanical subsystems (e.g., excavator cycle times, crusher throughput) remains an open challenge beyond this study&#x00027;s scope.</p>
<fig position="float" id="F13">
<label>Figure 13</label>
<caption><p>SoC time series for all trucks with (<bold>left</bold>) and without (<bold>right</bold>) the scheduling solution derived by the Hierarchical Framework coupled into DES.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0013.tif">
<alt-text content-type="machine-generated">Two line graphs compare the state of charge (SOC) over simulation time in minutes. Each graph shows multiple lines in different colors representing various data sets. The left graph depicts more overlapping patterns, while the right graph has distinct downward trends followed by mixed fluctuations. Both graphs range from zero to one for SOC and zero to one thousand minutes for time.</alt-text>
</graphic>
</fig>
<fig position="float" id="F14">
<label>Figure 14</label>
<caption><p>Total trips increase with the increase of the number of loading platforms.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fcomp-07-1729185-g0014.tif">
<alt-text content-type="machine-generated">Two bar graphs compare total trips with and without optimal scheduling. The left graph shows trips based on the number of loading platforms (2 to 4), indicating higher trips with optimal scheduling. The right graph examines unloading zones (1 to 3), also showing increased trips with optimal scheduling. Both graphs illustrate a clear advantage with optimization.</alt-text>
</graphic>
</fig>
</sec>
</sec>
<sec id="s5">
<label>5</label>
<title>Summary</title>
<p>This paper presents a hierarchical optimization model tailored for off-peak battery-swapping scheduling of electric trucks in open-pit mines, aiming to mitigate queuing inefficiencies caused by synchronized battery depletion. The proposed framework integrates an inner-layer Mixed-Integer Linear Programming (MILP) model and an outer-layer Bayesian Optimization (BO) layer. The MILP addresses station-level resource allocation by optimizing battery-swapping assignments and sequencing under capacity constraints, while the BO navigates the combinatorial complexity of scheduling decisions to proactively stagger battery-swapping demand across temporal and spatial dimensions. Through Discrete Event Simulation (DES) and case studies, the authors demonstrate significant reductions in queuing time (up to 94.5%) and improvements in haulage throughput (5.2-10.6%) compared to conventional strategies like First-Come-First-Served (FCFS). The framework effectively resolves synchronization-induced bottlenecks, particularly during initial battery-swapping cycles, by leveraging predictive state-of-charge monitoring and dynamic dispatching. Sensitivity analyses further reveal that optimized scheduling, rather than infrastructure scaling, offers a cost-effective solution for enhancing operational efficiency in resource-constrained environments.</p>
<p>Future research could extend this work by integrating real-time adaptive mechanisms to address dynamic uncertainties such as fluctuating energy demands, traffic disruptions, and equipment downtime. Enhancing the model&#x00027;s scalability for larger fleets or multi-mine networks, coupled with hybrid energy systems incorporating renewables, would broaden its applicability. Additionally, exploring synergies between battery-swapping schedules and mechanical subsystems (e.g., excavator loading rates, crusher throughput) could unlock holistic productivity gains. The integration of digital twins or edge computing for real-time decision-making, alongside lifecycle cost analysis of battery degradation and infrastructure investments, would further refine the framework&#x00027;s economic and environmental viability. Finally, extending the hierarchical model to multi-objective optimization&#x02014;balancing energy costs, carbon emissions, and workforce safety&#x02013;could align operational efficiency with broader sustainability goals in the mining sector.</p></sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>CM: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Validation, Visualization, Writing &#x02013; original draft, Writing &#x02013; review &#x00026; editing. YT: Funding acquisition, Supervision, Writing &#x02013; review &#x00026; editing. SX: Investigation, Writing &#x02013; review &#x00026; editing. XZh: Investigation, Writing &#x02013; review &#x00026; editing. XZe: Investigation, Writing &#x02013; review &#x00026; editing. LW: Investigation, Writing &#x02013; review &#x00026; editing.</p>
</sec>
<sec sec-type="COI-statement" id="conf1">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="ai-statement" id="s9">
<title>Generative AI statement</title>
<p>The author(s) declare that no Gen AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p></sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x00027;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Bai</surname> <given-names>R.</given-names></name> <name><surname>Fu</surname> <given-names>E.</given-names></name> <name><surname>Ma</surname> <given-names>L.</given-names></name> <name><surname>Zhao</surname> <given-names>H.</given-names></name> <name><surname>Chai</surname> <given-names>S.</given-names></name></person-group> (<year>2024</year>). <article-title>Collaborative mining technological system of safety-green-high efficiency-low carbon for open pit coal mine</article-title>. <source>J. China Coal Soc</source>. <volume>49</volume>, <fpage>298</fpage>&#x02013;<lpage>308</lpage>. doi: <pub-id pub-id-type="doi">10.13225/j.cnki.jccs.2023.1433</pub-id></mixed-citation>
</ref>
<ref id="B2">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Ban</surname> <given-names>M.</given-names></name> <name><surname>Zhang</surname> <given-names>Z.</given-names></name> <name><surname>Li</surname> <given-names>C.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Liu</surname> <given-names>Y.</given-names></name></person-group> (<year>2021</year>). <article-title>Optimal scheduling for electric vehicle battery swapping-charging system based on nanogrids</article-title>. <source>Int. J. Elect. Power Energy Syst</source>. <volume>130</volume>:<fpage>106967</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijepes.2021.106967</pub-id></mixed-citation>
</ref>
<ref id="B3">
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Banks</surname> <given-names>J.</given-names></name> <name><surname>Carson</surname> <given-names>J. S.</given-names></name> <name><surname>Nelson</surname> <given-names>B. L.</given-names></name> <name><surname>Nicol</surname> <given-names>D.M.</given-names></name></person-group> (<year>2013</year>). <source>Discrete-Event System Simulation, 5th Edn.</source> <publisher-loc>Essex</publisher-loc>: <publisher-name>Pearson Education Limited</publisher-name>.</mixed-citation>
</ref>
<ref id="B4">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Deng</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Z.</given-names></name> <name><surname>Yan</surname> <given-names>P.</given-names></name> <name><surname>Zhong</surname> <given-names>R.</given-names></name></person-group> (<year>2023</year>). <article-title>Battery swapping and management system design for electric trucks considering battery degradation</article-title>. <source>Transport. Res. Part D: Transp. Environm</source>. <volume>122</volume>:<fpage>103860</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.trd.2023.103860</pub-id></mixed-citation>
</ref>
<ref id="B5">
<mixed-citation publication-type="web"><person-group person-group-type="author"><collab>GlobalData</collab></person-group> (<year>2025</year>). <source>Global Electric Mining Truck Market Report</source>. Available online at: <ext-link ext-link-type="uri" xlink:href="https://www.globaldata.com/reportstore/global-electric-mining-truck-market/">https://www.globaldata.com/reportstore/global-electric-mining-truck-market/</ext-link> (Accessed May 15, 2025).</mixed-citation>
</ref>
<ref id="B6">
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>He</surname> <given-names>M.</given-names></name> <name><surname>Zhou</surname> <given-names>J.</given-names></name> <name><surname>Nie</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>&#x0201C;Ga-based resource transportation scheduling optimization of open-pit mine,&#x0201D;</article-title> in <source>2017 International Conference on Smart Grid and Electrical Automation (ICSGEA)</source> (<publisher-loc>Changsha</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>247</fpage>&#x02013;<lpage>250</lpage>.</mixed-citation>
</ref>
<ref id="B7">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huang</surname> <given-names>A.</given-names></name> <name><surname>Zhang</surname> <given-names>Y.</given-names></name> <name><surname>He</surname> <given-names>Z.</given-names></name> <name><surname>Hua</surname> <given-names>G.</given-names></name> <name><surname>Shi</surname> <given-names>X.</given-names></name></person-group> (<year>2021</year>). <article-title>Recharging and transportation scheduling for electric vehicle battery under the swapping mode</article-title>. <source>Adv. Prod. Eng. Managem</source>. <volume>16</volume>, <fpage>265</fpage>&#x02013;<lpage>388</lpage>. doi: <pub-id pub-id-type="doi">10.14743/apem2021.3.406</pub-id></mixed-citation>
</ref>
<ref id="B8">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Huayanca</surname> <given-names>D.</given-names></name> <name><surname>Bujaco</surname> <given-names>G.</given-names></name> <name><surname>Delgado</surname> <given-names>A.</given-names></name></person-group> (<year>2023</year>). <article-title>Application of discrete-event simulation for truck fleet estimation at an open-pit copper mine in Peru</article-title>. <source>Appl. Sci</source>. <volume>13</volume>:<fpage>4093</fpage>. doi: <pub-id pub-id-type="doi">10.3390/app13074093</pub-id></mixed-citation>
</ref>
<ref id="B9">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Jordehi</surname> <given-names>A. R.</given-names></name> <name><surname>Javadi</surname> <given-names>M. S.</given-names></name> <name><surname>Catalo</surname> <given-names>J. P. S.</given-names></name></person-group> (<year>2021</year>). <article-title>Optimal placement of battery swap stations in microgrids with micro pumped hydro storage systems, photovoltaic, wind and geothermal distributed generators</article-title>. <source>Int. J. Elect. Power Energy Syst</source>. <volume>125</volume>:<fpage>106483</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijepes.2020.106483</pub-id></mixed-citation>
</ref>
<ref id="B10">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Lindgren</surname> <given-names>L.</given-names></name> <name><surname>Grauers</surname> <given-names>A.</given-names></name> <name><surname>Ranggd</surname> <given-names>J.</given-names></name> <name><surname>Mki</surname> <given-names>R.</given-names></name></person-group> (<year>2022</year>). <article-title>Drive-cycle simulations of battery-electric large haul trucks for open-pit mining with electric roads</article-title>. <source>Energies</source>, <volume>15</volume>:<fpage>4871</fpage>. doi: <pub-id pub-id-type="doi">10.3390/en15134871</pub-id></mixed-citation>
</ref>
<ref id="B11">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Sun</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Borozan</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>G.</given-names></name> <name><surname>Qiu</surname> <given-names>J.</given-names></name> <name><surname>Strbac</surname> <given-names>G.</given-names></name></person-group> (<year>2024</year>). <article-title>Battery swapping dispatch for self-sustained highway energy system based on spatiotemporal deep-learning traffic flow prediction</article-title>. <source>IEEE Trans. Indust. Appl</source>. <volume>60</volume>, <fpage>1058</fpage>&#x02013;<lpage>1070</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TIA.2023.3321713</pub-id></mixed-citation>
</ref>
<ref id="B12">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Teng</surname> <given-names>S.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>L.</given-names></name> <name><surname>Ai</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>L.</given-names></name></person-group> (<year>2024</year>). <article-title>Scenario engineering for autonomous transportation: A new stage in open-pit mines</article-title>. <source>IEEE Trans. Intellig. Vehicl</source>. <volume>9</volume>, <fpage>4394</fpage>&#x02013;<lpage>4404</lpage>. doi: <pub-id pub-id-type="doi">10.1109/TIV.2024.3373495</pub-id></mixed-citation>
</ref>
<ref id="B13">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Vallera</surname> <given-names>A. M.</given-names></name> <name><surname>Nunes</surname> <given-names>P. M.</given-names></name> <name><surname>Brito</surname> <given-names>M. C.</given-names></name></person-group> (<year>2021</year>). <article-title>Why we need battery swapping technology</article-title>. <source>Energy Policy</source> <volume>157</volume>:<fpage>112481</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.enpol.2021.112481</pub-id></mixed-citation>
</ref>
<ref id="B14">
<mixed-citation publication-type="book"><person-group person-group-type="author"><collab>Vigerske S. and Gleixner, A..</collab></person-group> (<year>2016</year>). <article-title>&#x0201C;SCIP: global optimization of mixed-integer nonlinear programs in a branch-and-cut framework,&#x0201D;</article-title> in <source>Optimization Methods &#x00026; Software</source> (<publisher-loc>Abingdon</publisher-loc>: <publisher-name>Taylor &#x00026; Francis Group</publisher-name>), <fpage>1</fpage>&#x02013;<lpage>31</lpage>.</mixed-citation>
</ref>
<ref id="B15">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>F. Y.</given-names></name> <name><surname>Chen</surname> <given-names>Z.</given-names></name> <name><surname>Hu</surname> <given-names>Z.</given-names></name></person-group> (<year>2024</year>). <article-title>Comprehensive optimization of electrical heavy-duty truck battery swap stations with a soc-dependent charge scheduling method</article-title>. <source>Energy</source>, <volume>308</volume>:<fpage>132773</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.energy.2024.132773</pub-id></mixed-citation>
</ref>
<ref id="B16">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Wang</surname> <given-names>S.</given-names></name> <name><surname>Chen</surname> <given-names>A.</given-names></name> <name><surname>Wang</surname> <given-names>P.</given-names></name> <name><surname>Zhuge</surname> <given-names>C.</given-names></name></person-group> (<year>2023</year>). <article-title>Short-term electric vehicle battery swapping demand prediction: deep learning methods</article-title>. <source>Transport. Res. Part D: Transp. Environm</source>. <volume>119</volume>:<fpage>103746</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.trd.2023.103746</pub-id></mixed-citation>
</ref>
<ref id="B17">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>Y.</given-names></name> <name><surname>Zhou</surname> <given-names>W.</given-names></name> <name><surname>Luan</surname> <given-names>B.</given-names></name> <name><surname>Yang</surname> <given-names>K.</given-names></name> <name><surname>Yang</surname> <given-names>Y.</given-names></name></person-group> (<year>2024</year>). <article-title>Truck transportation scheduling for a new transport mode of battery-swapping trucks in open-pit mines</article-title>. <source>Appl. Sci</source>. <volume>14</volume>:<fpage>10185</fpage>. doi: <pub-id pub-id-type="doi">10.3390/app142210185</pub-id></mixed-citation>
</ref>
<ref id="B18">
<mixed-citation publication-type="book"><person-group person-group-type="author"><name><surname>Xu</surname> <given-names>T.</given-names></name> <name><surname>Shi</surname> <given-names>F.</given-names></name> <name><surname>Liu</surname> <given-names>W.</given-names></name></person-group> (<year>2019</year>). <article-title>&#x0201C;Research on open-pit mine vehicle scheduling problem with approximate dynamic programming,&#x0201D;</article-title> in <source>2019 IEEE International Conference on Industrial Cyber Physical Systems (ICPS)</source> (<publisher-loc>Taipei</publisher-loc>: <publisher-name>IEEE</publisher-name>), <fpage>571</fpage>&#x02013;<lpage>577</lpage>.</mixed-citation>
</ref>
<ref id="B19">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>X.</given-names></name> <name><surname>Yang</surname> <given-names>B.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name> <name><surname>Ma</surname> <given-names>K.</given-names></name> <name><surname>Xu</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2024</year>). <article-title>Online optimal scheduling for battery swapping charging systems with partial delivery</article-title>. <source>Electric Power Syst. Res</source>. <volume>235</volume>:<fpage>110629</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.epsr.2024.110629</pub-id></mixed-citation>
</ref>
<ref id="B20">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhan</surname> <given-names>W.</given-names></name> <name><surname>Wang</surname> <given-names>Z.</given-names></name> <name><surname>Zhang</surname> <given-names>L.</given-names></name> <name><surname>Liu</surname> <given-names>P.</given-names></name> <name><surname>Cui</surname> <given-names>D.</given-names></name> <name><surname>Dorrell</surname> <given-names>D. G.</given-names></name></person-group> (<year>2022</year>). <article-title>A review of siting, sizing, optimal scheduling, and cost-benefit analysis for battery swapping stations</article-title>. <source>Energy</source> <volume>258</volume>:<fpage>124723</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.energy.2022.124723</pub-id></mixed-citation>
</ref>
<ref id="B21">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>C.</given-names></name> <name><surname>Lu</surname> <given-names>X.</given-names></name> <name><surname>Chen</surname> <given-names>S.</given-names></name> <name><surname>Shi</surname> <given-names>M.</given-names></name> <name><surname>Sun</surname> <given-names>Y.</given-names></name> <name><surname>Wang</surname> <given-names>S.</given-names></name> <etal/></person-group>. (<year>2024</year>). <article-title>Synergies of variable renewable energy and electric vehicle battery swapping stations: case study for Beijing</article-title>. <source>eTransportation</source> <volume>22</volume>:<fpage>100363</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.etran.2024.100363</pub-id></mixed-citation>
</ref>
<ref id="B22">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><collab>Zhang X. and Wang, G..</collab></person-group> (<year>2016</year>). <article-title>&#x0201C;Optimal dispatch of electric vehicle batteries between battery swapping stations and charging stations,&#x0201D;</article-title> in <source>2016 IEEE Power</source> &#x00026; <italic>Energy Society General Meeting (PESGM)</italic> (Boston: IEEE).</mixed-citation>
</ref>
<ref id="B23">
<mixed-citation publication-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>F.</given-names></name> <name><surname>Li</surname> <given-names>L.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>K.</given-names></name> <name><surname>Lu</surname> <given-names>L.</given-names></name> <name><surname>Han</surname> <given-names>X.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Does the battery swapping energy supply mode have better economic potential for electric heavy-duty trucks</article-title>. <source>eTransportation</source> <volume>15</volume>:<fpage>100215</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.etran.2022.100215</pub-id></mixed-citation>
</ref>
</ref-list>
<fn-group>
<fn fn-type="custom" custom-type="edited-by" id="fn0001">
<p>Edited by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1052970/overview">Ningyi DAI</ext-link>, University of Macau, China</p></fn>
<fn fn-type="custom" custom-type="reviewed-by" id="fn0002">
<p>Reviewed by: <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3254479/overview">Jiang Dongchu</ext-link>, Hunan City University, China</p>
<p><ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3256606/overview">Yufeng Xiao</ext-link>, State Key Laboratory of Coal Resources and Safe Mining, China University of Mining and Technology, China</p></fn>
</fn-group>
</back>
</article>