CN106841821A - 一种薄膜覆盖海面等效介电常数计算方法 - Google Patents

Abstract

本发明公开了一种薄膜覆盖海面等效介电常数计算方法，将薄膜‑海水复合海面等效为介电常数为ε_e的单层反射介质，根据布拉格(Bragg)散射机制公式，求出仅与雷达入射角、反射面介电常数有关的极化比参数。真实的极化比可用空气‑薄膜‑海水多层介质散射理论公式计算，该极化比为雷达入射角和空气、薄膜、海水的相对介电常数的函数。由两种极化比相等所得公式，可用求解出等效介电常数，包括水平极化和垂直极化下的两种等效介电常数计算方法。本发明计算方法为污染物自动监测技术提供一种物理识别方法。

Description

一种薄膜覆盖海面等效介电常数计算方法

技术领域

本发明属于海洋技术领域，具体涉及一种基于双极化多层薄膜覆盖海面等效介电常数计算方法。

背景技术

为了防止、减少溢油灾害，对海面全天候、大范围、高分辨卫星遥感实时监测必不可少。对海面溢油监测主要采用卫星SAR遥感监测技术，但现有的溢油自动识别技术物理机制尚不完善，监测结果主要有两个方面不足，一是对于低风速区域海面和薄冰海面等非溢油海面误报为溢油污染，二是尚不能对海表面污染物进行区分，如矿物油，植物油，赤潮等污染物，形成当前污染物自动监测的技术瓶颈。我们提出一种基于海面多层电磁散射机制的薄膜覆盖海面等效介电常数计算方法，为污染物自动监测技术提供一种物理识别方法。

发明内容

本发明的目的在于克服现有技术的不足，提供一种基于双极化多层薄膜覆盖海面等效介电常数计算方法，计算结果用于区别不同海面污染物。

本发明采用的技术方案为：一种薄膜覆盖海面等效介电常数计算方法，将薄膜-海水复合海面等效为介电常数为ε_e的单层反射介质，根据布拉格(Bragg)散射机制公式，求出仅与雷达入射角、反射面介电常数有关的极化比参数。真实的极化比可用空气-薄膜-海水多层介质散射理论公式计算，该极化比为雷达入射角和空气、薄膜、海水的相对介电常数的函数。由两种极化比相等所得公式，可用求解出等效介电常数，包括水平极化和垂直极化下的两种等效介电常数计算方法。具体包括以下步骤：

步骤:1：考虑空气-薄膜-海水三层复合介质海面，它们的介电常数和磁导率分别为：ε₀,μ₀、ε₁,μ₁和ε₂,μ₂，三种介质均为非铁磁质，磁导率可近似为μ₁＝μ₂＝μ₀，薄膜厚度为h，波数为的电磁波，波长为λ，k_i入射角为θ_i入射到薄膜表面，透射角为当透射波透入海水中是透射角为在空气-薄膜界面即z＝0处，HH和VV极化反射系数和透射系数其中p＝HH,VV，分别为：

在薄膜-海水界面即z＝-h处：

其中由公式(1)-(8)，根据多层散射理论，在空气-薄膜界面的反射系数为：

步骤2：将薄膜-海水等效为一层等效介质，等效介电常数和等效磁导率为ε_e,μ_e，其中等效介质为非铁磁介质磁导率μ_e＝μ₀；在空气-等效介质界面即z＝0处，HH和VV极化等效反射系数分别为：

步骤3：水平极化和垂直极化介电常数分别为和他们分别由公式(12)和(13)计算。

有益效果：本发明提出的薄膜覆盖海面等效介电常数计算方法，是基于海面多层电磁散射机制的薄膜覆盖海面等效介电常数计算方法，为污染物自动监测技术提供一种物理识别方法。

附图说明

图1为薄膜覆盖海冰结构示意图；其中(a)部分为厚度为h的薄膜覆盖海面，(b)部分为与(a)部分等效的海面示意图；

图2为有效介电常数随薄膜相对介电常数的变化。

具体实施方式

下面结合附图和具体实施方式对本发明作进一步的说明。

本发明考虑水平极化(HH)和垂直极化(VV)两种情形，提出一种基于多层散射理论的薄膜覆盖海面等效介电常数计算方法，具体步骤主要包含：

步骤:1：考虑空气-薄膜-海水三层复合介质海面(见图1的a部分)，它们的介电常数和磁导率分别为：ε₀,μ₀、ε₁,μ₁和ε₂,μ₂，三种介质均为非铁磁质，磁导率可近似为μ₁＝μ₂＝μ₀，薄膜厚度为h，波数为的电磁波，波长为λ，k_i入射角为θ_i入射到薄膜表面，透射角为当透射波透入海水中是透射角为在空气-薄膜界面即z＝0处，HH和VV极化反射系数和透射系数其中p＝HH,VV，分别为：

在薄膜-海水界面即z＝-h处：

其中由公式(1)-(8)，根据多层散射理论，在空气-薄膜界面的反射系数为：

步骤2：将图1的a部分中薄膜-海水等效为一层等效介质(见图1的b部分)，等效介电常数和等效磁导率为ε_e,μ_e，其中等效介质为非铁磁介质磁导率μ_e＝μ₀；在空气-等效介质界面即z＝0处，HH和VV极化等效反射系数分别为：

步骤3：水平极化和垂直极化介电常数分别为和他们分别由公式(12)和(13)计算。

