- 153 Views
- Uploaded on

Download Presentation
## PowerPoint Slideshow about 'Stochastic Models for Operating Rooms Planning' - bradley-walter

**An Image/Link below is provided (as is) to download presentation**

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

Outline Outline Outline

for Operating Rooms Planning

Mehdi LAMIRI, Xiaolan XIE, Alexandre DOLGUI and Frédéric GRIMAUD

Centre Génie Industriel et Informatique

Centre Ingénierie et Santé

Outline

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

Problem description: Motivations

- Operating rooms represent one of the most expensive sector of the hospital
- Involves coordination of large number of resources
- Must deal with random demand for emergent surgery and unplanned activities
- Planning and scheduling operating rooms’ has become one of the major priorities of hospitals for reducing cost and improving service quality

Problem description

- How to plan elective cases when the operating rooms capacity is shared between two patients classes : elective and emergent patients
- Elective patients :

Electives cases can be delayed and planned for future dates

- Emergent patients :

Emergent cases arrive randomly and have to be performed in the day of arrival

Outline

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

TH

T1

H

1

2

Basic model: operating rooms capacity- We consider the planning of a set of elective surgery cases over an horizon of H periods under uncertain demand for emergency surgery
- Only aggregated capacity of the operating rooms is considered.

Tt : total operating rooms’ regular capacity

Ex: with a bloc of 5 operating rooms opened 10h each in day 1, T1 = 50 h

- Exceeding the regular capacity generates overtime costs (COt)

Basic model: Emergent patients

- In this work, the OR capacity needed for emergency cases for a period t is assumed to be a random variable ( wt ) based on:

- The distribution of the number of emergent patients in a given period estimated using information systems and / or by operating rooms’ manager

- The distribution of the OR time needed for emergent surgeries estimated from the historical data

TH

T1

Case1

Case 12

Case 5

Case 10

H

1

2

Basic model: Elective cases- At the beginning of the horizon, there are N requests for elective surgery

A plan that specifies the subset of elective cases to be performed in each period under the consideration of uncertain demand for emergency surgery

Basic model: Elective cases

- Each elective case i ( 1…N ) has the following characteristics :
- Operating Room Time needed for performing the case i : (pi)
- Estimated using information systems and/or surgeons’ expertises
- A release period (Bi)
- It represents hospitalisation date, date of medial test delivery
- A set of costs CEit ( t =Bi …H, H+1 )
- The CEit represents the cost of performing elective case i in period t
- CEi,H+1 : cost of not performing case i in the current plan

t

1

Bi

H

CEit

t

1

Bi

H

CEit

t

1

Bi

Li

H

Basic model: Elective cases related cost

- The cost structure is fairly general. It can represent many situations :
- Hospitalization costs / Penalties for waiting time
- Patient’s or surgeon’s preferences
- Eventual deadlines

Patient related cost

Overtime cost

Subject to:

overtime

Unplanned activities time

Planned activities time

Regular capacity

Basic model: Mathematical Model

- Decision:
- Assign case i to period t, Xit = 1
- Reject case i from plan, Xi,H+1 = 1

Basic model: Problem complexity

- The planning optimization problem is a stochastic combinatorial problem
- The stochastic planning problem is strongly NP-hard
- The resolution time increases exponentially as the size of the problem increases
- The problem is too difficult to be solved exactly within a reasonable amount of time

Outline

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

Monte Carlo optimization or Sample Average Algorithm

- Step 1. Generate randomly K different scenarios of emergency cases

Monte Carlo optimization or Sample Average Algorithm

- Step 2 : Estimate all performance functions according to these senarios

Monte Carlo optimization or Sample Average Algorithm

- Step 3 : Solving the optimisation problems using the estimated performance functions to obtain the sample optimum solution, also called Monte Carlo optimal solution

For our problem, the Monte Carlo optimization problem can be formulated as a mixed integer program

Monte Carlo optimization or Sample Average Algorithm

- Step 4 : Estimate the true criterion value for the Monte Carlo optimal solution

For our problem, it is estimated with a VERY large number K’ of senarios

