Next Article in Journal
Predicting PM2.5, PM10, SO2, NO2, NO and CO Air Pollutant Values with Linear Regression in R Language
Next Article in Special Issue
A Hybrid CNN-LSTM Random Forest Model for Dysgraphia Classification from Hand-Written Characters with Uniform/Normal Distribution
Previous Article in Journal
Complex Job Shop Simulation “CoJoSim”—A Reference Model for Simulating Semiconductor Manufacturing
Previous Article in Special Issue
DeepBreastCancerNet: A Novel Deep Learning Model for Breast Cancer Detection Using Ultrasound Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling of Multi-Level Planning of Shifting Bottleneck Resources Integrated with Downstream Wards in a Hospital

1
Department of Industrial Engineering, University of Engineering and Technology, Taxila 47050, Pakistan
2
Department of Computer Science, National Textile University, Faisalabad 37610, Pakistan
3
Department of Computer Science, Al Ain University, Abu Dhabi P.O. Box 112612, United Arab Emirates
4
Department of Software Engineering, Military College of Signals, National University of Science and Technology, Rawalpindi 04623, Pakistan
5
Department of Computer Science, College of Computer and Information Sciences, Princess Nourah Bint Abdulrahman University (PNU), Riyadh 11671, Saudi Arabia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(6), 3616; https://doi.org/10.3390/app13063616
Submission received: 8 February 2023 / Revised: 7 March 2023 / Accepted: 8 March 2023 / Published: 12 March 2023
(This article belongs to the Special Issue Artificial Intelligence and Robotics in Healthcare)

Abstract

:

Featured Application

The proposed model in this research is significant for hospitals for effective and efficient planning and scheduling of their resources. It has practical significance in a hospital where the whole system’s performance is affected by one bottleneck (i.e., either the operating room, ICU, or ward), and this bottleneck keeps changing in different planning horizons.

Abstract

Planning and scheduling critical resources in hospitals is significant for better service and profit generation. The current research investigates an integrated planning and scheduling problem at different levels of operating rooms, intensive care units, and wards. The theory of constraints is applied to make plans and schedules for operating rooms based on the capacity constraints of the operating room itself and downstream wards. A mixed integer linear programming model is developed considering shifting bottleneck resources among the operating room, intensive care unit, and hospital wards to maximize the utilization of resources at all levels of planning. Different sizes of planning and scheduling problems of the hospital, including small, medium, and large sizes, are created with variable arrivals and surgery durations and solved using a CPLEX solver for validating the developed models. Later, the application of the proposed models in the real world to develop planning systems for hospitals is discussed, and future extensions are suggested.

1. Introduction

