Abstract

Operating room scheduling is an important operational problem in most hospitals. In this paper, a novel mixed integer programming (MIP) model is presented for minimizing Cmax and operating room idle times in hospitals. Using this model, we can determine the allocation of resources including operating rooms, surgeons, and assistant surgeons to surgeries, moreover the sequence of surgeries within operating rooms and the start time of them. The main features of the model will include the chronologic curriculum plan for training residents and the real-life constraints to be observed in teaching hospitals. The proposed model is evaluated against some real-life problems, by comparing the schedule obtained from the model and the one currently developed by the hospital staff. Numerical results indicate the efficiency of the proposed model compared to the real-life hospital scheduling, and the gap evaluations for the instances show that the results are generally satisfactory.

1. Introduction

Health care expenditures comprise a meaningful portion of the Gross Domestic Product in both developed and developing countries. Expenditure on healthcare in the UK as a percentage of the UK Gross Domestic Product (GDP) was estimated to be 8.4% in 2007, from which the public share was 69% [1]. Also, according to the statics released by the WHO (World Health Organization), health care expenditures in 2007 in Iran as a developing country were estimated to be about 6.4% of its GDP, and the portion covered by the government was about 46.8%. This fact makes health systems an important research field for industrial engineering and operations research to improve their operational efficiency.

Operating rooms are simultaneously the largest cost centers and the greatest source of revenues for most hospitals. OR planning and scheduling is a key tool which can be useful to improve the productivity level of ORs and the related departments.

Basically, there are three OR scheduling strategies commonly employed:(1)blocked scheduling strategy,(2)open scheduling strategy, and(3)modified scheduling strategy.

Under a blocked scheduling strategy, individual surgeons or surgical groups are assigned times in a particular OR in a periodic (typically weekly or monthly) schedule. The planning within the framework of a blocked strategy consists of three stages. In the first stage, the OR capacity is divided among the surgeons, surgical groups, or departments on a strategic level. Then, a cyclic timetable called “Master Surgical Schedule” is constructed that defines the number and type of operating rooms available, the hours that ORs will be open, and the surgical groups or surgeon sessions for each OR [2]. The last stage which may be called “surgery process scheduling” splits into two subproblems called “advance scheduling” and “allocation scheduling” [3]. The first subproblem at a tactical level (one weak to one month) solves a planning phase by assigning an operating date to each patient over the time horizon. The second subproblem solves a scheduling phase which determines the sequence and resource allocation of cases in a given day [4, 5].

In the open strategy, the hospital does not hold operation rooms specific to a single surgeon although there is sometimes a weekly schedule for each surgeon. In this strategy, the intention is to accommodate all patients. The surgeons submit cases up until the day of the surgery, and all the cases are scheduled in ORs.

The modified strategy is similar to the blocked one except that certain slots in the master surgical schedule are left open for flexibility. In fact, this strategy is a mix of open and blocked strategies.

In this paper, we investigate the deterministic daily scheduling of surgical cases in teaching hospitals using the open scheduling strategy. A novel MIP model is presented with the objective of minimizing the total idle times of operating rooms and max which represent completion time of the last patient’s surgery. Using this model, we endeavor to allocate the required resources including operating rooms, surgeons, and assistant surgeons to surgeries and also determine the sequence of surgeries within operating rooms and the start time of each surgery. The major features of the model include the chronologic curriculum plan for training residents and the real-life constraints in teaching hospitals. We maintain that the problem defined here with its specific characteristics is not reported in the literature and that it finds real applications in teaching hospitals.

The rest of the paper is organized as follows. A brief review of previous work will first be presented in Section 2. Section 3 will provide a definition of the problem. The MIP model will be presented in Section 4. Section 5 is devoted to solving approach. Numerical results will be reported in Section 6. Conclusions will be drawn and presented in Section 7.

2. A Literature Review

Over the past 60 years, a large body of literature has evolved on the management of operating theaters. Magerlein and Martin [3] reviewed the literature on surgical demand scheduling and, as previously mentioned, focused on “surgery process scheduling” to differentiate between advance and allocation scheduling. Blake and Carter [6] extended this classification and added the domain of external resource scheduling which is defined as the process of identifying and reserving all resources external to the surgical suite necessary to ensure appropriate care for a patient before and after an instance of surgery. They further divided each domain into strategic, administrative, and operational levels [7]. Many authors (Magerlein and Martin, 1987; Kennedy [8]) classified the literature based on solution techniques [4], but Cardoen et al. [7] proposed a detailed classification based on 7 areas that are related to either the problem setting (e.g., performance measures or patient classes) or the technical features (e.g., solution technique or uncertainty incorporation) and analyzed the contributions in these areas.

