CN110609017A - 一种基于反射率光谱的溢油风化时间分析方法及应用 - Google Patents
一种基于反射率光谱的溢油风化时间分析方法及应用 Download PDFInfo
- Publication number
- CN110609017A CN110609017A CN201910867400.5A CN201910867400A CN110609017A CN 110609017 A CN110609017 A CN 110609017A CN 201910867400 A CN201910867400 A CN 201910867400A CN 110609017 A CN110609017 A CN 110609017A
- Authority
- CN
- China
- Prior art keywords
- oil
- reflectivity
- weathering
- film
- oil film
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/55—Specular reflectivity
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种基于反射率光谱的溢油风化时间分析方法及应用,所述方法根据不同风化时间的溢油油膜物化性质变化导致的反射率光谱数据的差异,进行分析和反演溢油油膜的风化时间,分别建立反射率特征峰的溢油油膜风化时间测量模型和实测反射率光谱与反射率光谱数据库中数据的整体相关性分析方法,既从反射率光谱微观特征,又从整体宏观相关性上提高了抗环境干扰能力,实现溢油油膜风化时间的反演和测量;本发明在确定溢油风化时间的同时,会根据溢油油膜反射率光谱的特征进行油种的判别,可以有针对性的开展溢油油膜风化时间反演,具有更强的针对性,提高了风化时间反演的精度。
Description
技术领域
本发明属于水面油膜检测分析技术领域,涉及一种基于反射率光谱的溢油风化时间分析方法及应用。
背景技术
溢油事故的发生会对海洋生态环境带来较大的损害,泄漏到海面的油膜会受到风、浪、光照等条件的复杂作用,油膜的物理化学性质会随着时间推移产生差异和变化。获取准确的海上溢油油膜的风化时间,对于推断溢油发生时间、油种鉴别、溯源和处置具有重要的意义。常规的分析检测方法主要有:同位素分析法,热红外分析法和三维荧光法等。同位素分析法的检测精度较高,但需要复杂的采样和样品处理过程,并且需要大型精密仪器在实验室进行分析,不适用野外现场的实时分析。热红外分析法根据油膜热辐射的不同来判断油膜的风化时间,但是易受到周围环境的影响,检测的误差较大。三维荧光法根据激光荧光光谱特征的差异反演油膜风化时间,但需要功率较大的激光器作为激发源,并且受到油种、油膜厚度和背景光的干扰。
光学遥感传感器和光谱分析方法的不断发展,为海洋环境现场原位测量提供了新的技术方法和解决方案。目前,光学遥感在海上溢油的探测的应用较为普遍,主要解决海上溢油目标的识别和分布问题,对于海上溢油油膜风化时间的光谱分析研究相对较少,并且缺少不同风化时间的溢油光谱数据和反演油膜风化时间的系统方法。所以,建立一种基于反射率光谱的溢油风化时间分析系统及方法,对于海洋溢油的鉴别和处置具有重要意义。
发明内容
鉴于已有技术存在的缺陷,本发明的目的是要提供一种基于反射率光谱的溢油风化时间分析方法及应用,该系统能够通过反射率光谱数据确定海上溢油的主要油种(柴油和原油),进一步根据经风化后油膜的反射率光谱数据确定溢油风化时间,以获取溢油发生的时间信息,并能够实时显示反演结果,为海上溢油的鉴别和处置提供有效的参考信息。
为了实现上述目的,本发明的技术方案:
一种基于反射率光谱的溢油风化时间分析方法,所述方法包括:
收集不同风化时间原油和柴油的反射率光谱数据,并存储;
实时采集经风化后海面溢油油膜的反射率光谱数据,作为油膜风化时间分析反演的实测反射率光谱数据;
对实测反射率光谱数据进行平滑处理,提取光谱数据的一阶导数信息;
根据提取的光谱数据的一阶导数信息,确定溢油油膜的种类为柴油或原油;
从反射率光谱数据中提取柴油油膜在波长为383nm处的反射率光谱特征波峰值,创建柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线,从反射率光谱数据中提取原油油膜在波长为760nm处的反射率光谱特征波峰值,创建原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线;
根据确认的溢油油膜种类,若为柴油,则提取383nm处溢油油膜的反射率光谱特征波峰值,与柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间,若为原油,则提取760nm处溢油油膜的反射率光谱特征波峰值,与原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;
将实测反射率光谱数据与各反射率光谱数据进行相关性分析,获得相关性参数最高的结果,进一步确定溢油油膜风化时间;
溢油油膜风化时间的确认首先根据特征峰高度与风化时间的拟合曲线确定一个初步的风化时间,但是不同环境下都有拟合函数,会对应不同环境下的多个风化时间,所以再进行整体光谱曲线的相关性分析,确定所处的唯一环境条件,从而确定具体的风化时间。
将溢油油膜风化时间输出并显示。
所述反射率光谱数据为不同风化时间、不同环境(自然光照强度、天气状况等)条件下的原油和柴油的350-1000nm波段反射率光谱数据,并根据不同环境条件进行分类存储,作为实测反射率光谱数据的参考和对比数据,并作为实测反射率光谱数据特征分析的基准和依据。
通过ASD地物光谱仪、高光谱仪或荧光光谱仪实时采集经风化后海面溢油油膜的反射率光谱数据,采集的波段范围涵盖350-1000nm。
所述提取光谱数据的一阶导数信息为提取350-450nm波段的光谱数据的一阶导数信息;
提取的350-450nm波段的光谱数据的一阶导数信息,若存在一阶导数的变化规律为由正变负,则确定溢油油膜的种类为柴油,若存在一阶导数始终为负数,则确定溢油油膜的种类为原油。
所述柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.056-0.014t+0.001t2
其中,h代表反射率光谱相对强度,t代表风化时间;
所述原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.018-0.000414t+0.00000412t2+0.0000000132t3
其中,h代表反射率光谱相对强度,t代表风化时间。
本发明还公开了一种基于反射率光谱的溢油风化时间分析系统,所述系统包括:
反射率光谱数据库,用于收集不同风化时间原油和柴油的反射率光谱数据,并存储;
反射率光谱数据输入单元,用于实时采集经风化后海面溢油油膜的反射率光谱数据,作为油膜风化时间分析反演的实测反射率光谱数据;
第一特征分析单元,与所述反射率光谱数据输入单元连接,用于对接收到的实测反射率光谱数据进行平滑处理,提取光谱数据的一阶导数信息;
判断油种模块,与所述第一特征分析单元连接,用于根据所述第一特征分析单元提取的光谱数据的一阶导数信息,确定溢油油膜的种类为柴油或原油;
第二特征分析单元,分别与所述反射率光谱数据库和所述判断油种模块连接,所述第二特征分析单元包括柴油风化光谱特征分析单元和原油风化光谱特征分析单元;所述柴油风化光谱特征分析单元,用于从反射率光谱数据中提取柴油油膜在波长为383nm处的反射率光谱特征波峰值,创建柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为柴油,则提取383nm处溢油油膜的反射率光谱特征波峰值,与柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;所述原油风化光谱特征分析单元,用于从反射率光谱数据中提取原油油膜在波长为760nm处的反射率光谱特征波峰值,创建原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为原油,则提取760nm处溢油油膜的反射率光谱特征波峰值,与原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;
光谱曲线相关性分析单元,分别与所述第二特征分析单元和所述反射率光谱数据库连接,用于将实测反射率光谱数据与各反射率光谱数据进行相关性分析,获得相关性参数最高的结果,进一步确定溢油油膜风化时间;
以及,结果输出单元,与所述光谱曲线相关性分析单元连接,用于将溢油油膜风化时间输出并显示。
所述反射率光谱数据输入单元为ASD地物光谱仪、高光谱仪或荧光光谱仪。
本发明还公开了一种终端,包括处理器、输入设备、输出设备和存储器,所述处理器、输入设备、输出设备和存储器相互连接,其中,所述存储器用于存储计算机程序,所述计算机程序包括程序指令,所述处理器被配置用于调用所述程序指令,执行如上述所述的方法。
本发明还公开了一种计算机可读存储介质,所述计算机存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被处理器执行时使所述处理器执行如上述所述的方法。
本发明的工作原理如下:
A、首先需要获取不同风化时间溢油的反射率光谱数据样本,对所采集的反射率光谱数据样本按照不同环境条件(自然光照强度、天气状况等)进行分类存储,作为实测反射率光谱数据的参考和对比数据,并作为实测反射率光谱数据特征分析的基准和依据。对实测采集的350-450nm处的反射率光谱数据进行平滑处理,然后进行一阶求导操作,获取一阶导数信息。
B、其次,根据一阶导数信息判断溢油油种,根据350-450nm实测反射率光谱的一阶导数信息,柴油的反射率光谱数据的一阶导数变化规律为由正变负,原油的反射率光谱数据的一阶导数变化规律为始终为负数。
C、再次分别对反射率光谱数据库中的光谱数据和实测反射率光谱进行特征分析,通过实测反射率光谱数据与反射率光谱数据库中的数据提取对比分析,确定溢油油膜的风化时间;具体的,由于油的风化过程会改变油品的密度、黏度、组份等物化性质,改变后的油膜会引起其对太阳光不同波长能量吸收、反射变化,即不同风化时间的油膜存在明显的反射率光谱特征差异,利用这一特征差异可以反演出溢油发生后的风化时间,从而推算出溢油发生时间。在建立不同风化时间柴油的特征反射率光谱曲线时,发现不同风化时间的柴油油膜在波长为383nm处,会呈现一个特征峰,并且特征峰高度随着风化时间增加而逐渐减小,并随着风化时间的增加,该特征峰会逐渐消失,因此其可以作为确定柴油油膜风化时间的特征波段;相应的发现不同风化时间的原油油膜在760nm处,会呈现一个特征峰,并且特征峰高度随着风化时间增加而逐渐减小,并随着风化时间的增加,该特征峰会逐渐消失,因此其可以作为确定原油油膜风化时间的特征波段。在所述反射率光谱数据库中根据存储的光谱数据,建立油膜风化时间与反射率光谱特征峰高度的函数关系,即分别建立柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线(函数),原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线(函数),以构建反演油膜风化时间的测量模型。
D、再次,为了增强抗干扰能力,尤其针对不同环境条件的适用能力,根据反射率光谱数据库中预存的不同环境条件下的反射率光谱数据,对实测反射率光谱数据进行相关性分析。由于不同环境条件对整个350-1000nm光谱的影响是较为一致的,所以对实测数据的整体与数据库中数据进行相关性分析,相关性系数越大的曲线,可以进一步确定溢油油膜风化时间。最终通过构建特征峰拟合函数确定的风化时间方法和整体数据相关性分析,可以确定精确的溢油油膜风化时间,并通过通讯技术将测量结果显示在不同的显示平台上。
与现有技术相比,本发明的有益效果是:
本发明提出一种基于反射率光谱的溢油风化时间分析方法及应用,其根据不同风化时间的溢油油膜物化性质变化导致的反射率光谱数据的差异,进行分析和反演溢油油膜的风化时间,分别建立反射率特征峰的溢油油膜风化时间测量模型和实测反射率光谱与反射率光谱数据库中数据的整体相关性分析方法,既从反射率光谱微观特征,又从整体宏观相关性上提高了抗环境干扰能力,实现溢油油膜风化时间的反演和测量;本发明在确定溢油风化时间的同时,会根据溢油油膜反射率光谱的特征进行油种的判别,可以有针对性的开展溢油油膜风化时间反演,具有更强的针对性,提高了风化时间反演的精度;对于海上未知源溢油油膜的检测,可以根据本发明确认油膜的风化时间确定溢油源,为溢油的鉴别和处置提供有效的参考信息。
基于上述理由本发明可在水面油膜检测分析等领域广泛推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做以简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明的具体实施方式中基于反射率光谱的溢油风化时间分析系统的结构框图;
图2是本发明的具体实施方式中反射率光谱数据库中存储的柴油反射率光谱曲线实例图;
图3是本发明的具体实施方式中反射率光谱数据库中存储的原油反射率光谱曲线实例图;
图4是本发明的具体实施方式中平滑处理后柴油反射率光谱曲线,340-450nm一阶导数提取示意图;
图5是本发明的具体实施方式中平滑处理后原油反射率光谱曲线,340-450nm一阶导数提取示意图;
图6是本发明的具体实施方式中柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线实例图;
图7是本发明的具体实施方式中原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线实例图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
海洋溢油发生后,海面油膜会受到风、浪、光照等环境条件的共同作用,油膜的厚度会随着扩散漂移运动发生变化,其成份、密度等物理化学性质也会随着风化时间的增加而改变。相比而言,原油由于其成份复杂并且黏度大,发生泄漏后会在海面存留较长时间,风化的周期也更长;相对而言,柴油的成份较为单一,黏度小,外观颜色也更浅,在发生泄漏后会迅速扩散,风化周期相对于原油就会很短。所以在进行油膜风化时间确认时,需要考虑不同油种的风化特点,即需要首先根据实测反射率光谱数据判断出溢油油种,进而有针对性的进行风化时间的反演。另外,不同油种在反射率光谱中具有不同的特征峰,柴油油膜在波长为383nm处会呈现一个特征峰,并且特征峰高度随着风化时间增加而逐渐减小;原油油膜在760nm处会呈现一个特征峰,并且特征峰高度随着风化时间增加而逐渐减小。通过建立油膜风化时间与反射率光谱特征峰高度的关系,即分别建立柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线(函数),原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线(函数),以构建反演油膜风化时间的测量模型。反射率光谱测量数据往往会受到测量环境(自然光照强度、天气状况等)影响,环境因素会对整个测量光谱数据造成影响,为了提高溢油油膜风化时间反演的精确度,对实测反射率光谱数据与不同环境下数据进行相关性分析,进一步确定油膜风化时间,提高分析方法和系统的抗干扰能力。
基于上述原理,如图1所示,本发明所述的基于反射率光谱的溢油风化时间分析系统,包括:反射率光谱数据库,其能够预存不同环境(自然光照强度、天气状况等)条件下的原油和柴油的350-1000nm波段反射率光谱数据;反射率光谱数据输入单元,实时采集经风化后海面溢油油膜的反射率光谱数据,波段范围涵盖350-1000nm;分别与所述反射率光谱数据库和所述反射率光谱数据输入单元连接的第一特征分析单元,对接收到的实测光谱数据进行平滑处理,提取350-450nm波段的光谱数据的一阶导数信息,如图4和图5所示,与所述第一特征分析单元连接的判断油种模块,其能够确定油膜的种类,即确认油膜为柴油还是原油;分别与所述反射率光谱数据库和所述判断油种模块连接的第二特征分析单元,其包括:柴油风化光谱特征分析单元和原油风化光谱特征分析单元。其能够提取溢油油膜的反射率光谱特征波峰值,并与反射率光谱数据库数据中的拟合曲线进行参照比对,初步确定溢油油膜风化时间;分别与所述第二特征分析单元和所述反射率光谱数据库连接的光谱曲线相关性分析单元,将实测反射率光谱数据与反射率光谱数据库中预存的各反射率光谱数据(曲线)进行相关性分析,获得相关性参数最高的结果,进一步确定溢油的风化时间;与所述光谱曲线相关性分析单元连接的结果输出单元,将溢油风化时间结果输出并显示。
所述反射率光谱数据库,由于环境光照、温度等条件的不同会对反射率光谱曲线产生不同的影响,反射率光谱会受到环境因素对其整体数据进行影响,为了更准确的反演和确定溢油油膜风化时间,需要存储不同环境(自然光照强度、天气状况等)条件下的柴油和原油的350-1000nm波段反射率光谱数据,其中存储的柴油和柴油不同风化时间的反射率光谱曲线实例图,如图2,3所示。与实测反射率光谱数据进行相关性分析,相关性系数越大的曲线,可以进一步确定溢油油膜风化时间。
所述的第二特征分析单元中,所述柴油风化光谱特征分析单元,用于从反射率光谱数据中提取柴油油膜在波长为383nm处的反射率光谱特征波峰值,创建柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为柴油,则提取383nm处溢油油膜的反射率光谱特征波峰值,与柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;所述原油风化光谱特征分析单元,用于从反射率光谱数据中提取原油油膜在波长为760nm处的反射率光谱特征波峰值,创建原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为原油,则提取760nm处溢油油膜的反射率光谱特征波峰值,与原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;
所述柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.056-0.014t+0.001t2
其中,h代表反射率光谱相对强度,t代表风化时间,如图6所示;
所述原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.018-0.000414t+0.00000412t2+0.0000000132t3
其中,h代表反射率光谱相对强度,t代表风化时间,如图7所示。
所述光谱曲线相关性分析单元,根据反射率光谱数据库中预存的不同环境条件下的反射率光谱数据,对实测反射率光谱进行相关性分析,相关性系数越大的曲线,可以进一步确定溢油油膜风化时间,将测量结果显示在不同的显示平台上。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种基于反射率光谱的溢油风化时间分析方法,其特征在于,所述方法包括:
收集不同风化时间原油和柴油的反射率光谱数据,并存储;
实时采集经风化后海面溢油油膜的反射率光谱数据,作为油膜风化时间分析反演的实测反射率光谱数据;
对实测反射率光谱数据进行平滑处理,提取光谱数据的一阶导数信息;
根据提取的光谱数据的一阶导数信息,确定溢油油膜的种类为柴油或原油;
从反射率光谱数据中提取柴油油膜在波长为383nm处的反射率光谱特征波峰值,创建柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线,从反射率光谱数据中提取原油油膜在波长为760nm处的反射率光谱特征波峰值,创建原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线;
根据确认的溢油油膜种类,若为柴油,则提取383nm处溢油油膜的反射率光谱特征波峰值,与柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间,若为原油,则提取760nm处溢油油膜的反射率光谱特征波峰值,与原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;
将实测反射率光谱数据与各反射率光谱数据进行相关性分析,获得相关性参数最高的结果,进一步确定溢油油膜风化时间;
将溢油油膜风化时间输出并显示。
2.根据权利要求1所述的方法,其特征在于:所述反射率光谱数据为不同风化时间、不同环境条件下的原油和柴油的350-1000nm波段反射率光谱数据,并根据不同环境条件进行分类存储。
3.根据权利要求1所述的方法,其特征在于:通过ASD地物光谱仪、高光谱仪或荧光光谱仪实时采集经风化后海面溢油油膜的反射率光谱数据,采集的波段范围涵盖350-1000nm。
4.根据权利要求1所述的方法,其特征在于:所述提取光谱数据的一阶导数信息为提取350-450nm波段的光谱数据的一阶导数信息;
提取的350-450nm波段的光谱数据的一阶导数信息,若存在一阶导数的变化规律为由正变负,则确定溢油油膜的种类为柴油,若存在一阶导数始终为负数,则确定溢油油膜的种类为原油。
5.根据权利要求1所述的方法,其特征在于:所述柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.056-0.014t+0.001t2
其中,h代表反射率光谱相对强度,t代表风化时间;
所述原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线满足以下公式:
h=0.018-0.000414t+0.00000412t2+0.0000000132t3
其中,h代表反射率光谱相对强度,t代表风化时间。
6.一种基于反射率光谱的溢油风化时间分析系统,其特征在于,所述系统包括:
反射率光谱数据库,用于收集不同风化时间原油和柴油的反射率光谱数据,并存储;
反射率光谱数据输入单元,用于实时采集经风化后海面溢油油膜的反射率光谱数据,作为油膜风化时间分析反演的实测反射率光谱数据;
第一特征分析单元,与所述反射率光谱数据输入单元连接,用于对接收到的实测反射率光谱数据进行平滑处理,提取光谱数据的一阶导数信息;
判断油种模块,与所述第一特征分析单元连接,用于根据所述第一特征分析单元提取的光谱数据的一阶导数信息,确定溢油油膜的种类为柴油或原油;
第二特征分析单元,分别与所述反射率光谱数据库和所述判断油种模块连接,所述第二特征分析单元包括柴油风化光谱特征分析单元和原油风化光谱特征分析单元;所述柴油风化光谱特征分析单元,用于从反射率光谱数据中提取柴油油膜在波长为383nm处的反射率光谱特征波峰值,创建柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为柴油,则提取383nm处溢油油膜的反射率光谱特征波峰值,与柴油油膜厚度在383nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;所述原油风化光谱特征分析单元,用于从反射率光谱数据中提取原油油膜在波长为760nm处的反射率光谱特征波峰值,创建原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线,若所述判断油种模块确认的溢油油膜种类为原油,则提取760nm处溢油油膜的反射率光谱特征波峰值,与原油油膜厚度在760nm处反射率光谱特征峰高度与风化时间的拟合曲线进行参照比对,初步确定溢油油膜风化时间;
光谱曲线相关性分析单元,分别与所述第二特征分析单元和所述反射率光谱数据库连接,用于将实测反射率光谱数据与各反射率光谱数据进行相关性分析,获得相关性参数最高的结果,进一步确定溢油油膜风化时间;
以及,结果输出单元,与所述光谱曲线相关性分析单元连接,用于将溢油油膜风化时间输出并显示。
7.根据权利要求6所述的系统,其特征在于,所述反射率光谱数据输入单元为ASD地物光谱仪、高光谱仪或荧光光谱仪。
8.一种终端,其特征在于,包括处理器、输入设备、输出设备和存储器,所述处理器、输入设备、输出设备和存储器相互连接,其中,所述存储器用于存储计算机程序,所述计算机程序包括程序指令,所述处理器被配置用于调用所述程序指令,执行如权利要求1-5任一权利要求书所述的方法。
9.一种计算机可读存储介质,其特征在于,所述计算机存储介质存储有计算机程序,所述计算机程序包括程序指令,所述程序指令当被处理器执行时使所述处理器执行如权利要求1-5任一权利要求书所述的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910867400.5A CN110609017A (zh) | 2019-09-12 | 2019-09-12 | 一种基于反射率光谱的溢油风化时间分析方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910867400.5A CN110609017A (zh) | 2019-09-12 | 2019-09-12 | 一种基于反射率光谱的溢油风化时间分析方法及应用 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110609017A true CN110609017A (zh) | 2019-12-24 |
Family
ID=68891201
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910867400.5A Pending CN110609017A (zh) | 2019-09-12 | 2019-09-12 | 一种基于反射率光谱的溢油风化时间分析方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110609017A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624164A (zh) * | 2020-03-19 | 2020-09-04 | 自然资源部第一海洋研究所 | 基于光谱标准差特征分析的海洋溢油油种高光谱识别方法 |
CN117029701B (zh) * | 2023-10-09 | 2023-12-15 | 交通运输部天津水运工程科学研究所 | 一种沿海水域非接触式溢油监测方法 |
CN117330550A (zh) * | 2023-10-10 | 2024-01-02 | 交通运输部天津水运工程科学研究所 | 基于遥感图像和光谱图像的溢油溯源方法、终端、系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102997856A (zh) * | 2012-12-12 | 2013-03-27 | 南京大学 | 一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法 |
KR20150022470A (ko) * | 2013-08-23 | 2015-03-04 | 안원배 | 차량용 엔진오일의 열화도 측정방법 및 상기 방법을 이용한 엔진오일의 교환시기 결정방법 |
CN106767454A (zh) * | 2016-12-02 | 2017-05-31 | 大连海事大学 | 一种基于光谱反射率特征的水面油膜厚度测量系统及方法 |
CN107356535A (zh) * | 2017-06-12 | 2017-11-17 | 湖北久之洋红外系统股份有限公司 | 一种基于光谱成像技术的海上溢油检测方法 |
-
2019
- 2019-09-12 CN CN201910867400.5A patent/CN110609017A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102997856A (zh) * | 2012-12-12 | 2013-03-27 | 南京大学 | 一种基于参数查找表的海洋溢油油膜厚度高光谱遥感估算方法 |
KR20150022470A (ko) * | 2013-08-23 | 2015-03-04 | 안원배 | 차량용 엔진오일의 열화도 측정방법 및 상기 방법을 이용한 엔진오일의 교환시기 결정방법 |
CN106767454A (zh) * | 2016-12-02 | 2017-05-31 | 大连海事大学 | 一种基于光谱反射率特征的水面油膜厚度测量系统及方法 |
CN107356535A (zh) * | 2017-06-12 | 2017-11-17 | 湖北久之洋红外系统股份有限公司 | 一种基于光谱成像技术的海上溢油检测方法 |
Non-Patent Citations (6)
Title |
---|
R.FERNÁNDEZ-VARELA ET AL.: "Screening the origin and weathering of oil slicks by attenuated total reflectance mid-IR spectrometry", 《TALANTA》 * |
YONGCHAO HOU ET AL.: "Design and Implementation of a Coastal-Mounted Sensor for Oil Film Detection on Seawater", 《SENSORS》 * |
任广波等: "海面溢油无人机高光谱遥感检测与厚度估算方法", 《海洋学报》 * |
刘丙新: "基于高光谱特征的水上油膜提取与分析研究", 《中国优秀博硕士学位论文全文数据库(博士) 工程科技Ⅰ辑》 * |
刘丙新等: "基于光谱反射率数据的水面油种鉴别研究", 《光谱学与光谱分析》 * |
尹晓楠等: "基于三维荧光光谱和主成分分析的溢油风化研究", 《湖南大学学报(自然科学版)》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111624164A (zh) * | 2020-03-19 | 2020-09-04 | 自然资源部第一海洋研究所 | 基于光谱标准差特征分析的海洋溢油油种高光谱识别方法 |
CN111624164B (zh) * | 2020-03-19 | 2023-06-02 | 自然资源部第一海洋研究所 | 基于光谱标准差特征分析的海洋溢油油种高光谱识别方法 |
CN117029701B (zh) * | 2023-10-09 | 2023-12-15 | 交通运输部天津水运工程科学研究所 | 一种沿海水域非接触式溢油监测方法 |
CN117330550A (zh) * | 2023-10-10 | 2024-01-02 | 交通运输部天津水运工程科学研究所 | 基于遥感图像和光谱图像的溢油溯源方法、终端、系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110609017A (zh) | 一种基于反射率光谱的溢油风化时间分析方法及应用 | |
CN105486655B (zh) | 基于红外光谱智能鉴定模型的土壤有机质快速检测方法 | |
CN101458213B (zh) | 海洋溢油的浓度辅助参量荧光光谱油种鉴别方法 | |
CN108801934A (zh) | 一种土壤有机碳含量高光谱预测模型的建模方法 | |
CN108534811B (zh) | 一种短腔光纤法珀传感器的腔长解调算法 | |
CN103884661A (zh) | 基于土壤可见-近红外光谱库的土壤全氮实时检测方法 | |
CN104007115B (zh) | 一种利用太赫兹时域光谱技术检测珠宝结构的方法及系统 | |
CN112014331A (zh) | 一种水体污染的检测方法、装置、设备以及存储介质 | |
CN111126497A (zh) | 变压器固体绝缘材料老化状态评估方法 | |
CN114739919A (zh) | 一种基于光谱反演分析的水质检测方法 | |
CN110907388A (zh) | 一种基于红外光谱分析的溢油种类鉴别的方法 | |
CN112504983A (zh) | 一种基于浊度色度补偿的硝酸盐浓度预测方法 | |
CN115905881A (zh) | 黄珍珠分类的方法以及装置、电子设备、存储介质 | |
Ruan et al. | A modified backward elimination approach for the rapid classification of Chinese ceramics using laser-induced breakdown spectroscopy and chemometrics | |
CN110823829A (zh) | 一种基于sg校准片的光谱数据补偿方法 | |
CN117491293A (zh) | 一种基于高光谱的高陡岸坡碳酸盐岩体溶蚀快速评估方法 | |
CN110887798B (zh) | 基于极端随机树的非线性全光谱水体浊度定量分析方法 | |
CN205506679U (zh) | 一种便携式近红外光谱检测系统 | |
CN105067556A (zh) | 一种快速检测生物液态样品的光学装置及光学方法 | |
Zhou et al. | Retrieving quantum yield of sun-induced chlorophyll fluorescence near surface from hyperspectral in-situ measurement in productive water | |
Xing et al. | Estimation of chlorophyll-a concentrations in the Pearl River Estuary using in situ hyperspectral data: A case study | |
Dmitriev et al. | Parameterization of light absorption by components of seawater in optically complex coastal waters of the Crimea Peninsula (Black Sea) | |
CN113092376A (zh) | 基于激光光谱技术的手持式白酒真伪检测装置及方法 | |
Zhan et al. | Estimation of Optical Properties using QAA-V6 model based on MODIS data | |
CN112229817A (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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191224 |
|
RJ01 | Rejection of invention patent application after publication |