CN113607658B - 一种基于油膜灰度值获取油膜衰减系数的方法 - Google Patents

一种基于油膜灰度值获取油膜衰减系数的方法 Download PDF

Info

Publication number
CN113607658B
CN113607658B CN202110898233.8A CN202110898233A CN113607658B CN 113607658 B CN113607658 B CN 113607658B CN 202110898233 A CN202110898233 A CN 202110898233A CN 113607658 B CN113607658 B CN 113607658B
Authority
CN
China
Prior art keywords
oil film
attenuation coefficient
oil
gray value
obtaining
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
CN202110898233.8A
Other languages
English (en)
Other versions
CN113607658A (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.)
Dalian Maritime University
Original Assignee
Dalian Maritime University
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 Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN202110898233.8A priority Critical patent/CN113607658B/zh
Publication of CN113607658A publication Critical patent/CN113607658A/zh
Application granted granted Critical
Publication of CN113607658B publication Critical patent/CN113607658B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/55Specular reflectivity
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/20Controlling water pollution; Waste water treatment
    • Y02A20/204Keeping clear the surface of open water from oil spills

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)
  • Spectroscopy & Molecular Physics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

本发明公开了一种基于油膜灰度值获取油膜衰减系数的方法,包括如下步骤:搭建海水溢油实验模型,计算每组数据的平均值得到海水与油膜的反射灰度值,根据Snell折射定律得到折射角和入射角的关系,再根据麦克斯韦方程组计算推导平行和垂直偏振光的菲涅尔反射公式,从而获得最终的菲涅尔反射率;将两个不同油膜厚度下的灰度值做比,根据菲涅尔反射率公式推导三层介质双光束干涉模型,获得灰度值的比与衰减系数之间的关系式,从而获得油膜衰减系数;该方法结合溢油模拟实验数据来分析溢油油膜的衰减特性,为提高水面溢油目标的反演精度提供理论和技术基础,为溢油污染区域的有效探测与识别奠定技术和数据基础。

Description