当薄膜厚度为0.01m，海水相对介电常数为80，模拟Radarsat2雷达波长为0.0595m，薄膜相对介电常数从1变为35，雷达入射角为θ_i＝30°。根据公式(12)和(13)，计算出的水平极化和垂直极化等效介电常数随薄膜的相对介电常数变化曲线如图2所示。

以上结合附图对本发明的实施方式做出详细说明，但本发明不局限于所描述的实施方式。对本领域的普通技术人员而言，在本发明的原理和技术思想的范围内，对这些实施方式进行多种变化、修改、替换和变形仍落入本发明的保护范围内。

Claims (1)

1.一种薄膜覆盖海面等效介电常数计算方法，其特征在于：包括以下步骤：

步骤:1：考虑空气-薄膜-海水三层复合介质海面，它们的介电常数和磁导率分别为：ε₀,μ₀、ε₁,μ₁和ε₂,μ₂，三种介质均为非铁磁质，磁导率可近似为μ₁＝μ₂＝μ₀，薄膜厚度为h，波数为的电磁波，波长为λ，k_i入射角为θ_i入射到薄膜表面，透射角为当透射波透入海水中是透射角为在空气-薄膜界面即z＝0处，HH和VV极化反射系数和透射系数其中p＝HH,VV，分别为：

$F_{HH}^{01}=\frac{k_{z0}-k_{z1}}{k_{z0}+k_{z1}}---\left(1\right)$

$F_{VV}^{01}=\frac{{\mathrm{\ϵ}}_1{k}_{z0}-{\mathrm{\ϵ}}_0{k}_{z1}}{{\mathrm{\ϵ}}_1{k}_{z0}+{\mathrm{\ϵ}}_0{k}_{z1}}---\left(2\right)$

$T_{HH}^{01}=\frac{2k_{z0}}{k_{z0}+k_{z1}}---\left(3\right)$

$T_{VV}^{01}=\frac{2{\mathrm{\ϵ}}_1{k}_{z0}}{{\mathrm{\ϵ}}_1{k}_{z0}+{\mathrm{\ϵ}}_0{k}_{z1}}---\left(4\right)$

在薄膜-海水界面即z＝-h处：

$F_{HH}^{12}=\frac{k_{z1}-k_{z2}}{k_{z1}+k_{z2}}---\left(5\right)$

$F_{VV}^{12}=\frac{{\mathrm{\ϵ}}_2{k}_{z1}-{\mathrm{\ϵ}}_1{k}_{z2}}{{\mathrm{\ϵ}}_2{k}_{z1}+{\mathrm{\ϵ}}_1{k}_{z2}}---\left(6\right)$

$T_{HH}^{12}=\frac{2k_{z1}}{k_{z1}+k_{z2}}---\left(7\right)$

$T_{VV}^{12}=\frac{2{\mathrm{\ϵ}}_2{k}_{z1}}{{\mathrm{\ϵ}}_2{k}_{z1}+{\mathrm{\ϵ}}_1{k}_{z2}}---\left(8\right)$

其中由公式(1)-(8)，根据多层散射理论，在空气-薄膜界面的反射系数为：

$R_p=\frac{F_p^{01}+F_p^{12}{e}^{j\left(2k_{z1}h\right)}}{1+F_p^{01}{F}_p^{12}{e}^{j\left(2k_{z1}h\right)}}---\left(9\right)$

步骤2：将薄膜-海水等效为一层等效介质，等效介电常数和等效磁导率为ε_e,μ_e，其中等效介质为非铁磁介质磁导率μ_e＝μ₀；在空气-等效介质界面即z＝0处，HH和VV极化等效反射系数分别为：

$F_{HH}^e=\frac{{\mathrm{cos\θ}}_i-{({\mathrm{\ϵ}}_e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}{{\mathrm{cos\θ}}_i+{({\mathrm{\ϵ}}_e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}---\left(10\right)$

$F_{VV}^e=\frac{{\mathrm{\ϵ}}_e{\mathrm{cos\θ}}_i-{({\mathrm{\ϵ}}_e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}{{\mathrm{\ϵ}}_e{\mathrm{cos\θ}}_i+{({\mathrm{\ϵ}}_e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}---\left(11\right)$

步骤3：水平极化和垂直极化介电常数分别为和他们分别由公式(12)和(13)计算。

$\frac{{\mathrm{cos\θ}}_i-{({\mathrm{\ϵ}}_{HH}^e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}{{\mathrm{cos\θ}}_i+{({\mathrm{\ϵ}}_{HH}^e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}=\frac{F_{HH}^{01}+F_{HH}^{12}{e}^{j\left(2k_{z1}h\right)}}{1+F_{HH}^{01}{F}_{HH}^{12}{e}^{j\left(2k_{z1}h\right)}}---\left(12\right)$

$\frac{{\mathrm{\ϵ}}_{VV}^e{\mathrm{cos\θ}}_i-{({\mathrm{\ϵ}}_{VV}^e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}{{\mathrm{\ϵ}}_{VV}^e{\mathrm{cos\θ}}_i+{({\mathrm{\ϵ}}_{VV}^e-{\sin }^2{\mathrm{\θ}}_i)}^{1/2}}=\frac{F_{VV}^{01}+F_{VV}^{12}{e}^{j\left(2k_{z1}h\right)}}{1+F_{VV}^{01}{F}_{VV}^{12}{e}^{j\left(2k_{z1}h\right)}}---\left(13\right)$