Deterministic and stochastic mathematical programming models, queuing models, simulation models, and heuristic approaches have all been widely used to investigate OR scheduling. In this paper, we will briefly focus on the studies that have a direct bearing on surgical process scheduling and will consider the deterministic programming models used for OR scheduling. A more comprehensive review can be found in [7].

Ogulata and Erol [9] modeled the hierarchical multiple criteria mathematical programming to assign surgeon groups and dates to patients. A goal programming has been proposed to solve this weekly programming. Vissers et al. [10] presented an MIP model for developing a master operating theater schedule for thoracic surgery. Chaabane et al. [2] introduced two planning methods to compare open and blocked strategies. A weekly surgical case planning has been introduced to assign OR, operation day, and start time to patients in the open strategy, and an MSS has been defined for the blocked strategy. Jebali et al. [11] split a surgical case scheduling into 2 stages. First, the surgical cases were assigned to OR and then sequenced based on two sequencing strategies. Testi et al. [12] used the blocked strategy to develop a three-phase hierarchical approach for the weekly planning of surgical cases. They also applied a simulation method to analyze different rules of surgical case sequencing such as LWT, LPT, and SPT. Fei et al. [13] focused on the planning stage and assigned a set of surgical cases to OR with the objective of minimizing total operating cost. They proposed a branch and price algorithm to solve this problem for the exact solution. This model is similar to a parallel machine scheduling problem. Pham and Klinkert [4] Proposed a novel extension of the job shop scheduling problem called multimode blocking job shop to construct a surgical case scheduling for elective and add-on cases. Riis and Burke [14] use a discrete time model for surgery scheduling problem. Their model is similar to the one-step job shop scheduling problem with multiple identical machines, and they present a metaheuristic local-search-based algorithm for solving it. Cardoen et al. [15] investigated the sequencing of outpatients in an ASC (Ambulatory Surgery Center) in which balancing the use of recovery beds and medical precaution constraints such as additional sterilizing of OR after surgery of infected patients and sterilizing of medical instruments were taken into account. A branch and price methodology was used to solve this NP-hard optimization problem. Fei et al. [16] employed the open strategy to design a weekly surgery schedule. They solved the problem in two phases to give the surgery date for each patient and to determine the sequence of operations in each operating room on each day. Roland et al. [5] introduced a RCPSP (Resource Constraint Project Scheduling Problem) to generate scheduling and planning or surgical cases while focusing as much attention on human resources as on economic factors. They proposed a heuristic solution procedure based on genetic algorithm for solving this problem.

Hanset et al. [17] studied a daily scheduling problem, including material and human constraints. Their problem is modeled and solved via a constraint programming approach. Hanset et al. [18] compared their constraint programming model, [17], to a mixed integer programming model previously published by Roland et al. [5] that was adapted for the need of their study.

In almost all the studies cited above, the surgeon for each case was known in advance. In this paper, the daily scheduling of surgical cases in teaching hospitals is investigated in which the training operations must be performed by postgraduate alumni, called “residents” or “fellows,” under the supervision of a relevant attending surgeon while they are also qualified to do certain operations independently.

Chronologic curriculum plan for training residents is addressed in resident scheduling problem (RSP) repeatedly, which is different from operating room scheduling problem. This problem involves assigning residents to day and night shifts over a given planning horizon subject to numerous working regulations and staffing requirements [19]. For example, Topaloglu and Ozkarahan [19] developed a mixed-integer programming model for scheduling residents’ duty hours considering the on-call night, dayoff, rest period, and total work-hour ACGME (Accreditation Council for Graduate Medical Education) regulations as well as the demand coverage requirements of the residency program. Sherali et al. [20] addressed the resident scheduling problem (RSP) at hospitals concerned with prescribed work nights for residents while considered departmental staffing and skill requirements as well as residents’ preferences.

Dexter et al. [21] developed a methodology to determine rotations consisting of combinations of specialties to be paired for purposes of trainee scheduling to reduce the incidence of daily assignments off rotation.

To the best of our knowledge among the body of work on OR planning and scheduling that has appeared in the literature, no report is available in the literature on operating room scheduling in teaching hospitals that takes into account real constraints such as the chronologic curriculum plan for training residents. Constraints related with allocation of operations based on different seniority levels of residents and fellowships and providing a good balance between total hours of operations that performed by them than the average in a month. Also most of these papers assume that surgeon of each patient is already known. So it seems, assigning surgeons, residents, and assistant surgeons to each instance of surgery in literature of OR planning and scheduling is new.