一种基于油膜灰度值获取油膜衰减系数的方法
技术领域
本发明涉及水面环境检测技术领域,尤其涉及一种基于油膜灰度值获取油膜衰减系数的方法。
背景技术
近年来,随着海洋石油勘探开发和海洋石油运输活动日益频繁,溢油事故频发,海洋溢油污染已成为海洋环境的最主要威胁之一。基于双光束干涉模型分析油膜对入射光的辐射传输过程,将均匀分布的浮油膜视为一个平行平板,建立大气-油膜-水面三层介质物理模型,建立油膜反射特性与光束入射角、传感器观测角、光谱波长之间的函数关系,分析光谱反射率的探测影响因素,具有明确的理论意义和良好的模拟精度。基于该理论可以建立浮油厚度与光谱反射率之间的定量关系,可以更好地解释入射光在海面浮油上的光学特性。光学遥感技术是检测海洋表层石油的最常用技术之一。这是因为油具有与水不同的光学特性,有不同的光谱对比度。因此,对被油污染的海水的光学特性的研究将为溢油的探测提供理论依据,对做出应急决策和保护海洋生态环境具有重要意义。光学遥感技术是检测海洋表层石油的最常用技术之一,由于油具有与水不同的光学特性,有不同的光谱对比度。因此,对被油污染的海水的光学特性的研究将为溢油的探测提供理论依据,对做出应急决策和保护海洋生态环境具有重要意义。由于实验仪器响应的存在,导致油膜衰减系数的最终结果不准确
对于传统介质的衰减系数主要是用椭偏仪测得的,如果样品是有机物薄膜,椭偏仪很难测出结果;另外椭偏仪测量得到的数据需要手工处理,建立模型,很复杂,非专业人员无法处理,暂时没有现成手段可以直接得油膜的衰减系数;测量衰减系数时不能同时测量透射光与反射光,所以得到的吸收系数存在一定误差。
发明内容
根据现有技术存在的问题,本发明公开了一种基于油膜灰度值获取油膜衰减系数的方法,具体包括如下步骤:
搭建海水溢油实验模型,设置不同的两个油膜厚度d1和d2,分别测量两个不同油膜厚度下的油膜表面灰度值;
采用光纤和准直镜头对同一测量目标进行重复测量剔除误差较大的数据,计算每组数据的平均值得到海水与油膜的反射灰度值;
将海水溢油实验模型推导至平静水面上,设空气折射率是n1,油膜折射率是n2,油膜厚度为d,海水折射率是n3,设两种介质的折射率分别为ni、nt,根据Snell折射定律得到折射角和入射角的关系:ni sinθi=nt sinθt,再根据麦克斯韦方程组计算推导平行和垂直偏振光的菲涅尔反射公式,从而获得最终的菲涅尔反射率;
将两个不同油膜厚度下的灰度值做比,根据菲涅尔反射率公式推导三层介质双光束干涉模型,获得灰度值的比与衰减系数之间的关系式,从而获得油膜衰减系数。
进一步的,所述油膜衰减系数a采用如下方式获取:
油膜表面反射率与仪器响应存在如下关系:
其中R1,R2为两个厚度下的反射率,N1,N2为两个厚度下的测量灰度值;
根据菲涅尔反射双光束干涉原理,油膜表面反射率与油膜的厚度和衰减系数存在如下关系:
其中E0为入射光电场强度,E1、E2为两道出射光的电场度,E为两束光干涉后的电场强度,Δ=2n2d,s=2d,r和t为反射和透射系数,通过菲涅尔反射公式求得;
将上述两式结合并令将化简得/>因此获得油膜衰减系数a。
由于采用了上述技术方案,本发明提供的一种基于油膜灰度值获取油膜衰减系数的方法,该方法具体是基于海水溢油实验,然后得到两个油膜厚度下对应的灰度曲线,利用于菲涅尔反射公式的物理模型并结合溢油模拟实验数据来分析溢油油膜的衰减特性,为提高水面溢油目标的反演精度提供理论和技术基础,为溢油污染区域的有效探测与识别奠定技术和数据基础,维护海洋安全、提升监测能力等提供理论与技术依据。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明方法的流程图;
图2为本发明中模拟溢油海表面的光路图。
具体实施方式
为使本发明的技术方案和优点更加清楚,下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚完整的描述:
如图1和图2所示的一种基于油膜灰度值获取油膜衰减系数的方法,具体包括如下步骤:
S1:搭建实验大致如图2所示,设光辐射照度Eλ,油膜厚度为d1时,光谱反射率为R1,仪器响应ρ,仪器灰度值为N1,则满足如下关系EλR1ρ=N1;同理,油膜厚度为d2时,光谱反射率为R2,仪器响应ρ,仪器灰度值为N2同样有如下关系EλR2ρ=N2。上述中反射率是唯一与衰减系数有关的量。
S2:不同厚度油膜反射灰度值数据测量时间选择在密闭的暗室下进行,目的是避免其他杂散光对实验造成干扰,且测量人员着深色服装;光谱测量设备为EDMUND OPTICS光谱仪,光谱范围为350-2500nm,数据采集前,先调试镜头到合适的焦距,目的是为了镜头获得最大范围进光量。测量时探头距离参考板20-50cm且垂直向下,视场角3°-25°;光源是内置原装进口高压短弧氙灯,最大特点是光谱范围宽,可覆盖紫外、可见和近红外,因其色温接近太阳,可用于模拟太阳;由于光源出来的光束是发散的,所以搭配一根长度大约为两米的光纤和准直镜头。准直光斑直径约为10mm,发散角0.2°,焦距40mm。对同一测量目标重复测量10次,剔除误差较大的数据后,计算每组数据的平均值,得到海水与油膜的反射灰度值。
S3:模型推导平静水面上,在发生溢油情况下,从上到下依次有空气、油膜和海水三层介质,由于不同介质所对应的不同折射率,所以会导致入射光在油膜层内折射并反射。假设空气折射率是n1;油膜折射率是n2,油膜厚度为d;海水折射率是n3
S4:已知两种介质的折射率分别为ni、nt首先根据Snell折射定律得到折射角和入射角的关系:nisinθi=ntsinθt,然后根据麦克斯韦方程组计算推导,平行和垂直偏振光的菲涅尔反射公式:
则最终的菲涅尔反射率为:
S5:设入射光电场强度为E0,波数为k=2π/λ,当入射光射入油膜层时,油膜上表面透射系数为t12,反射系数为r12,如图2所示,则入射光在油膜上表面的第一次反射的电场强度E1=r12E0;第一次透射入油膜的光电矢量强度为t12E0;油膜层内的透射光射入油膜与海水界面时,即射入油膜层下表面后光的电场强度为t12r23E0,油膜层内透射光经过油膜下表面反射后,从油膜层的上表面透射入空气油膜层中透射光经油膜下表面的反射,又经过油膜上表面透射至空气层,其电场强度为E2。考虑到E1和E2在传播路径上有长度为Δ的差别,两束光之间则存在相位延迟kΔ,若光束在油膜层传播过程中的消光系数为a,从油膜下表面反射回空气的光束在油膜中传播路径长度为s,则透射光电场强度为E2=t12r23t21eikΔ- asE0。平行平板干涉也就是反射光E1与油膜出射光E2,两者的电矢量叠加后产生双光束干涉。对于探测器,E1与E2叠加产生干涉,进入遥感器系统的电场E=E1+E2,为便于计算方便,取入射角为0°,则油膜反射率R可以计算得到,即,
其中Δ=2n2d,s=2d。
以上模型描述了油膜反射率与油膜厚度、入射光折射角、各界面反射率和透过率之间的理论关系,与灰度值得转换如下EλRρ=N,Eλ是入射光辐射照度,R是油膜光谱反射率,仪器响应ρ,仪器灰度值为N。
S6:对获得的灰度值左右作比,则可约去辐射照度与仪器响应参数,得到如下关系:将双光束干涉所得的油膜反射率R带入公式,令/>根据此关系得到油膜的衰减系数a。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,根据本发明的技术方案及其发明构思加以等同替换或改变,都应涵盖在本发明的保护范围之内。