- Step 1. Generate for each time period K samples of Wt
- Step 2. Formulate the Monte Carlo optimization problem (P')
- Step 3. Solve the mixed integer program (P') (with CPLEX in our case)
- Step 4. Evaluate the true criterion J(X*w,k) of the resulting Monte Carlo optimal solution X*w,k

- Theorem:

(a) limK Jk(X*w,K) = J(X*), where JK(X*w,K) is the “estimated” optimal criterion of the Monte Carlo optimization problem (P’) and J(X*) the “true” optimal criterion value

(b) N > 0, J(X*w,k) = J(X*), K > N, i.e. the Monte Carlo optimal solutionX*w,k becomes a true optimal solution.

Convergence : why it works well?

(Exponential convergence) :

P(X*w,k is not a true optimum) exp(- ck) where X*w,k isthe Monte Carlo optimal solution.

Note : The sample criterion Jw,k(X) converges to the true criterion Jw,k(X) at rate

(Optimal convergence rate) :

The convergence rate c is maximized by the common random variable scheme used in the Monte Carlo optimization

Remark: better convergence than independent evaluation of different solutions X.

Computation experiments: Comments

- Solutions provided by our optimization method are better than those of the deterministic method, even for small values of K (K=5)
- The proposed method achieves cost reduction of 4% with K=1000, comparing to the deterministic method.
- Disadvantage :
- The computation time increases beyond acceptable limit as the number of elective cases N > 60.

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

Model extensions: Multiple operating rooms

- Tts : regular capacity of OR-day (s, t)
- Wts : capacity needs for emergency cases in OR-day (s, t) (r.v.)
- CEits : cost of assigning case i to OR-day (s, t)
- COts : Overtime cost of OR-day (s, t)
- Decision variables:

Xits = 1 if case i is assigned to OR-day (s, t)

Model extensions: Overtime capacity and under utilization cost

- We introduce an additional penalty cost when the overtime capacity is exceeded

Operating Room related cost

ovetime cost

overtime capacity exceeded

under use

OR workload

regular capacity

overtime capacity

Overtime cost

Unplanned activities time

Planned activities time

Regular capacity

overtime

Extended model: Mathematical Model

Extended model: solution methods

- The Monte Carlo optimization method can be easily extended to solve the extended model
- The computation time quickly goes beyond acceptable limit as the number of cases increases
- Various methods such as Lagrangian relaxation and column generation have been tested.

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

- A “plan” is a possible assignment of patients to a particular OR-day
- p : plan for a particular OR-day is defined as follows
- aip = 1 if case i is in plan p
- btsp = 1 if plan p is assigned to OR-day (s, t)
- Cost of the plan :

Costs related to patients assigned to the paln

Overtime cost in the OR-day related to the plan

Column formulation for the planning problem

- : set of all possible plans
- Yp = 1, if plan p is selected and Yp = 0, otherwise

Master problem

Each patient is assigned at most to one selected plan

Subject to:

Each OR-day receives at most one plan

Relax the integrality constraints

Master Problem

Solve by Column Generation

Linear master problem (LMP)

Construct a “good” feasible solution

Optimal solution of the LMP

Near-optimal solution

Solving the linear master problem

simplex multipliers

pi , pt s

Pricing problem

minimizes reduced cost

Reduced Linear Master Problem

over Ω* Í Ω

min

min

st

st

reduced cost < 0

Y

add new column

N

STOP

- The pricing problem can be decomposed into H×M sub-problems
- One sub-problem for each OR-day

Simplex multipliers

Subject to:

Dynamic programming method

Relax the integrality constraints

Master Problem

Solve by Column Generation

Linear master problem (LMP)

Construct a “good” feasible solution

Optimal solution of the LMP

Near-optimal solution

Constructing a near optimal solution

- Step 1: Determine the corresponding patient assignment matrix (Xits) from the solution (Yp) of The Relaxed Master Problem.
- Step 2: Derive a feasible solution starting from (Xits)
- Step 3: Improve the solution obtained in Step 2

Derive a feasible solution

- Method I : Solving the integer master problem MP by restricting to generated columns
- Method II : Complete Reassignment

Fix assignment of cases in plans with Yp = 1

Reassign myopically but optimally all other cases one by one by taking into account scheduled cases.

- Method III : Progressive reassignment

Reassign each case to one OR-day by taking into account the current assignment (Xits) of all other cases, fractional or not.

Improvement of a feasible schedule

- Heuristic 1 : Local optimization of elective cases.

Reassign at each iteration the case that leads to largest improvement

- Heuristic 2 : Pair-wise exchange of elective cases (EX)
- Heuristic 3 : Period-based reoptimization (PB)

Re-optimize the planning of all cases assigned to a given OR-day (s, t) and all rejected cases.

Computation results

- The lower bound of the Column generation is very tight
- Solving the integer master problem with generated columns can be very poor and it very time consuming
- Progressive reassignment outperforms the complete reassignment as progressive reassignment preserves the solution structure of the column generation solution
- Numerical results show that the proposed solution methods are satisfying, and the best one can successfully solve the planning problem of about 240 interventions in 12 ORs within 10 minutes while promising the gap less than 0.6%.

- Motivation & Problem description
- Basic model
- Monte Carlo optimization method
- Model extensions
- A column generation approach
- Conclusions and perspectives

Conclusions and perspectives

- The proposed model can represent many real world constraints
- Monte Carlo simulation and MIP method provide good solutions
- Column generation is an efficient technique for providing provably good solutions in reasonable time for large problem.

Perspectives

- Make the stochastic model realistic enough to take into account random operating times, ...
- Develop exact algorithms able to solve problems with large size
- Test with field data

Download Presentation

Connecting to Server..