3. Problem Definition

The process of surgical scheduling consists of two steps. The first concerns a weekly plan to assign a specific date to each patient waiting for surgery while the second involves sequencing and scheduling surgical cases on a given date. In this paper, we consider a surgical case scheduling problem (the second step) in a teaching hospital. In such hospitals, a list is prepared of the patients waiting for operation on the following day. Both inpatients and outpatients are scheduled each day by the OR head nurse or manager who sequences the surgeries and assigns a surgeon, an operating room, and a start time to each surgery.

Patients are prioritized and sequenced by the head nurse based on both resource availability (the eligible surgeon, the required equipments and aids, and the appropriate OR) and patient (child or old) priority.

The training operations must be performed by postgraduate alumni (“residents” and “fellows”) under the supervision of relevant attending surgeons. The residents and fellows are usually grouped according to their seniority levels, and each of them is qualified to do certain operations based on their experience and qualifications acquired in the course of their curriculum plan. This plan determines the type of operations that can be performed by each resident or fellow group during each period of their education.

Certain operations also require an assistant surgeon. The distribution of surgeries among surgeon groups, residents, and fellows should be based on equal opportunity for all to acquire experience. The nontraining operations are performed by public surgeons.

Some surgeries require special equipments which are available in particular ORs. For each operation, a group of medical staff members collectively called the surgery aid group in this paper, consisting of a scrub nurse, a circular nurse, and an anesthetic technician, is required. Also several medical instruments are needed during the surgery. After each surgery, instruments possibly need to be sterilized for some periods and, hence, are not available for subsequent surgeries.

We should also deal with the likelihood of infection spread. Therefore, special sanitary procedures must be executed to avoid the transfer of infection from patient to patient. In particular, the operating room needs additional cleaning after operation on an infected patient.

In the next section, the objectives and constraints of the problem are identified in a mathematical formulation.

4. Mathematical Formulation of the Problem

This section describes an MIP model for determining the daily surgical scheduling of elective patients (inpatient and outpatient) that calls “HORS” (hospital operating room scheduling).

4.1. Notations

The following notations are used in this paper.

(i) Indexes
i: Index for surgery.r: Index for operating room.s, s′: Index for surgeon and assistant surgeon.p: Index for time slot (30-minute units).k: Index for medical instrument type.l: Index for group number.

(ii) Inputs
AtN(i, s): If operation i should be performed under the supervision of an attending surgeon s, it take 1; otherwise, 0.d(i): Duration of operation i. FAS(i, s): If surgeon s can be the assistant for operation i, it takes 1; otherwise, 0.FR(i, r): If surgery i can be performed in operating room r, it takes 1; otherwise, 0.FS(i, s): If surgery i can be performed by surgeon s it takes 1; otherwise, 0.Infect(i): If the patient of surgery i is infected, it takes 1; otherwise, 0.Ins(i, k): If surgery i needs instrument k, it takes 1; otherwise, 0.IP0(k): The initial inventory of instrument sets of type k. Nat(i): The number of required attending surgeons for the supervision of operation i. IsNA(i): If operation i needs an assistant surgeon, it takes 1; otherwise, 0.GN(s, l): If surgeon s to be of group number l, it takes 1; otherwise, 0.priority(i): If patient i be a priority patient, it takes 1; otherwise, 0.SAV(s, p): If surgeon s is available on block p, it takes 1; otherwise, 0.SGC(s): Surgeon group type: for residents it takes 1, for fellows it takes 2, and for attending surgeons it takes 3.TecNo(p): Number of technicians available on block p.TS(k): The duration of sterilization of instrument k. Kclean: The duration of additional cleaning of OR after the surgery of an infected patient.mean: Average total operation time assigned to group number l in one month.dev: Allowable deviation from mean .C0: Last time allowed for the operating room doing surgery.P0: Maximum time allowed to start of priority surgery.

(iii) Decision Variables
Xipsr: If surgery i starts in period p by surgeon s in room r, it takes 1; otherwise, 0.Yipsr: If surgery i occupies period p in room r by surgeon s, it takes 1; otherwise, 0.Aisp: If surgery i is assigned to assistant surgeon s in period p, it takes 1; otherwise, 0.APisp: If surgery i occupies period p by assistant surgeon s, it takes 1; otherwise, 0.De(i,k p): If instrument of type k is used in period p by surgery i, it takes 1; otherwise, 0.IP(k, p): The inventory level of instrument k in period p. Omax(r): The finish time of the last surgery of each operating room.Teta(r, p): The occupation status of room r in block p. If room r is occupied, it takes 1; otherwise, 0.Two(r): Idle time for each operating room.

