1.一种基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,所述方法包括以下步骤:1)获取历史数据信息,其中所述历史数据信息包括航季所有航班的计划信息以及前一年该航季每个月的历史运行数据信息;2)对得到的历史运行数据进行预处理,其中清洗掉取消航班的数据信息,并对清洗之后的航班数据进行数据筛选和历史数据匹配;3)构建数据驱动的优化模型,其中首先根据航班到时隙的分配定义0-1型整数决策变量,然后基于历史实际运行数据以最小化偏离航空公司请求时间和航班延误为目标建立最小化目标函数,随后确定约束条件,获得数据驱动的优化模型,其中,定义得到的决策变量为其中M表示所有航班号的集合,T表示所有时隙的集合,m∈M,t∈T,所建立的最小化目标函数为:其中f 1 (m,t)和f 2 (m,t)分别表示航班m分配到时隙t的计划延误代价函数和预期运行延误代价函数,并且其中f 1 (m,t)的计算公式为:f 1 (m,t)=c m d m |t-t m | (3)其中d m 表示航班m的实际执行的天数,c m 表示航班m延误的单位时间成本,t m 为航班m请求的时隙,t为待分配的时隙;f 2 (m,t)的计算公式为:f 2 (m,t)=c m d m E(g(m,t)),m∈M,t∈T (4)其中g(m,t)表示航班m在时隙t可能的延误时间,E(g(m,t))表示航班m在时隙t延误时间的期望值,c m 表示计划延误的单位时间成本,d m 表示航班m的实际执行的天数,约束条件包括如下:航班时隙唯一性约束:保证每个航班在同一天同一机场只能分配一个时隙航班时隙调整范围约束:保证航空公司可接受的分配时隙与请求时隙之差低于最大可接受值,即其中t max 为航班分配时隙与请求时隙的最大可接受值;容量约束:保证优化后的航班计划需要满足小时容量要求,为避免某些航班时隙过于密集,采用小时滚动容量限制;进场容量约束: 其中C a 为公布的最大进场容量,M a 为进场航班集合;离场容量约束: 其中C d 为公布的最大离场容量,M d 为离场航班集合;总容量约束: 其中C为机场运行的最大公布容量;最小过站时间约束:航空器滑行至停机位开启机门至航空器准备工作就绪关上机门之间的时间大于最小过站时间:其中 为航班m机型对应的最小过站时间,M 2 为过站航班的航班对集合,m a 为过站航班的进场航班,m d 为过站航班的离场航班;最大过站时间约束:为了不使航空公司在中转机场停留时间过长,有必要对中转航班的过站时间做如下最大过站时间限制:其中最大过站时间 为航班m机型对应的最大过站时间;4)根据建立的数据驱动的优化模型,利用切平面法求解,得到换季航班时隙的最优分配。
2.如权利要求1所述的基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,该航季的所有航班的计划信息包括航班计划中的航班号、机型、起航班场、计划起飞时间、落地机场、计划落地时间,而前一年该航季每个月的历史运行数据信息包括所有航班的运营日、航班号、计划起飞时间、计划落地时间、实际起飞时间、实际落地时间。
3.如权利要求1所述的基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,对清洗之后的航班数据进行数据筛选和历史数据匹配的步骤包括:对历史父权航班、历史非父权航班和新进入航班进行数据筛选,并对历史父权航班和历史非父权航班的历史运行数据进行匹配。
4.如权利要求1所述的基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,所述最大过站时间的确定是基于历史不同类型航班过站时间数据,利用百分位数法来确定,亦即假设所有过站航班可以分为Q类,记第q类航班的过站时间集合为 其中m q 表示第q类航班历史过站时间数据总个数,将集合S q 中的过站时间数据从小到大排序,排序后的集合记为 满足 其中 表示第k个小的过站时间;假设将p百分位值作为最大过站时间,根据百分位值的计算公式,得到第q类航班的最大过站时间为其中l为满足不等式l/m q <p<l+1/m q 的正整数。
5.如权利要求4所述的基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,步骤4)根据建立的数据驱动的优化模型,利用切平面法求解,得到换季航班时隙的最优分配包括:4.1)首先不考虑变量的取整约束来求解原整数规划对应的线性规划获得相应最优解;4.2)根据获得的最优解求切割方程;4.3)将所获得的切割方程添加到整数规划的约束条件中去,返回步骤4.1从而获得最优解。
6.如权利要求5所述的基于历史运行数据驱动的单机场航班换季时隙优化方法,其特征在于,所述步骤4)包括:4.1):不考虑变量的取整约束,求解原整数规划对应线性规划:设最优解为x * ;4.2):求一个切割方程:如果最优解的分量均为整数,则x * 为原整数规划的最优解;否则任选一个x * 中不是整数的分量,设其对应的基变量为x i ,x j 为非基变量,定义包含这个基变量的切割约束方程为:将b i 和a ij 都分解成整数部分N和非整数f之和,即:b i =N i +f i ,0<f i <1a ij =N ik +f ik ,0≤f ik ≤1 (10)其中N表示不超过b的最大整数,代入(9)式中得:提出变量为整数的条件: 这时式(11)中左边必须是整数,右边由于0<f i <1,不能为正,所以得到切割方程:4.3):将上述切割方程(12)添加到整数规划的约束条件中去,即对上述线性规划问题的可行域进行切割,然后返回步骤4.1。