Claims (1)

1.一种基于油膜灰度值获取油膜衰减系数的方法,其特征在于包括:
搭建海水溢油实验模型,设置不同的两个油膜厚度d1和d2,分别测量两个不同油膜厚度下的油膜表面灰度值;
采用光纤和准直镜头对同一测量目标进行重复测量剔除误差较大的数据,计算每组数据的平均值得到海水与油膜的反射灰度值;
将海水溢油实验模型推导至平静水面上,设空气折射率是n1,油膜折射率是n2,油膜厚度为d,海水折射率是n3,根据Snell折射定律得到折射角和入射角的关系再根据麦克斯韦方程组计算推导平行和垂直偏振光的菲涅尔反射公式,从而获得最终的菲涅尔反射率;
将两个不同油膜厚度下的灰度值做比,根据菲涅尔反射率公式推导三层介质双光束干涉模型,获得灰度值的比与衰减系数之间的关系式,从而获得油膜衰减系数;
所述油膜衰减系数a采用如下方式获取:
油膜表面反射率与仪器响应存在如下关系:
其中R1,R2为两个厚度下的反射率,N1,N2为两个厚度下的测量灰度值;
根据菲涅尔反射双光束干涉原理,油膜表面反射率与油膜的厚度和衰减系数存在如下关系:
其中E0为入射光电场强度,E1、E2为两道出射光的电场度,E为两束光干涉后的电场强度,Δ=2n2d,s=2d,r和t为反射和透射系数,通过菲涅尔反射公式求得;
将上述两式结合并令将化简得/>因此获得油膜衰减系数a。
CN202110898233.8A 2021-08-05 2021-08-05 一种基于油膜灰度值获取油膜衰减系数的方法 Active CN113607658B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110898233.8A CN113607658B (zh) 2021-08-05 2021-08-05 一种基于油膜灰度值获取油膜衰减系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110898233.8A CN113607658B (zh) 2021-08-05 2021-08-05 一种基于油膜灰度值获取油膜衰减系数的方法

Publications (2)

Publication Number Publication Date
CN113607658A CN113607658A (zh) 2021-11-05
CN113607658B true CN113607658B (zh) 2023-12-12

Family

ID=78307256

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110898233.8A Active CN113607658B (zh) 2021-08-05 2021-08-05 一种基于油膜灰度值获取油膜衰减系数的方法

Country Status (1)

