CN111257506A - 一种燃油硫含量的自动识别算法 - Google Patents
一种燃油硫含量的自动识别算法 Download PDFInfo
- Publication number
- CN111257506A CN111257506A CN202010041516.6A CN202010041516A CN111257506A CN 111257506 A CN111257506 A CN 111257506A CN 202010041516 A CN202010041516 A CN 202010041516A CN 111257506 A CN111257506 A CN 111257506A
- Authority
- CN
- China
- Prior art keywords
- time
- peak
- suspected
- follows
- sulfur content
- 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.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital
- G01N33/0067—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method, e.g. intermittent, or the display, e.g. digital by measuring the rate of variation of the concentration
Abstract
本发明公开了一种燃油硫含量的自动识别算法,首先对遥测数据对给定的嫌疑时刻进行判别,结合时间序列的导数对峰值的起始和结束时刻进行精准识别并对嫌疑峰值的真假进行判定,最后结合嗅探法原理根据尾气硫碳比推算燃油硫含量,解决了船舶遥测尾气SO2、CO2遥测数据人为识别嫌疑波峰存在工作量大,识别结果不精准且无量化依据的问题,能够自动识别使用硫含量超标燃油的嫌疑船舶,为海事监管提供依据和支撑。
Description
技术领域
本发明涉及大气环境保护检测技术领域,更具体的说是涉及一种燃油硫含量的自动识别算法。
背景技术
嗅探法能够通过船舶尾气中的SO2、CO2浓度有效推算船舶燃油中的硫含量,可应用于船舶燃油监管中。目前已形成了监测船舶尾气SO2、CO2浓度的船舶尾气遥测仪,通常布置在船舶航行的必经之路上,通过采集过往船舶尾气中的SO2、 CO2浓度推算燃油中的硫含量。
但是,船舶尾气SO2、CO2浓度遥测数据时间序列数据量大,背景噪声变化复杂且具有随机性,其嫌疑波峰具有显著的非同步性。目前形成了数据采集平台,需要人为对数据进行分析判别和计算,目前尚无基于船舶尾气遥测SO2、CO2浓度数据自动识别燃油硫含量的计算方法。
因此,提供一种可用于自动识别使用硫含量超标燃油的嫌疑船舶的燃油硫含量自动识别算法是本领域技术人员亟需解决的问题。
发明内容
有鉴于此,本发明提供了一种燃油硫含量的自动识别算法,能够用于自动识别使用硫含量超标燃油的嫌疑船舶,解决了人为识别嫌疑波峰工作量大,识别结果不精准且无量化依据的问题。
为了实现上述目的,本发明采用如下技术方案:
一种燃油硫含量的自动识别算法,包括以下顺序步骤:
1)获取监测点位时间T内的SO2及CO2浓度随时间变化的数据CS(tk)和CC(tk) 其中k=1,2,…,N;tk=kΔt为时间序列,Δt为采样时间间隔,N=T/ Δt为数据长度;确定时间尺度Tp=PΔt,P为时间尺度所对应的数据长度;给定的嫌疑时刻为tq,q为嫌疑时刻所对应的时间序号;
2)计算时间序列CS(tk)和CC(tk),其中k=1,2,…,N的一阶差分,计算公式如下:
3)计算嫌疑峰值的峰值时刻tp,在给定的嫌疑时刻tq附近邻域搜索SO2浓度的局部最大值,峰值所对应的时间点记为tp,p为嫌疑峰值时刻所对应的时间序号,计算公式如下:
其中,δqP={k∈Z|max(1,q-P/2)≤k≤min(N,q+P/2)}表示tk在时刻tq附近的时间邻域;
6)进一步判别嫌疑峰值的真假,嫌疑峰值的有效峰值时段为tp0至tp1,峰值时段长度为τp=tp1–tp0,SO2和CO2浓度背景值分别记为BS和BC,计算公式如下:
SO2和CO2浓度增幅分别记为ΔS和ΔC,计算公式如下:
当SO2和CO2浓度增幅小于阈值,即ΔS<ΔS0或ΔC<ΔC0时,该嫌疑峰值判定为假峰值;
7)计算每个真峰值的SO2和CO2浓度增量积分,记为IS和IC,计算公式如下:
8)计算真峰值的硫碳比,从而推算燃油硫含量S,计算公式如下:
进一步的,数据长度N以及时间尺度所对应的数据长度P均取值为偶数。
通过采取以上方案,本发明的有益效果是:
本发明首先对遥测数据对给定的嫌疑时刻进行判别,结合时间序列的导数对峰值的起始和结束时刻进行精准识别并对嫌疑峰值的真假进行判定,最后结合嗅探法原理根据尾气硫碳比推算燃油硫含量,解决了船舶遥测尾气SO2、CO2 遥测数据人为识别嫌疑波峰存在工作量大,识别结果不精准且无量化依据的问题,能够自动识别使用硫含量超标燃油的嫌疑船舶,为海事监管提供依据和支撑。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据提供的附图获得其他的附图。
图1附图为本发明提供的船舶尾气遥测仪采集的SO2、CO2浓度的原始数据;
图2附图为本发明提供的出现的嫌疑峰值识别结果的曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例公开了一种燃油硫含量的自动识别算法,包括以下顺序步骤:
1)获取监测点位时间T内的SO2及CO2浓度随时间变化的数据CS(tk)和CC(tk) 其中k=1,2,…,N;tk=kΔt为时间序列,Δt为采样时间间隔,N=T/ Δt为数据长度;确定时间尺度Tp=PΔt,P为时间尺度所对应的数据长度;给定的嫌疑时刻为tq,q为嫌疑时刻所对应的时间序号;
2)计算时间序列CS(tk)和CC(tk),其中k=1,2,…,N的一阶差分,计算公式如下:
3)计算嫌疑峰值的峰值时刻tp,在给定的嫌疑时刻tq附近邻域搜索SO2浓度的局部最大值,峰值所对应的时间点记为tp,p为嫌疑峰值时刻所对应的时间序号,计算公式如下:
其中,δqP={k∈Z|max(1,q-P/2)≤k≤min(N,q+P/2)}表示tk在时刻tq附近的时间邻域;
6)判别嫌疑峰值的真假,嫌疑峰值的有效峰值时段为tp0至tp1,峰值时段长度为τp=tp1–tp0,SO2和CO2浓度背景值分别记为BS和BC,计算公式如下:
SO2和CO2浓度增幅分别记为ΔS和ΔC,计算公式如下:
当SO2和CO2浓度增幅小于阈值,即ΔS<ΔS0或ΔC<ΔC0时,该嫌疑峰值判定为假峰值;
7)计算每个真峰值的SO2和CO2浓度增量积分,记为IS和IC,计算公式如下:
8)计算真峰值的硫碳比,从而推算燃油硫含量S,计算公式如下:
本发明首先对遥测数据对给定的嫌疑时刻进行判别,结合时间序列的导数对峰值的起始和结束时刻进行精准识别并对嫌疑峰值的真假进行判定,最后结合嗅探法原理根据尾气硫碳比推算燃油硫含量,解决了船舶遥测尾气SO2、CO2遥测数据人为识别嫌疑波峰存在工作量大,识别结果不精准且无量化依据的问题,能够自动识别使用硫含量超标燃油的嫌疑船舶,为海事监管提供依据和支撑。
具体的,数据长度N以及时间尺度所对应的数据长度P均取值为偶数。
应用案例:
如图1、2所示,对苏通大桥船舶尾气遥测仪采集的船舶尾气中的SO2、CO2浓度时程进行处理,得到原始数据,其中,Δt=5s,对12:22出现的嫌疑峰值所对应的船舶燃油含量进行识别,判别阈值α取0.3;
首先对时间序列求导数,找到峰值时刻为12:22:05,再在其附近的时间林宇内找到峰值的起始时刻12:21:05,峰值的结束时刻12:23:05,得到峰值持续长度为120s,计算SO2和CO2浓度背景值分别为BS=-0.79,BC=485.65,计算SO2和 CO2浓度浓度增量积分IS=33.35,IC=1.84;最后根据硫碳比推算燃油硫含量S= 4.2%,远超过限值0.5%,海事部门登船检查燃油硫含量约为4%,与推算结果吻合,说明方法具有较高精度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
对所公开的实施例的上述说明,使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的,本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下,在其它实施例中实现。因此,本发明将不会被限制于本文所示的这些实施例,而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (2)
1.一种燃油硫含量的自动识别算法,其特征在于,包括以下顺序步骤:
1)获取监测点位时间T内的SO2及CO2浓度随时间变化的数据CS(tk)和CC(tk)其中k=1,2,…,N;tk=kΔt为时间序列,Δt为采样时间间隔,N=T/Δt为数据长度;确定时间尺度Tp=PΔt,P为时间尺度所对应的数据长度;给定的嫌疑时刻为tq,q为嫌疑时刻所对应的时间序号;
2)计算时间序列CS(tk)和CC(tk),其中k=1,2,…,N的一阶差分,计算公式如下:
3)计算嫌疑峰值的峰值时刻tp,在给定的嫌疑时刻tq附近邻域搜索SO2浓度的局部最大值,峰值所对应的时间点记为tp,p为嫌疑峰值时刻所对应的时间序号,计算公式如下:
其中,δqP={k∈Z|max(1,q-P/2)≤k≤min(N,q+P/2)}表示tk在时刻tq附近的时间邻域;
6)进一步判别嫌疑峰值的真假,嫌疑峰值的有效峰值时段为tp0至tp1,峰值时段长度为τp=tp1–tp0,SO2和CO2浓度背景值分别记为BS和BC,计算公式如下:
SO2和CO2浓度增幅分别记为ΔS和ΔC,计算公式如下:
当SO2和CO2浓度增幅小于阈值,即ΔS<ΔS0或ΔC<ΔC0时,该嫌疑峰值判定为假峰值;
7)计算每个真峰值的SO2和CO2浓度增量积分,记为IS和IC,计算公式如下:
8)计算真峰值的硫碳比,从而推算燃油硫含量S,计算公式如下:
2.根据权利要求1所述的一种燃油硫含量的自动识别算法,其特征在于,数据长度N以及时间尺度所对应的数据长度P均取值为偶数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010041516.6A CN111257506B (zh) | 2020-01-15 | 2020-01-15 | 一种燃油硫含量的自动识别算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010041516.6A CN111257506B (zh) | 2020-01-15 | 2020-01-15 | 一种燃油硫含量的自动识别算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111257506A true CN111257506A (zh) | 2020-06-09 |
CN111257506B CN111257506B (zh) | 2021-02-02 |
Family
ID=70947025
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010041516.6A Active CN111257506B (zh) | 2020-01-15 | 2020-01-15 | 一种燃油硫含量的自动识别算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111257506B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112763465A (zh) * | 2020-12-24 | 2021-05-07 | 交通运输部天津水运工程科学研究所 | 一种抗干扰的船用柴油硫含量超标嗅探识别方法 |
CN113281264A (zh) * | 2021-05-08 | 2021-08-20 | 安徽大学 | 一种基于卫星遥感确定环境大气污染点源一氧化碳排放通量的方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105428A (zh) * | 2013-01-14 | 2013-05-15 | 深圳市清时捷科技有限公司 | 一种判断水中重金属是否超限的方法 |
CN104307101A (zh) * | 2014-10-13 | 2015-01-28 | 中国人民解放军第三军医大学第三附属医院 | 用于产生驱动胃肠起搏系统的刺激电流的方法及装置 |
CN104382569A (zh) * | 2014-12-08 | 2015-03-04 | 天津工业大学 | 光纤传感智能服装及其心音、血压参数的处理方法 |
CN104597159A (zh) * | 2015-01-20 | 2015-05-06 | 力合科技(湖南)股份有限公司 | 基于色谱法的污染预警方法及装置 |
CN106248780A (zh) * | 2016-09-12 | 2016-12-21 | 安徽理工大学 | 一种反向脉冲hadamard离子迁移谱 |
CN107167833A (zh) * | 2017-05-10 | 2017-09-15 | 上海市计量测试技术研究院 | 一种γ谱假峰甄别方法、存储介质和系统 |
CN107179310A (zh) * | 2017-06-01 | 2017-09-19 | 温州大学 | 基于鲁棒噪声方差估计的拉曼光谱特征峰识别方法 |
CN107589100A (zh) * | 2017-09-08 | 2018-01-16 | 交通运输部天津水运工程科学研究所 | 一种船用燃油硫含量嗅探估算法 |
-
2020
- 2020-01-15 CN CN202010041516.6A patent/CN111257506B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103105428A (zh) * | 2013-01-14 | 2013-05-15 | 深圳市清时捷科技有限公司 | 一种判断水中重金属是否超限的方法 |
CN104307101A (zh) * | 2014-10-13 | 2015-01-28 | 中国人民解放军第三军医大学第三附属医院 | 用于产生驱动胃肠起搏系统的刺激电流的方法及装置 |
CN104382569A (zh) * | 2014-12-08 | 2015-03-04 | 天津工业大学 | 光纤传感智能服装及其心音、血压参数的处理方法 |
CN104597159A (zh) * | 2015-01-20 | 2015-05-06 | 力合科技(湖南)股份有限公司 | 基于色谱法的污染预警方法及装置 |
CN106248780A (zh) * | 2016-09-12 | 2016-12-21 | 安徽理工大学 | 一种反向脉冲hadamard离子迁移谱 |
CN107167833A (zh) * | 2017-05-10 | 2017-09-15 | 上海市计量测试技术研究院 | 一种γ谱假峰甄别方法、存储介质和系统 |
CN107179310A (zh) * | 2017-06-01 | 2017-09-19 | 温州大学 | 基于鲁棒噪声方差估计的拉曼光谱特征峰识别方法 |
CN107589100A (zh) * | 2017-09-08 | 2018-01-16 | 交通运输部天津水运工程科学研究所 | 一种船用燃油硫含量嗅探估算法 |
Non-Patent Citations (6)
Title |
---|
化学工业部人事教育司,化学工业部教育培训中心组织: "《化工分析仪表 2》", 31 March 1997, 化学工业出版社 * |
徐东群等: "《空气污染对人群健康影响数据清洗及评价方法》", 31 December 2016, 湖北科学技术出版社 * |
汪鲁才等: "水质检测中色谱谱峰寻峰算法的设计与实现", 《计算机工程与应用》 * |
胡健波等: "嗅探法: 通过监测尾气估算船用燃油硫含量的技术", 《水道港口》 * |
胡健波等: "通过尾气估算船用燃油硫含量的岸基嗅探法实验", 《水道港口》 * |
胡健波等: "通过尾气估算船用燃油硫含量的跟船嗅探法实验", 《水道港口》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112763465A (zh) * | 2020-12-24 | 2021-05-07 | 交通运输部天津水运工程科学研究所 | 一种抗干扰的船用柴油硫含量超标嗅探识别方法 |
CN112763465B (zh) * | 2020-12-24 | 2022-09-09 | 交通运输部天津水运工程科学研究所 | 一种抗干扰的船用柴油硫含量超标嗅探识别方法 |
CN113281264A (zh) * | 2021-05-08 | 2021-08-20 | 安徽大学 | 一种基于卫星遥感确定环境大气污染点源一氧化碳排放通量的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111257506B (zh) | 2021-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111257506B (zh) | 一种燃油硫含量的自动识别算法 | |
CN106354995A (zh) | 一种基于拉格朗日插值与时间序列的预测方法 | |
US9239372B2 (en) | Extracting spectral features from a signal in a multiplicative and additive noise environment | |
CN111475544A (zh) | 船舶轨迹数据中离群点的检测方法和装置 | |
RU2004128249A (ru) | Способ и устройство для обнаружения механических воздействий импульсного типа на компонент установки | |
CN111080976A (zh) | 一种温度变化场景下实时监测天然气泄漏的方法及装置 | |
CN113887324A (zh) | 基于卫星遥感数据的火点检测方法 | |
CN117436005A (zh) | 一种环境空气自动监测过程中异常数据处理方法 | |
CN111308015B (zh) | 一种船舶尾气遥测数据嫌疑波峰自动识别方法 | |
JP2007333515A (ja) | 目標検出装置 | |
CN117476140A (zh) | 基于物联网多维数据分析的大气环境污染治理方法 | |
CN116359876B (zh) | 一种雷达信号处理方法 | |
CN110400446B (zh) | 一种用于游泳池溺水检测的方法 | |
CN112217584A (zh) | 高动态ads-b电文的低误报率报头检测方法 | |
CN108737319B (zh) | 一种目标ofdm信号的实时检测方法及装置 | |
CN116255896A (zh) | 基于卫星的滑坡预警方法、装置和电子设备 | |
CN113534096A (zh) | 基于样条函数的LiDAR信号几何特征提取方法及系统 | |
CN114567340A (zh) | 基于相关算法的ads-b信号接收基带处理方法及系统 | |
Gangeskar | An adaptive method for estimation of wave height based on statistics of sea surface images | |
CN109945075B (zh) | 一种供水管道泄漏程度检测方法及装置 | |
KR20220153229A (ko) | 선박 테일 가스의 원격 측정 데이터의 의심 파봉 자동 식별 방법 | |
CN111401120A (zh) | 一种河道船流量分析系统及方法 | |
CN114305442B (zh) | 基于滑动窗编码的房颤发生起止点的检测方法 | |
NL2028253B1 (en) | Method and automatic recognition algorithm for sulfur content in fuel | |
CN115563889B (zh) | 一种目标对象的抛锚预测方法、装置及设备 |
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 |