Providing high-quality and affordable health care is one of the greatest challenges faced by the world’s nations. In developed countries, the aging populations are straining the healthcare budget, while many third-world countries face constant hardships in the form of low life expectancy, high infant mortality due to harsh environmental conditions, and scarce medical resources. The healthcare system must become more efficient at delivering care and preventing disease to solve these human and economic health problems [1]. A very important way to improve the efficiency of a healthcare system is by improving the planning and scheduling procedures in healthcare facilities for the efficient use of the critical resources of hospitals, including operating rooms, intensive care units (ICU), wards, patients, doctors, nurses, etc. Specifically, the operating room is considered one of the most critical units of hospitals in terms of cost and profit generation.
Hence, the current study is focused on the planning and scheduling problems of operating rooms.
Planning and scheduling problems in healthcare facilities are divided into three different levels including strategic, tactical, and operational level planning [2]. Strategic level planning is also called higher-level planning and is based on forecast and surgical demands, where the operating room time is allocated to different surgical specialties. Tactical level planning is also termed medium-level planning, and it involves the division of allocated time for different surgical specialties and preparing a master surgery schedule. Operational level planning is also termed lower-level planning, where the daily plans of the operating room are made and patients are sequenced within the allocated operating rooms. Some articles investigate decision-making regarding operating room planning and scheduling in hospitals at four hierarchical decision levels: strategic, tactical, offline operational, and online operational. The offline operational level deals with assigning patients to dates and sequencing the patients in the operating rooms and involves monitoring and controlling the schedule execution [3,4].
The literature contains studies of the different levels of planning and scheduling problems to improve the performance of healthcare facilities [5,6,7]. For example, McRae and Brunner [8] studied case mix planning decisions for higher-level planning and proposed a mixed integer linear programming model to solve this problem. Freeman et al. [9] proposed an iterative approach for case mix planning decisions. They developed mathematical models to generate a pool of candidate solutions and used a simulation method to evaluate the effectiveness of each candidate solution for higher-level and lower-level performance measures, i.e., the overutilization of operating rooms and utilization of downstream wards. The case mix planning problem was also investigated by Yahia et al. [10]. They developed a stochastic model to find the optimal case mix of patients considering surgery duration, length of stay, and arrival time as uncertain variables.
Medium-level planning problems have also been studied in the literature. For example, Santos and Marques [11] developed a master surgery schedule and proposed a stochastic programming model to integrate the medium-level decision-making of operating rooms with downstream departments such as intensive care units and wards. Spratt and Kozan [12] studied the master surgery scheduling problem on a medium level along with the surgical case assignment problem. They proposed a mixed integer nonlinear programming model and hybrid metaheuristics to solve the problem. Additionally, Fügener et al. [13] studied the medium-level planning problem to calculate the demand distribution of downstream departments for a given master surgery schedule in a hospital.
Lower-level planning has also been investigated in the literature. For example, Pham and Klinkert [14] studied operation room planning and the scheduling problem as a generalized job shop scheduling problem. They developed a mixed integer linear programming model to minimize the makespan of surgeries at a lower-level planning horizon. Oliveira et al. [15] studied the problem of elective patient scheduling at a lower-level planning horizon to maximize operating room utilization based on patient prioritization. They developed an integer programming model and tested it on case hospital data. Younespour et al. [16] considered the scheduling problem of patients at a lower level to minimize the overtime costs, makespan, and completion time costs of surgeons considering parallel surgeries. They proposed a mixed integer programming model and constraint programming models for the scheduling problem.
All the planning level decisions are interrelated, and the solution of one level influences the solution of other levels. In the literature, all these planning levels are dealt with in a hierarchical manner where the output solution of one planning level is moved to the next lower planning level as input. All three decision levels are linked to each other so that decisions made at the strategic level influence the quality of decisions made at the tactical level. Both the higher-level decisions influence further operational-level decisions. Such interdependence of these decision levels provides the basis for research to investigate all these levels concurrently and solve the problem of multi-level planning and scheduling.
Some researchers have focused on multi-level planning in healthcare facilities. For example, Ma and Demeulemeester [17] investigated all levels of decision-making, such as case mix planning, master surgery scheduling, and patient sequencing, at the operational level in a hierarchical manner. All three planning levels were integrated iteratively in three phases with different objective functions for each phase. In addition, Fügener [18] integrated higher-level planning and medium-level planning decisions. They proposed a method to calculate the distribution of patient demand in downstream departments and a model to generate a master surgery schedule with an assumption of fixed capacities. Guido and Conforti [19] studied planning and scheduling problems at medium and lower levels. They proposed an algorithm for allocating operating rooms to specialty groups for a pre-determined period and scheduling patients in the allocated operating rooms. They developed a cyclic master surgery schedule where the patients were allocated to operating rooms based on the developed master surgery schedule.
Reviewing the literature on operating room planning and scheduling shows that the importance of all levels of decision-making can be inferred. In the literature, however, when efforts to integrate any of two or three levels are made, and integration occurs in a hierarchical manner, it proceeds that the results of one planning level are used to make decisions at the lower planning level. Hierarchically, the quality of decisions made at the lower level is greatly determined and influenced by the higher-level decisions. All levels’ performance can be increased if all decision-making levels are integrated. This integration can help to improve the robustness and flexibility of operating room schedules [3]. To increase the performance of the healthcare system, the integration of all three decision-making levels is important, and, hence, the present study is focused on integrating such levels.
In addition, most of the literature has considered the capacity of constraint resources of the hospital during planning and scheduling decisions. For example, the planning and scheduling of the operating theatre complex have gained much attention, as it is considered one of the hospital’s critical resources and its performance has a significant impact on other departments [20,21]. An in-depth review of the literature on operating theatre planning and scheduling is available in many articles [3,22,23]. The literature on operating theatre planning and scheduling has covered many perspectives to improve the operations of the healthcare facility, considering maximizing the utilization of operating theatre [9,15,24,25] and minimizing the patient waiting time [26], surgeon overtime [16], the cost of operating theatre [27], the makespan [28], the overtime [29], etc. Roshanaei and Naderi [30] integrated the problem of patient allocation to a day in an operating room and sequenced the allocated patients within the operating rooms to maximize the total scheduled surgical time. They developed mixed integer and constraint programming models and various bender decomposition algorithms to solve the developed models to optimality. The authors of [31] studied the operating room planning and scheduling problem integrated with downstream wards. They proposed a two-stage artificial bee colony algorithm for solving the problem. In [32], researchers developed a robust optimization model that combines staffing and scheduling decisions to minimize the impact of variations in surgery duration, staff availability, and emergency arrivals.
In most of the literature, the planning and scheduling of surgeries focus on the optimal use of the available operating theatre capacity. However, the performance of healthcare facilities is constrained not only by the capacity of the operation theatre complex but also by the capacity constraint of other critical resources of hospitals [33]. For example, in hospitals, scheduled surgeries can be canceled in large numbers due to the unavailability of beds for post-operation recovery [34]. Furthermore, in most hospitals, the scheduling systems consider the available capacity of beds, while most surgery planning systems used in hospitals consider the capacity limit of the operating theatre. Since there is an interconnection between the operation theatre complex, the available number of beds in wards, and the available number of beds in the ICU, it is, therefore, significant to study the capacity consideration of different interlinked resources, including the capacity of the operating theatre, the number of available beds in the ICU, and the number of available beds in the wards, etc. Therefore, the capacity consideration of critical resources without consideration of their interconnection with other resources may lead to their suboptimal use [35]. Thus, independent optimization of the operating theatre’s resources or the ward’s optimization may not lead to global optimization of the healthcare system. Therefore, some researchers have addressed the operating theatre planning problem with other critical resources, such as the ICU and wards, at a different level of the planning problem [11,36,37].
Testi et al. [38] developed a three-phase hierarchical approach for the weekly scheduling of operating rooms where the optimal case mix is identified in the first phase to maximize the overall benefit. In the second phase, they developed a master surgery schedule. In the third phase, they performed a simulation to evaluate the operational performance of the master surgery schedule and the utilization of downstream departments. Chow et al. [39] proposed a simulation and mixed integer programming model to integrate and improve surgical scheduling and predict ward and leveling bed occupancy, respectively. Fügener et al. [13] investigated a medium-level planning problem to calculate the demand distribution of downstream departments for a given master surgery schedule. They developed a method to optimize the master surgery schedule and analyzed the impact of resulting block allocation on the bed requirements in the ICU and general wards. Fügener [18] integrated the strategic and tactical master surgery scheduling while considering the impact on downstream resources such as the ICU and general patient wards. The study aimed to increase hospital earnings by optimizing the master surgery schedule. Freeman et al. [9] studied the case mix planning problem with consideration of the resources of downstream wards. They developed mathematical models to generate a pool of solutions with different case mix plans and employed a simulation to evaluate each solution for the overutilization of the operating theatre and variability in bed usage in downstream wards. In their multi-phase solution approach, they developed four mixed integer programming models for case mix planning, block allocation, and master surgery scheduling, and used simulations to assess the quality of the generated solution. In another article, Santos and Marques [11] developed a master surgery schedule and proposed a stochastic programming model to integrate the medium-level decision-making of operating rooms with downstream departments such as intensive care units and wards. In their approach, they estimated the bed requirements at the operational level and developed master surgery schedules based on these estimates of bed requirements.
In most of the literature, the planning and scheduling problems did not consider the interconnection of the critical resources with the other interlinked wards. For a feasible and realistic solution to the health care system, the interconnection of the operating room with the ICU and ward cannot be ignored, and the optimization of the operating theatre alone may result in the underutilization or congestion of other resources and may lead to infeasible schedules. Hence, it is proposed in the present study that the optimal utilization of operating rooms integrated with the ICU and ward can lead to feasible and optimal solutions and enhance the performance of the healthcare facility as a whole because the operating rooms also pace the activities of the linked units. Optimizing the utilization of the operating room alone without considering the capacity constraints of downstream units may result, in some cases, in the underutilization of these units and, in other cases, congestion in these downstream units, which ultimately leads to early discharge or even surgery cancellations. For the balanced use of all the hospital resources, it is necessary to consider the planning and scheduling of operating rooms in combination with interconnected downstream wards. This work, therefore, aims to integrate higher-level, medium-level, and operational-level planning and scheduling along with the integration of the planning and scheduling of operating rooms with downstream units such as beds in the ICU and ward. Such a problem of the simultaneous consideration of all decision levels and critical resources of the healthcare system is new and, to the best of the authors’ knowledge, has not been addressed so far in the literature.
In the literature, the theory of constraints was used to improve the overall system performance while focusing on only the critical resource [40,41]. The theory of constraints applies the drum buffer rope method to identify the resource with limiting capacity, called the drum, and a rope mechanism to provide the planning information to upstream resources of the drum for effective planning and execution. To fully utilize the drum resource, buffers in the form of time and material are provided to prevent the drum from starvation. The drum buffer rope method has been applied successfully in different planning and scheduling problems [42,43,44,45]. For example, Ronen et al. [46] applied the drum buffer rope concept for aircraft scheduling, Gilland [43] used the DBR method in a serial production line for production planning and control, Pegels and Watrous [47] applied the DBR method for the optimization of assembly shops, Sirikrai and Yenradee [45] used the DBR method on production planning and control, and Georgiadis and Politou [42] used the DBR method to solve production planning problems; they introduced a dynamic DBR approach for a time buffer for production planning and control in two machine-capacitated flow shops. Saif et al. [44] recently applied DBR concepts for multi-level integrated production planning and scheduling problems in a flow shop considering industry 4.0 concepts. However, in the literature for planning and scheduling on multiple planning levels in the healthcare industry, the DBR method of the theory of constraints is rarely found; for example, [45,48] used the theory of constraints and identified the bottlenecks among the human resources, such as anesthesiologists, doctors, and nurses, which has limited its application in operating room planning. In the literature, no study is reported on the method in which planning problems either on a single decision level or in an integrated manner are solved using the theory of constraints. The current research used the theory of constraints method to solve multi-level planning and scheduling problems in hospitals, integrating all decision levels and critical resources of the healthcare system. To the best of the authors’ knowledge, the theory of constraints has not yet been applied to the multi-level planning and scheduling of critical resources of hospitals while considering the capacity constraints of downstream wards. It is considered for the first time in the literature in the current study.
In the literature, mathematical programming, particularly mixed integer linear programming, has been widely used for planning and scheduling in healthcare facilities [3]. In mathematical programming, minimizing, or maximizing objectives are formulated subject to the constraints related to the considered problems. The developed model is then solved using standard software such as CPLEX [49,50] or by developing exact solution methods [17,51]. Further, heuristic methods, metaheuristics [12,52], and mathematical modeling and analytical procedures such as the Markov Decision process and queuing theory have also been used for the planning and scheduling of operating rooms [53,54,55,56]. Careful analysis of the research methodologies employed by various researchers reveals that most studies focus on obtaining the optimal solution to the formulated problems using heuristics or exact solutions. In the case mix level, the authors developed mixed integer programming and stochastic models in a study and used a sample average approximation and simulation to solve the models [8]. Many others developed mathematical models and heuristics to solve a case mix planning problem [9,10]. On the tactical level, mathematical modeling and exact method or heuristics are used to solve the developed models [11,12,13]. On the operational level, mathematical modeling and heuristics are employed to solve operating room allocation and sequencing problems [14,15,16]. On integrated levels, the methodology employed is similar [17,18]. For example, Roshanaei and Naderi [30] solved the integrated problem of patient allocation to the operating room, day, and sequencing of patients using mixed integer programming and constraint programming and developed a heuristic algorithm to solve the developed models. However, heuristics and metaheuristics give near-optimal solutions to the considered problem and use the constraints presented in the proposed model of the problem. The exact methods for linear programming models provide an accurate solution to the problem, which is significant for validating new mathematical models. Since the current problem is new in the literature and our study proposes a mathematical model for it, the present study uses CPLEX to solve the problems used for its validation. The CPLEX solver can be used to code and solve integer programming, mixed integer programming, multi-objective optimization, and quadratic programming problems. However, for large-sized problems, it requires more computational time, data handling from various data sources, and a greater number of parameters to define the problem.
The objective of the current study is to develop and solve a mathematical model using the theory of constraints for operating room planning and scheduling considering all the decision levels, such as the strategic, tactical, and operational levels, and integrating them with downstream wards.
The current study contributes to the planning and scheduling literature in the following ways.
  • The current research is new in integrating all the planning levels of the hospital, considering the higher-level, medium-level, and lower-level planning considering constraints of the interlinked resources, including the operation theatre, ICU, and wards;
  • The current research is new to applying the theory of constraints for multi-level planning in hospitals;
  • The current research proposes a new mixed integer linear programming model for multi-level planning and scheduling in hospitals considering the theory of constraints concept;
  • The current research develops a new mixed integer linear programming model considering the capacity constraints of the operating room, ICU, and wards.
The rest of the paper is organized as follows: Section 2 presents the problem description and mathematical model. Section 3 and Section 4 present the solution approach, computational experiments, and results. Finally, Section 5 shows the conclusion of the research with important findings and highlights the limitations and future extensions of the work.

2. Problem Description