Country Link
CN (1) CN113607658B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117006956B (zh) * 2023-10-07 2023-12-22 广东海洋大学 一种液体表面的油膜厚度测量方法及系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000028484A (ja) * 1998-07-08 2000-01-28 Toppan Printing Co Ltd 薄膜の光学定数測定方法及びこの装置
RU141304U1 (ru) * 2014-01-31 2014-05-27 Дмитрий Владимирович Кизеветтер Устройство для оперативного контроля качества технического масла
CN104614324A (zh) * 2015-03-04 2015-05-13 哈尔滨工业大学 一种微藻培养基光学常数及微藻光谱衰减系数的联合测量方法
CN104897550A (zh) * 2015-06-17 2015-09-09 大连理工大学 一种确认碳纤维增强复合材料孔隙率与超声衰减系数关系的方法
RU2014152767A (ru) * 2014-12-24 2016-07-20 Федеральное государственное бюджетное учреждение науки "Морской гидрофизический институт РАН" (ФГБУН "МГИ РАН") Способ определения спектрального показателя ослабления направленного света в морской воде "in situ"
CN105823787A (zh) * 2016-03-15 2016-08-03 南昌航空大学 激光快速成形tc4材料射线衰减系数与锻造tc4射线衰减系数差异比较测量工艺
CN106767454A (zh) * 2016-12-02 2017-05-31 大连海事大学 一种基于光谱反射率特征的水面油膜厚度测量系统及方法
CN110160452A (zh) * 2019-06-19 2019-08-23 哈尔滨工业大学(威海) 基于激光拉曼及激光荧光的水面油膜厚度的测量方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000028484A (ja) * 1998-07-08 2000-01-28 Toppan Printing Co Ltd 薄膜の光学定数測定方法及びこの装置
RU141304U1 (ru) * 2014-01-31 2014-05-27 Дмитрий Владимирович Кизеветтер Устройство для оперативного контроля качества технического масла
RU2014152767A (ru) * 2014-12-24 2016-07-20 Федеральное государственное бюджетное учреждение науки "Морской гидрофизический институт РАН" (ФГБУН "МГИ РАН") Способ определения спектрального показателя ослабления направленного света в морской воде "in situ"
CN104614324A (zh) * 2015-03-04 2015-05-13 哈尔滨工业大学 一种微藻培养基光学常数及微藻光谱衰减系数的联合测量方法
CN104897550A (zh) * 2015-06-17 2015-09-09 大连理工大学 一种确认碳纤维增强复合材料孔隙率与超声衰减系数关系的方法
CN105823787A (zh) * 2016-03-15 2016-08-03 南昌航空大学 激光快速成形tc4材料射线衰减系数与锻造tc4射线衰减系数差异比较测量工艺
CN106767454A (zh) * 2016-12-02 2017-05-31 大连海事大学 一种基于光谱反射率特征的水面油膜厚度测量系统及方法
CN110160452A (zh) * 2019-06-19 2019-08-23 哈尔滨工业大学(威海) 基于激光拉曼及激光荧光的水面油膜厚度的测量方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Oil-film extinction coefficient inversion based on thickness difference;ZHENDUO ZHANG 等;Optics Express;第30卷(第17期);第30368-30378页 *
基于激光拉曼光谱的水面油膜厚度测量方法研究;蔡宗岐 等;光谱学与光谱分析(第06期);第7-10页 *
用于海上油膜厚度测量的激光光强调节装置设计;陈曦 等;天津大学学报(自然科学与工程技术版);第52卷(第6期);第576-584页 *

Also Published As

Publication number Publication date
CN113607658A (zh) 2021-11-05

Similar Documents

Publication Publication Date Title
Gilerson et al. Fluorescence component in the reflectance spectra from coastal waters. Dependence on water composition
Wettle et al. Assessing the effect of hydrocarbon oil type and thickness on a remote sensing signal: A sensitivity study based on the optical properties of two different oil types and the HYMAP and Quickbird sensors
Ma et al. Monte Carlo simulation of spectral reflectance and BRDF of the bubble layer in the upper ocean
CN105891152A (zh) 一种大量程折射率测量的方法
Deng et al. Pure water absorption coefficient measurement after eliminating the impact of suspended substance in spectrum from 400 nm to 900 nm
CN113607658B (zh) 一种基于油膜灰度值获取油膜衰减系数的方法
Chang et al. Toward closure of upwelling radiance in coastal waters
Badewien et al. Towards continuous long-term measurements of suspended particulate matter (SPM) in turbid coastal waters
Bernacki et al. Modeling liquid organic thin films on substrates
CN105352915A (zh) 一种折射率二维分布的动态测量方法
CN102243065B (zh) 一种基于背向补偿的透明基底薄膜厚度测量系统
Bahrami et al. Dual-wavelength spectroscopy of a metallic-grating-coupled surface plasmon resonance biosensor
Bulgarelli et al. Measured and modeled radiometric quantities in coastal waters: toward a closure
Froidefond et al. Introducing a mini-catamaran to perform reflectance measurements above and below the water surface
Tanaka et al. Alternative measuring method for water-leaving radiance using a radiance sensor with a domed cover
RU2605640C2 (ru) СПОСОБ ОПРЕДЕЛЕНИЯ СПЕКТРАЛЬНОГО ПОКАЗАТЕЛЯ ОСЛАБЛЕНИЯ НАПРАВЛЕННОГО СВЕТА В МОРСКОЙ ВОДЕ "in situ"
CN113433077B (zh) 一种基于油膜相对厚度获取油膜衰减系数的方法
RU2300077C1 (ru) Дистанционный способ измерения толщины толстых пленок нефтепродуктов на поверхности воды
Pavlova et al. Deck spectroradiometer for measuring remote sensing reflectance
Shang et al. Impact of ship on radiometric measurements in the field: a reappraisal via Monte Carlo simulations
Sun et al. Effects of water salinity on the multi-angular polarimetric properties of light reflected from smooth water surfaces
CN116843634A (zh) 一种基于改进双光束模型与偏振图像结合的油厚解算方法
Otremba The influence of an oil-film covered sea surface on the reflection and upward transmission of light
Nikonenko et al. Measurement of the thickness of thin polymer layers by infrared frustrated total internal reflection spectroscopy
Shvets A method for measuring the phase of the reflection coefficient in the visible range of the spectrum

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