CN109738035A - 基于ads-b航迹数据的航空器油耗计算方法 - Google Patents
基于ads-b航迹数据的航空器油耗计算方法 Download PDFInfo
- Publication number
- CN109738035A CN109738035A CN201910057201.8A CN201910057201A CN109738035A CN 109738035 A CN109738035 A CN 109738035A CN 201910057201 A CN201910057201 A CN 201910057201A CN 109738035 A CN109738035 A CN 109738035A
- Authority
- CN
- China
- Prior art keywords
- stage
- lto
- flight
- aircraft
- tim
- 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.)
- Pending
Links
Landscapes
- Traffic Control Systems (AREA)
Abstract
本发明涉及一种基于ADS‑B航迹数据的航空器油耗计算方法包括:获取数据信息,存入数据库;通过航班时刻信息和ADS‑B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;依据数据信息确定每个阶段的燃油流量;依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。可计算航班完整运行过程的总油耗,并且可根据阶段划分算出各个阶段内的油耗,以便于阶段性分析与统计方面的需要。
Description
技术领域
本发明涉及航空器运行状态识别与性能分析领域,具体涉及一种基于ADS-B 航迹数据的航空器油耗计算方法。
背景技术
随着全球空中交通活动保持持续快速增长,航空运输已经成为最便捷的交通出行方式之一,但同时也对全球和地区的环境造成了严重影响,如果地球生态系统和气候系统突破了某个的临界点,则会导致不可逆转的气候变化,例如高温、雾霾、大雾和强雷暴等极端天气的频发。而这些极端天气又会造成大面积航班延误和取消,反过来又会影响和制约到航空运输业的发展。因此,航空运行节能减排成为推进绿色民航发展的重要内容。同时,降低航空器排放量也意味着减少燃油消耗,对于航空公司来说降低运营成本是至关重要的。另外各航空公司由于商业秘密的考虑,不会向外提供航班的载重与油耗等数据,因此在进行航班油耗分析与统计之前需要有一个油耗估算方法。
当前国内外对油耗估算方法主要有三类:(1)基于航班计划数据建立航空器性能仿真计算平台;(2)仅针对机场地区使用国际民航组织(ICAO, International CivilAviation Organization)的起飞着陆循环(LTO,Landing and Take-off Cycle)方法;(3)基于航空器QAR(QUick Access Record)监视数据建立神经网络模型对空中油耗进行估算。这些方法主要有两个方面的限制:
仿真模型中假设航空器按照理想航迹进行飞行,而航空器实际飞行轨迹与理想航迹之间存在一定的差距;
LTO方法仅适用于机场上空3000英尺范围内的油耗计算,而对3000英尺以上的航路飞行则无法适用;
神经网络方法估算的精度取决于学习库的范围,因此需要建立完整而庞大的学习库才能获得最佳效果。
如何解决上述问题,是目前亟待解决的。
发明内容
本发明的目的是提供一种基于ADS-B航迹数据的航空器油耗计算方法。
为了解决上述技术问题,本发明提供了一种基于ADS-B航迹数据的航空器油耗计算方法包括:
获取数据信息,存入数据库;
通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;
依据数据信息确定每个阶段的燃油流量;
依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;
依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;
根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。
本发明的有益效果是,本发明提供了一种基于ADS-B航迹数据的航空器油耗计算方法包括:获取数据信息,存入数据库;通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;依据数据信息确定每个阶段的燃油流量;依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。可计算航班完整运行过程的总油耗,并且可根据阶段划分算出各个阶段内的油耗,以便于阶段性分析与统计方面的需要。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1是本发明所提供的基于ADS-B航迹数据的航空器油耗计算方法的流程图。
图2是本发明所提供的航班历史时刻与ADS-B飞行轨迹数据的主要字段示意图。
图3是本发明所提供的航班完整运行阶段划分示意图。
图4是本发明所提供的航路上地速GS、风速WS和真空速TAS矢量关系示意图。
具体实施方式
现在结合附图对本发明作进一步详细的说明。这些附图均为简化的示意图,仅以示意方式说明本发明的基本结构,因此其仅显示与本发明有关的构成。
实施例1
请参阅图1为本实施例1提供了一种基于ADS-B航迹数据的航空器油耗计算方法。基于ADS-B航迹数据的航空器油耗计算方法可计算航班完整运行过程的总油耗,并且可根据阶段划分算出各个阶段内的油耗,以便于阶段性分析与统计方面的需要。基于ADS-B航迹数据的航空器油耗计算方法包括:
S110:获取数据信息,存入数据库;
S120:通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;
S130:依据数据信息确定每个阶段的燃油流量;
S140:依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;
S150:依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;
S160:根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。
在本实施例中,步骤S110包括:
准备静态数据,存入数据库;
其中,静态数据包括机场信息资料、BADA机型性能数据、ICAO EEDB发动机排放数据和高空航路气象观测资料,机场信息资料提供机场的地理位置信息和标高,BADA提供机型性能数据,EEDB提供机型发动机在7%、30%、85%和100%四种推力调定值的燃油流量、HC排放指数、CO排放指数和NOx排放指数;
建立BADA机型与EEDB发动机型号的关联表;如下表给出了几种常见机型与发动机型号的关联信息。
准备动态数据,存入数据库;
其中,动态数据包括航班历史时刻信息和ADS-B飞行轨迹数据。这两种数据通过航班唯一标识码FlightID进行关联,即通过航班FlightID对该航班航空器的ADS-B飞行轨迹数据进行读取。动态数据的具体字段内容如图2所示。
对该航班的ADS-B飞行轨迹数据进行去噪处理,以保证后续步骤的稳定执行。
请参阅图3,在本实施例中,步骤S120包括:
起飞机场DEP与目标机场ARR的上空3000英尺,气压高度分别为:AGLDEP,3000和AGLARR,3000,其中,AGL_DEP,3000意思是起飞机场上空3000英尺所对应的气压高度值,同理,AGL_ARR,3000意思是目标机场上空3000英尺所对应的气压高度值范围之内为LTO阶段,其他区域为CCD,即:
航班的ADS-B飞行轨迹记录点{A0,A1,A2,A3,…An}归属阶段LTODEP、LTOARR和CCD,其中,LTODEP、LTOARR分别为起飞机场范围的LTO阶段和目标机场的LTO 阶段划分的具体过程为:
从A0到An进行增序循环,对每一个轨迹点Ai,如果满足以下条件,则Ai归属为起飞机场的LTODEP阶段:
其中,MTime_i表示第i个ADS-B记录点上的监视时刻,HP_i表示第i个 ADS-B记录点上的气压高度,ATOT表示实际起飞时间;
从An到A0进行降序循环,对每一个轨迹点Ai,如果满足以下条件,则Ai归属为目标机场的LTOARR阶段:
其中,ALDT表示实际落地时间;
不属于LTODEP和LTOARR阶段的轨迹点,则归属于CCD阶段。
在本实施例中,步骤S130包括:
获取数据库中航班基本信息、航班关键时刻、航空器基本信息和航班ADS-B 飞行轨迹数据;
根据航空器机型与EEDB发动机型号关联信息将该航班所执飞航空器的机型与发动机型号进行关联,以确定该航空器在滑行模式、起飞模式、爬出模式和进近着陆模式下的燃油流量。
在本实施例中,S140包括:
计算LTO阶段内的爬出模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn;
依据模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn以及LTO阶段的燃油流量计算LTO阶段的燃油消耗量。
其中,所述计算LTO阶段内的爬出模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn的方法包括:
依据最后一个起飞机场LTO阶段的轨迹点ADEP,Last和第一个CCD阶段的轨迹点ACCD,First,使用线性插值方法计算出航空器在起飞机场上空3000英尺高度时的时间MTimeDEP,AGL3000,然后结合ATOT和起飞所花时间TIMTakeoff,计算出航空器在起飞机场LTO阶段内的爬出过程所花时间TIMClimhout,即:
TIMClimbout=MTimeDEP,AGL3000-ATOT-TIMTakeoff;
依据第一个目标机场LTO阶段的轨迹点AARR,First和最后一个CCD阶段的轨迹点ACCD,Last,使用线性插值方法计算出航空器在目标机场上空3000英尺高度时的时间MTimeARR,AGL3000,然后结合ALDT计算出航空器在目标机场LTO阶段内的进近过程所花时间TIMApproach,即:
TIMApproach=ALDT-MTimeARR,3000;
滑出模式和滑入模式时间可直接由AOBT、ATOT、ALDT和AIBT得到,即:
TIMTaxiOut=ATOT-AOBT,TIMTaxiIn=AIBT-ALDT;
其中,AOBT,表示实际撤轮档时间,ATOT,表示实际起飞时间,ALDT,表示实际落地时间,AIBT,表示实际挡轮档时间,依据航班历史时刻信息得出。
其中,所述依据模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn以及LTO阶段的燃油流量计算LTO阶段的燃油消耗量的,即:
依据ICAO EEDB中的滑行模式、起飞模式、爬出模式和进近着陆模式中发动机燃油消耗率FFTaxi、FFTakeoff、FFClimbout和FFApproach与计算得出的TIMTaxiOut、TIMTaxiIn、 TIMClimbout和TTMApproach,以及起飞所用时间TTMTakeoff,计算得到航空器在这些过程中的燃油消耗量TFTaxi、TFTakeoff、TFClimbout和TFApproach。最后可得到LTO阶段燃油消耗量TFDEP、TFARR和TFLTO,即:
TFTaxiOut=FFTaxi·TIMTaxiOut,TFTaxiIn=FFTaxi·TIMTaxiIn
TFTakeoff=FFTakeoff·TIMTakeoff,TFClimbout=FFClimbout·TIMClimbout
TFApproach=FFApproach·TIMApproach
TFDEP=TFTaxiOut+TFClimbout+TFTakeoff,TFARR=TFApproach+TFTaxiIn
TFLTO=TFDEP+TFARR。
在本实施例中,步骤S150包括:
根据CCD阶段内的航空器ADS-B飞行轨迹数据中的飞行高度及其变化对记录点Ai进行飞行状态识别,飞行状态分为飞行过程PHSi和气动构型CFGi;
依据航路高空气象观测资料数据中读取对应时段内航空器ADS-B飞行轨迹中每一个记录点Ai所在位置和所在飞行高度上的风速WSi和风向WDi结合航空器的地速GSi和磁航向MHi,结合航空器的地速GSi和磁航向MHi计算出真空速TASi;
对CCD阶段内每一个ADS-B飞行轨迹记录点Ai,根据飞行状态飞行过程PHSi、气动构型CFGi和真空速TASi,利用BADA性能模型(AircraftPermance Model, APM)和对应机型数据计算出航空器在记录点Ai时的升力系数CL,i、阻力系数CD,i、阻力Di和推力Ti,最后由Ti结合PHSi和CFGi由发动机性能模型得到燃油流量FFi;
计算Ai到Ai+1航空器的耗油量为TFi,再求出整个CCD阶段内航空器的燃油消耗量TFCCD。
其中,根据CCD阶段内的航空器ADS-B飞行轨迹数据中的飞行高度及其变化对记录点Ai进行飞行状态识别,飞行状态分为飞行过程PHSi和气动构型CFGi,即:
其中飞行过程PHSi分为上升CLIMB、下降DESCENT和平飞LEVEL三种,气动构型CFGi分为光滑构型CLEAN和进近构型APPROACH;
对记录点Ai,其飞行过程PHSi和气动构型CFGi的识别过程为:
其中,Ri,i+1为Ai到Ai+1的高度变化率,HPi+1和HPi为Ai和Ai+1时的飞行高度, RCmin和RDmin分别为最小上升率和最小下降率设定值,且均为正数。
请参阅图4,在本实施例中,依据航路高空气象观测资料数据中读取对应时段内航空器ADS-B飞行轨迹中每一个记录点Ai所在位置和所在飞行高度上的风速WSi和风向WDi(只取水平风,忽略垂直风),结合航空器的地速GSi和磁航向 MHi,结合航空器的地速GSi和磁航向MHi计算出真空速TASi的具体计算过程为:
TAi=|WDi-MHi|,DAi=|GAi-MHi|
WAi=180-DAi-(180-TAi)=TAi-DAi
并且,如果TAi≥180,则TAi=|WDi-MHi|-180。
在本实施例中,所述对CCD阶段内每一个ADS-B飞行轨迹记录点Ai,根据飞行状态飞行过程PHSi、气动构型CFGi和真空速TASi,利用BADA性能模型 (AircraftPermanceModel,APM)和对应机型数据计算出航空器在记录点Ai 时的升力系数CL,i、阻力系数CD,i、阻力Di和推力Ti,最后由Ti结合PHSi和CFGi由发动机性能模型得到燃油流量FFi的方法,即:
FFCR,i=Cfcr·FFnom,i;
其中Ci,i+1为Ai到Ai+1的TAS变化率,CD0,AP,CD2,AP,CD0,AR和CD2,CR为航空器机型阻力系数参数,mi为Ai上航空器重量,ρi为Ai高度上的大气密度,S为航空器机翼总面积,g为重力加速度,Cf1、Cf2、Cf3、Cf4和Cfcr为燃油流量系数。
在本实施例中,步骤S160包括:
据出的LTO阶段耗油量和CCD阶段油耗量,可计算出航班整个运行过程的耗油量TF:TF=TFLTO+TFCCD。
在本实施例中,将每个阶段的燃油流量替换为HC排放指数、CO排放指数和 NOx排放指数,即可计算整个运行过程中的HC、CO、NO2的排放含量;便于为大规模航班的污染物排放的统计分析提供数据支持和技术参考。
通过上述实施方式对2018年8月1日南京禄口机场飞往中国香港国际机场的某航班燃油消耗进行了计算,与实际燃油消耗总量相比,燃油计算结果的相对误差仅为1.11%,通过实验证明,本发明具备较高的精度,可为航班污燃油消耗与染物排放方面相关研究提供参考。
综上所述,本发明提供了一种基于ADS-B航迹数据的航空器油耗计算方法包括:获取数据信息,存入数据库;通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;依据数据信息确定每个阶段的燃油流量;依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。可计算航班完整运行过程的总油耗,并且可根据阶段划分算出各个阶段内的油耗,以便于阶段性分析与统计方面的需要。
以上述依据本发明的理想实施例为启示,通过上述的说明内容,相关工作人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容,必须要根据权利要求范围来确定其技术性范围。
Claims (10)
1.一种基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,包括:
获取数据信息,存入数据库;
通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段;
依据数据信息确定每个阶段的燃油流量;
依据LTO阶段的燃油流量计算LTO阶段燃油消耗量;
依据CCD阶段的燃油流量计算CCD阶段燃油消耗量;
根据LTO阶段和CCD阶段的燃油流量,计算出航班整个运行过程的燃油消耗量TF。
2.如权利要求1所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述获取数据信息,存入数据库的方法包括:
准备静态数据,存入数据库;
其中,静态数据包括机场信息资料、BADA机型性能数据、ICAO EEDB发动机排放数据和高空航路气象观测资料,机场信息资料提供机场的地理位置信息和标高,BADA提供机型性能数据,EEDB提供机型发动机在7%、30%、85%和100%四种推力调定值的燃油流量、HC排放指数、CO排放指数和NOx排放指数;
建立BADA机型与EEDB发动机型号的关联表;
准备动态数据,存入数据库;
其中,动态数据包括航班历史时刻信息和ADS-B飞行轨迹数据。
3.如权利要求2所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述通过航班时刻信息和ADS-B飞行轨迹记录点的高度信息将航班运行分为两个阶段:LTO阶段和CCD阶段的方法包括:
起飞机场DEP与目标机场ARR的上空3000英尺,气压高度分别为:AGLDEP,3000和AGLARR,3000,其中,AGL_DEP,3000意思是起飞机场上空3000英尺所对应的气压高度值,同理,AGL_ARR,3000意思是目标机场上空3000英尺所对应的气压高度值范围之内为LTO阶段,其他区域为CCD,即:
航班的ADS-B飞行轨迹记录点{A0,A1,A2,A3,…An}归属阶段LTODEP、LTOARR和CCD,其中,LTODEP、LTOARR分别为起飞机场范围的LTO阶段和目标机场的LTO阶段划分的具体过程为:
从A0到An进行增序循环,对每一个轨迹点Ai,如果满足以下条件,则Ai归属为起飞机场的LTODEP阶段:
其中,MTime_i表示第i个ADS-B记录点上的监视时刻,HP_i表示第i个ADS-B记录点上的气压高度,ATOT表示实际起飞时间;
从An到A0进行降序循环,对每一个轨迹点Ai,如果满足以下条件,则Ai归属为目标机场的LTOARR阶段:
其中,ALDT表示实际落地时间;
不属于LTODEP和LTOARR阶段的轨迹点,则归属于CCD阶段。
4.如权利要求3所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述依据数据信息确定每个阶段的燃油流量的方法包括:
获取数据库中航班基本信息、航班关键时刻、航空器基本信息和航班ADS-B飞行轨迹数据;
根据航空器机型与EEDB发动机型号关联信息将该航班所执飞航空器的机型与发动机型号进行关联,以确定该航空器在滑行模式、起飞模式、爬出模式和进近着陆模式下的燃油流量。
5.如权利要求4所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述依据LTO阶段的燃油流量计算LTO阶段燃油消耗量的方法包括:
计算LTO阶段内的爬出模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn;
依据模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn以及LTO阶段的燃油流量计算LTO阶段的燃油消耗量。
6.如权利要求5所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述计算LTO阶段内的爬出模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn的方法包括:
依据最后一个起飞机场LTO阶段的轨迹点ADEP,Last和第一个CCD阶段的轨迹点ACCD,First,使用线性插值方法计算出航空器在起飞机场上空3000英尺高度时的时间MTimeDEP,AGL3000,然后结合ATOT和起飞所花时间TIMTakeoff,计算出航空器在起飞机场LTO阶段内的爬出过程所花时间TIMClimbout,即:
TIMClimbout=MTimeDEP,AGL3000-ATOT-TIMTakeoff;
依据第一个目标机场LTO阶段的轨迹点AARR,First和最后一个CCD阶段的轨迹点ACCD,Last,使用线性插值方法计算出航空器在目标机场上空3000英尺高度时的时间MTimeARR,AGL3000,然后结合ALDT计算出航空器在目标机场LTO阶段内的进近过程所花时间TIMApproach,即:
TIMApproach=ALDT-MTimeARR,3000;
滑出模式和滑入模式时间可直接由AOBT、ATOT、ALDT和AIBT得到,即:
TIMTaxiOut=ATOT-AOBT,TIMTaxiIn=AIBT-ALDT;
其中,AOBT,表示实际撤轮档时间,ATOT,表示实际起飞时间,ALDT,表示实际落地时间,AIBT,表示实际挡轮档时间,依据航班历史时刻信息得出。
7.如权利要求6所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述依据模式时间TIMClimbout、进近着陆模式时间TIMApproach、滑出模式时间TIMTaxiOut和滑入模式时间TIMTaxiIn以及LTO阶段的燃油流量计算LTO阶段的燃油消耗量的,即
依据ICAO EEDB中的滑行模式、起飞模式、爬出模式和进近着陆模式中发动机燃油消耗率FFTaxi、FFTakeoff、FFClimbout和FFApproach与计算得出的TIMTaxiOut、TIMTaxiIn、TIMClimbout和TIMApproach,以及起飞所用时间TIMTakeoff,计算得到航空器在这些过程中的燃油消耗量TFTaxi、TFTakeoff、TFClimbout和TFApproach;最后可得到LTO阶段燃油消耗量TFDEP、TFARR和TFLT0,即:
TFTaxxiOut=FFTaxi·TIMTaxiOut,TFTaxiIn=FFTaxi·TIMTaxiIn
TFTakeoff=FFTakeoff·TIMTakeoff,TFClimbout=FFClimbout·TIMClimbout
TFApproach=FFApproach·TIMApproach
TFDEP=TFTaxiOut+TFClimbout+TFTakeoff,TFARR=TFApproach+TFTaxiIn
TFLTO=TFDEP+TFARR。
8.如权利要求7所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述依据CCD阶段的燃油流量计算CCD阶段燃油消耗量的方法包括:
根据CCD阶段内的航空器ADS-B飞行轨迹数据中的飞行高度及其变化对记录点Ai进行飞行状态识别,飞行状态分为飞行过程PHSi和气动构型CFGi;
依据航路高空气象观测资料数据中读取对应时段内航空器ADS-B飞行轨迹中每一个记录点Ai所在位置和所在飞行高度上的风速WSi和风向WDi结合航空器的地速GSi和磁航向MHi,结合航空器的地速GSi和磁航向MHi计算出真空速TASi;
对CCD阶段内每一个ADS-B飞行轨迹记录点Ai,根据飞行状态飞行过程PHSi、气动构型CFGi和真空速TASi,利用BADA性能模型(AircraftPermance Model,APM)和对应机型数据计算出航空器在记录点Ai时的升力系数CL,i、阻力系数CD,i、阻力Di和推力Ti,最后由Ti结合PHSi和CFGi由发动机性能模型得到燃油流量FFi;
计算Ai到Ai+1航空器的耗油量为TFi,再求出整个CCD阶段内航空器的燃油消耗量TFCCD。
9.如权利要求8所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,根据CCD阶段内的航空器ADS-B飞行轨迹数据中的飞行高度及其变化对记录点Ai进行飞行状态识别,飞行状态分为飞行过程PHSi和气动构型CFGi,即:
其中飞行过程PHSi分为上升CLIMB、下降DESCENT和平飞LEVEL三种,气动构型CFGi分为光滑构型CLEAN和进近构型APPROACH;
对记录点Ai,其飞行过程PHSi和气动构型CFGi的识别过程为:
其中,Ri,i+1为Ai到Ai+1的高度变化率,HPi+1和HPi为Ai和Ai+1时的飞行高度,RCmin和RDmin分别为最小上升率和最小下降率设定值,且均为正数。
10.如权利要求9所述的基于ADS-B航迹数据的航空器油耗计算方法,其特征在于,所述对CCD阶段内每一个ADS-B飞行轨迹记录点Ai,根据飞行状态飞行过程PHSi、气动构型CFGi和真空速TASi,利用BADA性能模型(AircraftPermance Model,APM)和对应机型数据计算出航空器在记录点Ai时的升力系数CL,i、阻力系数CD,i、阻力Di和推力Ti,最后由Ti结合PHSi和CFGi由发动机性能模型得到燃油流量FFi的方法,即:
FFCR,i=Cfcr·FFnom,i;
其中Ci,i+1为Ai到Ai+1的TAS变化率,CD0,AP,CD2,AP,CD0,AR和CD2,CR为航空器机型阻力系数参数,mi为Ai上航空器重量,ρi为Ai高度上的大气密度,S为航空器机翼总面积,g为重力加速度,Cf1、Cf2、Cf3、Cf4和Cfcr为燃油流量系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057201.8A CN109738035A (zh) | 2019-01-21 | 2019-01-21 | 基于ads-b航迹数据的航空器油耗计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910057201.8A CN109738035A (zh) | 2019-01-21 | 2019-01-21 | 基于ads-b航迹数据的航空器油耗计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109738035A true CN109738035A (zh) | 2019-05-10 |
Family
ID=66365583
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910057201.8A Pending CN109738035A (zh) | 2019-01-21 | 2019-01-21 | 基于ads-b航迹数据的航空器油耗计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109738035A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110503349A (zh) * | 2019-09-10 | 2019-11-26 | 深圳市瑞达飞行科技有限公司 | 航班油耗数据的处理方法、处理装置和电子设备 |
CN111507642A (zh) * | 2020-04-28 | 2020-08-07 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种基于航迹运行的经济性评估方法及装置 |
CN111581801A (zh) * | 2020-04-29 | 2020-08-25 | 南京航空航天大学 | 一种航空器巡航温室效应计算方法 |
CN111696389A (zh) * | 2020-05-29 | 2020-09-22 | 航科院中宇(北京)新技术发展有限公司 | 一种基于航班飞行计划的航空器燃油估计方法及系统 |
CN113094938A (zh) * | 2021-06-10 | 2021-07-09 | 北京航空航天大学 | 一种面向海上搜救任务仿真的直升机油耗模型构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655215A (zh) * | 2013-08-12 | 2015-05-27 | 株式会社堀场制作所 | 燃料消耗量计算装置、方法和测量装置以及排气测量装置 |
CN106768123A (zh) * | 2016-11-30 | 2017-05-31 | 中国直升机设计研究所 | 一种无人直升机燃油预估方法 |
US20170152776A1 (en) * | 2015-11-30 | 2017-06-01 | United Technologies Corporation | Method of Real-Time Oil Consumption Detection |
CN107016880A (zh) * | 2017-04-10 | 2017-08-04 | 南京航空航天大学 | 一种基于模拟退火算法的连续下降进近油耗优化方法 |
-
2019
- 2019-01-21 CN CN201910057201.8A patent/CN109738035A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104655215A (zh) * | 2013-08-12 | 2015-05-27 | 株式会社堀场制作所 | 燃料消耗量计算装置、方法和测量装置以及排气测量装置 |
US20170152776A1 (en) * | 2015-11-30 | 2017-06-01 | United Technologies Corporation | Method of Real-Time Oil Consumption Detection |
CN106768123A (zh) * | 2016-11-30 | 2017-05-31 | 中国直升机设计研究所 | 一种无人直升机燃油预估方法 |
CN107016880A (zh) * | 2017-04-10 | 2017-08-04 | 南京航空航天大学 | 一种基于模拟退火算法的连续下降进近油耗优化方法 |
Non-Patent Citations (2)
Title |
---|
李杰等: "《首都国际机场航空器排放清单的计算分析》", 《中国环境科学》 * |
王蕾: "《基于航班轨迹的空中交通运行燃油效率分析》", 《中国优秀硕士论文全文数据库工程科技Ⅱ辑》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110503349A (zh) * | 2019-09-10 | 2019-11-26 | 深圳市瑞达飞行科技有限公司 | 航班油耗数据的处理方法、处理装置和电子设备 |
CN111507642A (zh) * | 2020-04-28 | 2020-08-07 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种基于航迹运行的经济性评估方法及装置 |
CN111507642B (zh) * | 2020-04-28 | 2023-08-29 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种基于航迹运行的经济性评估方法及装置 |
CN111581801A (zh) * | 2020-04-29 | 2020-08-25 | 南京航空航天大学 | 一种航空器巡航温室效应计算方法 |
CN111581801B (zh) * | 2020-04-29 | 2024-05-24 | 南京航空航天大学 | 一种航空器巡航温室效应计算方法 |
CN111696389A (zh) * | 2020-05-29 | 2020-09-22 | 航科院中宇(北京)新技术发展有限公司 | 一种基于航班飞行计划的航空器燃油估计方法及系统 |
CN113094938A (zh) * | 2021-06-10 | 2021-07-09 | 北京航空航天大学 | 一种面向海上搜救任务仿真的直升机油耗模型构建方法 |
CN113094938B (zh) * | 2021-06-10 | 2021-08-13 | 北京航空航天大学 | 一种面向海上搜救任务仿真的直升机油耗模型构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110930770B (zh) | 一种基于管制意图和飞机性能模型的四维航迹预测方法 | |
CN109738035A (zh) | 基于ads-b航迹数据的航空器油耗计算方法 | |
CN110046735B (zh) | 基于飞行数据分析的飞机离场燃油消耗评估方法 | |
Nuic et al. | BADA: An advanced aircraft performance model for present and future ATM systems | |
Pagoni et al. | Calculation of aircraft fuel consumption and CO2 emissions based on path profile estimation by clustering and registration | |
Lovegren et al. | Estimation of potential aircraft fuel burn reduction in cruise via speed and altitude optimization strategies | |
CN111177851A (zh) | 一种无人机运行安全风险评估中对地风险的评估方法 | |
CN109493644A (zh) | 一种基于历史航迹数据挖掘的四维航迹推测方法 | |
CN103294062A (zh) | 用于沿着飞行路径驾驶飞行器的方法 | |
CN106768123A (zh) | 一种无人直升机燃油预估方法 | |
CN107123315A (zh) | 一种考虑环境影响的终端区进离场航线优化方法 | |
CN109785462A (zh) | 航空器油耗计算系统 | |
Lee et al. | Closed-form takeoff weight estimation model for air transportation simulation | |
CN111192481B (zh) | 一种基于碰撞风险的进离场程序无人机管控区边界确定方法 | |
WO2009101988A1 (ja) | 乱気流予測システムおよび乱気流予測方法 | |
CN109598815A (zh) | 一种机载燃油系统耗油估计及健康监测方法 | |
Holzäpfel et al. | Assessment of wake-vortex encounter probabilities for crosswind departure scenarios | |
Murrieta-Mendoza et al. | Performance DataBase creation using a level D simulator for Cessna Citation X aircraft in cruise regime | |
Fernandes de Oliveira et al. | Benefits of optimal flight planning on noise and emissions abatement at the Frankfurt airport | |
CN106648673A (zh) | 一种连续下降运行程序的分析和设计方法 | |
CN106651014A (zh) | 一种运输机飞行轨迹优化方法 | |
Jiang et al. | Research on method of trajectory prediction in aircraft flight based on aircraft performance and historical track data | |
Bronsvoort et al. | Enhanced descent wind forecast for aircraft | |
Andreeva-Mor et al. | Scheduling of arrival aircraft based on minimum fuel burn descents | |
Pawlak et al. | Emission of selected exhaust compounds in jet engines of a jet aircraft in cruise phase |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20190510 |
|
WD01 | Invention patent application deemed withdrawn after publication |