In a healthcare facility, elective and emergency patients flow through different departments. The general flow of patients moving through critical units of a hospital is illustrated in Figure 1. Figure 1 indicates that elective patients are added to the waiting list while emergency patients are sent to the planning process and assigned to the required resources. It further explains the flow of patients from one resource such as the operating room to the other resources such as the ICU and ward. For capacity planning during higher-level planning decisions, the forecasted patients are considered. It is considered that the patients can directly come for treatment in the operating room (OR), intensive care units (ICU), or wards.
It can be seen from Figure 1 that the patients arriving in the OR have the possibility to move into the ICU or the wards or can be discharged after treatment in the OR; the patients arriving in the ICU can move into a ward or can be discharged after treatment in the ICU. Current research considered the multi-level planning of patients on critical resources. The process considered in the current research for the multi-level planning of critical resources is presented in Figure 2.
Figure 2 illustrates that the higher-level planning involves the planning of patients on each critical resource of the hospital. At a higher level, the patients are allocated to the planning horizons based on their arrival, considering the infinite capacity of all critical resources. At the medium level, the patients from higher levels are allocated to their required critical resources, including the operating room, ICU, and ward. Then, according to the theory of constraints, the bottleneck resource is identified. The bottleneck resource is identified based on patients’ workload on all resources and utilization of the resources. After identifying the bottleneck resources, at the lower level, the respective planning and scheduling model is performed, and a plan is released which is communicated to the other upstream departments and resources. All other resources follow the released plan and scheduling of the bottleneck resources. Figure 2 further explains the patients for each resource; for example, the patients in the ICU include the patients from the operating room and randomly arriving patients. The mathematical relations for the multi-level planning and scheduling of patients on each critical resource are presented in this section.

2.1. Higher-Level Planning

At a higher planning level, the randomly arriving patients are inserted on different days of the planning horizon. The surgery patients and patients needing the ICU and ward are inserted into the planning horizon as explained in this section.

2.1.1. Allocation of Patients to Operating Room

On the higher level, the available information of arriving surgery patients, such as arrival date, due date, preoperative time, surgery duration, etc., is used to insert the patients into different days of the planning horizon. The capacity of all resources is considered infinite at a higher level, and all patients are placed in suitable planning horizons. Equation (1) explains the objective function of higher-level planning.
o b j = maximize ( p = 1 P H p d τ )   0 < d D     τ
d = 1 D H p d τ = 1   p P   τ
H p d τ = 1   p P   D D p d E P S T p   τ
EPST p = A p + t p p r e o p   p P
The objective function (see Equation (1)) is aimed to maximize the number of patients inserted into each day of the planning horizon. Equation (2) gives the constraint that each patient can be inserted for one day of the planning horizon only. Equation (3) gives the constraint that a patient can be assigned to a planning horizon only if the patient’s earliest possible start time (EPST) is the day d of the planning horizon. It is worth mentioning that ( D D p d E P S T p ) in Equation (3) means that a patient can be inserted into a day equal to EPST or greater than EPST in the planning horizon. It also states that the patient’s due date is considered while assigning a patient to a day. This equation has practical significance when there is a larger number of patients and limited capacity. Equation (4) shows the method to calculate the EPST of a patient. The earliest possible start time of a patient is calculated from the arrival time of the patient and after the time required for the preoperative activities of the patient. For the patients arriving at the ICU and ward, the preoperative time is zero, and patients are inserted into the planning horizon based on their arrival date using Equation (4).

2.1.2. Allocation of Patients to the ICU

The patients who require the ICU in a planning horizon are already current patients receiving care in the ICU, patients from operating rooms, and randomly arriving patients during a planning horizon. The number of patients who need ICU after surgery is determined using Equation (5), where N P d τ o i is the number of surgery patients who require the ICU on day d of the planning horizon τ . The total number of patients who require the ICU on day d in a planning horizon τ is calculated using Equation (6), where N P d τ i exp is the expected number of randomly arriving patients in the ICU while N P d τ i o c c represents the number of patients already occupying a bed in the ICU.
N P d τ o i = e = 1 E Pr N P e × N P e   0 < d D ,   τ
N P d τ i = N P d τ o i + N P d τ i exp + N P d τ i o c c   0 < d D ,   τ
The surgery patients based on their date of surgery and the randomly arriving patients based on their date of arrival are inserted into the days of the planning horizon.

2.1.3. Allocation of Patients to the Ward

The patients in the ward are those patients who are already receiving care in the ward, patients who will arrive at the ward after surgery, and randomly arriving patients. The number of patients who require a ward bed after surgery is determined using Equation (7), where N P d τ o w is the number of patients who require a ward bed after surgery. The total number of patients that require a ward in a day of the planning horizon τ is calculated using Equation (8), where N P d τ w exp is the expected number of patients randomly arriving at the ward, N P d τ w o c c is the number of patients already occupying beds in the ward, and N P d τ i w is the number of patients that require a ward bed after the ICU.
N P d τ o w = e = 1 E Pr N P e × N P e   0 < d D ,   τ
N P d τ w = N P d τ o w + N P d τ i w + N P d τ w exp + N P d τ o c c   0 < d D   τ
The output of higher-level planning is a list of patients for the operating room, ICU, and ward for each day d of the planning horizon τ .

2.2. Medium-Level Planning

