CN107341620B - 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法 - Google Patents

基于bada燃油消耗率的短期天气下进场航班延误成本计算方法 Download PDF

Info

Publication number
CN107341620B
CN107341620B CN201710614866.5A CN201710614866A CN107341620B CN 107341620 B CN107341620 B CN 107341620B CN 201710614866 A CN201710614866 A CN 201710614866A CN 107341620 B CN107341620 B CN 107341620B
Authority
CN
China
Prior art keywords
flight
approach
cost
fuel consumption
delay
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710614866.5A
Other languages
English (en)
Other versions
CN107341620A (zh
Inventor
张明
刘凯
孔祥鲁
黄倩文
喻珏
张一帆
仇志峰
吴翰林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201710614866.5A priority Critical patent/CN107341620B/zh
Publication of CN107341620A publication Critical patent/CN107341620A/zh
Application granted granted Critical
Publication of CN107341620B publication Critical patent/CN107341620B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/40Business processes related to the transportation industry

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Game Theory and Decision Science (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Traffic Control Systems (AREA)
  • Navigation (AREA)

Abstract

本发明公开一种基于BADA燃油消耗率的短期天气下进场航班延误成本计算方法,即,短期天气下进场航班延误总成本包括所有航空器执飞进场和进近程序因航班延误产生的时间成本和燃油成本之和,其中,进场及进近航段包括:进场航段、起始进近航段、中间进近航段和最后进近航段;燃油成本是基于BADA数据库中相关机型发动机的燃油消耗率计算得到的因航班延误产生的燃油消耗成本;时间成本包括因航班延误产生的航空器相关的机组小时费成本、维修费成本、飞机和发动机折旧费成本以及其它费用成本。本方法从航空器的发动机耗油率角度可分析不同进场运行航段不同机型的耗油状况,从而为准确航空器耗油成本计算提供更为科学和准确的依据。

Description

基于BADA燃油消耗率的短期天气下进场航班延误成本计算 方法
技术领域
本发明属于民航技术领域,具体涉及一种短期天气下进场航班延误成本计算方法。
背景技术
终端空域存在短期天气变化,引起终端空域容量和运行效率急剧下降,航班延误成本的急剧上升,影响航班的正常进场,甚至有可能引发飞行安全事故。恶劣天气变化是航班延误的重要原因。相关数据表明,中国2015年客运航空公司共执行航班337.3万班次,平均航班正常率仅为68.33%,天气原因引起的航班延误所占比例29.53%;2015年,美国国家空域系统由天气原因导致的航班延误比例高达53.1%。欧洲的单一天空计划(SESAR)和美国的下一代航空运输系统(NEXTGEN)都将天气变化下的空中交通延误分析作为重要的研究议题。
航班延误成本分析是飞行运行决策需要考虑的重要方面。对于短期天气变化下的终端区某进场航班来说,进场航路通常被气象云团堵塞而必须实施改航策略,航空器在原来既定航路基础上多飞行了一定距离;短期天气条件变化引起的机场跑道容量减小,进场航班服务率降低,会在实施进近程序之前进行空中盘旋等待。以上改航策略和空中等待程序的实施将导致航班的燃油消耗量和进场时间增多,引起燃油成本和时间成本的增大,即航班延误成本增大,给航空公司带来直接经济损失。因此,获取准确的航班延误运行成本应当是除了飞行安全外,进场过程中改航和等待运行决策需要考虑的重要依据。
现有技术中已从多方面研究了航班延误对航空公司飞行运行成本和旅客出行成本的影响,但在空中交通流量管理预战术阶段,终端空域的航班由于受到天气影响产生的改航和等待等飞行预案,需要结合原飞行计划和天气雷达图,建立新的飞行改航和等待航迹,并考虑不同机型和航段燃油消耗的差异,才能准确获取航空器进场的航班延误运行成本。而目前研究工作较少从该角度讨论终端空域恶劣天气影响下航班改航和空中等待带来的延误时间成本和燃油成本。
发明内容
为精确计算终端区内进场航班的延误成本,本发明针对短期天气影响下的终端区改航策略和空中等待程序造成的航班延误成本问题,将航班的延误时间成本和燃油成本进行定量分析,提出一种短期天气下进场航班延误运行成本预测方法。
本发明所公开一种基于BADA燃油消耗率的短期天气下进场航班延误成本计算方法,即,短期天气下进场航班延误总成本包括所有航空器执飞进场和进近程序因航班延误产生的时间成本和燃油成本之和,其中,进场及进近航段包括:进场航段、起始进近航段、中间进近航段和最后进近航段;燃油成本是基于BADA数据库中相关机型发动机的燃油消耗率计算得到的因航班延误产生的燃油消耗成本;时间成本包括因航班延误产生的航空器相关的机组小时费成本、维修费成本、飞机和发动机折旧费成本以及其它费用成本。
进一步的,因航班延误产生的燃油成本
Figure BDA0001360347100000021
为:
Figure BDA0001360347100000022
式中,f表示进场航班,Ce表示单位燃油价格,Ef表示燃油消耗量,
Figure BDA0001360347100000023
表示航班计划进场燃油消耗量,
Figure BDA0001360347100000024
表示航班实际进场燃油消耗量,
Figure BDA0001360347100000025
表示航班空中等待燃油消耗量;
Figure BDA0001360347100000026
式中,
Figure BDA0001360347100000027
表示航班空中等待燃油消耗量,Ne表示发动机数量,
Figure BDA0001360347100000028
表示航班空中等待时间,fcr表示航班空中等待阶段的燃油消耗率。
进一步的,利用BADA数据库,并结合终端空域恶劣天气影响下的航空器改航策略和空中等待程序,将航空器进场和进近飞行程序划分为若干个平飞和下降航段,将航空器执飞进场和进近程序所需的燃油消耗量结合实际的航空器改航情况进行平飞航段和下降航段的叠加组合,并建立航空器执飞进场程序的燃油消耗量计算模型,基于该燃油消耗量计算模型计算得到的所有航空器执飞进场和进近程序因航班延误产生的燃油消耗量,其中,燃油消耗量计算模型是:
Figure BDA0001360347100000029
式中,
Figure BDA0001360347100000031
表示燃油消耗量,p表示机型,f表示进场航班,r表示执行的进场程序,a表示平飞航段数量,b表示下降航段数量,Ne表示发动机数量,
Figure BDA0001360347100000032
表示第i个平飞航段的飞行速度,
Figure BDA0001360347100000033
表示第i平飞航段的飞行距离,
Figure BDA0001360347100000034
表示第i平飞航段的燃油消耗率,
Figure BDA0001360347100000035
表示第j个下降航段的飞行速度,
Figure BDA0001360347100000036
表示第j个下降航段的飞行距离,
Figure BDA0001360347100000037
表示第j个下降航段的燃油消耗率。
进一步的,因航班延误产生的时间成本为:
Figure BDA0001360347100000038
式中,
Figure BDA0001360347100000039
表示因航班延误产生的燃油成本,f表示航班类型,Ct表示单位时间延误成本,Tf表示延误时间,
Figure BDA00013603471000000310
表示机组小时费,
Figure BDA00013603471000000311
表示维修费,
Figure BDA00013603471000000312
表示飞机和发动机折旧费,
Figure BDA00013603471000000313
表示其它费用;
Figure BDA00013603471000000314
表示实际进场时间,
Figure BDA00013603471000000315
表示计划航班执行进场程序的时间,
Figure BDA00013603471000000316
表示航班空中等待时间。
进一步的,机组小时费包括机长、副驾、安全员、乘务长和乘务员的小时工资;飞机和发动机的折旧费按照飞机价格除以飞机使用小时数得到。
进一步的,其它费用包括根据延误时间的长短酌情给予旅客和/或机组小组的食宿费用。
本发明还公开一种基于BADA燃油消耗率的短期天气下进场航班延误成本计算方法,包括以下步骤:
步骤1,建立终端空域进场航路结构,获取终端空域范围内的气象雷达回波图像及进场航班计划;
步骤2,基于图像颜色与降水等级、回波强度的关系,对气象雷达回波图像进行颜色识别与轮廓提取;
步骤3,确定受天气影响的进场航路,规划与天气状况对应的改航航路;
步骤4,确定航班的执飞机型,并查阅BADA数据库关于航空器燃油消耗率的性能参数,并利用BADA数据库关于航空器进场平飞和下降阶段的燃油消耗率模型,将进场航路划分为若干个平飞和下降航段,建立进场航班的燃油消耗计算模型;
步骤5,确定进场航路的航路点、航段距离高度以及航空器速度,分析航空器在各个航段的飞行姿态,并计算计划进场航路和改航航路各个航段的航空器燃油消耗率;
步骤6,根据航班进场计划航路与改航航路的航路数据,计算因航班延误产生的延误时间及其产生的时间成本;
步骤7,利用燃油消耗计算模型计算航班执行计划进场航路和改航航路的燃油消耗量,根据计划和改航的航空器燃油消耗量计算因航班延误产生的燃油成本;
步骤8,计算航班延误总成本,即所有航空器执飞进场和进近程序时因航班延误产生的时间成本和燃油成本之和。
进一步的,各航班有且仅能由一种机型p,执行进场程序r降落于机场跑道,延误总成本CDelay的计算公式为:
Figure BDA0001360347100000041
Figure BDA0001360347100000042
Figure BDA0001360347100000043
式中,F表示终端区进场航班的集合,
Figure BDA0001360347100000044
P表示终端区进场航班的航空器机型集合,
Figure BDA0001360347100000045
R表示终端区进场进近程序的集合,
Figure BDA0001360347100000046
Figure BDA0001360347100000047
表示航班f的状态变量,为0-1离散变量;
Figure BDA0001360347100000048
表示航班f的延误成本;
Figure BDA0001360347100000049
表示航班f的延误时间成本;
Figure BDA00013603471000000410
表示航班f的延误燃油成本。
进一步的,因航班延误产生的燃油成本为:
Figure BDA00013603471000000411
式中,
Figure BDA00013603471000000412
表示因航班延误产生的燃油成本,f表示进场航班,Ce表示单位燃油价格,Ef表示燃油消耗量,
Figure BDA00013603471000000413
表示航班计划进场燃油消耗量,
Figure BDA00013603471000000414
表示航班实际进场燃油消耗量,
Figure BDA0001360347100000051
表示航班空中等待燃油消耗量;
其中,航班空中等待燃油消耗量
Figure BDA0001360347100000052
Figure BDA0001360347100000053
式中,
Figure BDA0001360347100000054
表示航班空中等待燃油消耗量,Ne表示发动机数量,
Figure BDA0001360347100000055
表示航班空中等待时间,fcr表示航班空中等待阶段的燃油消耗率;
其中,航班计划进场和航班实际进场的燃油消耗量采用的燃油消耗计算模型计算公式如下:
Figure BDA0001360347100000056
式中,
Figure BDA0001360347100000057
表示燃油消耗量,p表示机型,f表示进场航班,r表示执行的进场程序,a表示平飞航段数量,b表示下降航段数量,Ne表示发动机数量,
Figure BDA0001360347100000058
表示第i个平飞航段的飞行速度,
Figure BDA0001360347100000059
表示第i平飞航段的飞行距离,
Figure BDA00013603471000000510
表示第i平飞航段的燃油消耗率,
Figure BDA00013603471000000511
表示第j个下降航段的飞行速度,
Figure BDA00013603471000000512
表示第j个下降航段的飞行距离,
Figure BDA00013603471000000513
表示第j个下降航段的燃油消耗率。
进一步的,因航班延误产生的时间成本为:
Figure BDA00013603471000000514
式中,
Figure BDA00013603471000000515
表示因航班延误产生的时间成本,f表示航班类型,Ct表示单位时间延误成本,Tf表示延误时间,
Figure BDA00013603471000000516
表示机组小时费,
Figure BDA00013603471000000517
表示维修费,
Figure BDA00013603471000000518
表示飞机和发动机折旧费,
Figure BDA00013603471000000519
表示其它费用;
Figure BDA00013603471000000520
表示实际进场时间,
Figure BDA00013603471000000521
表示计划航班执行进场程序的时间,
Figure BDA00013603471000000522
表示航班空中等待时间。
进一步的,机组小时费包括机长、副驾、安全员、乘务长和乘务员的小时工资;飞机和发动机的折旧费按照飞机价格除以飞机使用小时数得到。
进一步的,其它费用包括根据延误时间的长短酌情给予旅客和/或机组小组的食宿费用。
本发明提出的基于BADA数据库(Base of Aircraft Data)燃油消耗率的短期天气下进场航班延误成本计算方法,具有以下有益效果:
基于BADA的燃油消耗动力学模型确定的航班延误成本,可为飞行计划制定和航班高效运行提供依据,有助于提高航空公司的航班运行控制决策的科学性。
基于BADA燃油消耗的短期天气影响下的延误成本分析,可作为空中交通流量管理系统的组成,为空管部门在保证航班运行安全性的前提下,实施进场航班改航和排序调度的决策预案提供科学依据。
本发明提出的进场航班运行成本不仅仅包含燃油成本,还有与航空器运行时间相关成本,对于时间成本进行分类,并考虑航班的改航和空中等待的飞行差异,在燃油成本计算模型中分类计算,从而得到航班延误总运行成本。
相对于以往的研究文献对航空器燃油消耗成本多以历史数据驱动建立统计学的数学模型去分析挖掘出不同的航班延误成本影响因素不同,本发明从航空器的发动机耗油率角度可分析不同进场运行航段不同机型的耗油状况,从而为准确航空器耗油成本计算提供更为科学和准确的依据。
附图说明
图1终端空域进离场航路和改航路径
具体实施方式
终端区进场程序由一系列不同飞行状态组合构成,典型的进场及进近航段包括:进场航段、起始进近航段、中间进近航段和最后进近航段。
实施例中公开一种基于BADA燃油消耗率的短期天气下进场航班延误成本计算方法,其包括以下步骤:
步骤1,建立终端空域进场航路结构,获取终端空域范围内的气象雷达回波图像及进场航班计划;
步骤2,利用MATLAB的图像处理工具箱对气象雷达回波图像进行颜色识别(识别天气等级1及以上)与轮廓提取,其图像颜色与降水等级、回波强度的关系如表1所示;
表1危险天气等级
Figure BDA0001360347100000071
步骤3,确定受天气影响的进场航路,规划与天气状况对应的改航航路;
实施步骤1~3可得到终端空域进离场航路、和受天气影响的进场航路及其改航航路,如图1。
如图1所示,AND方向进场航路为BK-BAVIK-IGLIT-XSY-PVG,航空器受到恶劣天气影响执行实时规划的进场改航航路AND-DADAT-BELOP-A-B-C-XSY-PVG。
步骤4,确定航班的执飞机型,查阅BADA数据库(Base of Aircraft Data)中关于航空器燃油消耗率的性能参数,并利用BADA数据库关于航空器进场平飞和下降阶段的燃油消耗率模型,将进场航路划分为若干个平飞和下降航段,建立进场航班的燃油消耗计算模型。
燃油消耗模型:终端空域进场进近航段由平飞阶段和下降阶段组成,对于某一进场航班f,机型为p,发动机数量为Ne,执行的进场程序r可划分为a个平飞航段、b个下降航段。假设下降、平飞过程中的燃油消耗率恒定,则对应于某一机型的航空器执行进场程序r的燃油消耗量
Figure BDA0001360347100000072
表示为:
Figure BDA0001360347100000073
式中,
Figure BDA0001360347100000074
—平飞阶段第i个航段的长度,nm;
Figure BDA0001360347100000075
—平飞阶段第i个航段的平均真空速,knot;
Figure BDA0001360347100000081
—下降阶段第j个航段的长度,nm;
Figure BDA0001360347100000082
—下降阶段第j个航段的平均真空速,knot;
Figure BDA0001360347100000083
表示第i平飞航段的燃油消耗率;
Figure BDA0001360347100000084
表示第j个下降航段的燃油消耗率。
Figure BDA0001360347100000085
表示计划燃油消耗量,
Figure BDA0001360347100000086
表示实际燃油消耗量
根据航班计划,AND方向进场航班执飞机型为A320-212、A320-212、A321、A330-301、A340-313、B737-800、B757-200、B767-300ER、B777-300,其性能参数如表2。
表2机型性能参数
Figure BDA0001360347100000087
表1中关于航空器燃油消耗率的性能参数中:CTc,1表示第一最大爬升推力系数(单位:牛顿),CTc,2表示第二最大爬升推力系数(单位:英尺),CTc,3表示第三最大爬升推力系数(单位:1/英尺2),CTc,4表示第一推力温度系数(单位:C);CTc,5表示第二推力温度系数(单位:1/C),Cf1表示第一推力燃油消耗系数(单位:公斤/分钟╳千牛),Cf2表示第二推力燃油消耗系数(单位:节),Cf3表示第一下降燃油流量系数(单位:公斤/分钟),Cf4表示第二下降燃油流量系数(单位:英尺),Cfcr表示巡航燃油流量校正系数(单位:无量纲)。
步骤5,确定进场航路的航路点、航段距离高度以及航空器速度,分析航空器在各个航段的飞行姿态,并分别计算计划进场航路和改航航路中各个航段的航空器燃油消耗率;其中飞行姿态包括下降匀速,平飞匀速,下降减速,下降减速,等待匀速5种姿态。
计算得到表2中所列的A320-212、A321、A330-301、A340-313、B737-800、B757-200、B767-300ER、B777-300机型在计划进场航路BK-BAVIK-IGLIT-XSY-PVG和改航航路AND-DADAT-BELOP-A-B-C-XSY-PVG各阶段的燃油消耗率,如表3和表4所示。
表3 AND方向计划进场航路数据
Figure BDA0001360347100000091
表4 AND方向进场改航航路数据
Figure BDA0001360347100000101
步骤6,根据航班进场计划航路与改航航路的航路数据,计算航班延误时间及其延误成本;
进场时间包括航班计划进场时间、航班实际进场时间和航班实际进场时间,分别对应的计算模型如下:
(1)航班计划进场时间
根据目的地机场的标准进场程序,计划航班由终端空域走廊口到降落至目的机场所执行的进场航路是固定的,进场程序r可划分为a个平飞航段、b个下降航段,计划航班执行进场程序r的时间
Figure BDA0001360347100000102
为:
Figure BDA0001360347100000103
式中,
Figure BDA0001360347100000104
表示第i个平飞航段的飞行速度,
Figure BDA0001360347100000105
表示第i平飞航段的飞行距离,
Figure BDA0001360347100000106
表示第j个下降航段的飞行速度,
Figure BDA0001360347100000107
表示第j个下降航段的飞行距离,
(2)航班实际进场时间
由于终端区短期天气条件变化,航班可能实施改航策略引起进场程序发生变化,改航之后的进场程序r′可划分为个c个平飞航段、d个下降航段,相应的实际进场时间
Figure BDA0001360347100000111
为:
Figure BDA0001360347100000112
式中,
Figure BDA0001360347100000113
表示第k个平飞航段的飞行速度,
Figure BDA0001360347100000114
表示第k平飞航段的飞行距离,
Figure BDA0001360347100000115
表示第l个下降航段的飞行速度,
Figure BDA0001360347100000116
表示第l个下降航段的飞行距离,
(3)航班空中等待时间
由于短期天气条件的变化导致跑道进场容量减小,进场服务率降低,使得进场航班不得不在起始进近定位点IAF处进行盘旋等待,假设航班空中等待时间为
Figure BDA0001360347100000117
即航班进场平均延误时间。本文仅讨论近距平行跑道一起一降模式仪表飞行规则下的平均进场延误时间。其他进场运行模式的航班空中等待时间可类推。
变量定义:
T:给定一天中的时间范围(00:00~24:00);
s:讨论的短期天气影响下的时间片长短,取整数,s∈T;
t:第t个时间片,t=1,2,…,24;
Figure BDA0001360347100000118
第t个时间片的机场进场容量;
Figure BDA0001360347100000119
第t个时间片的进场航班需求;
第t个时间片的服务率;
Figure BDA00013603471000001110
第t个时间片的到达率;
Figure BDA00013603471000001111
定义ρ为服务机构的利用率,表示在相同时间片内航班到达率与航班服务率的比值,故进场跑道的利用率分别为
Figure BDA00013603471000001112
第t个时间片的进场航班的平均等待时间,即进场航班延误时间Wt a
Figure BDA0001360347100000121
则s时间片内的进场航班等待时间
Figure BDA0001360347100000122
经计算得到,航班执行计划航路与改航航路所产生的延误时间为19.7min。
相应的,根据对上述三种进场时间的计算,得到因航班延误产生的时间成本计算模型:
对于航班f的延误时间成本
Figure BDA0001360347100000123
等于单位时间延误成本Ct乘以延误时间Tf。单位时间的延误成本包括机组小时费
Figure BDA0001360347100000124
维修费
Figure BDA0001360347100000125
飞机和发动机折旧费
Figure BDA0001360347100000126
其它费用
Figure BDA0001360347100000127
得到航班f的延误时间成本
Figure BDA0001360347100000128
为:
Figure BDA0001360347100000129
式中,Ct表示延误成本,Tf表示延误时间,
Figure BDA00013603471000001210
表示机组小时费,
Figure BDA00013603471000001211
表示维修费,
Figure BDA00013603471000001212
表示飞机和发动机折旧费,
Figure BDA00013603471000001213
表示旅客食宿费用,
Figure BDA00013603471000001214
表示航班执行进场程序的时间,
Figure BDA00013603471000001215
表示实际进场时间,
Figure BDA00013603471000001216
表示空中等待时间。
本实施例中,空中等待时间为0,故:
Figure BDA00013603471000001217
其中,机组小时费包括机长、副驾、安全员、乘务长和乘务员的小时工资;飞机和发动机的折旧费按照飞机价格除以飞机使用小时数得到,飞机使用小时根据飞机寿命15年,每年飞行320天,每天飞行10小时来计算,已知A320-212、A321、A330-301、A340-313、B737-800、B757-200、B767-300ER、B777-300机型价格分别为0.99亿、1.16亿、2.59亿、3.56亿、0.96亿、1.42、1.97亿、3.39亿美元,不同机型的单位时间成本如表5所示,以下价格单位均为人民币元。
表5不同机型的单位时间成本(元/小时)
Figure BDA0001360347100000131
经计算,航班延误时间成本分别为6630、7461、15583、20797、6488、8999、12646、19991元。
步骤7,利用燃油消耗模型计算航班执行计划进场航路和改航航路的燃油消耗量,根据计划与改航的航空器燃油消耗量计算航班延误的燃油成本,如表6所示;
表6 AND进场航班燃油消耗量
Figure BDA0001360347100000132
因航班延误产生的燃油成本计算模型为:
对于航班f的延误燃油成本
Figure BDA0001360347100000133
等于单位燃油价格Ce乘以燃油消耗量Ef,即:
Figure BDA0001360347100000134
式中,f表示进场航班,Ce表示单位燃油价格,Ef表示燃油消耗量,
Figure BDA0001360347100000135
表示航班计划进场燃油消耗量,
Figure BDA0001360347100000136
表示航班实际进场燃油消耗量,
Figure BDA0001360347100000137
表示航班空中等待燃油消耗量;
Figure BDA0001360347100000138
式中,
Figure BDA0001360347100000139
表示航班空中等待燃油消耗量,Ne表示发动机数量,
Figure BDA00013603471000001310
表示航班空中等待时间,fcr表示航班空中等待阶段的燃油消耗率。
本实施例中,空中等待时间为0,故:
Figure BDA0001360347100000141
假设单位燃油成本Ce为3000元/吨,根据延误燃油成本计算模型计算得到各种机型的航班延误燃油成本分别为:3480、4395、8610、17755、3783、5612、8036、8871元。
步骤8,计算航班延误总成本,计算模型如下:
终端空域进场航班的延误总成本CDelay为所有进场航班延误成本之和,包括每个航班的延误时间成本和延误燃油成本;而且航班f有且仅能由一种机型p,执行进场程序r降落于机场跑道。
Figure BDA0001360347100000142
Figure BDA0001360347100000143
Figure BDA0001360347100000144
式中,F表示终端区进场航班的集合,
Figure BDA0001360347100000145
P表示终端区进场航班的航空器机型集合,
Figure BDA0001360347100000146
R表示终端区进场进近程序的集合,
Figure BDA0001360347100000147
Figure BDA0001360347100000148
表示航班f的状态变量,为0-1离散变量;
Figure BDA0001360347100000149
表示航班f的延误成本;
Figure BDA00013603471000001410
表示航班f的延误时间成本;
Figure BDA00013603471000001411
表示航班f的延误燃油成本。
经计算,因航班延误带来的时间成本、燃油成本以及总成本如表7所示。
表7航班延误时间成本、燃油成本以及总成本
Figure BDA00013603471000001412
尽管以上结合附图对本发明的实施方案进行了描述,但本发明并不局限于上述的具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的、指导性的,而不是限制性的。本领域的普通技术人员在本说明书的启示下,在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式,这些均属于本发明保护之列。

Claims (1)

1.一种基于BADA燃油消耗率的短期天气下进场航班延误成本计算方法,其特征在于,包括以下步骤:
步骤1,建立终端空域进场航路结构,获取终端空域范围内的气象雷达回波图像及进场航班计划;
步骤2,基于图像颜色与降水等级、回波强度的关系,对气象雷达回波图像进行颜色识别与轮廓提取;
步骤3,确定受天气影响的计划进场航路,规划与天气状况对应的改航航路;
步骤4,确定航班的执飞机型,并查阅BADA数据库关于航空器燃油消耗率的性能参数,并利用BADA数据库关于航空器进场平飞和下降阶段的燃油消耗率模型,将计划进场航路划分为若干个平飞和下降航段,建立进场航班的燃油消耗计算模型;
步骤5,确定进场航路的航路点、航段距离高度以及航空器速度,分析航空器在各个航段的飞行姿态,并计算计划进场航路和改航航路各个航段的航空器燃油消耗量;
步骤6,根据航班计划进场航路与改航航路的航路数据,计算因航班延误产生的延误时间及因航班延误产生的时间成本;
步骤7,利用燃油消耗计算模型计算航班执行计划进场航路和改航航路的燃油消耗量,根据所述航空器燃油消耗量计算因航班延误产生的燃油成本;
步骤8,计算短期天气下进场航班延误总成本,即所有航空器执飞进场和进近程序时因航班延误产生的时间成本和燃油成本之和;
其中:各终端区进场航班有且仅能由一种航空器机型p,执行进场程序r降落于机场跑道,进场航班延误总成本的计算公式为:
Figure FDA0002940393000000011
Figure FDA0002940393000000012
Figure FDA0002940393000000013
式中,CDelay表示延误总成本,F表示终端区进场航班f的集合,P表示终端区进场航班的航空器机型p的集合,R表示终端区执行进场程序r的集合,
Figure FDA0002940393000000021
表示终端区进场航班f的状态变量,为0-1离散变量,
Figure FDA0002940393000000022
表示终端区进场航班f的延误成本,
Figure FDA0002940393000000023
表示终端区进场航班f的延误时间成本,
Figure FDA0002940393000000024
表示终端区进场航班f的延误燃油成本;
终端区进场航班f的延误燃油成本
Figure FDA0002940393000000025
的计算公式如下:
Figure FDA0002940393000000026
式中,Ce表示单位燃油价格,Ef表示终端区进场航班延误产生的燃油消耗量,
Figure FDA0002940393000000027
表示航班计划进场航路的燃油消耗量,
Figure FDA0002940393000000028
表示航班改航航路的燃油消耗量,
Figure FDA0002940393000000029
表示航班空中等待的燃油消耗量;
航班空中等待的燃油消耗量
Figure FDA00029403930000000210
的计算公式如下:
Figure FDA00029403930000000211
式中,Ne表示发动机数量,
Figure FDA00029403930000000212
表示航班空中等待时间,fcr表示航班空中等待阶段的燃油消耗率;
其中:航班执行计划进场航路和航班改航航路各个航段的航空器燃油消耗量
Figure FDA00029403930000000213
的计算公式如下:
Figure FDA00029403930000000214
式中,a表示平飞航段数量,b表示下降航段数量,
Figure FDA00029403930000000215
表示第i个平飞航段的飞行速度,
Figure FDA00029403930000000216
表示第i平飞航段的飞行距离,
Figure FDA00029403930000000217
表示第i平飞航段的燃油消耗率,
Figure FDA00029403930000000218
表示第j个下降航段的飞行速度,
Figure FDA00029403930000000219
表示第j个下降航段的飞行距离,
Figure FDA00029403930000000220
表示第j个下降航段的燃油消耗率;
其中:因航班延误产生的时间成本
Figure FDA00029403930000000221
为:
Figure FDA0002940393000000031
式中,Ct表示单位时间延误成本,Tf表示延误时间,
Figure FDA0002940393000000032
表示机组小时费,
Figure FDA0002940393000000033
表示维修费,
Figure FDA0002940393000000034
表示飞机和发动机折旧费,
Figure FDA0002940393000000035
表示食宿费用,
Figure FDA0002940393000000036
表示航班实际进场时间,
Figure FDA0002940393000000037
表示航班执行进场程序的时间,
Figure FDA0002940393000000038
表示航班空中等待时间。
CN201710614866.5A 2017-07-26 2017-07-26 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法 Active CN107341620B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710614866.5A CN107341620B (zh) 2017-07-26 2017-07-26 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710614866.5A CN107341620B (zh) 2017-07-26 2017-07-26 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法

Publications (2)

Publication Number Publication Date
CN107341620A CN107341620A (zh) 2017-11-10
CN107341620B true CN107341620B (zh) 2021-04-23

Family

ID=60216137

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710614866.5A Active CN107341620B (zh) 2017-07-26 2017-07-26 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法

Country Status (1)

Country Link
CN (1) CN107341620B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111143972A (zh) * 2019-12-05 2020-05-12 南京航空航天大学 考虑运行环境限制的航空器场面滑行油耗和排放评估方法
CN111507642B (zh) * 2020-04-28 2023-08-29 中国商用飞机有限责任公司北京民用飞机技术研究中心 一种基于航迹运行的经济性评估方法及装置
CN111696389B (zh) * 2020-05-29 2021-07-02 航科院中宇(北京)新技术发展有限公司 一种基于航班飞行计划的航空器燃油估计方法及系统
US11694556B2 (en) 2022-04-11 2023-07-04 The 28Th Research Institute Of China Electronics Technology Group Corporation Time-space conversion method of flight sequencing information
CN114446094B (zh) * 2022-04-11 2022-06-17 中国电子科技集团公司第二十八研究所 一种航班排序信息的时空转换方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692315A (zh) * 2009-09-25 2010-04-07 民航总局空管局技术中心 基于实时雷达数据的高精度4d飞机航迹分析方法
CN103810907A (zh) * 2012-11-07 2014-05-21 霍尼韦尔国际公司 用于在恶劣天气下着陆期间增强飞行员决策的系统和方法
CN106768123A (zh) * 2016-11-30 2017-05-31 中国直升机设计研究所 一种无人直升机燃油预估方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8645009B2 (en) * 2012-02-23 2014-02-04 Ge Aviation Systems Llc Method for flying an aircraft along a flight path

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101692315A (zh) * 2009-09-25 2010-04-07 民航总局空管局技术中心 基于实时雷达数据的高精度4d飞机航迹分析方法
CN103810907A (zh) * 2012-11-07 2014-05-21 霍尼韦尔国际公司 用于在恶劣天气下着陆期间增强飞行员决策的系统和方法
CN106768123A (zh) * 2016-11-30 2017-05-31 中国直升机设计研究所 一种无人直升机燃油预估方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于空管数据的飞行器终端区域能耗分析;刘钊;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20170315;第三章3.1节,第四章4.2节及第五章5.1节至5.3节 *
恶劣天气下终端区航迹智能优化;杨超;《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》;20170315;第二章2.1节至2.5节,第三章3.1节至3.3节以及第四章4.1节至4.3节 *

Also Published As

Publication number Publication date
CN107341620A (zh) 2017-11-10

Similar Documents

Publication Publication Date Title
CN107341620B (zh) 基于bada燃油消耗率的短期天气下进场航班延误成本计算方法
CN110689764B (zh) 一种基于动态模拟的飞机离场放行排序方法
CN110335506A (zh) 一种终端区短时空域利用率计算方法
Çeçen et al. Aircraft sequencing and scheduling in TMAs under wind direction uncertainties
CN110264787A (zh) 一种航班航路飞行时间可靠性评价方法及系统
Carr et al. Airline arrival prioritization in sequencing and scheduling
Gilbert Historical development of the air traffic control system
Bell Operational research into air traffic control
CN109902924A (zh) 一种军用机场单跑道训练容量评估方法
Malik et al. Runway scheduling for Charlotte Douglas international airport
Kim et al. Runway assignment by minimizing emissions in terminal airspace
Hartjes et al. Evaluation of intermediate stop operations in long-haul flights
Kaplan et al. A multi-objective nonlinear integer programming model for mixed runway operations within the TMAs
Andreeva-Mori et al. Joint Strategic and Tactical Air Traffic Management of Arrival Flights at Tokyo International Airport for Reduced Climate Impact
DeArmon et al. Metering during severe en route weather via Advanced Flight-Specific Trajectories (AFST)
Sprong Measuring the effects of RNAV departure procedures on airport efficiency
Putra Review Study on Runway Capacity Parameters and Improvement
Jiawei Airport Intelligent Dispatching Model Based on Double Queue and Tabu Search
Wieland Safety and performance implications of wake vortex buffers
Donahue The US Air Transportation System: A Bold Vision for Change
Rome et al. Ripple delay and its mitigation
Price Air traffic control: Its effect on fuel conservation
Wang et al. Research on Identification of Risk Factors in Civil Airports
Yıldırım PERFORMANCE CALCULATION: MIDDLE-CLASS AIRCRAFT
FEDERAL AVIATION ADMINISTRATION WASHINGTON DC OFFICE OF AVIATION POLICY A Proposed Aviation Energy Conservation Program for the National Aviation System. Volume I. The Short Run, 1977-1978.

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant