The mathematical model for airline network seat inventory control problem is usually investigated to maximize the total revenue under some constraints such as capacities and demands. This paper presents a chance-constrained programming model based on the uncertainty theory for network revenue management, in which the fares and the demands are both uncertain variables rather than random variables. The uncertain programming model can be transformed into a deterministic form by taking expected value on objective function and confidence level on the constraint functions. Based on the strategy of nested booking limits, a solution method of booking control is developed to solve the problem. Finally, this paper gives a numerical example to show that the method is practical and efficient.
1. Introduction
After the deregulations in the airline industry, the revenue management techniques have become indispensable for airline seat inventory control. A central problem in airline revenue management is determining optimal decision rules for sequentially accepting or denying itinerary requests. So it is necessary to develop mathematical models to determine complex booking control strategies.
The optimization methods on seat inventory control problem with multiple-fare classes can be separated into the single-leg optimization method and the network optimization method. Based on the optimized frequency, each optimization method can be sorted out the static method and the dynamic method. In brief, the dynamic nature of the arrivals of the requests over time is not explicitly considered in the static method, whereas the mutability in the demand is taken into account in the dynamic method as the end of the reservation period approaches and the seat capacity diminishes.
The single-leg optimization method firstly appeared in the research of Littlewood [1]. He studied a seat inventory control problem with two-fare classes on a single leg and proposed a marginal seat revenue rule applied into a two-price, single-leg model. Belobaba [2, 3] extended this idea to a multiclass problem and introduced the expected marginal seat revenue heuristic for the general approach. Wollmer [4], Brumelle and McGill [5], and Robinson [6] further studied the single-leg problem with multiple-fare classes. They developed algorithms to find the optimal booking control policy under the assumption that the probability distributions of the demands for different fare classes were known. Lee and Hersh [7] developed a discrete-time dynamic programming model to find an optimal booking control policy without requiring any assumptions about the arrival mode for the manifold booking classes. Liang [8] proposed a continuous-time, stochastic, dynamic programming model and showed that a threshold control policy was optimal. Feng and Xiao [9] presented a stochastic control model to dynamically tackle with seat inventory control problem.
For the network optimization method, Glover et al. [10] initially described a minimum cost network flow formulation with deterministic demand without focusing on the stochastic elements. After that, a solution method for the sequential allocation of seats under the assumption of stochastic demand was provided by Wang [11]. Wollmer [12] proposed a linear programming model that considered stochastic demand. Dror et al. [13] proposed a similar deterministic network minimum cost flow formulation that allowed for cancellations as deterministic losses on arcs in the network. Curry [14] developed a combined mathematical model for a multiclass seat inventory control problem. Williamson [15] studied two network-based mathematical programming models. The first model incorporated probabilistic demand and the second model simplified the problem by substituting stochastic demand by its expectation. Wong et al. [16] applied nesting techniques into a multiclass seat inventory control problem. de Boer et al. [17] proposed stochastic linear programming for network revenue management and developed the nesting technique of Williamsion. Bertsimas and De Boer [18] and van Ryzin and Vulcano [19] used simulation-based optimization methods that also investigated nesting over the network. Cooper and Homem-De-Mello [20] proposed a decomposition method combining mathematical programming methods and Markov decision process. Recently, İlker Birbil et al. [21] proposed a framework for solving airline revenue management problems on large networks.
With the fast development of civil aviation industry, many airlines often create new routes. Due to lacking reliable data and accurate information, the methods in the literature above become invalid for these new routes. On the other hand, when unconventional sudden events such as war, atrocious weather, and earthquake, happen, the cumulative data of computer reservation system is no longer trustworthy. There are some limitations when the traditional stochastic models above deal with such problems at this situation. In the two cases, we have to invite some experts to evaluate their degree of belief that each event will occur. However, humans tend to overweigh unlikely events (Kahnema and Tversky [22]); thus, the degree of belief may have a much larger range than the real frequency. In this situation, if we insist on dealing with the degree of belief using the probability theory, some counterintuitive results will be obtained (Liu [23]). In revenue management of the above two cases, as we stated before, the domain experts invited are likely to overrate the market demand on the new routes and underestimate the market demand under the circumstances of unconventional sudden events. If the belief degree of the market demand is treated as probability, we have no choice but to increase the capacity on the new routes and reduce the capacity under the circumstances of unconventional sudden events. This will cause great losses in revenue for airlines. This conclusion seems unacceptable and then the belief degree cannot be treated as probability.
In order to deal with the experts' degree of belief, the uncertainty theory was founded by Liu [24] and refined by Liu [25] in 2013. Many researchers have contributed to this area. The uncertainty theory has been applied to uncertain programming, uncertain risk analysis, uncertain game, uncertain inference, uncertain logic, uncertain finance, and uncertain optimal control (Liu [25]). Nowadays, the uncertainty theory has become a branch of axiomatic mathematics to model human uncertainty (Liu [26]).
Depending on the analysis as mentioned above, we think that it is necessary to apply the uncertainty theory as a basic approach to model the uncertainty in the revenue management of the two cases above. In this paper, we propose the chance-constrained programming model based on the uncertainty theory to deal with the uncertain factors. The rest of this paper is structured as follows. In Section 2, some basic concepts and properties in uncertainty theory used throughout this paper are introduced. In Section 3, an uncertain programming model is constructed. According to inverse uncertainty distribution, the model can be transformed to its deterministic form. In Section 4, we present a solution method of booking control on the basis of the strategy of nested booking limits. After that, a numerical example is given in Section 5. At last, a brief summary is presented in Section 6.
2. Preliminaries
In this section, some basic definitions and arithmetic operations of uncertainty theory needed throughout this paper are presented.
Definition 1 (Liu [24]).
Let Γ be a nonempty set and L a σ-algebra over Γ. Each element Λ∈L is called an event. The set function M is called an uncertain measure if it satisfies the following four axioms:
Axiom 3 (Self-Duality). M{Λ}+M{Λc}=1 for any event Λ;
Axiom 4 (Countable Subadditivity). For every countable sequence of events {Λi}, we have(1)M{⋃i=1∞Λi}≤∑i=1∞M{Λi}.
Definition 2 (Liu [24]).
Let Γ be a nonempty set, L a σ-algebra over Γ, and M an uncertain measure. Then the triple (Γ,L,M) is called on uncertainty space.
Definition 3 (Liu [24]).
An uncertain variable ξ is a measurable function from an uncertainty space (Γ,L,M) to the set of real numbers; that is, for any Borel set B of real numbers, the set
(2){ξ∈B}={γ∈Γ∣ξ(γ)∈B}
is an event.
For a sequence of uncertainty variables ξ1,ξ2,…,ξn and a measurable function f, Liu [24] proved that ξ=f(ξ1,ξ2,…,ξn) defined as ξ(γ)=f(ξ1(γ),ξ2(γ),…,ξn(γ)),∀γ∈Γ is also an uncertain variable. In order to describe an uncertain variable, a concept of uncertainty distribution is introduced as follows.
Definition 4 (Liu [24]).
The uncertainty distribution Φ of an uncertain variable ξ is defined by
(3)Φ(x)=M{ξ≤x},
for any real number x.
Peng and Iwamura [27] proved that a function Φ:R→[0,1] is an uncertainty distribution if and only if it is a monotone increasing function except for Φ(x)≡0 or Φ(x)≡1. The inverse function Φ-1 is called the inverse uncertainty distribution of ξ. Inverse uncertainty distribution is an important tool in the operation of uncertain variables.
Theorem 5 (Liu [24]).
Let ξ1,ξ2,…,ξn be independent uncertain variables with regular uncertainty distributions Φ1,Φ2,…,Φn, respectively. If f(x1,x2,…,xn) is strictly increasing with respect to x1,x2,…,xm and strictly decreasing with respect to xm+1,xm+2,…,xn, then
(4)ξ=f(ξ1,ξ2,…,ξn)
is an uncertain variable with inverse uncertainty distribution
(5)Ψ-1(α)=f(Φ1-1(α),…,Φm-1(α),Φm+1-1(1-α),…,mmmmmlmlΦn-1(1-α)).
Expected value is the average of an uncertain variable in the sense of uncertain measure. It is an important index to rank uncertain variables.
Definition 6 (Liu [24]).
Let ξ be an uncertain variable. Then the expected value of ξ is defined by
(6)E[ξ]=∫0∞M{ξ≥r}dr-∫-∞0M{ξ≤r}dr,
provided that at least one of the two integrals is finite.
In order to calculate the expected value via inverse uncertainty distribution, Liu and Ha [28] proved that
(7)E[ξ]=∫01f(Φ1-1(α)),…,Φm-1(α),Φm+1-1(1-α),…,mmmmmlΦn-1(1-α)dα
under the condition described in Theorem 5. Generally, the expected value operator E has no linearity property for arbitrary uncertain variables. But, for independent uncertain variables ξ and η with finite expected values, we have
(8)E[aξ+bη]=aE[ξ]+bE[η],
for any real numbers a and b.
3. Uncertain Programming Model for Multiple-Leg Network Seat Inventory Control3.1. Problem Description
Because an airline wants to maximize revenue from the whole network, the researchers on this field focus on the network-based models now. Airlines usually provide thousands of such combinations of origin, destination and fare class (ODF). Therefore determining a comprehensive booking control strategy for the entire network is crucially important.
The objective of network seat inventory control is to maximize the airline's expected revenue from its supply of ODF combinations. Each ODF in the network is constitutive of one or more flight legs. The limited capacity on each flight leg has to be made full use of in the most profitable way. This can be achieved by limiting the number of seats available to the less lucrative classes. So the problem is to allocate all seats of each flight leg to the related ODF in the most profitable way. Due to its economic importance in the airline, the problem has been extensively studied.
In this paper, the network seat inventory control problem will be modeled by the chance-constrained programming based on the uncertainty theory in which the fare and the demand of each ODF are assumed to be uncertain variables with given uncertainty distributions.
3.2. Model Development
At first, we introduce the following notations to represent the mathematical formulation throughout the remainder of this paper:
: the number of seats reserved for each separate ODF;
: the total number of flight legs in the ODF network;
: the set of ODF combinations available on flight leg;
: the seat capacity on leg l;
: the deterministic aggregated demand for each ODF;
: the fare required for an ODF.
In order to facilitate the analysis, we make some reasonable assumptions as follows.
The flight market demand exceeds its capacity supply.
Overbooking is not considered by the model discussed here.
Next, based on the analysis of the decision making process, the general problem is formulated as follows [15]:
(9)max∑ODFfODFxODFs.t.ll∑ODF∈SlxODF≤Cl∀flightlegsl=1,…,N,lmmvlilxODF≤DODF∀ODF,lmmmllxODF≥0integer∀ODF.
In the above model, the quantities fODF and DODF are all assumed to be crisp numbers. However, when there are new routes created by the airlines or the emergency takes place sometimes, the quantities generally are not fixed but obtained from experience evaluation or expert knowledge. In this case, we may assume the quantities are uncertain variables. Then the model (9) is only a conceptual model rather than a mathematical model because there does not exist a natural ordership in an uncertain world. Here we take expected value criterion on the objective function and confidence level on the constraint functions (Liu [25]). Then the model (9) turns into the following mathematical model:
(10)maxE[∑ODFfODFxODF]s.t.l∑ODF∈SlxODF≤Cl∀flightlegsl=1,…,N,mmmllmmlM{xODF≤DODF}≥βODF∀ODF,mlmmlxODF≥0integer∀ODF,
where βODF are some predetermined confidence levels for all ODF.
In practical applications, the uncertainty distributions of uncertain variables fODF and DODF and the confidence levels βODF are determined by linear interpolation method, the principle of least squares, the method of moments, and the Delphi method from expert's experimental data (Liu [25]). How do we obtain expert's experimental data? Liu [25] proposed a questionnaire survey for collecting expert’s experimental data. In this paper, we assume that the uncertainty distributions of uncertain variables fODF and DODF and the confidence levels βODF have been determined.
In order to solve model (10), firstly, we introduce two corollaries which were from the uncertainty theory (Liu [25]).
Corollary 7.
Assume the objective function f(x,ξ1,ξ2,…,ξn) is strictly increasing with respect to ξ1,ξ2,…,ξm and strictly decreasing with respect to ξm+1,ξm+2,…,ξn. If ξ1,ξ2,…,ξn are independent uncertain variables with uncertainty distributions Φ1,Φ2,…,Φn, respectively, then the expected objective function E[f(x,ξ1,ξ2,…,ξn)] is equal to
(11)∫01f(x,Φ1-1,…,Φm-1,Φm+1-1,…,Φn-1)dα.
Corollary 8.
Assume the constraint function g(x,ξ1,ξ2,…,ξn) is strictly increasing with respect to ξ1,ξ2,…,ξk and strictly decreasing with respect to ξk+1,ξmk2,…,ξn. If ξ1,ξ2,…,ξn are independent uncertain variables with uncertainty distributions Φ1,Φ2,…,Φn, respectively, then the chance constrain
(12)M{g(x,ξ1,ξ2,…,ξn)≤0}≥α
holds if and only if
(13)g(x,Φ1-1(α),…,Φk-1(α),Φk+1-1(1-α),…,Φn-1(1-α))≤0.
Secondly, the next theorem shows that the model (10) is equivalent to a deterministic model, for which many efficient algorithms have been designed.
Theorem 9.
Assume that fODF and DODF are independent uncertain variables with uncertainty distributions ϕODF and ψODF. Then the model (10) is equivalent to the following model:
(14)max∑ODFxODF∫01ϕODF-1(α)dαs.t.0000∑ODF∈SlxODF≤Cl∀flightlegsl=1,…,N,ll00000mmxODF≤ψODF-1(1-βODF)∀ODF,ll00000mmxODF≥0integer∀ODF.
Proof.
The function ∑ODFfODFxODF is strictly increasing with respect to fODF and fODF are independent uncertain variables with uncertainty distributions ϕODF, respectively. By using Corollary 7, we obtain
(15)E[∑ODFfODFxODF]=∑ODFxODF∫01ϕ-1ODF(α)dα.
Since that
(16)M{xODF≤DODF}≥βODF
is equivalent to
(17)M{-DODF+xODF≤0}≥βODF
and the function -DODF+xODF is strictly decreasing with respect to DODF with uncertainty distribution ψODF, it follows from Corollary 8 that we have
(18)-ψODF-1(1-βODF)+xODF≤0,
that is,
(19)xODF≤ψODF-1(1-βODF).
The theorem is thus verified.
4. Solution Method of Booking Control4.1. The Strategy of Nested Booking Limits
In the model (14), we use f¯ODF to denote ∫01ϕODF-1(α)dα and D¯ODF to denote ψODF-1(1-βODF). Then the relaxation of the model (14) is
(20)max∑ODFf¯ODFxODFs.t.mmllll∑ODF∈SlxODF≤Cl∀flightlegsl=1,…,N,mmmmmllllixODF≤D¯ODF∀ODF,mmlmmmlmxODF≥0∀ODF.
The whole booking time should be partitioned into a few time periods of reservation; for example, a day is a time period of reservation. In order to facilitate the analysis, we describe the model (20) in the form of the matrix and the vector considering the time period of reservation. For this, we introduce the following notations:
I: the total number of flight legs in the ODF network;
i: index for set of flight legs;
J: the total number of the ODF;
j: index for set of the ODF;
t: index for the time period of reservation;
Ct=[C1t,C2t,…,Cit,…,CIt]T: the seat capacity on each flight leg in the reservation time period t;
Ft=[f-1t,f-2t,…,f-jt,…f-Jt]T: the corresponding f¯ODF for each separate ODF in the reservation time period t;
Dt=[D-1t,D-2t,…,D-jt,…,D-Jt]T: the corresponding D¯ODF for each separate ODF in the reservation time period t;
Xt=[x1t,x2t,…,xjt,…,xJt]T: the number of seats reserved for each separate ODF in the reservation time period t;
A=(aij)I×J: the matrix denotation of flight legs that each ODF travels, where if the ODF j travels the flight leg i, then aij=1, otherwise aij=0;
Aj: the jth column of the matrix A, denoting the flight legs that the ODF j travels.
The model (20) is described as follows in the form of the matrix and the vector:
(21)LP(Ct,Dt)=max(Ft)TXtmls.t.mmmmvAXt≤Ct,mmlmmv0≤Xt≤Dt.
The dual problem of the model above can be described as follows:
(22)DLP(Ct,Dt)=min[(Pt)TCt+(Qt)TDt]mls.t.mmmmmATPt+Qt≥Ft,mmmmlvPt,Qt≥0,
where Pt=[P1t,P2t,…,Pit,…,PIt]T denotes the shadow prices corresponding to each flight leg and Qt=[Q1t,Q2t,…,Qjt,…,QJt]T denotes the shadow prices corresponding to each D¯ODF.
Bid price control method is one of the prevalent methods of network seat inventory control. The bid price of each ODF is equal to the sum of shadow prices of the flight legs that the ODF crosses. A booking request for a passenger from the ODF is rejected if the bid price of the ODF exceeds the fare for the ODF and is accepted otherwise. Although bid price control method has been used in the actual operations of the airlines, it has a few shortcomings as follows.
Each ODF’s contribution to network revenue is not considered in the bid price control method.
When calculating shadow prices using the relevant models, the solution of the model may be degenerate solution. This will cause the multiple bid prices of an ODF.
The fares of most of the passengers on the flight just exceed the bid prices so that airlines suffer losses.
For this, we present a nesting control method based on the network contribution value for the above uncertain programming model.
First, we define an ODF’s net contribution value to network revenue in the reservation time period t as the expected fare for the ODF in the reservation time period t minus the opportunity cost of the ODF in the reservation time period t, that is,
(23)NCVjt=f-jt-OCjt,
where NCVjt denotes the net contribution value of the ODF j to network revenue in the reservation time period t and OCjt denotes the opportunity cost of the ODF j in the reservation time period t.
The opportunity cost OCjt of the ODF j in the reservation time period t is calculated based on the DLP model as follows:
(24)OCjt=(Pt)TAj
for 1≤j≤J.
Now the opportunity cost of the ODF in the reservation time period t is known, so we figure up the net contribution value of the ODF to network revenue in the reservation time period t. However, the solution of the model DLP may be degenerate solution and this phenomenon will cause multiple net contribution values of the ODF to network revenue. In this case, we will use the following method to calculate the net contribution value of the ODF to network revenue in the reservation time period t. If Ct-Aj≥0, then NCVjt=f-jt-[LP(Ct,Dt)-LP(Ct-Aj,Dt)]. Otherwise,NCVjt=-∞ where ∞ denotes a large enough number.
Finally, we rank the ODF on the basis of their net contribution value to network revenue to determine the nesting level. If some of the ODFs have the same net contribution value to network revenue, we can rank the ODF on the basis of their expected fare.
4.2. The Algorithm for Nested ODF-Based Booking Control
Every time a booking request arrives for any ODF in the network, a quick decision should be made whether or not to accept the request. We have to specify a booking control strategy for the decision. We propose the algorithm for nested ODF-based booking control.
The notations used in the following algorithm are given as below:
Xt=[xt1,xt2,…,xtj,…,xtJ]T: the number of seats reserved for each separate ODF in the reservation time period t after ranking the elements of Xt according to the nesting level;
PLt=[PLt1,PLt2,…,PLtj,…,PLtJ]T: seat protect level for each ranked ODF in the reservation time period t;
Mk(1≤k≤j): the set of flight legs that the kth ODF travels;
bt=[bt1,bt2,…,btj,…,btJ]: the number of booking requests for the ranked ODF that have already been accepted in the reservation time period t;
BLtj: seat booking limit for the jth ODF.
The heuristic algorithm for nested ODF-based booking control is as follows.
Step 1.
Calculate each ODF's net contribution value to network revenue in the reservation time period t and determine the nesting level. By solving the model above, we can obtain Xt=[x1t,x2t,…,xjt,…,xJt]T.
Step 2.
Determine seat protect level for each ODF. According to ranking results of Step 1, we rank the elements of Xt and obtain Xt=[xt1,xt2,…,xtj,…,xtJ]T. Seat protect level is determined as PLt=[PLt1,PLt2,…,PLtj,…,PLtJ]=[xt1,xt2,…,xtj,…,xtJ]T.
Step 3.
At the beginning of the reservation time period t, it is obvious that bt=0. Determine seat booking limit for the ranked ODF as follows.
For the first ODF,
(25)BLt1=min{Nit∣1≤i≤I,i∈M1};
for the jth ODF (2≤j≤J), on the flight leg i(1≤i≤I), the total number of the seats protected for the lower level ODF than the jth ODF is Πit=∑k=1j-1PRkt, where the kth ODF travels the flight leg i and the number of the seats available on the flight leg i is (Cit-Πit). So we have
(26)BLtj=min{Cit-Πit∣1≤i≤I,i∈Mj}.
Step 4.
If BLtj>btj, then accept the booking request, let btj=btj+1; if BLtj=btj, then decline the request.
Step 1 determines the nesting level for each separate ODF; Step 2 determines seat protect level for each ranked ODF; Step 3 determines seat booking limit for the ranked ODF; Step 4 develops the standard of accepting or rejecting the booking requests.
When entering the next reservation time period, let
(27)Ct=Ct-btT[A1A2⋯Aj⋯AJ]T,t=t-1,bt=0,
and then go into the algorithm above.
5. Numerical Experiment
In this section, we apply the model and the algorithm of the former two sections to airline seat inventory control and give an optimal policy.
An airline prepares to open a new route, which is from A to C through B. So there are two legs and three segments. Every segment has two fare classes with Y and T, so there are six ODFs. These are shown in Figure 1.
Schematic.
That is,
(28)N=2,l=1,2,S1={ABY,ABT,ACY,ACT},S2={BCY,BCT,ACY,ACT}.
Furthermore, there are 140 seats available in the flight and the booking period is partitioned into 11 time intervals.
In order to simulate the booking process, some simulated reservation data about each time interval are given in Table 1.
Simulated reservation data.
ODF
t = 10
t = 9
t = 8
t = 7
t = 6
t = 5
t = 4
t = 3
t = 2
t = 1
t = 0
ABY
3
3
3
5
4
1
3
7
8
11
8
ABT
13
5
8
7
1
2
6
2
5
5
2
BCY
0
2
1
5
4
8
6
4
7
7
4
BCT
15
5
11
9
9
8
2
3
4
2
0
ACY
2
0
0
1
2
9
1
7
5
3
5
ACT
10
2
5
6
2
4
7
4
3
1
1
When t=10, there is C1=C2=140.
According to experts' experience, the demands for ABY, ABT, BCY, BCT, ACY, and ACT follow a normal uncertainty distribution N(ei,σi), i=1,2,3,4,5,6, respectively. The fare fi for ABY, ABT, BCY, BCT, ACY, and ACT follows a normal uncertainty distribution N(ei′,σi′), i=1,2,3,4,5,6, respectively. Tables 2 and 3 give the value of the quantities.
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(42, 2)
(66, 2)
(41, 3)
(71, 4)
(35, 1)
(50, 3)
Parameters of normal distribution of fare.
N(ei′,σi′)
1
2
3
4
5
6
(1000, 80)
(800, 40)
(400, 20)
(320, 10)
(1200, 90)
(960, 50)
Note that the normal uncertain variable N(e,σ) has an expected value e and an inverse uncertainty distribution
(29)Φ-1(x)=e+3σπln1-xx.
Furthermore, β=0.9.
According to Section 3 and Section 4, the result of the seat inventory control for t=10 is obtained in Table 4.
The seat inventory control for t = 10.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 10
ABY
1000
200
42
140
3
3
ABT
800
0
63
63
13
13
BCY
400
80
41
105
0
0
BCT
320
0
64
64
15
15
ACY
1200
80
35
98
2
2
ACT
960
−160
10
0
10
0
When t=9,8,7,6,5,4,3,2,1,0, the results of the seat inventory control are listed in Appendix.
Finally, the result of simulation is listed as Table 5.
The seat inventory control result.
The combination of the ODF
The total number of accepted seats
ABY
52
ABT
54
BCY
48
BCT
58
ACY
30
ACT
4
The total expected revenue
172800
The numerical experiment shows that the computation time of the algorithm is 0.4–0.8 seconds.
We use bid price control method for the above simulation data and the result is listed as Table 6.
The seat inventory control result of the bid price control method.
The combination of the ODF
The total number of accepted seats
ABY
48
ABT
54
BCY
45
BCT
57
ACY
30
ACT
8
The total expected revenue
171120
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(40, 2)
(56, 2)
(39, 3)
(61, 4)
(34, 1)
(42, 3)
The seat inventory control for t = 9.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 9
ABY
1000
200
40
122
3
3
ABT
800
0
48
48
5
5
BCY
400
80
39
89
2
2
BCT
320
0
50
50
5
5
ACY
1200
80
34
82
0
0
ACT
960
−160
0
0
2
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(38, 2)
(47, 2)
(36, 3)
(51, 4)
(32, 1)
(36, 3)
The seat inventory control for t = 8.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 8
ABY
1000
200
38
114
3
3
ABT
800
0
44
44
8
8
BCY
400
80
36
84
1
1
BCT
320
0
48
48
11
11
ACY
1200
80
32
76
0
0
ACT
960
−160
0
0
5
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(36, 2)
(39, 2)
(33, 3)
(42, 4)
(30, 1)
(30, 3)
The seat inventory control for t = 7.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 7
ABY
1000
200
36
103
5
5
ABT
800
0
37
37
7
7
BCY
400
80
33
74
5
5
BCT
320
0
41
41
9
9
ACY
1200
80
30
67
1
1
ACT
960
−160
0
0
6
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(33, 2)
(32, 2)
(30, 3)
(34, 4)
(28, 1)
(24, 3)
The seat inventory control for t = 6.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 6
ABY
1000
200
23
90
4
4
ABT
800
0
29
29
1
1
BCY
400
80
30
61
4
4
BCT
320
0
31
31
9
9
ACY
1200
80
28
57
2
2
ACT
960
−160
0
0
2
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(30, 2)
(25, 2)
(27, 3)
(26, 4)
(25, 1)
(19, 3)
The seat inventory control for t = 5.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 5
ABY
1000
200
30
83
1
1
ABT
800
0
25
28
2
2
BCY
400
80
27
49
8
8
BCT
320
0
19
19
8
8
ACY
1200
80
25
53
9
9
ACT
960
−160
3
3
4
3
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(27, 2)
(20, 2)
(23, 3)
(20, 4)
(22, 1)
(15, 3)
The seat inventory control for t = 4.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 4
ABY
1000
200
27
68
3
3
ABT
800
0
19
19
6
6
BCY
400
80
23
24
6
6
BCT
320
0
1
1
2
1
ACY
1200
80
22
41
1
1
ACT
960
−160
0
0
7
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(23, 2)
(15, 2)
(19, 3)
(14, 4)
(19, 1)
(11, 3)
The seat inventory control for t = 3.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 3
ABY
1000
440
23
58
7
7
ABT
800
240
15
16
2
2
BCY
400
0
18
18
4
4
BCT
320
−80
0
0
3
0
ACY
1200
240
19
35
7
7
ACT
960
0
1
1
4
1
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(19, 2)
(11, 2)
(15, 3)
(10, 4)
(16, 1)
(8, 3)
The seat inventory control for t = 2.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 2
ABY
1000
200
19
41
8
8
ABT
800
0
6
6
5
5
BCY
400
0
10
10
7
7
BCT
320
−80
0
0
4
0
ACY
1200
0
16
22
5
5
ACT
960
−240
0
0
3
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(14, 2)
(8, 2)
(10, 3)
(6, 4)
(12, 1)
(5, 3)
The seat inventory control for t = 1.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 1
ABY
1000
200
14
23
11
11
ABT
800
0
5
5
5
5
BCY
400
0
10
10
7
7
BCT
320
−80
0
0
2
0
ACY
1200
0
4
9
3
3
ACT
960
−240
0
0
1
0
Parameters of normal distribution of demand.
N(ei,σi)
1
2
3
4
5
6
(9, 2)
(3, 2)
(5, 3)
(2, 1)
(7, 1)
(2, 1)
The seat inventory control for t = 0.
ODF
fjt-
NCVjt
xjt
BLtj
Simulated reservation data
The number of accepted seats
t = 0
ABY
1000
0
4
4
8
4
ABT
800
200
0
0
2
0
BCY
400
0
4
4
4
4
BCT
320
−80
0
0
0
0
ACY
1200
−200
0
0
5
0
ACT
960
−440
0
0
1
0
Comparing Tables 5 and 6, we can see that the total expected revenue using the bid price control method is 171120 RMB, while the total expected revenue using the proposed method in this paper is 172800 RMB. Numerical simulation results show that the proposed method in this paper is effective for improving the airline's revenue.
6. Conclusions
To consider network revenue management problem under conditions of new routes and unconventional sudden events, we established an uncertain programming model. Based on the strategy of nested booking limits, a heuristic algorithm for booking control was developed. Numerical test was performed to evaluate the model and the solution algorithm. The test results show that the model and the solution are all effective. There are some suggestions for future research:
the impact of the nesting heuristics on total revenue;
more complex hub-spoke network;
considering dynamic factors such as the arrival order of the requests;
integrated uncertain and stochastic model.
Appendix
When t=9, there is C1=122, C2=123, β=0.8. (See Tables 7, 3, and 8.)
When t=8, there is C1=114, C2=116, β=0.8. (See Tables 9, 3, and 10.)
When t=7, there is C1=103, C2=104, β=0.8. (See Tables 11, 3, and 12.)
When t=6, there is C1=90, C2=89, β=0.8. (See Tables 13, 3, and 14.)
When t=5, there is C1=83, C2=74, β=0.8. (See Tables 15, 3, and 16.)
When t=4, there is C1=68, C2=46, β=0.8. (See Tables 17, 3, and 18.)
When t=3, there is C1=58, C2=38, β=0.8. (See Tables 19, 3, and 20.)
When t=2, there is C1=41, C2=26, β=0.8. (See Tables 21, 3, and 22.)
When t=1, there is C1=23, C2=14, β=0.8. (See Tables 23, 3, and 24.)
When t=0, there is C1=4, C2=4, β=0.8. (See Tables 25, 3, and 26.)
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This work has been financed by the Fundamental Research Funds for the Central Universities (no. 3122014 K007).
LittlewoodK.Forecasting and control of passengers19721295117BelobabaP. P.1987Cambridge, Mass, USAFlight Transportation Laboratory, Massachusetts Institute of TechnologyBelobabaP. P.Application of a probabilistic decision model to airline seat inventory control19893721831972-s2.0-0024629453WollmerR. D.An airline seat management model for a single leg when lower fare classes book first1992402637BrumelleS. L.McGillJ. I.Airline seat allocation with multiple nested fare classes199341127137RobinsonL. W.Optimal and approximate control policies for airline booking with sequential nonmonotonic fare classes199543252263LeeT. C.HershM.Model for dynamic airline seat inventory control with multiple seat bookings19932732522652-s2.0-0027647530LiangY.Solution to the continuous time dynamic yield management model19993311171232-s2.0-0000492371FengY.XiaoB.A dynamic airline seat inventory control model and its optimal policy20014969389492-s2.0-0035521463GloverF.GloverR.LorenzoJ.McMillanC.The passenger mix problem in the scheduled airlines1982127379WangK.Optimal seat allocations for multi-leg flights with multiple fare types198323225246WollmerR. D.A hub-spoke seat management model1986Long Beach, Calif, USAMcDonnell Douglas CorporationDrorM.TrudeauP.LadanyS. P.Network models for seat allocation on flights19882242392502-s2.0-0011407208CurryR. E.Optimal airline seat allocation with fare classes nested by origins and destinations19902431932042-s2.0-0025477551WilliamsonE. L.1992Cambridge, Mass, USAMassachusetts Institute of TechnologyWongJ. T.KoppelmanF. S.DaskinM. S.Flexible assignment approach to itinerary seat allocation199327133482-s2.0-0006857511de BoerS. V.FrelingR.PiersmaN.Mathematical programming for network revenue management revisited2002137172922-s2.0-003711675410.1016/S0377-2217(01)00096-0BertsimasD.De BoerS.Simulation-based booking limits for airline revenue management2005531901062-s2.0-1464442900210.1287/opre.1040.0164van RyzinG.VulcanoG.Simulation-based optimization of virtual nesting controls for network revenue management20085648658802-s2.0-6134910379510.1287/opre.1080.0550CooperW. L.Homem-De-melloT.Some decomposition methods for revenue management20074133323532-s2.0-7044962352410.1287/trsc.1060.0184İlker BirbilŞ.FrenkJ. B. G.GromichoJ. A. S.ZhangS.A network airline revenue management framework based on decomposition by origins and destinations201310.1287/trsc.2013.0469KahnemaD.TverskyA.Prospect theory: an analysis of decision under risk197947263292LiuB.Why is there a need for uncertainty theory?2012613102-s2.0-84863012623LiuB.20072ndBerlin, GermanySpringerMR2683154LiuB.20134thBerlin, GermanySpringerLiuB.2010Berlin, GermanySpringerPengZ.IwamuraK.A sufficient and necessary condition of uncertainty distribution201013327728510.1080/09720502.2010.10700701MR2721956ZBL1229.28029LiuY.HaM.Expected value of function of uncertain variables20104181186