Based on the lists of patients from the higher-level planning for each operating room, ICU, and ward, the bottleneck is identified from these three units based on the load of each planning horizon. The percentage utilization of all three units, such as the operating rooms, ICU, and ward, is calculated to identify the system’s bottleneck. The unit with the highest percentage utilization is identified as the bottleneck. The workload of surgery patients in operating rooms can be calculated using Equation (9). The total workload of all operating rooms is calculated by adding the workload of all patients on the day of the planning horizon as given by Equation (9).
TWL o d τ = p = 1 P ( ( t p p r e + t p s e t + t p s u r + t p c l ) × H p d τ )   0 < d D   τ
UTL o d τ % = TWL o d τ A t o d τ × N o d τ × η o a v g   o O d τ   0 < d D   τ
UTL i d τ % = N i d τ w l N i d τ × η i a v g   i I d τ   0 < d D   τ
UTL w d τ % = N w d τ w l N w d τ × η w a v g   w W d τ   0 < d D ,   τ
N i d τ w l = N P d τ i   0 < d D   τ
N w d τ w l = N P d τ w   0 < d D   τ
The percentage utilization of operating rooms is calculated as the ratio of the total workload of operating rooms to the total available capacity of all operating rooms, and is given by Equation (10), where N o d τ is the total number of operating rooms on the day d of the planning horizon τ and η o a v g is the average efficiency of the operating rooms. Similarly, the percentage utilization of the ICU and the downstream ward is calculated in Equation (11) and Equation (12), respectively. Where N i d τ and N w d τ are the total number of beds in the ICU and ward, respectively, on the day d of the planning horizon τ ; N i d τ w l and N w d τ w l are the number of occupied beds in the ICU and downstream ward, respectively; and η i a v g and η w a v g are the average efficiency of the ICU and ward, respectively. Equations (13) and (14) provide the method to calculate workload in the ICU and ward, respectively.
After calculating the percentage utilization of the operating room, ICU, and ward, the bottleneck is identified. The resource with the highest percentage utilization is identified as a bottleneck on that day of the planning horizon. The following is the condition for a resource to be considered a bottleneck. Consider that R is a set of critical resources consisting of the OT, ICU, and ward, and R* is the resource with maximum utilization and is the bottleneck resource as given in Equation (15)
R * = r e s o u r c e     w i t h     max { U T L % } r d τ   r R = [ O T , I C U , W a r d ]
The condition for any resource to be bottlenecked is represented by the expression given below:
if { U T L % r > U T L % r , w h e r e   r R ,   r [ R r ]       r = R *       f o r   a n y   r r = R U T L % r < U T L % r , w h e r e   r R ,   r [ R r ]       r R *                                   r r = ϕ
After identifying the bottleneck, the respective planning and scheduling model is triggered, as explained here.

2.2.1. Operating Room Is the Bottleneck Resource (R*)

As a result of the allocation of patients and utilization of the operating room, ICU, and ward, if the operating room is identified as the bottleneck, the model for the operating room planning operates. Figure 3 explains the procedure of multi-level integrated planning when the operating room is identified as the bottleneck.
After identifying the bottleneck, the patients are allocated to the operating rooms considering the available capacity of operating rooms on each day d of the medium-level planning horizon τ . It is worth mentioning that if the available capacity of the operating rooms is less than the required capacity of the operating room, the patients are moved to the next day by considering the patients’ due dates. The patients are then allocated to the operating rooms to maximize the workload of the operating rooms such that the patients of the same specialty are allocated to the same operating room to minimize the operating room turnover time when a patient of one specialty is operated in the same operating room after a patient of a different specialty. The objective function of medium-level planning when the operating room is identified as bottleneck resource R* is given by Equation (16).
o b j = maximize ( W L ) o d τ   o O   o < d D   τ
W L o d τ = p = 1 P ( ( t p p r e + t p s e t + t p s u r + t p c l ) × X p o d τ )   o O   0 < d D ,   τ
W L o d τ A t o d τ   o O   0 < d D   τ
d = 1 D o = 1 O X p o d τ = 1   p P   τ
s = 1 S U s p o d τ = 1   p P   o O   0 < d D   t
p = 1 P U s p o d τ N s d τ a l l   s S   0 < d D ,   τ
U s p o d τ = 1   p P e   s S e   0 < d D ,   τ
o = 1 O U s p o d τ = 1   s S   0 < d D ,   τ
The objective function of the medium-level planning of operating rooms is to maximize the workload of all operating rooms. Equation (17) provides the expression to calculate the workload of each operating room. The objective function of the medium level is significant to increase the utilization of the operating rooms. Greater operating room utilization leads to less idle time and hence reduced costs. Equation (18) puts an upper bound on an operating room’s workload such that the operating room’s workload should be less than or equal to the available capacity on the day d of the planning horizon τ . Equation (19) ensures that each patient is allocated to one operating room and only one day of the planning horizon. Equation (20) ensures that each patient is assigned to one surgeon only. Equation (21) puts an upper bound to the number of surgeries a surgeon performs on a day d of the planning horizon τ . Equation (22) states that only a surgeon of the same specialty can be assigned for the surgery of a patient from a specialty. Equation (23) states that each surgeon can be assigned to one operating room on a day d of the planning horizon τ .
The output of medium-level planning is the allocation of patients to operating rooms and surgeons’ assignment for patients’ surgery. As a result of the allocation of surgery patients to operating rooms, the corresponding number of patients that require beds in the ICU and ward is calculated using Equations (7) and (8), respectively. The number of beds required by the surgery patients in the ICU and ward is fixed, and the information is sent to both the ICU and ward to freeze the capacity for patients in the operating room, which is the bottleneck in the considered planning horizon.
The medium-level plan based on the allocation of patients to operating rooms is released for downstream departments, such as the ICU and ward, to reserve capacity for planned patients and to upstream departments for the preparation of the surgery of planned patients before the actual start of surgery.

2.2.2. ICU Is the Bottleneck Resource (R*)

The planning is performed based on the capacity of the ICU when it is identified as the bottleneck resource. Figure 4 represents when the intensive care unit becomes the system’s bottleneck. In this case, the ICU’s capacity is fixed, and this capacity information is sent through a rope to the operating room. Operating room planning and scheduling now take the capacity information of the ICU, and the patient allocation to the operating rooms is performed based on the capacity of the bottleneck of the ICU.
o b j = maximize p = 1 P Y p i d τ   0 < d D   τ
i = 1 I d = 1 D Y p i d τ = 1   p P   t
p = 1 P Y p i d τ N i d τ   0 < d D   t
N i d τ o c c = N i ( d 1 ) τ o c c N i ( d 1 ) τ d i s c + N i d τ o + N i d t exp   i I   0 < d D   τ
N i d τ o c c N i d τ   i I d τ   0 < d D   τ
t p i d i s c = t p i a r r + L O S p i   i I   p P
The following are the objective function and constraints of medium-level planning when the ICU becomes a bottleneck. The objective function given by Equation (24) aims to maximize the number of patients requiring a bed in the ICU.
Equation (25) ensures that on a given day d of the planning horizon τ , a patient can be assigned to one bed only in ICU. Equation (26) ensures that the number of patients requiring the ICU is less than the total number of ICU beds. Equation (27) provides the method of calculating the number of occupied beds in the ICU, where N i d τ d i s c is the number of beds released after the discharge of patients from the ICU. N i d τ o and N i d τ exp are the number of beds required in the ICU for the patients from the operating room and other expected patients. Similarly, based on the information on the available capacity of the ICU, Equation (28) puts an upper bound on the number of beds occupied in the ICU that should be less than or equal to the total number of beds in the ICU. Where N i d τ is the total number of beds in the ICU on the day d of the planning horizon τ . N i d τ o c c is the number of occupied beds in the ICU on d of planning horizon τ . Equation (29) provides the patient’s discharge time from the ICU. Based on the available capacity of the ICU, the ICU patients’ schedule is released to the operating room, and operating room planning is performed based on the released schedule of the ICU. The information on the bed occupancy in the ICU is then sent to the operating rooms to plan surgeries as per the available capacity of the ICU.

2.2.3. Ward Is the Bottleneck Resource (R*)

In the integrated planning of the operating room with the ICU and ward, when the percentage utilization of beds in the ward is maximized, the ward becomes the bottleneck. Figure 5 explains the procedure for multi-level planning when the ward becomes the system’s bottleneck. In such a case, the capacity information of the ward beds is provided through the rope mechanism to the ICU and operating rooms, and some portion of the capacity allocation of the operating room becomes fixed. The operating room allocation of patients is then performed based on this capacity constraint of ward beds. The objective function given by Equation (30) is aimed at maximizing the number of patients that require a bed in the ward.
o b j = maximize p = 1 P Z p w d τ   0 < d D   τ
Z p w d τ = 1   p P , 0 < d D   τ
p = 1 P d = 1 D Z p w d τ N w d τ   p P   τ
N w d τ o c c = N w ( d 1 ) τ o c c N w ( d 1 ) τ d i s c + N w d τ o + N w d t i c u   i I   0 < d D   τ
N w d τ o c c N w d τ   i W d τ   0 < d D   τ
t p w d i s c = t p w a r r + L O S p w   w W   p P
Equation (31) ensures that on a given day d of the planning horizon τ , a patient can be assigned to one bed only in the ward. Equation (32) ensures that the number of patients requiring a bed in the ward is less than the total number of beds in the ward. Equation (33) provides the method to calculate the number of occupied beds in wards, where N w d τ d i s c is the number of beds in wards from where patients are discharged on the day d . N w d τ o and N w d τ i c u are the number of beds required in wards for the patients from the operating room and ICU, respectively.
Equation (34) puts an upper bound on the wards’ workload such that the number of beds occupied should be less than or equal to the total number of beds in the ward. Where N w d τ is the total number of beds in the ward on the day d of the planning horizon τ . N w d τ o c c is the number of occupied beds in downstream wards on d of planning horizon τ . Equation (35) provides the patient’s discharge time from the ward. The information on the bed occupancy in the ward is then sent to the operating rooms and ICU to plan surgeries and to make a schedule of the ICU as per the available capacity of the ward beds.

2.3. Lower-Level Planning

At the lower level of planning, the sequencing of patients is performed. Based on the load at the medium-level planning, any unit such as the operating room, ICU, and ward can become the bottleneck. The lower-level planning is performed based on the bottleneck resource at the medium-level planning.

2.3.1. Lower-Level Planning when the Operating Room Is the Bottleneck Resource (R*)

When the operating room becomes the bottleneck, the lower-level planning model for the operating room is triggered. The objective function of operating room planning at the lower level is given by Equation (36).
( o b j ) o d τ = minimize ( m a k e s p a n ) o d τ   o O d τ   0 < d D ,   τ
( m a k e s p a n ) o d τ = p = 1 P t p c o m p × X p o d τ   o O   0 < d D   τ
t p c o m p ( t s p a r r + t p s e t + t p s u r + t p c l ) × W p o d τ   p P
t p c o m p ( t s p a r r + t p p s d s + t p s u r + t p c l ) × S p p o d τ   p P
p = 1 P S p p o d τ = 1 V p o d t   o O   0 < d D   τ
p = 1 P S p p o d τ = 1 W p o d t   o O   0 < d D   τ
p = 1 O V p o d τ = 1   o O , 0 < d D   τ
p = 1 P W p o d τ = 1   o O   0 < d D   τ
t p c o m p t p c o m p + t p p s d s + t p s u r + t p c l M ( 1 S p p o d τ )   p , p P   0 < d D   τ
( m a k e s p a n ) o d τ 0   o O   0 < d D ,   τ
S t p i D t p o   p P o i
S t p w D t p w   p P o w
The objective function given by Equation (36) aims to minimize the makespan of the operating room. Equation (37) represents the method of calculating the makespan of an operating room. The objective function considered at lower-level planning is significant to minimize the completion time of all surgeries in the allocated operating rooms. The minimum makespan leaves time to accommodate more surgeries, thus leading to increased utilization of operating rooms and cost savings. Equation (38) represents the method of calculating the completion time of the surgery of a patient if they were operated on at the first position in the operating room o. Equation (39) represents the method of calculating the completion time of the surgery of patient p’ if they were operated on immediately after patient p in the operating room o. Equations (38) and (39) also state that each patient can be operated on after the arrival of the surgeon for the patient p. Equation (40) states that each patient must precede another patient unless they are operated on last. Equation (41) states that each patient must be preceded by another patient unless they are operated on first. Equations (42) and (43) impose the constraint that only one patient can be operated on as first and last in the operating room o. Equation (44) imposes that each patient’s surgery must be started before the surgery of the previous patient is completed. Equation (45) states that the makespan of an operating room must be non-negative. Equations (46) and (47) provide constraints for the ICU and ward, respectively, when the operating room is the bottleneck resource (R*). These constraints limit the starting time of surgery patients in the ICU and ward. The output of the lower-level plan is the optimal schedule of patients in the allocated operating rooms. The optimal schedule is released to the upstream units of the operating room as well as to the ICU and ward to accommodate the patients as per the capacity of the operating room, which is bottleneck resource (R*) on day d of the planning horizon τ .

2.3.2. Lower-Level Planning When the Intensive Care Unit Is the Bottleneck Resource (R*)

As a result of medium-level planning, if the ICU is identified as the bottleneck, the lower-level planning of the ICU is performed, the schedule is released for ICU, and information is sent through the rope mechanism to the operating room to follow the release schedule of the ICU. Operating room planning and scheduling now take the capacity information of the ICU, and patient allocation and sequencing to the operating rooms are performed based on the capacity of the bottleneck of the ICU.
The patients arriving at the ICU are admitted ICU based on their priority as per the procedure given below, where Q is the set of patients arriving at the ICU. α q is used to represent the relative priorities of the ICU patients in set Q. A patient q is assigned to the ICU if α q is the highest.
if { α q > α q   q Q ,   q [ Q q ]                   Assign ICU α q < α q   q Q ,   q [ Q q ]           Do not assign ICU
α = c p L O S p i
D t p o S t p i   p P o i
S t p w D t p i   p P i w
The relative priority of a patient to the ICU is calculated by Equation (48), where c p is the criticality weight of the patient assigned by the doctors and L O S p i is the length of the patient’s stay in the ICU. Equations (49) and (50) represent the extra constraints for the operating room and ward when the ICU is identified as the bottleneck.
Based on the available capacity of the ICU, the ICU patients’ schedule is released to the operating room, and operating room planning is performed based on the released schedule of the ICU.

2.3.3. Lower-Level Planning When the Ward Is the Bottleneck Resource (R*)

In the integrated planning of operating rooms with the ICU and ward, when the percentage utilization of beds in wards is maximized, the ward becomes the bottleneck. Figure 5 explains the procedure when the ward becomes the system’s bottleneck. In such a case, the capacity information of ward beds is provided through the rope mechanism to the ICU and operating rooms, and some portion of the capacity allocation of the operating room becomes fixed. The operating room allocation of patients is then performed based on this capacity constraint of ward beds. The following procedure is adopted for lower-level planning when the ward becomes the bottleneck. The patients arriving for a bed in the ward are allocated a ward bed based on their relative priority as given below.
if { α q > α q   q Q ,   q [ Q q ]                   Assign ward α q < α q   q Q ,   q [ Q q ]           Do not assign ward
α = c p L O S p w
D t p o S t p w   p P o w
D t p i S t p w   p P i w
The relative priority of a patient for the ward beds is calculated by Equation (51), where c p is the criticality weight of the patient assigned by the doctors and L O S p w is the length of stay of the patient in the ward. Equations (52) and (53) represent the extra constraints for the operating room and ICU when the ward is identified as the bottleneck.
The optimal sequence of patients in the operating room when the ward is the bottleneck is determined based on the release schedule of the ward, and the optimal schedule of operating rooms is released. Finally, Equations (54) and (55) give the binary domain decision variables and non-negativity conditions of parameters, respectively.
H p d τ , X p o d τ , Y p i d τ , Z p w d τ , W p o d τ , R p o d τ , U s p o d τ , S p p o d τ { 0 , 1 }   p P , o O , 0 < d < D ,   τ
W L o d τ 0

3. Solution Methods

In hospitals, any resource among the operating rooms, ICUs, and wards can become bottlenecks in the considered planning horizon. Random samples of different problem sizes with varying lengths of stay in the operating room and arrival patterns of patients for the ICU and ward were created to solve the developed mathematical models. Three different sets of problems were generated considering the operating room, ICU, and ward as the bottleneck resource. The surgery duration data from a real hospital was collected, and more samples were created using a lognormal distribution with the parameters of the real data. The data collected from the hospital contained the data of surgery patients for one year. Table 1 represents the characteristics of the three different-sized problems.

4. Computational Experiments

Various experiments were designed to evaluate the performance of the developed method and validate the models. The three major problem categories were small, medium, and large-sized problems. These problems were further divided based on the patients’ length of stay in the operating room and the identification of the bottleneck resource. Table 2 represents the experiments conducted to evaluate the developed method and models. A total of 270 experiments were conducted to test the developed models.
The numerical experiments in this section were carried out on a laptop Intel® Core™ i3-5005U CPU @ 2.00 GHz with a memory of 8.0 GB. The developed models were implemented in the IBM ILOG CPLEX 12.10 optimization studio.

4.1. The Data

For the arrival of elective surgery patients, it is assumed that the elective patients follow a discrete uniform distribution derived from historical data. The data collected from the case hospital consisted of 10,664 patients in operating rooms. Out of these patients, 95% require a ward after surgery and 5% require the ICU after surgery. For each type of patient, the underlying distribution is derived, and the number of patients for each type for the planning horizon is randomly generated using Equation (56) based on the relative frequencies of each type of patient in the historical data where N P e is the total number of surgery patients of specialty, e N P is the total number of surgery patients for all specialties in the considered planning horizon, and f e is the frequency of patients for each specialty in the historical data.
N P e = f e × N P   e E
Patient arrivals in the ICU and ward are assumed to follow a Poisson distribution with a known mean [1]. Patients’ length of stay in the operating room, ICU, and ward follows lognormal distribution [9,36]. The binomial distribution is used to find the number of patients visiting the ICU and ward after surgery. The probability that out of total surgery P patients who had surgery on a given day, the p patients in the ICU or ward can be calculated using the binomial distribution [57]. The formula for binomial distribution is given by Equation (57), where p is the probability that a patient will go to the ICU or ward after surgery. Using Equation (57), the number of patients from each specialty who will move to the ICU after surgery is determined.
Pr N P e = ( N P N P e ) p N P e q N P N P e
It is considered that 98% of the patients might follow one of the following paths after surgery. About 95% of patients are admitted to the ward before discharge and around 5% go to the ICU and ward before discharge. Only about 1% are discharged directly from the operating room [13,37]. After generating the number of patients arriving and length of stay data in the operating room, ICU, and ward, the elective surgery patient list is generated by assigning random arrivals and due dates to the patients within the planning horizon. Moreover, the preoperative time, operating room setup time, and cleaning time after the surgery of patients are known parameters and are generated randomly. The patients are inserted into the operating rooms, ICUs, and wards in the higher-level planning horizon. For experimentation using patients’ arrival information, it is assumed that 95% of patients operated in operating rooms are moved to the ward for post-operative recovery, and 5% of patients are admitted to ICU after surgery. Further, for experimentation, the efficiency of all the resources is considered to be 100% while performing the experiments. In practice, the efficiency of a resource can be limited due to the expertise level of deputed staff and the efficiency of the equipment in the operating room, ICU, and ward. Regarding the arrival of patients into the ICU and ward, it is considered that patient arrival follows a Poisson distribution. The arrival patterns considered for small, medium, and large problems are represented in Table 3 with the mean ( λ ) values for the ICU and ward.

4.2. Results

To validate our developed method, computational experiments are performed with different problem sizes and surgery durations. The arrival pattern for the ICU and ward is kept variable to generate scenarios of bottleneck resources among the operating room, ICU, and ward. As a first step, the patients are inserted into higher-level planning horizons based on their planned start time for the operating room and their arrival time for both the ICU and ward. After the higher-level of planning, the patients are allocated to different resources based on their requirements during their stay in the hospital. After the allocation of patients to all resources, based on the workload and utilization of the resources such as the operating room, ICU, and ward, the bottleneck resource with the highest utilization is identified. After identifying the bottleneck resource, the lower-level planning model for the bottleneck resource is run, and the schedule obtained is sent through a rope to be followed by the upstream departments. Table 4 represents the summary of results with average values for all samples based on the percentage utilization when the operating room, ICU, and ward are identified as a bottleneck.
The resource with the highest percentage utilization is the bottleneck resource. The percentage utilization is related to the idle time as given in Equation (58) where Idle% is the percentage idleness of the resource. As the percentage utilization increases, the idle time decreases, which in turn leads to cost savings.
I d l e % = 1 U T L %
When the operating room is the bottleneck resource, the operating room’s lower-level plan is operated, patients are allocated to operating rooms, and the optimal sequence of operating rooms is identified to minimize the makespan for all operating rooms. Table 5 represents the summary results of the lower-level plan when the operating room is the bottleneck resource.
When the optimal sequence of patients in operating rooms is obtained, this sequence is followed by the ICU and ward such that the arrival times of patients in the ICU and ward are scheduled after the discharge times of patients from the operating room as represented by Equations (46) and (47). Table 6 represents the summary results of the lower-level plan when the ICU is identified as the bottleneck resource. When ICU becomes the bottleneck resource, the lower-level plan for the ICU is made, and it is followed by both the operating room and ward such that for the operating room, the discharge time of patients from the operating room is planned before the arrival time of the patients in the ICU as explained in Equation (49).
In the same way, the patient’s arrival time in the ward is planned after the patient’s discharge from the ICU as given by Equation (50). Table 7 represents the summary results of the lower-level plan when the ward is the bottleneck resource.
When the ward becomes the bottleneck resource, the lower-level plan for the ward is made, and it is followed by both the operating room and ICU such that for the operating room, the discharge time of patients from the operating room is planned before the arrival time of the patients in the ward, as explained in Equation (52). In the same way, the patient’s discharge time from ICU is planned before the patient’s arrival time in the ward as given by Equation (53).

CPU Time

The average CPU time in seconds for each category of problem for one planning horizon is shown in Figure 6. The CPU time is proportional to the size of the problem in cases when the operating room and ward are the bottlenecks. While in the scenario where the ICU is the bottleneck, the trend in the CPU time for the medium-sized problem is random. It is proposed that the CPU time can be reduced using heuristic algorithms without compromising the solution’s optimality.

5. Conclusions, Limitations, and Future Research

The planning and scheduling of hospital resources are significant for better service and profit generation. The current research proposed an integrated planning and scheduling problem at different levels for operating rooms, intensive care units, and recovery wards. The theory of constraints was applied, and a mixed integer linear programming model was developed for the integrated planning and scheduling of the operating room, considering the constraints of the downstream ward. Different models for planning and scheduling were made considering the shifting of the bottleneck resources of hospitals. The integrated models for planning at different levels when any one resource such as the operating room, ICU, or ward becomes a bottleneck were proposed for the first time in the literature. The proposed model for integrated planning and scheduling was solved using the CPLEX solver on different sizes of test problems. The small, medium, and large-sized test problems were made for each scenario where the operating room, ICU, or ward became the bottlenecked resource.
The proposed model is significant for hospitals for effective and efficient planning and scheduling of their resources. This approach has practical significance in a hospital where the whole system’s performance is affected by one bottleneck which keeps changing in different planning horizons. Planning based on the capacity of the bottleneck resource subordinates the rest of the resources to the bottleneck resource to achieve optimal hospital performance. The theory of constraints is significant to identify the bottleneck resource among all the critical resources of the hospital. Therefore, the mathematical models developed using the theory of constraints in the current study can be used to develop hospital planning systems in which planning is based on the capacity of bottleneck resources.
The current study has its limitation in using the CPLEX solver, which provides an exact method to solve the mathematical models. In the future, constraint-based heuristics can be developed to reduce the CPU time of the different sizes of the considered problems, particularly for larger-sized instances. In addition, the current study uses deterministic surgery durations and lengths of stay in downstream units. Further extension of the work includes the use of data mining and machine learning methods to predict patient arrivals and length of stay in the operating room and downstream units and using the predicted data as the inputs in the models developed in the current study.

Author Contributions

Conceptualization and supervision, S.U.; methodology, A.T.; validation, A.T. and T.M.; formal analysis, A.T.; writing—original draft preparation, A.T.; writing—review and editing, A.T., T.M., Y.Y.G., B.L. and H.A.; funding acquisition, H.A. All authors have read and agreed to the published version of the manuscript.

Funding

Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2023R54), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

The study did not involve humans.

Data Availability Statement

The data used to support the findings of the study are available from the first author upon request.

Acknowledgments

The authors are grateful to the Department of Industrial Engineering, University of Engineering and Technology, Taxila.

Conflicts of Interest

The authors declare that they have no conflict of interest to report regarding the present study.

Notations

We describe the notations used in the current study as under:
C Index used to represent the operating room
O Set of operating rooms { o , o + 1 , , O }
i Index used to represent beds in the intensive care unit
I Set of beds in the intensive care unit { i , i + 1 , , I }
w Index used to represent beds in the wards
W Set of beds in wards { w , w + 1 , , W }
p Index used to represent a patient
P Set of patients { p , p + 1 , , P }
e Index used to represent a surgical specialty
E Set of surgical specialties { e , e + 1 , , E }
s Index to represent a surgeon
S Set of all surgeons { s , s + 1 , , S }
τ Planning horizon, a week
d Index used to represent the days in a planning horizon τ
D Set of days in the planning horizon τ
t p p r e o p Time taken by the patient p for preoperative tests and stay in the hospital etc.
t o p s e t Set up time of an operating room o for the surgery of patient p
t p s u r Surgery time of a patient p
t p p r e Preparation time of a patient p for surgery
t o p c l Cleaning time of operating room o after surgery of patient p
t p p S D S Sequence-dependent setup time when p P e , p P e
t s p a r r Arrival time of a surgeon s for surgery of patient p
A p Arrival time of patient p
D D p Due date of patient p
E P S T p Earliest possible start time of surgery for the patient p
t p i a r r Arrival time of patient p in ICU
t p i d i s c Discharge time of patient p from ICU
t p w a r r Arrival time of patient p in the ward
t p w d i s c Discharge time of patient p from the ward
M A big positive number
L O S p i Length of stay of the patient p in ICU on bed i
L O S p w Length of stay of the patient p in the Ward on bed w
W L o d τ Workload of the operating room o on the day d of the planning horizon τ
T W L o d τ The total workload of all operating rooms on the day d of the planning horizon τ
A t o d τ Available capacity of the operating room o on the day d of the planning horizon τ
N s d τ a l l The upper limit on the number of surgeries to be performed by a surgeon s on the day d of the planning horizon τ
Y p i d τ =1, if a patient p on the day d of planning horizon τ requires a bed i in ICU; 0 otherwise
Z p w d τ =1, if a patient p on the day d of planning horizon τ requires a bed w in the wards; 0 otherwise
Decision variables
H p d τ =1, if a patient p is assigned to a day d of planning horizon τ ; 0 otherwise
X p o d τ =1, if a patient p is assigned to an operating room o on the day d of planning horizon τ ; 0 otherwise
W p o d τ =1, if a patient p is operated on as first in the operating room o on the day d of planning horizon τ ; 0 otherwise
V p o d τ =1, if a patient p is operated on as last in the operating room o on the day d of planning horizon τ ; 0 otherwise
S p p o d τ =1 if a patient p’ is operated on immediately after patient p in the operating room o on the day d of planning horizon τ ; 0 otherwise
U s p o d τ =1 if the surgeon s S e is assigned to perform surgery on the patient p P e on the day d of planning horizon τ ; 0 otherwise

References

  1. Hillier, F.S.; Price, C.C. International Series in Operations Research & Management Science; Springer: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  2. Samudra, M.; Van Riet, C.; Demeulemeester, E.; Cardoen, B.; Vansteenkiste, N.; Rademakers, F.E. Scheduling operating rooms: Achievements, challenges and pitfalls. J. Sched. 2016, 19, 493–525. [Google Scholar] [CrossRef] [Green Version]
  3. Wang, L.; Demeulemeester, E.; Vansteenkiste, N.; Rademakers, F.E. Operating Room Planning and Scheduling for Outpatients and Inpatients: A Review and Future Research; KU Leuven, Faculteit Economie en Bedrijfswetenschappen: Leuven, Belgium, 2021. [Google Scholar] [CrossRef]
  4. Hans, E.W.; Vanberkel, P.T. Operating theatre planning and scheduling. In Handbook of Healthcare System Scheduling; Springer: Berlin/Heidelberg, Germany, 2012; pp. 105–130. [Google Scholar] [CrossRef]
  5. Akbarzadeh, B.; Moslehi, G.; Reisi-Nafchi, M.; Maenhout, B. The re-planning and scheduling of surgical cases in the operating room department after block release time with resource rescheduling. Eur. J. Oper. Res. 2019, 278, 596–614. [Google Scholar] [CrossRef]
  6. Bargetto, R.; Garaix, T.; Xie, X. A branch-and-price-and-cut algorithm for operating room scheduling under human resource constraints. Comput. Oper. Res. 2023, 152, 106136. [Google Scholar] [CrossRef]
  7. Aringhieri, R.; Duma, D.; Landa, P.; Mancini, S. Combining workload balance and patient priority maximisation in operating room planning through hierarchical multi-objective optimisation. Eur. J. Oper. Res. 2022, 298, 627–643. [Google Scholar] [CrossRef]
  8. McRae, S.; Brunner, J.O. Assessing the impact of uncertainty and the level of aggregation in case mix planning. Omega 2019, 97, 102086. [Google Scholar] [CrossRef]
  9. Freeman, N.; Zhao, M.; Melouk, S. An iterative approach for case mix planning under uncertainty. Omega 2018, 76, 160–173. [Google Scholar] [CrossRef]
  10. Yahia, Z.; Eltawil, A.B.; Harraz, N.A. The operating room case-mix problem under uncertainty and nurses capacity constraints. Health Care Manag. Sci. 2016, 19, 383–394. [Google Scholar] [CrossRef]
  11. Santos, D.; Marques, I. Designing master surgery schedules with downstream unit integration via stochastic programming. Eur. J. Oper. Res. 2021, 299, 834–852. [Google Scholar] [CrossRef]
  12. Spratt, B.; Kozan, E. Waiting list management through master surgical schedules: A case study. Oper. Res. Health Care 2016, 10, 49–64. [Google Scholar] [CrossRef]
  13. Fügener, A.; Hans, E.W.; Kolisch, R.; Kortbeek, N.; Vanberkel, P.T. Master surgery scheduling with consideration of multiple downstream units. Eur. J. Oper. Res. 2014, 239, 227–236. [Google Scholar] [CrossRef]
  14. Pham, D.-N.; Klinkert, A. Surgical case scheduling as a generalized job shop scheduling problem. Eur. J. Oper. Res. 2008, 185, 1011–1025. [Google Scholar] [CrossRef]
  15. Oliveira, M.; Bélanger, V.; Marques, I.; Ruiz, A. Assessing the impact of patient prioritization on operating room schedules. Oper. Res. Health Care 2020, 24, 100232. [Google Scholar] [CrossRef]
  16. Younespour, M.; Atighehchian, A.; Kianfar, K.; Esfahani, E.T. Using mixed integer programming and constraint programming for operating rooms scheduling with modified block strategy. Oper. Res. Health Care 2019, 23, 100220. [Google Scholar] [CrossRef]
  17. Ma, G.; Demeulemeester, E. A multilevel integrative approach to hospital case mix and capacity planning. Comput. Oper. Res. 2013, 40, 2198–2207. [Google Scholar] [CrossRef] [Green Version]
  18. Fügener, A. An integrated strategic and tactical master surgery scheduling approach with stochastic resource demand. J. Bus. Logist. 2015, 36, 374–387. [Google Scholar] [CrossRef]
  19. Guido, R.; Conforti, D. A hybrid genetic approach for solving an integrated multi-objective operating room planning and scheduling problem. Comput. Oper. Res. 2017, 87, 270–282. [Google Scholar] [CrossRef]
  20. Beliën, J.; Demeulemeester, E.; Cardoen, B. Visualizing the demand for various resources as a function of the master surgery schedule: A case study. J. Med. Syst. 2006, 30, 343–350. [Google Scholar] [CrossRef] [Green Version]
  21. Macario, A.; Vitez, T.; Dunn, B.; McDonald, T. Where are the costs in perioperative care?: Analysis of hospital costs and charges for inpatient surgical care. Anesthesiol. J. Am. Soc. Anesthesiol. 1995, 83, 1138–1144. [Google Scholar] [CrossRef]
  22. Cardoen, B.; Demeulemeester, E.; Beliën, J. Operating room planning and scheduling: A literature review. Eur. J. Oper. Res. 2010, 201, 921–932. [Google Scholar] [CrossRef] [Green Version]
  23. Guerriero, F.; Guido, R. Operational research in the management of the operating theatre: A survey. Health Care Manag. Sci. 2011, 14, 89–114. [Google Scholar] [CrossRef]
  24. Roshanaei, V.; Booth, K.E.; Aleman, D.M.; Urbach, D.R.; Beck, J.C. Branch-and-check methods for multi-level operating room planning and scheduling. Int. J. Prod. Econ. 2019, 220, 107433. [Google Scholar] [CrossRef]
  25. Schneider, A.T.; van Essen, J.T.; Carlier, M.; Hans, E.W. Scheduling surgery groups considering multiple downstream resources. Eur. J. Oper. Res. 2019, 282, 741–752. [Google Scholar] [CrossRef]
  26. Kamran, M.A.; Karimi, B.; Dellaert, N.; Demeulemeester, E. Adaptive operating rooms planning and scheduling: A rolling horizon approach. Oper. Res. Health Care 2019, 22, 100200. [Google Scholar] [CrossRef]
  27. Vancroonenburg, W.; De Causmaecker, P.; Berghe, G.V. Chance-constrained admission scheduling of elective surgical patients in a dynamic, uncertain setting. Oper. Res. Health Care 2019, 22, 100196. [Google Scholar] [CrossRef]
  28. Burdett, R.L.; Kozan, E. An integrated approach for scheduling health care activities in a hospital. Eur. J. Oper. Res. 2018, 264, 756–773. [Google Scholar] [CrossRef] [Green Version]
  29. Kroer, L.R.; Foverskov, K.; Vilhelmsen, C.; Hansen, A.S.; Larsen, J. Planning and scheduling operating rooms for elective and emergency surgeries with uncertain duration. Oper. Res. Health Care 2018, 19, 107–119. [Google Scholar] [CrossRef] [Green Version]
  30. Roshanaei, V.; Naderi, B. Solving integrated operating room planning and scheduling: Logic-based Benders decomposition versus Branch-Price-and-Cut. Eur. J. Oper. Res. 2020, 293, 65–78. [Google Scholar] [CrossRef]
  31. Tayyab, A.; Saif, U. A Two-Stage Genetic Artificial Bee Colony Algorithm for Solving Integrated Operating Room Planning and Scheduling Problem With Capacity Constraints of Downstream Wards. IEEE Access 2022, 10, 131109–131127. [Google Scholar] [CrossRef]
  32. Breuer, D.J.; Lahrichi, N.; Clark, D.E.; Benneyan, J.C. Robust combined operating room planning and personnel scheduling under uncertainty. Oper. Res. Health Care 2020, 27, 100276. [Google Scholar] [CrossRef]
  33. Adan, I.; Bekkers, J.; Dellaert, N.; Vissers, J.; Yu, X. Patient mix optimisation and stochastic resource requirements: A case study in cardiothoracic surgery planning. Health Care Manag. Sci. 2009, 12, 129–141. [Google Scholar] [CrossRef]
  34. Kheiri, A.; Lewis, R.; Thompson, J.; Harper, P. Constructing operating theatre schedules using partitioned graph colouring techniques. Health Syst. 2021, 10, 286–297. [Google Scholar] [CrossRef]
  35. Smith-Daniels, V.L.; Schweikhart, S.B.; Smith-Daniels, D.E. Capacity management in health care services: Review and future research directions. Decis. Sci. 1988, 19, 889–919. [Google Scholar] [CrossRef]
  36. Jebali, A.; Diabat, A. A stochastic model for operating room planning under capacity constraints. Int. J. Prod. Res. 2015, 53, 7252–7270. [Google Scholar] [CrossRef]
  37. Moosavi, A.; Ebrahimnejad, S. Robust operating room planning considering upstream and downstream units: A new two-stage heuristic algorithm. Comput. Ind. Eng. 2020, 143, 106387. [Google Scholar] [CrossRef]
  38. Testi, A.; Tanfani, E.; Torre, G. A three-phase approach for operating theatre schedules. Health Care Manag. Sci. 2007, 10, 163–172. [Google Scholar] [CrossRef]
  39. Chow, V.S.; Puterman, M.L.; Salehirad, N.; Huang, W.; Atkins, D. Reducing surgical ward congestion through improved surgical scheduling and uncapacitated simulation. Prod. Oper. Manag. 2011, 20, 418–430. [Google Scholar] [CrossRef]
  40. Goldratt, E.M. Essays on the Theory of Constraints; North River Press: Great Barrington, MA, USA, 1990. [Google Scholar]
  41. Goldratt, E.M.; Cox, J. The Goal: A Process of Ongoing Improvement; Routledge: London, UK, 2016. [Google Scholar] [CrossRef]
  42. Georgiadis, P.; Politou, A. Dynamic Drum-Buffer-Rope approach for production planning and control in capacitated flow-shop manufacturing systems. Comput. Ind. Eng. 2013, 65, 689–703. [Google Scholar] [CrossRef]
  43. Gilland, W.G. A simulation study comparing performance of CONWIP and bottleneck-based release rules. Prod. Plan. Control 2002, 13, 211–219. [Google Scholar] [CrossRef]
  44. Saif, U.; Guan, Z.; Wang, C.; He, C.; Yue, L.; Mirza, J. Drum buffer rope-based heuristic for multi-level rolling horizon planning in mixed model production. Int. J. Prod. Res. 2019, 57, 3864–3891. [Google Scholar] [CrossRef]
  45. Sirikrai, V.; Yenradee, P. Modified drum–buffer–rope scheduling mechanism for a non-identical parallel machine flow shop with processing-time variation. Int. J. Prod. Res. 2006, 44, 3509–3531. [Google Scholar] [CrossRef]
  46. Ronen, B.; Gur, R.; Pass, S. Focused management in military organizations: An avenue for future industrial engineering. Comput. Ind. Eng. 1994, 27, 543–544. [Google Scholar] [CrossRef]
  47. Pegels, C.C.; Watrous, C. Application of the theory of constraints to a bottleneck operation in a manufacturing plant. J. Manuf. Technol. Manag. 2005, 16, 302–311. [Google Scholar] [CrossRef]
  48. Huang, W.-T.; Chen, P.-S.; Liu, J.J.; Chen, Y.-R.; Chen, Y.-H. Dynamic configuration scheduling problem for stochastic medical resources. J. Biomed. Inform. 2018, 80, 96–105. [Google Scholar] [CrossRef]
  49. Hachicha, H.K.; Mansour, F.Z. Two-MILP models for scheduling elective surgeries within a private healthcare facility. Health Care Manag. Sci. 2018, 21, 376–392. [Google Scholar] [CrossRef]
  50. Anjomshoa, H.; Dumitrescu, I.; Lustig, I.; Smith, O.J. An exact approach for tactical planning and patient selection for elective surgeries. Eur. J. Oper. Res. 2018, 268, 728–739. [Google Scholar] [CrossRef]
  51. Nemati, S.; Shylo, O.V.; Prokopyev, O.A.; Schaefer, A.J. The surgical patient routing problem: A central planner approach. Inf. J. Comput. 2016, 28, 657–673. [Google Scholar] [CrossRef] [Green Version]
  52. Marques, I.; Captivo, M.E.; Pato, M.V. A bicriteria heuristic for an elective surgery scheduling problem. Health Care Manag. Sci. 2015, 18, 251–266. [Google Scholar] [CrossRef]
  53. Liu, N.; Truong, V.A.; Wang, X.; Anderson, B.R. Integrated Scheduling and Capacity Planning with Considerations for Patients’ Length-of-Stays. Prod. Oper. Manag. 2019, 28, 1735–1756. [Google Scholar] [CrossRef]
  54. Siqueira, C.L.; Arruda, E.F.; Bahiense, L.; Bahr, G.L.; Motta, G.R. Long-term integrated surgery room optimization and recovery ward planning, with a case study in the Brazilian National Institute of Traumatology and Orthopedics (INTO). Eur. J. Oper. Res. 2018, 264, 870–883. [Google Scholar] [CrossRef]
  55. Barz, C.; Rajaram, K. Elective patient admission and scheduling under multiple resource constraints. Prod. Oper. Manag. 2015, 24, 1907–1930. [Google Scholar] [CrossRef] [Green Version]
  56. González, P.; Herrero, C. Optimal sharing of surgical costs in the presence of queues. Math. Methods Oper. Res. 2004, 59, 435–446. [Google Scholar] [CrossRef] [Green Version]
  57. Vanberkel, P.T.; Boucherie, R.J.; Hans, E.W.; Hurink, J.L.; Van Lent, W.A.; Van Harten, W.H. An exact approach for relating recovering surgical patient workload to the master surgical schedule. J. Oper. Res. Soc. 2011, 62, 1851–1860. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Patient flow considered in this paper.
Figure 1. Patient flow considered in this paper.
Applsci 13 03616 g001
Figure 2. Multi-level planning considered in the current paper.
Figure 2. Multi-level planning considered in the current paper.
Applsci 13 03616 g002
Figure 3. Multi-level integrated planning when the operating room is the bottleneck.
Figure 3. Multi-level integrated planning when the operating room is the bottleneck.
Applsci 13 03616 g003
Figure 4. Multi-level integrated planning when the ICU is the bottleneck.
Figure 4. Multi-level integrated planning when the ICU is the bottleneck.
Applsci 13 03616 g004
Figure 5. Multi-level integrated planning when the ward is the bottleneck.
Figure 5. Multi-level integrated planning when the ward is the bottleneck.
Applsci 13 03616 g005
Figure 6. CPU time for each category of problem: (a) operating room is the bottleneck, (b) ICU is the bottleneck, and (c) ward is the bottleneck.
Figure 6. CPU time for each category of problem: (a) operating room is the bottleneck, (b) ICU is the bottleneck, and (c) ward is the bottleneck.
Applsci 13 03616 g006
Table 1. Characteristics of the different-sized problems.
Table 1. Characteristics of the different-sized problems.
Problem SizeNo. of ORNo. of Beds in ICUNo. of Beds in Ward
Small3–47–1025–35
Medium6–811–1545–55
Large10–1216–1880–100
Table 2. Matrix of experiments.
Table 2. Matrix of experiments.
Sr. No.Problem SizeSurgery DurationBottleneck ResourceNumber of Instances
1SmallSmallOR10
2SmallMediumOR10
3SmallLargeOR10
4SmallSmallICU10
5SmallMediumICU10
6SmallLargeICU10
7SmallSmallWard10
8SmallMediumWard10
9SmallLargeWard10
10MediumSmallOR10
11MediumMediumOR10
12MediumLargeOR10
13MediumSmallICU10
14MediumMediumICU10
15MediumLargeICU10
16MediumSmallWard10
17MediumMediumWard10
18MediumLargeWard10
19LargeSmallOR10
20LargeMediumOR10
21LargeLargeOR10
22LargeSmallICU10
23LargeMediumICU10
24LargeLargeICU10
25LargeSmallWard10
26LargeMediumWard10
27LargeLargeWard10
Table 3. Value of λ for different problem sizes.
Table 3. Value of λ for different problem sizes.
Problem SizePoisson λ for ICU Poisson λ for ward
Small2–35–6
Medium5–710–12
Large10–1215–17
Table 4. Summary of results.
Table 4. Summary of results.
Problem SizeSurgery DurationNo. of Patients in the ORWorkload of ORUtilization of ORNo. of Patients in the ICUUtilization of ICUNo. of Patients in the WardUtilization of Ward
The operating room becomes a Bottleneck
SmallSmall16169593%560%2272%
SmallMedium11145284%559%1066%
Small Large9141376%558%1861%
MediumSmall33354292%869%4182%
MediumMedium22301783%868%3673%
MediumLarge21319084%972%3671%
LargeSmall43490587%1171%6168%
LargeMedium35506788%1272%5966%
LargeLarge28518890%1170%6775%
ICU becomes Bottleneck
SmallSmall12139779%787%2171%
SmallMedium9122471%790%1963%
SmallLarge9130274%788%2169%
MediumSmall26293376%1192%3775%
MediumMedium20288276%1191%3468%
MediumLarge17315582%1191%3468%
LargeSmall33384368%1591%6876%
LargeMedium28394771%1595%6370%
LargeLarge22402370%1591%6370%
Ward becomes Bottleneck
SmallSmall14155085%670%2894%
SmallMedium10136071%559%2892%
SmallLarge7131368%568%2791%
MediumSmall23261771%866%4590%
MediumMedium20290176%971%4387%
MediumLarge16287775%870%4386%
LargeSmall34391973%1276%8291%
LargeMedium28401774%1172%8291%
LargeLarge23415174%1277%8594%
Table 5. Summary results of the lower-level plan when the operating room becomes the bottleneck.
Table 5. Summary results of the lower-level plan when the operating room becomes the bottleneck.
PSSDOperating RoomICUWard
PPMakespanPOROB EPPORPICUOBEP
SmallSmall16134712215043
SmallMedium11124212210154
Small Large913860228154
Medium Small332858215310010
Medium Medium222624144211411
Medium Large212620135201411
LargeSmall433930218412216
LargeMedium3547102193311114
LargeLarge2845211282722315
PS = Problem Size; SD = Surgery Duration; PP = Patients Planned; POR = Patient from OR; OB = Occupied Beds; EP = Expected Patients; and PICU = Patients from ICU.
Table 6. Summary results of the lower-level plan when the ICU becomes the bottleneck.
Table 6. Summary results of the lower-level plan when the ICU becomes the bottleneck.
PSSDOperating RoomICUWard
PPMakespanPOROB EPPORPICUOBEP
SmallSmall12106313413125
SmallMedium9111905210144
Small Large9122105210146
Medium Small26254215525247
Medium Medium202572155192211
Medium Large172923146161512
LargeSmall3333072583221915
LargeMedium2836211682632113
LargeLarge2236851592142414
PS = Problem Size; SD = Surgery Duration; PP = Patients Planned; POR = Patient from OR; OB = Occupied Beds; EP = Expected Patients; and PICU = Patients from ICU.
Table 7. Summary results of the lower-level plan when the ward becomes the bottleneck.
Table 7. Summary results of the lower-level plan when the ward becomes the bottleneck.
PSSDOperating RoomICUWard
PPMakespanPOROBEPPORPICUOBEP
SmallSmall14122013213186
SmallMedium10121203292125
SmallLarge7123803271154
MediumSmall2322731252221111
MediumMedium202557135192139
MediumLarge1626301351521610
LargeSmall3431562193223513
LargeMedium2836911192753615
LargeLarge23352411102164315
PS = Problem Size; SD = Surgery Duration; PP = Patients Planned; POR = Patient from OR; OB = Occupied Beds; EP = Expected Patients; and PICU = Patients from ICU.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Tayyab, A.; Ullah, S.; Mahmood, T.; Ghadi, Y.Y.; Latif, B.; Aljuaid, H. Modeling of Multi-Level Planning of Shifting Bottleneck Resources Integrated with Downstream Wards in a Hospital. Appl. Sci. 2023, 13, 3616. https://doi.org/10.3390/app13063616

AMA Style

Tayyab A, Ullah S, Mahmood T, Ghadi YY, Latif B, Aljuaid H. Modeling of Multi-Level Planning of Shifting Bottleneck Resources Integrated with Downstream Wards in a Hospital. Applied Sciences. 2023; 13(6):3616. https://doi.org/10.3390/app13063616

Chicago/Turabian Style

Tayyab, Aisha, Saif Ullah, Toqeer Mahmood, Yazeed Yasin Ghadi, Bushra Latif, and Hanan Aljuaid. 2023. "Modeling of Multi-Level Planning of Shifting Bottleneck Resources Integrated with Downstream Wards in a Hospital" Applied Sciences 13, no. 6: 3616. https://doi.org/10.3390/app13063616

APA Style

Tayyab, A., Ullah, S., Mahmood, T., Ghadi, Y. Y., Latif, B., & Aljuaid, H. (2023). Modeling of Multi-Level Planning of Shifting Bottleneck Resources Integrated with Downstream Wards in a Hospital. Applied Sciences, 13(6), 3616. https://doi.org/10.3390/app13063616

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

Article Metrics

Back to TopTop