4.2. The Mathematical Model

This model aims at minimizing operating room idle time and max, where max is, as mentioned earlier, the completion time of the latest patient’s surgery in operating room. The most important objective of decision maker is minimizing the operating room idle time. But considering this objective alone in the model makes some operating rooms remain empty while there may be some operating rooms which are occupied up to a maximum possible. Therefore, minimizing max is considered as the primary objective, and minimizing operating room idle time is the second level objective. Here, is defined max and is the total idle time of operating rooms. Problem-solving process will be explained in a next section.

The MIP model is then formulated as follows: The constraints are where

Equations (4.1) and (4.2) define (max) and (total idle time of operating rooms), respectively. Equation (4.3) states that each surgery should be assigned to just one operating room, one surgeon, and one time block as a start block. The variable is defined in (4.4). At most, one surgery can be operated in each operating room in each time block (see (4.5)). Equation (4.6) states that the number of simultaneous operations performed by each surgeon cannot exceed one. Equation (4.7) specifies that the operation should not be assigned to the surgeon that is not qualified for it. Constraint (4.8) prohibits assigning a surgical case to a surgeon who is not available in that period. Constraints (4.9) to (4.13) relate to assigning assistant surgeons to operations. Constraints (4.14) and (4.15) are related to the availability of attending surgeons during the operations. Equation (4.16) defines the balanced distribution of operations between several surgeon groups. Equation (4.17) forbids an operation in a room where the required equipment is not available. Constraint (4.18) ensures that the required technicians for all the surgical cases performed simultaneously cannot exceed the total number of technicians available in that period. Constraints (4.19) to (4.21) state that the number of required instruments cannot exceed the available capacity. We consider medical precaution constraints for infected patients also considered by [16]. Constraint (4.22) states that, when the surgery of an infected patient is scheduled for a specific period in an OR, the noninfected patient cannot be assigned to that OR until additional cleaning are performed. Constraint (4.23) determines the finish time of the last surgery of each operating room. The occupation of rth room in pth time block is defined by the variable Teta in (4.24). Constraint (4.25) is used to define the total idle time (TWO) for each operating room. Equation (4.26) states that no surgery can start too late but must be planned to be completed within the operation room opening hours. Equations (4.27) and (4.28) ensure that end time of last operation should not exceed form allowable time and, finally, (4.29) defines priority constraint.

5. Solving Approach

As previously mentioned, this model aims at minimizing operating room idle time and max. Initially the problem was solved using the Lexicograph method. Due to the high time to solve, some changes were made in solving method. Below, and are max and idle time of operating rooms, respectively, and shows constraints of the model. Solving process is described below.

Step 1. Choose an appropriate upper bound for .

Step 2. Solve the following problem:

Step 3. If the problem is feasible, deduce the upper bound and go to Step 2.

Step 4. If the problem has no feasible solution, then exit with the current feasible solution as the optimum solution.

Flow diagram of this method is shown in Figure 1.

6. Computational Experiments

We did not come across any problem as defined in this paper. Although there are different works cited in the literature, the differences exist in the constraints considered. In this paper, the real constraints of teaching hospitals have been duly taken into account.

To evaluate the proposed mathematical model, several real problems were collected from the Hasheminejad Kidney Center (HKC), an academic hospital, in Iran. The collected data consisted of lists of patients for each day, start and end times of surgeries, details of the patients such as infection, the surgeon assigned, and the operating room assigned. To collect some required data that is not registered in the HIS (hospital information system), we have designed an ICR (intelligent character recognition) form that was completed by operating room personnel for 6 months and selected the most complete forms that were filled out within a one-month period. The designed form is shown in the appendix section.

There are four different surgical services at HKC: urologic, endoscopic, laparoscopic, and vascular. The operating theatre is composed of 7 operating rooms. All the operating rooms are multifunctional, but endoscopic operations can be performed in only three. There are four resident groups of urology according to their seniority levels and two fellowship groups. Each resident group consists of three residents with the same skill, and each fellowship group consists of two fellows. There is a trainee rotation that determines the days each resident or fellow should serve in the operating theatre.

