1. National Center of Technology Innovation for Digital Construction, Huazhong University of Science and Technology, Wuhan 430074, China; School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China
2. National Center of Technology Innovation for Digital Construction, Huazhong University of Science and Technology, Wuhan 430074, China; School of Civil and Hydraulic Engineering, Huazhong University of Science and Technology, Wuhan 430074, China; Institute of Artificial Intelligence, Huazhong University of Science and Technology, Wuhan 430074, China
3. Shantui Construction Machinery Co., Ltd, Jining 272000, China
charleschou@163.com
Show less
History+
Received
Accepted
Published
2024-05-07
2024-09-23
2025-03-15
Issue Date
Revised Date
2024-10-30
2024-09-19
PDF
(6341KB)
Abstract
Large-scale machinery operated in a coordinated manner in earthworks for mining constitutes high safety risks. Efficient scheduling of such machinery, factoring in safety constraints, could save time and significantly improve the overall safety. This paper develops a model of automated equipment scheduling in mining earthworks and presents a scheduling algorithm based on deep reinforcement learning with spatio-temporal safety constraints. The algorithm not only performed well on safety parameters, but also outperformed randomized instances of various sizes set against real mining applications. Further, the study reveals that responsiveness to spatio-temporal safety constraints noticeably increases as the scheduling size increases. This method provides important noticeable improvements to safe automated scheduling in mining.
Mining projects are often characterized by complex and hazardous terrain conditions, involving a significant number of construction machines, frequently of considerable size. These operations may sometimes neglect established safety regulations. For instance, between 2008 and 2021, there were a total of 528 mining-related fatalities reported in the United States, with 83 deaths specifically attributed to safety rule violations. This makes such violations the leading contributing factor to accidents (Quintero 2023). Mine earthwork represents a critical aspect of mining engineering and requires the coordination of numerous bulldozers and trucks. However, challenges such as blind spots (Gauthier et al., 2022; Imam et al., 2023) and driver fatigue (Inayah et al., 2023) make it unreliable to depend on operators to consistently adhere to safety regulations over extended durations. Consequently, scheduling mining earthmoving machinery with a focus on safety regulations presents a viable solution.
In contrast, extensive research has been devoted to automated equipment scheduling in earthworks for construction projects. Early methods were based on databases and progress mapping (Hong et al., 2023; Lee et al., 2023; Lee et al., 2022; Li and Lu 2019; Stipanovic et al., 2021). With advancements in artificial intelligence, genetic algorithms, and deep reinforcement learning (DRL), these techniques are gradually being integrated into automated scheduling approaches (Choi and Han 2023; Deng et al., 2020; Hwang et al., 2014; Liu et al., 2022; Shehadeh et al., 2022; Woo et al., 2018). Furthermore, the focus of scheduling has shifted from single-objective to multi-objective optimization (Hwang et al., 2014; Liu et al., 2022; Shehadeh et al., 2022). Unfortunately, a critical aspect, safety, is often overlooked. This concern is especially pressing in mining engineering, where large-scale machinery introduces significant safety risks. Therefore, it is crucial to prioritize scheduling methods that effectively incorporate safety considerations.
Operations management incorporating safety considerations has been researched. A common method for safety risk prediction is to build this prediction using early information modeling (Feng et al., 2023; Liang and Liu, 2022; Zhou et al., 2024) and more advanced artificial intelligence prediction (Ayhan and Tokdemir 2020; Gondia et al., 2023; Hallowell et al., 2013; Kulikova and Balovtsev 2020; Liu et al., 2020; Lu et al., 2023; Ou et al., 2021; Sadeghi et al., 2020). Real time safety control, for instance, utilization of safety rules or real time sensing (Deepak and Mahesh 2023; Lee D et al., 2023; Saurin et al., 2008), has been of recent interest. Nevertheless, though these are management approaches at the macro level including a system of safety indicators routine assessment of safety risks and work control perception system, the perception of safety work of subordinates reduces to the passive fulfillment of the requirements of safety procedures. Furthermore, safety constraints are an appropriate method for use in project scheduling. Shop operation scheduling has been a focus of current research rather than safety scheduling in earthwork projects.
Artificial intelligence technology has been developed widely and a great number of intelligent methods and tools have been used in the engineering field for industrial automation (Lu et al., 2024; Qin et al., 2024; You et al., 2023; Zhou et al., 2023). Consequently, the focus of this research is on mine earthworks and a method for the modeling of the automated scheduling of machinery is proposed. In addition, it presents operational constraint for mine earthworks incorporating safety. This research then proceeds to develop an automated machinery scheduling method that incorporates spatio-temporal safety constraints. This paper presents an optimization based on the DRL method which optimizes the machinery scheduling scheme in conjunction with safety constraints. This research proposes a method with high value in reducing safety accidents in mine earthworks.
2 Related work
2.1 Planning and management for mine construction
Mining engineering is a highly complex field that has attracted extensive scholarly attention, particularly in the areas of planning and management. The primary focus of mining engineering lies in mining operations, making the planning and management of these operations the most researched topic within the discipline. The principal aim is to maximize profitability, which is achieved through the optimization of mining planning, taking into account capital flow and mining costs (Guo et al., 2021; Paithankar & Chatterjee, 2019). However, profitability represents just one dimension of mining operations. Other crucial aspects of mine planning development include the socio-economic context (Liu et al., 2023) and the production lifespan of the mine (Xu et al., 2023). Furthermore, various in-mine factors, such as ventilation and cooling requirements, play a significant role in creating a safe and productive working environment for operators (Ogunmodede et al., 2022). Additionally, warehouse capacity considerations are vital (Flores-Fonseca et al., 2022). Beyond macro-level planning, substantial research exists on micro-level operations management, which centers on short-term resource allocation constraints (Joshi et al., 2022; Manríquez et al., 2023). Nonetheless, effective short-term planning often necessitates optimization across multiple objectives to achieve satisfactory mining operation planning solutions (Both et al., 2020; Paithankar et al., 2020).
In addition to mining operations, the transshipment of materials within the mine requires the coordinated efforts of various machines, which has become a prominent area of study. Trucks serve as the primary transportation machinery in mining, making the optimization of truck management a critical research focus. Successful truck management involves considerations of fleet size, fleet allocation, and other relevant factors (Mohtasham et al., 2022). Given the typically large size of mining trucks, reducing idle time (Mirzaei-Nasirabad et al., 2023) and minimizing energy consumption (Zhang et al., 2022) are essential objectives. Conveyor transportation is also a widely used method in mines, and its effective planning and management can significantly lower electricity costs (Mu et al., 2020). Moreover, the advent of intelligent technology has paved the way for the use of robotic transloading systems in mining operations (Yan et al., 2023). Intelligent methods are increasingly being employed for scheduling and managing transloading operations (Huo et al., 2023; Wang et al., 2023). Additionally, some researchers have explored the scheduling and management of resources such as water (Bo et al., 2022) and electricity (Nikolaev et al., 2023).
While significant research has focused on planning and managing mining and transshipment operations, there has been relatively less exploration of mine earthworks. However, mine earthwork operations involve a substantial number of bulldozers and trucks working in concert, necessitating careful planning and management in this domain.
2.2 Automated equipment scheduling for earthwork
Earthwork is vital to construction projects, characterized by an extensive construction cycle and the need for coordination among a large array of engineering machinery. In recent years, the advancement of automated construction has elevated the automation of earthwork equipment scheduling into a prominent area of research. Effective scheduling prioritizes resources, prompting some researchers to devise scheduling methods based on resource databases (Lee C et al., 2023) or the status of assets (Stipanovic et al., 2021) involved in earthwork. However, these methods often require considerable time and financial investment to establish the necessary databases. As a result, some researchers have investigated the automation of scheduling by directly extracting information from progress diagrams (Hong et al., 2023; Lee et al., 2022; Li and Lu, 2019). Furthermore, simulation has emerged as an effective approach for enhancing scheduling decisions by developing simulation models through the use of sensors or drones (Rachmawati et al., 2023; Zamani et al., 2023).
With the rise of artificial intelligence, intelligent algorithms, such as genetic algorithms, are becoming increasingly prominent in the automatic scheduling of earthwork equipment. For example, earthwork can be conceptualized as a path assignment problem, solvable through genetic algorithms (Deng et al., 2020). A particularly notable advantage of genetic algorithms is their capability for multi-objective optimization, which is essential for the successful execution of earthwork projects (Hwang et al., 2014; Liu et al., 2022; Shehadeh et al., 2022). Consequently, many researchers advocate for the application of genetic algorithms in the scheduling of earthworks from a multi-objective optimization standpoint. Additionally, DRL—combining deep learning capabilities with reinforcement learning—has proven to be an effective strategy for addressing scheduling challenges. Some researchers have begun implementing DRL in earthwork scheduling, yielding significant results (Choi and Han, 2023; Woo et al., 2018).
While numerous methods for scheduling automated equipment for earthworks exist, many primarily focus on cost, time, and similar considerations, often overlooking safety. In practice, safety is significant due to the diverse range of construction machinery involved in earthwork operations. This concern is heightened in mining contexts, where the number of pieces of equipment is substantially greater, underscoring the necessity of incorporating safety as a critical factor in scheduling.
2.3 Project management and scheduling considering safety
In construction projects, inherent safety risks must be continuously addressed. Consequently, an increasing number of scholars are investigating strategies to integrate safety into operations management. One practical and efficient method is to predict safety risks prior to the commencement of construction. Building information modeling (BIM) is employed not only during the design phase but also plays a crucial role in forecasting safety risks. The synergy between BIM and the Internet of Things (IoTs) presents significant potential for effectively anticipating risks throughout the construction process (Feng et al., 2023; Liang and Liu, 2022). Furthermore, researchers have begun to consolidate essential safety risk assessment principles and indicators through an analysis of existing literature standards and accident case studies (Hallowell et al., 2013; Kulikova and Balovtsev, 2020; Liu et al., 2020). With advancements in artificial intelligence, machine learning, and neural network-based methods have become mainstream for safety risk prediction, yielding promising results (Ayhan and Tokdemir, 2020; Gondia et al., 2023; Lu et al., 2023; Ou et al., 2021; Sadeghi et al., 2020).
However, implementing a safety risk prediction approach in operations management is not a one-time solution. In practical engineering scenarios, it is often crucial to actively manage safety risks in real-time during construction. For instance, a systems engineering approach can be employed to analyze the construction process and identify critical safety management principles (Saurin et al., 2008). Establishing a comprehensive safety framework based on safety indicators is another effective strategy (Deepak and Mahesh, 2023). Additionally, the integration of safety regulations with sensor technologies, such as cameras, can significantly enhance real-time monitoring and safety control, thereby facilitating intelligent safety management (Lee D et al., 2023).
For project scheduling challenges, incorporating constraints to ensure safety is an effective strategy. In the context of shop scheduling, it is crucial to consider ergonomic risks while optimizing energy consumption, thereby safeguarding the well-being of shop workers through adherence to safety goals (Destouet et al., 2024). Worker fatigue is also a significant safety concern (Aribi et al., 2023). Beyond physical collisions, environmental factors such as dust and noise pose additional threats to worker safety, and these issues have been the subject of various research studies (Fu et al., 2020). Furthermore, the reliability of safety can be quantitatively represented as safety time (Yang and Xu, 2021) and safety distance (Zhou and Liao, 2020). However, current scheduling methodologies that account for safety primarily focus on the scheduling of shop jobs within factories, with limited research dedicated to the safe scheduling of earthwork operations at construction sites. Earthwork tasks entail the coordinated operation of numerous large machinery units; therefore, incorporating safety considerations into scheduling constraints is both important and valuable.
3 Methodology
As illustrated in Fig.1, this research proposes an automated equipment scheduling method tailored specifically for mine earthworks. Initially, the study models the scheduling issues associated with mine earthworks and introduces operational constraints that emphasize safety considerations. Subsequently, spatio-temporal safety constraints are formulated and integrated into a DRL algorithm. This approach allows for a more scientific methodology in scheduling machinery for mine earthworks, taking potential risks into account to enhance overall construction safety.
3.1 Modeling of scheduling for automated mining earthwork equipment
3.1.1 Problem description
Mining earthwork is a vital component of mining engineering, particularly in coal mines and other mineral resource sites. These operations produce significant amounts of waste soil that must be transported to designated discharge locations for dumping, followed by subsequent earthwork activities performed by bulldozers. It is important to note that mining earthwork differs from conventional earthwork practices. The detailed process is illustrated in Fig.2 and consists of four key steps: a) trucks enter the site to transport large quantities of soil to the discharge point for unloading; b) trucks adjust their positions accordingly for this unloading; c) upon completion of the unloading, trucks exit the site while bulldozers enter the area to carry out bulldozing operations at the unloading point; d) after finishing their tasks, bulldozers leave the discharge yard in preparation for the next work cycle. Specifically, as shown in Fig.3, each machine operates along roughly parallel paths known as discharge lines. The end of these lines is where machinery pushes and fills soil—this location is referred to as the soil discharge point. In actual mine earthwork scenarios, continuous filling leads to an ongoing extension of the discharge line as it gradually advances toward new areas needing attention.
Given the collaboration of multiple large-scale machines in mine earthmoving operations, it is essential to implement a scientifically designed automated scheduling system. Such a system not only enhances standardization and efficiency within the workflow but also improves overall productivity. The machinery involved in mine earthmoving operations consists of various bulldozers and trucks. The primary work area is the soil discharge field, which contains multiple soil discharge points where machinery unloads and performs soil pushing operations in a specific sequence. The scheduling of all jobs and machines adheres to the following rules:
Rule 1: A machine can work only on one line at one time step.
Rule 2: No more than one machine can work on a single line at the same time
Rule 3: Once a machine starts a job, it cannot be interrupted.
Rule 4: Operations within the same job have a sequential order.
Thus, the scheduling of machinery for mine earth removal operations can be regarded as a problem analogous to the Flexible Job Shop Scheduling Problem (FJSP).
As illustrated in Fig.4, this paper proposes a model for scheduling earthmoving operations in mining within the framework of FJSP. In this model, the discharge field is segmented into parallel discharge lines, through which the bulldozers and trucks involved in the operations access the earth displacement area. Each line culminates at a discharge point, which serves as the principal site for the machinery’s tasks. At each discharge point, the machinery is required to complete a series of operations. Consequently, the FJSP model for mine earthmoving operations emphasizes the scientific scheduling of multiple bulldozers and trucks at each discharge point.
3.1.2 Operational constraints considering safety
The engineering machinery utilized in mine earthmoving operations is often large, leading to numerous blind spots for operators. This increases the risk of accidents, as well as potential casualties and property damage. Therefore, while optimizing the scheduling of machinery in earthwork operations, it is essential to ensure that adequate safety distances and durations are maintained between large machines. Fig.5 illustrates the proposed operational constraints derived from the modeling of mine earthmoving operations, with an emphasis on safety considerations.
The safety constraints proposed in this paper include two aspects. The first aspect pertains to the X-axis constraint, which governs the independent functioning of parallel row lines. Although machinery may operate concurrently in adjacent row lines, their close proximity significantly elevates the risk of accidents. Consequently, the X-axis safety constraint highlights the importance of maintaining a safe distance between machinery in parallel rows, thereby prohibiting simultaneous work in neighboring rows.
The second aspect involves the Y-axis constraint, which focuses on the sequence of operations within a row. Multiple pieces of equipment enter and exit a row of lines, and failing to maintain a safe time interval between them may result in safety incidents. As such, the Y-axis safety constraint is designed to ensure a safe exit time for successive operations performed by machines in the same row of soil lines. This requires maintaining a safe time interval between the start time of the latter machine’s work and the end time of the former machine’s work.
3.2 A deep reinforcement learning algorithm with spatio-temporal safety constraints
3.2.1 Markov decision modeling
For the solution to the scheduling problem in mine earthwork, we utilize a Markov decision process. The agent selects an action at each moment based on the current state and transitions to the next state while receiving an instantaneous reward. This instantaneous reward accumulates progressively as the scheduling process continues. The primary objective of scheduling is to identify a sequence of actions that maximizes the cumulative reward.
In the context of earthwork within a mine, the agent includes the overall scheduling solution, making choices about scheduling actions based on the current state St at each time step. As illustrated in Fig.6, the scheduling actions involve two key considerations. First, the agent must decide which job to perform, determining which operation at which discharge point will be executed. Second, the agent needs to select which machine will be employed to carry out the job. Depending on the selected action, the agent receives a reward value and updates the state to (Tab.1).
3.2.2 Model architecture
To address the Markov decision problem related to mine earthworks, this paper employs a DRL algorithm. DRL represents an advanced methodology that integrates deep learning with reinforcement learning, making it particularly effective for resolving complex scheduling challenges.
The algorithm employed in this research is grounded in the multi-action DRL framework (Lei et al., 2022), which serves as an end-to-end solution for efficiently and swiftly addressing scheduling problems. This framework adeptly integrates graph neural networks with multi-Proximal Policy Optimization (multi-PPO), yielding outstanding performance across various scheduling scenarios. The structure of the algorithm is illustrated in Fig.7. Each job contains specific structural information, such as priority constraints and the set of compatible machines for each operation. This structural information can be extracted and leveraged effectively. The Graph Isomorphism Network (GIN), a variant of Graph Neural Networks (GNNs), demonstrates even greater efficacy (Wang and Zhang, 2022). The GIN encodes the state nodes guiding the job selection strategy, which are then processed through a multilayer perceptron (MLP). There is no inherent structural information between the machine's state data; therefore, it can be directly encoded using a fully connected layer. Subsequently, both the job selection and machine selection outputs are linked in series with MLP-based decoders, utilizing an identical network structure while maintaining separate parameters. During the decoding stage, each decoder computes the probability distribution for the job operation space or the machine operation space, respectively, producing both the job operation score and the machine operation score.
3.2.3 Spatio-temporal safety constraints
To integrate safety considerations into job constraints, spatio-temporal safety constraints are implemented within the context of DRL. As illustrated in Fig.8, this research employs spatio-temporal safety constraints to refine certain parameters of the Markov decision process. In this process, the state matrix defined in this study incorporates StartTime and Dur as correction parameters, in addition to the three standard values of , , and . Specifically, at time step , the StartTime and Dur of state are transferred to two matrices, Safety_jz[] and Dur_jz[], for safety assessments. The new state ’, generated after action selection, does not immediately enter the next decision-making phase but must first undergo the safety assessment dictated by the spatio-temporal safety constraints. Only upon passing this assessment can it be deemed the new state . If it fails, the StartTime is adjusted to ensure compliance with the safety criteria.
Pseudo-code for spatio-temporal safety constraints is presented in Algorithm 1. Initially, this research establishes two matrices, Safety_jz[] and Dur_jz[] to log the StartTime and duration information throughout the scheduling process. The index value i in both matrices signifies the number of epochs, while the dimensions of the matrices are defined by the number of jobs and operations, represented by the row and column indices, respectively. It is crucial to note that these matrices are updated with a delay; focusing solely on historical states when making decisions aligns with the principles of a Markov decision process. These two matrices serve primarily for comparing safety decisions in subsequent analyses.
Specifically, the code includes five conditional statements, each addressing different safety constraints. The first if statement pertains to Y-axis constraints, stipulating that, within the same row of soil lines, the time interval between the start time of the subsequent machine and the end time of the preceding machine must exceed the operational duration of the latter. This requirement ensures that the former machine can safely exit after completing its task. Regarding X-axis constraints, this research mandates that machinery operating in adjacent rows of soil lines cannot function simultaneously, as indicated by the second through fifth if statements. In instances where these regulations are breached, the StartTime will be reverted to a point that complies with the specified safety constraints.
3.3 Optimization strategies for deep reinforcement learning
3.3.1 Instances generation
Publicly available data sets are insufficient in providing scheduling examples for mining earthwork operations. To affirm the scientific validity of the proposed method, this paper outlines a data pipeline for automatic instance generation rooted in engineering practices. The generation of instances adheres to the following three principles:
Rule 1: The total number of machines available for dispatch is contingent upon the size of the schedule; however, there is variability in the number of bulldozers and trucks assigned to each schedule.
Rule 2: Variability in machinery efficiency arises from factors such as the machinery’s performance and the operator’s proficiency.
Rule 3: The first operation in any row must be a truck unloading dirt, while the final operation must be a bulldozer pushing dirt.
Algorithm 2 illustrates the process in detail. Initially, the method employs two random arrays to generate the foundational data for the examples. The ‘time0’ array is a three-dimensional array with values predominantly ranging between 30 and 60, signifying that the majority of equipment operates within this time frame. Conversely, the ‘time1’ array is similar to ‘time0’, but consists of a single column with values ranging between 20 and 30. This indicates that a small number of devices are permitted to perform exceptionally, operating between 20 and 30. These two arrays are subsequently concatenated and reshaped to form an array named “times,” with dimensions n_j * n_m * n_m.
Building upon the time array and the established construction rules, various factors must be considered and evaluated. First, it is essential to determine the respective number of bulldozers and trucks using a random number generation method described in this paper. A random value, denoted as num_kk and falling between the maximum and minimum number of bulldozers, is selected to represent the number of bulldozers for this example, while the number of trucks is calculated as n_m - num_kk. Following this, specific construction rules must be established. It is evident that unloading of trucks should be the first operation in any job, while bulldozer operations should occur last. Therefore, in the generation of arithmetic examples, this research designates the bulldozer working time for the first operation and the truck working time for the final operation of each job as −999, indicating that these operations are unavailable. Moreover, to illustrate a defined construction rhythm, a random interpolation method is utilized to render certain operations unavailable. Consequently, an example representing mine earthmoving operations is generated.
3.3.2 Optimization strategies
Leveraging the instance generation method discussed in Section 3.3.1, the training and optimization of the DRL algorithm can proceed efficiently. Since the instance generation module operates independently, it can be invoked multiple times within each epoch to produce varying instances for model learning. Specifically, an epoch comprises several batches, with each batch generating and inputting a fixed number of instances for computation. The results are then aggregated by epoch to produce reward and loss values for model iteration and optimization.
In this study, the model aims to allocate an appropriate workload to the machinery while optimizing the solution. The primary optimization objective is to minimize the total scheduling time. The following formula represents the total makespan:
where is the value of makespan; where denotes the completion time of job i; and n× m denotes the size of the FJSP instance.
The scheduling problem addressed in this research focuses on strategies for both job selection and machine selection. To resolve this issue, a multi-PPO algorithm is employed, which utilizes a multi-subject-critic architecture to learn two sub-policies, using Proximal Policy Optimization (PPO) as the policy optimization method (Wu et al., 2022). The multi-PPO architecture consists of two component networks, with each network responsible for a singular policy. Within this architecture, a single critic network with parameters is utilized, serving as an estimator to learn the approximate value state function. The specific computation is depicted below:
where is the empirical tuple; is the estimator used to compute the variance-approximate minus dominance function; and γ is the discount factor.
4 Model evaluation and implementation
4.1 Training setup
4.1.1 Environment setting
This research utilized a computing platform featuring an 11th Gen Intel(R) Core(TM) i7-11800 CPU and an RTX 3090 (24GB) GPU. The algorithms were implemented using Python 3.8 and configured with essential libraries, including PyTorch 1.8.1 and CUDA 11.1.
4.1.2 Input conditions
To comprehensively evaluate the performance of the proposed methodology, three scales with varying instance sizes were designed. The instance sizes are determined by the number of available machines. The small-sized instances comprise five available machines, featuring a minimum of 2 bulldozers and a maximum of 3. The specific sizes include (J× M): 5 × 5, 5 × 10, and 5 × 20. The medium-sized instances consist of 10 available machines, with a minimum of 4 bulldozers and a maximum of 6. The designated sizes include (J× M): 10 × 10 and 10 × 20. The large-sized instances feature 20 available machines, with a minimum of 10 bulldozers and a maximum of 15, represented by the specific scale of (J× M): 20 × 20.
4.2 Scheduling results of generated instances
4.2.1 Analysis of scheduling results
This section analyzes the training and scheduling results at different scales, as presented in Tab.2 and Fig.8. The model demonstrates faster convergence during small-scale mechanical training, achieving convergence around the 1,000th batch. Due to the smaller scheduling scale, model training initiates with a lower loss. In medium-scale training, although convergence is still satisfactory, it is slightly slower compared to small-scale training. In large-scale training, the model initially experiences a higher loss due to the increased scheduling size. However, as training progresses, the loss curve for large-scale scheduling displays good convergence trends.
As shown in Tab.2, the incremental increase in scheduling size results in only a modest rise in total scheduling time. The implementation of the surface security rule has a minimal effect on the overall scheduling time, contributing only a negligible increase. This indicates that our proposed method can effectively adapt and converge in the presence of spatio-temporal safety constraints while preserving strong learning capabilities.
Moreover, as illustrated in Fig.9, although scaling up the scheduling leads to a proportional increase in the loss associated with the training start time, this does not significantly affect the model’s training process. In fact, the DRL model demonstrates rapid and consistent learning and convergence, even when applied at varying scheduling scales.
4.2.2 Analysis of spatio-temporal analytic map
To demonstrate the performance of the proposed method in this study, spatio-temporal analysis plots at various scales are presented in this section to visualize the specific effects of scheduling. It is important to note that the algorithms employed in this research, including those that generate the spatio-temporal analysis graphs, are automatically produced without any human intervention.
The spatio-temporal analysis plot is a three-dimensional representation, where the x and y axes correspond to jobs and their respective scheduling times, while the chosen machinery is depicted as distinct colored orbs. Notably, the z-axis in this context illustrates the trend of the total safety response value (SafetySum) in the scheduling results, representing the delays incurred due to the influence of spatio-temporal safety constraints at each time step.
In Fig.10, the spatio-temporal analysis diagram is presented at a scale of 5 × 10. Due to the smaller-scale scheduling, the spatio-temporal relationships are distinctly visible, with all five machines being relatively evenly engaged in the work. Furthermore, as time progresses, it is clear that the SafetySum value gradually increases, indicating that the spatio-temporal safety constraints consistently affect work progress, with later scheduling times experiencing a higher effect from these constraints.
Transitioning to Fig.11, which displays the spatio-temporal analysis graph at a scale of 10 × 10, the overall trend becomes more evident as the scheduling scale increases. The graph exhibits a consistent upward trend from the bottom left to the top right, reinforcing the notion that spatio-temporal safety constraints exert a more pronounced effect as time progresses.
Fig.12 presents the spatio-temporal analysis graph at a scale of 20 × 20. The graph reveals a nearly uniform plate-like pattern and shows a discernible upward trend from the bottom left to the top right. This indicates that the influence of spatio-temporal safety constraints gradually intensifies over time, following a linear and uniform trajectory. Consequently, it is evident that as scheduling progresses, the interactions between machinery become more complex, leading to an increased occurrence of safety rule violations. This highlights the importance of implementing spatio-temporal safety constraints.
4.2.3 Analysis of response to spatio-temporal safety constraints
To further evaluate the significance of spatio-temporal safety constraints, this study conducts a statistical analysis of their responses at various scales, based on the scheduling results. It is crucial to note that the response values in this context do not represent instantaneous reactions at individual time steps; rather, they reflect cumulative responses over time.
As illustrated in Tab.3 and Fig.13, in small-scale scheduling, both safetyXsum and safetyYsum show an upward trend, albeit with certain irregularities. Furthermore, as time advances, the slopes of safetyYsum exhibit an increasing trajectory, while the slope of safetyXsum gradually diminishes, potentially indicating a shift in the prioritization of X and Y responses. Fig.14 and Fig.15 present scheduling results at medium and large scales, respectively, and their fundamental trends correspond with those observed in small-scale scheduling. However, it becomes apparent that with the increasing scale of scheduling, the curves of both safetyXsum and safetyYsum tend to exhibit smoother changes. SafetyYsum demonstrates an increasing slope, whereas safetyXsum exhibits a decreasing slope.
To mitigate the influence of response priority, this study aggregates safetyXsum and safetyYsum to further analyze the combined response, referred to as safetySum. As illustrated in Fig.16, safetySum consistently increases with an upward trend as time progresses, aligning with the findings presented in Section 4.2.2. Notably, safetySum experiences a significant surge at higher multiples as the scheduling scale expands. This indicates that instances of safety rule violations become more frequent as scheduling scales increase, emphasizing the importance of employing scheduling methods that account for spatio-temporal safety constraints on a larger scale.
To gain a deeper understanding of the effect of spatio-temporal safety constraints on scheduling results, this study examined the changes in results during model training iterations. As shown in Fig.17, without these constraints, scheduling results exhibit minimal variation throughout the iterations. Additionally, the makespan of scheduling results without spatio-temporal safety constraints for epochs 20 and 100 reveals a difference of just under 100. In contrast, when spatio-temporal safety constraints are incorporated, the model demonstrates significant fluctuations during the iterations. For instance, at epoch 20, Machine 1 has yet to be utilized, but this situation improves markedly by epoch 60. This suggests that the presence of spatio-temporal safety constraints has a substantial effect on the scheduling of a significant portion of machinery. However, as the iterations progress, improved convergence is observed by epoch 100, indicating that the spatio-temporal safety constraints among the machines are being satisfactorily met. Ultimately, the makespan difference between epochs 20 and 100 for scheduling results with spatio-temporal safety constraints exceeds 300, highlighting the robust iterative capabilities of the model proposed in this research.
4.3 Scheduling results for multiple earthwork sites
In practical engineering contexts, there are scenarios in which multiple earthmoving yards are dispatched uniformly. This research provides further investigation into this operational scenario. The working conditions for multiple rows of earthmoving yards are largely similar to those for single rows, with the primary difference being that the X-axis safety constraints do not need to be considered when machinery is situated in different rows of earthmoving yards.
In this research, model training and testing were conducted with various numbers of earth discharge fields. It is important to note that the number of discharge yards functions as an intrinsic variable in the model training, which facilitates the direct application of scheduling across different numbers of discharge yards using a single training process. The training loss curves are illustrated in Fig.18. A scale of 10 × 20 was employed for training, and the model exhibited good convergence, showing minimal variations in the loss curves when compared to a single row of soil fields.
Following the training, scheduling tests were conducted using the obtained weights to evaluate instances with varying numbers of earth displacement fields, as presented in Tab.4. It is significant to emphasize that, in the instances illustrated in Tab.4, the machinery’s working time for each job remains constant; the only variation lies in the division of the earth displacement yard. It is clear that the total scheduling time shows a decreasing trend as the number of discharge yards increases. This trend can be attributed to the fact that a greater number of earthmoving fields allows for more rows to operate without needing to comply with X-axis safety constraints. Consequently, the time delays arising from these constraints are diminished, resulting in a reduction in total scheduling time. Additionally, Fig.19 displays the scheduling results for different numbers of earth rows, which are consistent with the findings in Tab.4. While the overall structure of the four scheduling result graphs remains similar, slight variations at certain node locations arise due to differing safety constraints.
5 Discussion
5.1 Performance comparison with different parameters
The model's performance can be influenced by its parameters. This section presents comparative results of the model's performance utilizing different parameter combinations. Two key parameters of interest are the learning rate and the size of the experience pool. An excessively large learning rate may hinder model convergence, while a rate that is too small can extend the convergence time. The second parameter, the experience pool size, must also be carefully considered; a pool that is too small may result in poor model performance, whereas an overly large pool can lead to an excess of experience, potentially causing early termination of model training.
As indicated in Tab.5, this study conducted performance comparisons across three different learning rates and three different empirical pool sizes. At a learning rate of 0.01, the higher learning rate resulted in model makespan values exceeding 900, suggesting that the model did not converge as effectively as it could have. Furthermore, the increased learning rate led to the premature convergence of the experience pools at both 3200 and 6400, indicating an underutilization of these pools. Conversely, at a learning rate of 0.0001, the 6400 experience pool achieved a makespan of 875.87, while the other experience pools performed inadequately, implying that a lower learning rate resulted in less efficient learning. Lastly, under a learning rate of 0.001, the 3200 experience pool demonstrated optimal performance and fulfilled training requirements effectively, indicating that it fully utilized the experience pool.
5.2 Comparison with other models
To further validate the superiority of the proposed method in this study, classical scheduling algorithms were selected for comparison. The First-In-First-Out (FIFO) scheduling algorithm organizes tasks according to their arrival order, queuing them on a first-come-first-served basis. The Longest Waiting Time with Known Release (LWKR) algorithm is typically employed for real-time tasks and aims to minimize task waiting time. The Most Operations Remaining Next (MOPNR) algorithm selects the job with the most remaining operations, intending to balance system load and reduce job waiting time.
As demonstrated in Tab.6 and Fig.20, the method proposed in this paper significantly outperforms the aforementioned classical scheduling algorithms regarding makespan. It is essential to note that a 10x10 arithmetic case at a medium scale was chosen to evaluate the performance of the different methods in this study. The makespan for the scheduling using DRL was 939, whereas the other scheduling algorithms exceeded 1000, with MOPNR reaching 1745. As illustrated in Fig.18, classical scheduling algorithms exhibit varying degrees of concentration in scheduling due to their lack of advantages in learning and iteration. Consequently, certain high-performing machines are assigned a disproportionate amount of work. Specifically, Fig.20 shows that LWKR and MOPNR demonstrate significant scheduling breaks attributable to waiting times. In contrast, DRL capitalizes on continuous learning and iteration to allocate work in a scientifically and rational manner within complex machinery, resulting in enhanced scheduling results.
5.3 Testing and evaluation in a real mine earthworks project
The method described in this research has been implemented in the earthworks of an active mine located in Hulun Buir, Inner Mongolia, China, where lignite is the primary mineral output. As illustrated in Fig.21, the mine spans a vast area with proven reserves of 4.166 billion tons, and it has the capacity to produce up to 35 million tons of coal annually. Currently, the mine features a large centralized earth removal site, employing two bulldozers and four trucks for its earth removal operations. This research conducts field tests using the proposed method within the operational context of this earth removal site.
In the field tests, this study examines the impact of spatio-temporal safety constraints on scheduling outcomes, as depicted in Fig.22. The primary goal of these constraints is to ensure that machinery maintains a safe distance and time during operations. Consequently, during the initial phases of scheduling—particularly at the start of each job—the influence of safety constraints is minimal. Most operations demonstrate consistency in machinery selection; however, variations in time distribution are observed.
As the time steps progress, the effects of temporal safety constraints become increasingly apparent. In this section, we quantify both the frequency and total working duration for each type of machinery involved in scheduling, as presented in Tab.7. The number of bulldozers was sufficiently low that safety restraints had negligible effects on their utilization. Conversely, trucks displayed significant changes in distribution patterns: without safety constraints, Machine 6 performed the majority of tasks; however, with the implementation of these constraints, previously underutilized machines, most notably Machine 3, were incorporated into the schedule, resulting in marked increases in both usage frequency and total operational duration.
It is evident that, in the absence of safety constraints, efforts to minimize makespan can lead to an uneven allocation that favors specific machines. In contrast, the implementation of spatio-temporal safety constraints necessitates broader participation from various types of machinery to uphold operational safety standards. Thus, the spatio-temporal safety constraints proposed in this research not only enhance overall safety but also contribute to a more equitable and scientifically grounded scheduling framework.
5.4 Limitations
Although this paper provides a scientific discussion on the safety of scheduling mine earthworks, it does have some limitations. The focus is primarily on the static scheduling problem, which considers safety during the initial conditions. The method proposed herein can mitigate the risk of accidents caused by inadequate safe distance and safe time in scheduling. However, it lacks sufficient consideration of environmental factors such as visibility and ground conditions, as well as dynamic factors like personnel interference. Future research could develop a more comprehensive dynamic scheduling technique building on the findings presented in this paper.
Additionally, this paper employs an idealized model of parallel rows of soil lines for mine earthworks. In actual construction scenarios, machinery may operate along intersecting paths. Thus, further exploration of the safe scheduling problem under complex path conditions is warranted based on the research conducted in this study.
6 Conclusions
Numerous collaborative machinery operations engaged in mine earthmoving may and obviously involve a heavy-weight portion of safety risk. This research addresses the scheduling of machinery for mine earthwork operations and provides an automated scheduling method based on DRL with spatio-temporal safety constraints. The proposed method yields promising results across various scales and has undergone real-world validation in an authentic mining environment. Thus:
1) This study proposes a model for the mine earthwork scheduling problem that prioritizes safety for the first time-a reference for safe construction practices in mining.
2) The presented spatio-temporal safety constraints are now an extension of and merged with DRL algorithms, allowing commendable performance.
Additionally, spatio-temporal safety constraints introduced in this paper enhance existing DRL algorithms by imposing restrictions on state updates during each iteration. Before long, improved structures and computational logic will be developed for such algorithms to offer fine scheduling solutions with better performance for given engineering constraints.
AribiDDrissO BEl HaouziH B (2023). Multi-objective optimization of the dynamic and flexible job shop scheduling problem under workers fatigue constraints. Proceedings of ICAART (3), 301–308
[2]
Ayhan B U,Tokdemir O B, (2020). Accident analysis for construction safety using latent class clustering and artificial neural networks. Journal of Construction Engineering and Management, 146( 3): 04019114–04019114
[3]
Bo L,Li Z H,Liu Y,Yue Y G,Zhang Z H,Wang Y Y, (2022). Research on multi-level scheduling of mine water reuse based on improved whale optimization algorithm. Sensors, 22( 14): 5164–5164
[4]
BothCDimitrakopoulosR J OEngineering (2020). Joint stochastic short-term production scheduling and fleet management optimization for mining complexes. 21(4): 1717–1743
[5]
Choi G,Han S G, (2023). Reinforcement learning-based dynamic planning of cut and fill operations for earthwork optimization. Automation in Construction, 156: 105111
[6]
Deepak M D,Mahesh G, (2023). A framework for enhancing construction safety through knowledge-based safety culture indicators. International Journal of Construction Management, 23( 12): 2039–2047
[7]
DengNLiXSu Y J E S W oC (2020). Optimization of earthwork allocation path as vehicle route problem based on genetic algorithm. E3S Web of Conferences, 165(1): 04057
[8]
Destouet C,Tlahig H,Bettayeb B,Mazari B J C, (2024). Multi-objective sustainable flexible job shop scheduling problem: Balancing economic, ecological, and social criteria. Computers & Industrial Engineering, 195: 110419
[9]
FengZShaYLiuM (2023). Application of BIM in the design and construction of fabricated buildings. Computer-Aided Design and Applications
[10]
Flores-Fonseca C,Linfati R,Escobar J W, (2022). Exact algorithms for production planning in mining considering the use of stockpiles and sequencing of power shovels in open-pit mines. Operations Research, 22( 3): 2529–2553
[11]
Fu Y M,Li Z H,Chen N J,Qu C W, (2020). A discrete multi-objective rider optimization algorithm for hybrid flowshop scheduling problem considering makespan, Noise and Dust Pollution. IEEE Access: Practical Innovations, Open Solutions, 8: 88527–88546
[12]
Gauthier S,Leduc M,Perfetto S J,Godwin A, (2022). Use of virtual reality to increase awareness of line-of-sight hazards around industrial equipment. Safety, 8( 3): 52–52
[13]
Gondia A,Moussa A,Ezzeldin M,El-Dakhakhni W, (2023). Machine learning-based construction site dynamic risk models. Technological Forecasting and Social Change, 189: 122347
[14]
Guo H Q,Nguyen H,Vu D A,Bui X N, (2021). Forecasting mining capital cost for open-pit mining projects based on artificial neural network approach. Resources Policy, 74: 101474
[15]
Hallowell M R,Hinze J W,Baud K C,Wehle A, (2013). Proactive construction safety control: Measuring, monitoring, and responding to safety leading indicators. Journal of Construction Engineering and Management, 139( 10): 04013010–04013010
[16]
Hong Y,Xie H Y,Agapaki E,Brilakis I, (2023). Graph-based automated construction scheduling without the use of BIM. Journal of Construction Engineering and Management, 149( 2): 05022020–05022020
[17]
Huo D,Sari Y A,Kealey R,Zhang Q, (2023). Reinforcement learning-based fleet dispatching for greenhouse gas emission reduction in open-pit mining operations. Resources, Conservation and Recycling, 188: 106664
[18]
Hwang S I,Son J H,Lee S H, (2014). development of scheduling model for earth work using genetic algorithm. KSCE Journal of Civil Engineering, 18( 6): 1618–1624
[19]
Imam M,Baina K,Tabii Y,Ressami E,Adlaoui Y,Benzakour I,Abdelwahed E, (2023). The future of mine safety: A comprehensive review of anti-collision systems based on computer vision in underground mines. Sensors, 23( 9): 4294–4294
[20]
InayahZCahyonoY HBudionoNDP(2023). Work shift system and heavy equipment operator fatigue in mine. 6(2): 223–229
[21]
JoshiDPaithankarAChatterjeeSEqueenuddinS M J M (2022). Integrated parametric graph closure and branch-and-cut algorithm for open pit mine scheduling under uncertainty. 2(1): 32–51
[22]
KulikovaE YBalovtsevS V (2020). Risk control system for the construction of urban underground structures. In: International Conference on Construction, Architecture and Technosphere Safety, Iccats 2020, 962
[23]
LeeCChaMKimHChoH (2023). Detailed activity-based earthwork scheduling model to aid during the planning stage of road construction. Applied Sciences, 13(23)
[24]
Lee D,Khan N,Park C, (2023). Rigorous analysis of safety rules for vision intelligence-based monitoring at construction jobsites. International Journal of Construction Management, 23( 10): 1768–1778
[25]
Lee J H,Park S M,Kim S K,Kang L S, (2022). Linear 4D system using schedule-location charts for infrastructure projects. Automation in Construction, 141: 104413
[26]
Lei K,Guo P,Zhao W C,Wang Y,Qian L M,Meng X Y,Tang L S, (2022). A multi-action deep reinforcement learning framework for flexible Job-shop scheduling problem. Expert Systems with Applications, 205: 117796
[27]
Li D S,Lu M, (2019). Classical planning model-based approach to automating construction planning on earthwork projects. Computer-Aided Civil and Infrastructure Engineering, 34( 4): 299–315
[28]
Liang Y,Liu Q X, (2022). Early warning and real-time control of construction safety risk of underground engineering based on building information modeling and internet of things. Neural Computing & Applications, 34( 5): 3433–3442
[29]
Liu G W,Guo W Q,Fu E S,Yang C D,Li J M, (2023). Dynamic optimization of open-pit coal mine production scheduling based on ARIMA and fuzzy structured element. Frontiers of Earth Science, 10: 1040464
[30]
LiuPXieMBianJLiHSongL (2020). A hybrid PSO–SVM model based on safety risk prediction for the design process in metro station construction. International Journal of Environmental Research and Public Health, 17(5): 1714
[31]
Liu Y,You K,Jiang Y T,Wu Z G,Liu Z Y,Peng G,Zhou C, (2022). Multi-objective optimal scheduling of automated construction equipment using non-dominated sorting genetic algorithm (NSGA-III). Automation in Construction, 143: 104587
[32]
Lu Y A,Qin W B,Zhou C,Liu Z H, (2023). Automated detection of dangerous work zone for crawler crane guided by UAV images via Swin Transformer. Automation in Construction, 147: 104744
[33]
Lu Y A,You K,Zhou C,Chen J X,Wu Z A,Jiang Y T,Huang C, (2024). Video surveillance-based multi-task learning with swin transformer for earthwork activity classification. Engineering Applications of Artificial Intelligence, 131: 107814
[34]
Manríquez F,González H,Morales N, (2023). Short-term open-pit production scheduling optimizing multiple objectives accounting for shovel allocation in stockpiles. Optimization and Engineering, 24( 1): 681–707
[35]
Mirzaei-Nasirabad H,Mohtasham M,Askari-Nasab H,Alizadeh B, (2023). An optimization model for the real-time truck dispatching problem in open-pit mining operations. Optimization and Engineering, 24( 4): 2449–2473
[36]
MohtashamMMirzaei-NasirabadHAskari-NasabHAlizadehB J I J o MReclamationEnvironment (2022). Multi-stage optimization framework for the real-time truck decision problem in open-pit mines: A case study on Sungun copper mine. 36(7): 461–491
[37]
Mu Y F,Yao T A,Jia H J,Yu X D,Zhao B,Zhang X S,Ni C W,Du L J, (2020). Optimal scheduling method for belt conveyor system in coal mine considering silo virtual energy storage. Applied Energy, 275: 115368
[38]
Nikolaev A V,Vöth S,Kychkin A V, (2023). Application of the cybernetic approach to price-dependent demand response for underground mining enterprise electricity consumption. Journal of Mining Institute, 261: 403–414
[39]
Ogunmodede O,Lamas P,Brickey A,Bogin G Jr,Newman A, (2022). Underground production scheduling with ventilation and refrigeration considerations. Optimization and Engineering, 23( 3): 1677–1705
[40]
Ou G Z,Jiao Y Y,Zhang G H,Zou J P,Tan F,Zhang W S, (2021). Collapse risk assessment of deep-buried tunnel during construction and its application. Tunnelling and Underground Space Technology, 115: 104019
[41]
Paithankar A,Chatterjee S, (2019). Open pit mine production schedule optimization using a hybrid of maximum-flow and genetic algorithms. Applied Soft Computing, 81: 105507
[42]
Paithankar A,Chatterjee S,Goodfellow R,Asad M W A, (2020). Simultaneous stochastic optimization of production sequence and dynamic cut-off grades in an open pit mining operation. Resources Policy, 66: 101634
[43]
Qin W B,Chen E J,Wang F,Liu W L,Zhou C, (2024). Data-driven models in reliability analysis for tunnel structure: A systematic review. Tunnelling and Underground Space Technology, 152: 105928
[44]
QuinteroM F (2023). Analysis of the root causes of fatal injuries in the United States surface mines between 2008 and 2021
[45]
Rachmawati T S N,Park H C,Kim S, (2022). A scenario-based simulation model for earthwork cost management using unmanned aerial vehicle technology. Sustainability, 15( 1): 503–503
[46]
Sadeghi H,Mohandes S R,Hosseini M R,Banihashemi S,Mahdiyar A,Abdullah A, (2020). Developing an ensemble predictive safety risk assessment model: Case of malaysian construction projects. International Journal of Environmental Research and Public Health, 17( 22): 8395–8395
[47]
Saurin T A,Formoso C T,Cambraia F B, (2008). An analysis of construction safety best practices from a cognitive systems engineering perspective. Safety Science, 46( 8): 1169–1183
[48]
Shehadeh A,Alshboul O,Tatari O,Alzubaidi M A,Salama A H E S, (2022). Selection of heavy machinery for earthwork activities: A multi-objective optimization approach using a genetic algorithm. Alexandria Engineering Journal, 61( 10): 7555–7569
[49]
StipanovicIBukhshZ ARealeCGavinK (2021). A Multiobjective decision-making model for risk-based maintenance scheduling of railway earthworks. Applied Sciences, 11(3)
[50]
WangG RYanX XKouZ MDengH SWangK S (2023). Research on operation conflict of auxiliary transport locomotive in complex mine based on extended petri net. Machines, 11(5)
[51]
WangX YZhangM H (2022). How powerful are spectral graph neural networks. International Conference on Machine Learning, Vol 162
[52]
WooSYeonJJiMMoonI CParkJ (2018). Deep reinforcement learning with fully convolutional neural network to solve an earthwork scheduling problem. In: 2018 Ieee International Conference on Systems, Man, and Cybernetics (Smc), 4230–4236
[53]
Wu W J,Sun D H,Jin K Q,Sun Y,Si P B, (2022). Proximal policy optimization-based committee selection algorithm in blockchain-enabled mobile edge computing systems. China Communications, 19( 6): 50–65
[54]
Xu X C,Gu X W,Wang Q,Zhao Y Q,Kong W Y,Zhu Z G,Wang F D, (2023). Enumeration optimization of open pit production scheduling based on mobile capacity search domain. Scientific Reports, 13( 1): 91–91
[55]
Yan X X,Wang G R,Jiang K S,Kou Z M,Wang K S,Zhang L X, (2023). Multi-objective scheduling strategy of mine transportation robot based on three-dimensional loading constraint. Minerals, 13( 3): 431–431
[56]
Yang S, Xu Z Intelligent scheduling for permutation flow shop with dynamic job arrival via deep reinforcement learning. Proceedings of 2021 IEEE 5th Advanced Information Technology, Electronic and Automation Control Conference (IAEAC), IEEE, 2672–2677
[57]
You K,Zhou C,Ding L Y, (2023). Deep learning technology for construction machinery and robotics. Automation in Construction, 150: 104852
[58]
ZamaniVTaghaddosHGholipour Y J A iC (2023). Sensor-fused simulation-based decision support system for earthmoving projects. Automation in Construction, 156, 105117
[59]
Zhang X T,Guo A,Ai Y F,Tian B,Chen L, (2022). Real-time scheduling of autonomous mining trucks via flow allocation-accelerated tabu search. IEEE Transactions on Intelligent Vehicles, 7( 3): 466–479
[60]
Zhou B H,Liao X M, (2020). Particle filter and Levy flight-based decomposed multi-objective evolution hybridized particle swarm for flexible job shop greening scheduling with crane transportation. Applied Soft Computing, 91: 106217
[61]
Zhou C,Gao Y Y,Chen E J,Ding L Y,Qin W B, (2023). Deep learning technologies for shield tunneling: Challenges and opportunities. Automation in Construction, 154: 104982
[62]
Zhou C,Gao Y Y,Zhou Y,She W,Shi Y S,Ding L Y,Miao C W, (2024). Properties and Characteristics of Regolith-Based Materials for Extraterrestrial Construction. Engineering, 37: 159–181
RIGHTS & PERMISSIONS
Higher Education Press
AI Summary 中Eng×
Note: Please be aware that the following content is generated by artificial intelligence. This website is not responsible for any consequences arising from the use of this content.