The number of technician groups varies on different work shifts. There are 35 different types of important instruments used in the operating theatre (e.g., laparoscopic trolley, C-Arm, etc.), each case requiring a specific set. Furthermore, surgery durations, required instruments for each surgery, instrument inventory, and the corresponding sterilization times are known. Also the duration of setup time is non-sequence-dependent and is considered in the processing time of the operation.

Time data are expressed in 30-minute units, and the scheduling horizon comprises 25 time units.

The model was implemented in GAMS and solved using CPLEX 12.1.0. The allowed computation time for the solution procedure was limited to 400 seconds per instance on a 2, 66 GHz Pentium 4 PC with the Windows 7 operating system.

The model was evaluated by solving several real instances at HKC, and the gap from optimal solution was reported.

The evaluation was based on 30 real test problems with 18 to 45 patients and 8 to 13 surgeons. To study the performance of the model, the MIP gap was evaluated. Table 1 shows the resulting MIP gap for each test instance defined as the absolute gap between the values obtained for the objective function and the lower bound given by CPLEX, the number of discrete variables, constraints, patients, and surgeons.

In addition, to analyze the surgical case scheduling resulting from the proposed model, the final solution for one month of chosen instances was explored and the actual scheduling of hospital compared with the proposed scheduling.

Since number and type of surgery varies in different days of a week, 7 patterns to choose samples in different days were considered. Training operations are not performed on Thursday and Friday. Also transplant operations with their specific constraints are performed on Sunday and Tuesday. Thereby the problem instances are not chosen from these days of the week (because in this model, focus is on training operations and constraints of transplant operations are not considered).

Because of the complexities of this problem, no predictive scheduling considering all surgeries and all constraints can be done by the OR head nurse or manager of the operating theatre. So the comparison of the obtained schedule from HORS and the predictive schedule of the hospital is not possible. In HKC, after completing each surgery, the OR head nurse schedules the next case considering the actual system status at that moment. Hence, an online evaluation is proposed that includes the actual system status. In this validation method, the proposed model is updated each hour periodically. In addition the rescheduling is done after the end of each operation or emergency patient entry. Thereby, the Hybrid Rescheduling Policy is used.

The gap evaluations for the real instances show that the results are in general satisfactory. The smallest and greatest absolute gaps were 0 and 5 in the allowed computation time and 73% of tested cases reached optimal solutions in the specified time; even the time needed to find optimal solutions for the all of the test instances was lower than 1700 s (about 30 min) and could, hence, be still considered as reasonable. Figure 2 shows the frequency of the results in the range, verifying the efficiency of the proposed approach.

The result of comparison between actual scheduling and proposed model (HORS) based on the value of operating room idle time and max are presented in Table 2.

From the results presented in Table 2, the solutions provided by “HORS” scheduling are better than those of the HKC scheduling. Table 3 shows the mean of obtained result.

Averagely, idle time operating room in HKC was 323 minutes per day and max was 18:30 while these results for proposed model, respectively, were 22.5 minutes per day and 16:22.

Figures 3 and 4 show comparison of “operating room idle time” and max between the Hasheminejad scheduling (HKC) and the proposed model (HORS), respectively.

7. Conclusions and Recommendations

In this paper, the second step of the surgical process scheduling was investigated for scheduling elective patients waiting surgery. We formulated a mixed integer problem to assign a set of operations to some resources consisting of human resources (surgeon, assistant surgeon) and applied resource (OR) and sequenced them simultaneously. Taking into account the training plan of residents and fellows, real constraints in teaching hospitals and assignment of assistant surgeons to each operation form the main features of the proposed model.

The proposed method was evaluated and verified through solving several real instances at a teaching hospital. Also, the final solution for one month of chosen instances was explored and the actual scheduling of hospital compared with the proposed scheduling.

Numerical results indicated the efficiency of the proposed model compared to the actual hospital scheduling, and the gap evaluations for the real instances showed that the results were generally satisfactory.

The proposed model can be extended by including other resources such as ICU bed or ward bed in this model. We are working in a deterministic context; for future work, we will consider uncertainty related to surgery time and emergency patient arrival. Another suggestion for future work is to use the constraint programming, as a viable and powerful method for modeling and solving many combinatorial problems, to solve the described problem.

Appendix

The Designed ICR (Intelligent Character Recognition) Form

For more details, see Figure 5.

Acknowledgments

This research was conducted in part using data and resources from the Hasheminejad Kidney Center (HKC), a tertiary hospital of the Tehran University of Medical Sciences (TUMS). The authors are grateful to the contribution of Dr. Etemadian in providing access to key research materials. The contributions of HKC’s operating rooms staff are gratefully acknowledged.