CN105891434A - 一种海水盐度的在线检测方法及其装置 - Google Patents

Abstract

本发明提供了海水盐度的在线检测方法及其装置，用扫频同步信号控制扫频激光光源，使其所输出调频光波的波长是一个周期性的锯齿波信号，将该调频光波分成两路，分别送至海水中的折射率探头和温度探头；该折射率探头为干涉仪结构，其返回的干涉光强信号的频率值与海水折射率相关，通过对该干涉光强信号做离散傅里叶变换，计算出海水折射率；该温度探头内置光纤Bragg光栅，其反射谱Bragg波长与海水温度相关，对上述扫频同步信号和光纤Bragg光栅反射光强信号进行同步的离散抽样，按照光栅温度传感解调算法，计算出海水温度值；根据所得的折射率、温度值及调频光波的平均波长，通过求解经验方程，得到被测海水的盐度值，从而实现海水盐度的在线检测。

Description

一种海水盐度的在线检测方法及其装置

技术领域

本发明涉及海水盐度检测技术，特别涉及一种海水盐度在线检测的方法及装置。

背景技术

海水盐度是海水中含盐量的一个标度，盐度是反映海水的物理过程和化学过程的基本参数。海水盐度的实时在线检测，在海洋环境保护、海洋科学、海洋工程和军事国防等领域有重要意义。近年来，关于海水盐度的检测技术得到了的广泛重视。

中国发明专利“棱镜模型多次折射的海水盐度检测装置”(专利申请号：201010603445.0)，利用激光器输出650nm波长的红光，照射参考液体和被测液体，并经棱镜多次折射后，红光被位置敏感器件PSD接收，根据PSD的位置变化信息计算出液体的盐度值。由于PSD是光电器件，需要提供电源才能工作，属于海水测量端(俗称湿端)带电的方法，不适合长时间的海上在线监测。

中国发明专利“高精度海水盐度测量仪”(专利申请号：201210244182.8)，分别将被测海水和标准海水盛入两个不同的恒温槽中。通过正弦波发生器、高精确标准电阻及电压变换器等部件测量出海水样品和标准海水的电导率，根据海水样品与标准海水的电导率之比，换算出海水的盐度值。由于该方法需要将正弦波电压加到恒温槽中的电极之上，且恒温槽本身要加电工作，也是海水测量端(湿端)带电的方法。仅适用实验室内检测，不适用海上在线监测。

中国发明专利“海水盐度的检测装置与方法”(专利申请号：201410425894.9)，其利用微纳光纤的强倏逝场与海水相互作用，在微纳光纤环形腔中产生谐振，以谐振波长按公式计算出海水的盐度值。该专利技术的检测原理先进，有好的应用前景。但需要用到光谱分析仪，以扫描方式获得被测海水的光谱图，读出被测海水的谐振峰波长。因此，该专利技术实施成本高，设备体积较大。

中国发明专利“利用布里渊频移和线宽同步反演海水温度和盐度的方法”(专利申请号：201410386610.X)。其仅给出了在一定海水温度和盐度条件下，布里渊散射频移和线宽的计算公式与反演方法，并未涉及海水温度和盐度测量方法与装置设计。

综上所述，现有的海水盐度检测技术存在多种技术缺陷，尤其不适用长时间的海上在线监测。

发明内容

为了解决现有技术中问题，本发明提供了一种海水盐度的在线检测方法及其装置，用扫频同步信号控制扫频激光光源，使其所输出调频光波的波长是一个周期性的锯齿波信号，将该调频光波分成两路，分别送至海水中的折射率探头和温度探头；该折射率探头为干涉仪结构，其返回的干涉光强信号的频率值与海水折射率相关，通过对该干涉光强信号做离散傅里叶变换，计算出海水折射率；该温度探头内置光纤Bragg光栅，其反射谱Bragg波长与海水温度相关，对上述扫频同步信号和光纤Bragg光栅反射光强信号进行同步的离散抽样，按照光栅温度传感解调算法，计算出海水温度值；根据所得的折射率、温度值及调频光波的平均波长，通过求解经验方程，得到被测海水的盐度值，从而实现海水盐度的在线检测。光纤Bragg光栅的英文：Fiber Bragg Grating，文献中常将光纤Bragg光栅，简称为FBG。

本发明的有益效果是：

一、采用波长以周期性的锯齿波变化的调频光波，能有效避免相位衰落现象。该折射率检测探头为干涉仪结构，干涉仪的传感臂由被测海水样本与反射镜构成，其参考臂由已知折射率的介质与反射镜构成，且传感臂与参考臂的长度相等，进入折射率检测探头的调频光波，在干涉仪的作用下，形成传感光和参考光，二者因干涉效应，产生干涉光，并由折射率检测探头返回。

根据双光束相干原理，该折射率检测探头返回的干涉光强度信号I(t)为：

$I\left(t\right)=I_S+I_R+2\sqrt{I_S{I}_R}cos\left(\mathrm{\Δ}\mathrm{\φ}\right)---\left(1\right)$

其中，t为时间变量，I_S、I_R分别为干涉仪的传感光强和参考光强，I_S≥0，I_R≥0，Δφ为干涉仪的传感光与参考光之间的相位差。因为传感臂与参考臂的长度相等，故相位差Δφ为：

$\mathrm{\Δ}\mathrm{\φ}=\frac{4\mathrm{\π}l(n_S-n_R)}{\mathrm{\λ}\left(t\right)}---\left(2\right)$

其中，l为传感臂和参考臂的长度，n_S、n_R分别为海水的待测折射率、参考臂介质的已知折射率，λ(t)为激光光源输出光波的波长。

如果激光光源输出固定波长的光波，即λ(t)＝λ₀，λ₀为常数，根据式(1)、(2)，则，折射率检测探头返回的干涉光强度I(t)为：

$I\left(t\right)=I_S+I_R+2\sqrt{I_S{I}_R}cos\[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}\]---\left(3\right)$

其中，忽略I_S、I_R、n_R随时间的变化，则干涉光强度I(t)对时间的微分

$\frac{dI\left(t\right)}{dt}=-2\sqrt{I_S{I}_R}\left\{\sin \[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}\]\right\}\frac{4\mathrm{\π}l}{{\mathrm{\λ}}_0}\frac{{\mathrm{dn}}_S}{dt}---\left(4\right)$

其中，为海水折射率对时间的微分。

由(4)式可知，

当时，其中m为整数，有下式成立：

$\frac{dI\left(t\right)}{dt}=0$

以上意味着，当激光光源输出光波为固定波长λ₀，且当时，即使海水折射率n_S的随时间随变化，即干涉光强度I(t)对时间的微分也为零，I(t)完全不随海水折射率n_S而变化。此时，折射率检测探头工作在最不灵敏的区域，出现所谓的相位衰落现象。因此，必须采用技术手段，才能避免相位衰落现象的影响，准确地实现物理量的传感。

采用扫频激光光源输出调频光波，波长λ(t)是一个扫频周期为T的锯齿波信号，具体形式如下：

$\mathrm{\λ}\left(t\right)=\left\{\begin{array}{cc}{\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}(t-nT)& nT\≤t\≤nT+\frac{T}{2}\\ {\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}(nT+T-t)& nT+\frac{T}{2}\≤t\≤(n+1)T\end{array}\right.---\left(5\right)$

其中，t为时间变量，n为非负整数，λ₀为波长λ(t)在每个扫频周期T的起始或终止时刻，即t＝nT或t＝(n+1)T对应的波长值，且λ₀为扫频过程中的最小波长，Δλ＞0，Δλ为波长的扫频变化范围，λ₀+Δλ为λ(t)在每个扫频周期的中间时刻，即对应的波长值，且λ₀+Δλ为扫频过程中的最大波长。

$\frac{1}{\mathrm{\λ}\left(t\right)}=\left\{\begin{array}{cc}\frac{1}{{\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}(t-nT)}& nT\≤t\≤nT+\frac{T}{2}\\ \frac{1}{{\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}(nT+T-t)}& nT+\frac{T}{2}\≤t\≤(n+1)T\end{array}\right.---\left(6\right)$

$\frac{1}{\mathrm{\λ}\left(t\right)}=\{\begin{array}{cc}\frac{1}{{\mathrm{\λ}}_0}\[\frac{1}{1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(t-nT\right)}\]& nT\≤t\≤nT+\frac{T}{2}\\ \frac{1}{{\mathrm{\λ}}_0}\[\frac{1}{1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(nT+T-t\right)}\]& nT+\frac{T}{2}\≤t\≤\left(n+1\right)T\end{array}---\left(7\right)$

如果Δλ＜＜λ₀，利用级数展开，忽略高阶小量：

$\frac{1}{\mathrm{\λ}\left(t\right)}=\{\begin{array}{cc}\frac{1}{{\mathrm{\λ}}_0}\[1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(t-nT\right)\]& nT\≤t\≤nT+\frac{T}{2}\\ \frac{1}{{\mathrm{\λ}}_0}\[1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(nT+T-t\right)\]& nT+\frac{T}{2}\≤t\≤\left(n+1\right)T\end{array}---\left(8\right)$

将式(8)代入式(2)：

$\mathrm{\Δ}\mathrm{\φ}=\left\{\begin{array}{cc}\frac{4\mathrm{\π}l\left(n_S-n_R\right)}{{\mathrm{\λ}}_0}\[1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(t-nT\right)\]& nT\≤t\≤nT+\frac{T}{2}\\ \frac{4\mathrm{\π}l\left(n_S-n_R\right)}{{\mathrm{\λ}}_0}\[1-\frac{2\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{T\λ}}_0}\left(nT+T-t\right)\]& nT+\frac{T}{2}\≤t\≤\left(n+1\right)T\end{array}\right.---\left(9\right)$

将式(9)代入式(1)：

$I\left(t\right)=\left\{\begin{array}{cc}{I}_S+I_R+2\sqrt{I_S{I}_R}\cos \[\frac{4\mathrm{\π}l\left(n_S-n_R\right)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}\left(t-nT\right)}{{\mathrm{T\λ}}_0})\]& nT\≤t\≤nT+\frac{T}{2}\\ I_S+I_R+2\sqrt{I_S{I}_R}\cos \[\frac{4\mathrm{\π}l\left(n_S-n_R\right)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}\left(nT-t\right)}{{\mathrm{T\λ}}_0})\]& nT+\frac{T}{2}\≤t\≤\left(n+1\right)T\end{array}\right.---\left(10\right)$

如果忽略I_S、I_R在一个扫频周期T内的变化，将一个扫频周期T内的I_S、I_R视为常数，则以上干涉光强度信号I(t)为直流分量I_DC与交流分量I_AC(t)之和：、I(t)＝I_DC+I_AC(t) (11)

其中，干涉光强度I(t)的直流分量I_DC：

I_DC＝I_S+I_R (12)

干涉光强度I(t)的交流分量I_AC(t)为：

$I_{AC}\left(t\right)=\left\{\begin{array}{cc}2\sqrt{I_S{I}_R}cos\[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(t-nT)}{{\mathrm{T\λ}}_0})\]& nT\≤t\≤nT+\frac{T}{2}\\ 2\sqrt{I_S{I}_R}\cos \[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(nT-t)}{{\mathrm{T\λ}}_0})\]& nT+\frac{T}{2}\≤t\≤(n+1)T\end{array}\right.---\left(13\right)$

上式的λ₀、Δλ、l及n_R为常数。海水折射率n_S一般变化较慢，n_S在一个扫频周期T内的也可视为常数，则一个扫频周期T内的交流分量I_AC(t)为单频信号，其频率值ω_s

${\mathrm{\ω}}_s=\frac{8\mathrm{\π}(n_S-n_R)l}{{\mathrm{\λ}}_0}\frac{\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{\λ}}_{0}T}---\left(14\right)$

可见，一个扫频周期T内的交流分量I_AC(t)为频率值ω_s与海水折射率n_S相关。

利用式(14)，对其变换，得到下式：

进而有：

$n_S=n_R+\frac{{\mathrm{\ω}}_s{{\mathrm{T\λ}}_0}^2}{8\mathrm{\π}\mathrm{\Δ}\mathrm{\λ}l}---\left(15\right)$

通过对干涉光强度信号I(t)做离散傅里叶变换，得到频率值ω_s，利用式(15)，可计算出海水折射率n_S。因此，采用扫频激光光源输出调频光波，可有效避免相位衰落现象的影响，准确地实现海水折射率的检测。

另外，以周期性的锯齿波变化的波长，是连续变化的，相对周期性的方波形式，其激光光源输出变化更平滑，性能更稳定。

二、采用波长以周期性的锯齿波变化的调频光波，方便光纤Bragg光栅实现温度传感解调。

本温度探头内置光纤Bragg光栅，其反射光谱的Bragg波长与海水温度值相关，对上述扫频同步信号和光纤Bragg光栅反射光强度信号进行同步的离散抽样，得到当前温度状态下的反射光谱，由该反射光谱中的最大值及所对应的扫频同步信号电压值，得到相应的Bragg波长，再根据光纤Bragg光栅温度传感器的特性参数，计算出当前的海水温度值。因此，采用波长以周期性的锯齿波变化的调频光波，方便光纤Bragg光栅实现温度传感解调。

三、通过求解经验方程，计算海水的盐度值。

根据所得的折射率、温度值及调频光波的平均波长，通过求解经验方程，得到被测海水的盐度值，从而实现海水盐度的在线检测。参考文献：Xiaohong Quan and EdwardS.Fry.Empirical equation for the index of refraction of seawater，APPLIEDOPTICS[J].1995.Vol.34,No.18:3477-3480。经验方程具体形式如下：

$n_S=n_0+(n_1+n_2{T}_S+n_3{T_S}^2)S+n_4{T_S}^2+\frac{n_5+n_{6}S+n_7{T}_S}{\stackrel{\‾}{\mathrm{\λ}}}+\frac{n_8}{{\stackrel{\‾}{\mathrm{\λ}}}^2}+\frac{n_9}{{\stackrel{\‾}{\mathrm{\λ}}}^3}---\left(16\right)$

其中，各系数为：

n₀＝1.31405，n₁＝1.779×10^-4，n₂＝-1.05×10^-6，n₃＝1.6×10^-8，n₄＝-2.02×10^-6，n₅＝15.868，n₆＝0.01155，n₇＝-0.00423，n₈＝-4382，n₉＝1.1455×10⁶。

以上经验方程将海水的折射率视为海水盐度值、温度值及光波的平均波长的函数，参考文献表明，其计算结果与实验数据吻合度好，准确度可靠。四、本发明属于海水测量端(湿端)完全不带电的方法。

本发明的一种海水盐度在线检测装置，其包括水上平台、海水折射率检测探头和海水温度检测探头。该水上平台需供电工作，而放置在海水中的两种检测探头，均由光学无源器件组成，无需供电。其与水上平台没有任何电气连接，仅通过两根光纤相连，属于海水测量端(湿端)完全不带电的方法。具有防漏电、防腐蚀、安装简单、维护方便的特点。更适合野外环境下的海水盐度在线检测。

附图说明

图1是本发明海水盐度在线检测装置结构示意图。

具体实施方式

下面结合附图对本发明做进一步说明。

实施方式1：

一种海水盐度在线检测方法，包括以下步骤：

步骤1：控制扫频同步信号源，使其输出的扫频同步信号V(t)是扫频周期为T的锯齿波电压信号，其具体形式如下：

$V\left(t\right)=\left\{\begin{array}{cc}{V}_0+\frac{2(V_P-V_0)}{T}(t-nT)& nT\≤t\≤nT+\frac{T}{2}\\ V_P+\frac{2(V_P-V_0)}{T}(nT+\frac{T}{2}-t)& nT+\frac{T}{2}\≤t\≤nT+T\end{array}\right.$

其中，t为时间变量，n为非负整数，n＝0 1 2,......，V₀为V(t)在每个扫频周期T的起始或终止时刻，即t＝nT或t＝(n+1)T对应的电压值，而V_P为V(t)在每个扫频周期的中间时刻，即对应的电压值，且V₀≥0，V_P＞V₀；

步骤2：在所述的扫频同步信号源的作用下，扫频激光光源所输出的调频光波，其波长λ(t)随扫频同步信号V(t)而线性变化，波长λ(t)也是一个扫频周期为T的锯齿波信号，具体形式如下：

$\mathrm{\λ}\left(t\right)=\left\{\begin{array}{cc}{\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}\left(t-nT\right)& nT\≤t\≤nT+\frac{T}{2}\\ {\mathrm{\λ}}_0+\frac{2\mathrm{\Δ}\mathrm{\λ}}{T}\left(nT+T-t\right)& nT+\frac{T}{2}\≤t\≤\left(n+1\right)T\end{array}\right.$

其中，λ₀为波长λ(t)在每个扫频周期T的起始或终止时刻，即t＝nT或t＝(n+1)T对应的波长值，且λ₀为扫频过程中的最小波长，Δλ＞0，Δλ为波长的扫频变化范围，而λ₀+Δλ为λ(t)在每个扫频周期的中间时刻，即对应的波长值，且λ₀+Δλ为扫频过程中的最大波长；

步骤3：扫频激光光源输出的调频光波被光波导器件分成两路光，分别进入放置在被测海水中的折射率检测探头和温度检测探头，该折射率检测探头为干涉仪结构，干涉仪的传感臂主要由被测海水样本与第一反射镜构成，干涉仪的参考臂主要由已知折射率的介质与第二反射镜构成，且传感臂与参考臂的长度相等，二者均为l，进入折射率检测探头的光，在干涉仪的作用下，形成传感光和参考光，二者因干涉效应，产生干涉光，该干涉光由折射率检测探头返回；所述的温度检测探头内置光纤Bragg光栅温度传感器，进入温度检测探头的调频光波，在该光纤Bragg光栅温度传感器的作用下，形成反射光及反射光谱，因该反射光谱的Bragg波长与海水温度值相关，将该反射光称为温度反射光；

所述的温度检测探头内置的光纤Bragg光栅温度传感器，其反射谱Bragg波长随海水温度的变化而移动，且在海水温度的整个变化范围内，该温度检测探头的Bragg波长的移动区间不超过扫频激光光源的扫频波长范围，即该Bragg波长在λ₀至λ₀+Δλ之间左右移动；

步骤4：采用光电探测器，测量上述干涉光强度信号I(t)和温度反射光强度信号G(t)，将整个测量时间，分成若干个长度相等的测量时间段，该时间段的长度取为2个扫频周期，即为2T，在每个时间段的第一个扫频周期T内进行干涉光强度I(t)的测量，第二个扫频周期T内进行温度反射光强度G(t)的测量；

其中的干涉光强度I(t)：

$I\left(t\right)=\left\{\begin{array}{cc}{I}_S+I_R+2\sqrt{I_S{I}_R}cos\[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(t-nT)}{{\mathrm{T\λ}}_0})\]& nT\≤t\≤nT+\frac{T}{2}\\ I_S+I_R+2\sqrt{I_S{I}_R}\cos \[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(nT-t)}{{\mathrm{T\λ}}_0})\]& nT+\frac{T}{2}\≤t\≤(n+1)T\end{array}\right.$

其中，I_S、I_R分别为干涉仪的传感光强和参考光强，如果忽略I_S、I_R在一个扫频周期T内的变化，将一个扫频周期T内的I_S、I_R视为常数，则以上干涉光强度I(t)为直流分量I_DC与交流分量I_AC(t)之和，即：

I(t)＝I_DC+I_AC(t)

其中，干涉光强度I(t)的直流分量I_DC：

I_DC＝I_S+I_R

干涉光强度I(t)的交流分量I_AC(t)为：

$I_{AC}\left(t\right)=\left\{\begin{array}{cc}2\sqrt{I_S{I}_R}cos\[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(t-nT)}{{\mathrm{T\λ}}_0})\]& nT\≤t\≤nT+\frac{T}{2}\\ 2\sqrt{I_S{I}_R}\cos \[\frac{4\mathrm{\π}l(n_S-n_R)}{{\mathrm{\λ}}_0}(1-\frac{2\mathrm{\Δ}\mathrm{\λ}(nT-t)}{{\mathrm{T\λ}}_0})\]& nT+\frac{T}{2}\≤t\≤(n+1)T\end{array}\right.$

上式的λ₀、Δλ、l及n_R为常数，海水折射率n_S一般变化较慢，n_S在一个扫频周期T内的也可视为常数，则该交流分量I_AC(t)为单频信号，其频率值ω_s

${\mathrm{\ω}}_s=\frac{8\mathrm{\π}(n_S-n_R)l}{{\mathrm{\λ}}_0}\frac{\mathrm{\Δ}\mathrm{\λ}}{{\mathrm{\λ}}_{0}T}.$

步骤5：先令n＝0，设置当前时间段的起始时刻为t＝nT；

步骤6：从t＝nT时刻开始，用A/D转换器，对当前时间段第一个扫频周期T内的干涉光强度I(t)信号，以采样周期为T₁，进行离散抽样，得到第一个扫频周期T内，即t＝nT至t＝(n+1)T时间的抽样信号序列I(m)，该序列的长度为L₁，m为序列序号且是非负整数，m＝0 12......L₁-1，将I(m)保存在信号处理单元中；

其中的采样周期T应满足采样定理的要求，即满足以下条件：

$T_1<\frac{2\mathrm{\&pi;}}{2{\mathrm{\&omega;}}_s}=\frac{{{\mathrm{\&lambda;}}_0}^{2}T}{8(n_S-n_R)l\mathrm{\&Delta;}\mathrm{\&lambda;}},$

且抽样信号序列I(m)的长度L₁为：

$L_1=\frac{T}{T_1}.$

步骤7：从t＝(n+1)T开始，用A/D转换器，对当前时间段第二个扫频周期T内的扫频同步信号V(t)和温度反射光强度G(t)，以采样周期为T₂，进行离散抽样，得到第二个扫频周期T内，即t＝(n+1)T至t＝(n+2)T时间的扫频同步信号序列V(h)、温度反射光强信号序列G(h)，序列V(h)和G(h)的长度均为L₂，h为序列序号且为非负整数，h＝0 1 2......L₂-1，将V(h)和G(h)保存在信号处理单元中；

其中的采样周期T₂为：

扫频同步信号序列V(h)和温度反射光强信号序列G(h)长度均为L₂，且L₂为：

步骤8：信号处理单元对当前时间段第一个扫频周期T内的干涉光强度抽样信号序列I(m)，做快速傅立叶变化FFT，得到干涉光强度I(t)在该时间内的频谱分布，由该频谱计算出以上时间内干涉光强度I(t)交流分量I_AC(t)的频率值ω_s；

步骤9：根据干涉光强度交流分量I_AC(t)的频率值ω_s与海水折射率n_S之间的关系，按下式计算出当前时间段第一个扫频周期T内，即t＝nT至t＝(n+1)T内的被测海水样本折射率n_S：

$n_S=n_R+\frac{{\mathrm{\&omega;}}_s{{\mathrm{T\&lambda;}}_0}^2}{8\mathrm{\&pi;}\mathrm{\&Delta;}\mathrm{\&lambda;}l}$

其中，l为折射率检测探头传感臂和参考臂的长度，n_R为该探头参考臂介质的已知折射率；

步骤10：利用当前时间段第二个扫频周期T内的扫频同步信号序列V(h)和温度反射光强度信号序列G(h)，按照光纤Bragg光栅温度传感解调算法，计算得到第二个扫频周期T内，即t＝(n+1)T至t＝(n+2)T内的被测海水温度值T_S；

其中的光纤Bragg光栅温度传感解调算法，其包括以下步骤：

第一：由温度反射光强度信号序列G(h)，找到温度反射光强的最大值G(h_M)及所对应的序号h_M；

第二：由以上温度反射光强的最大值G(h_M)所对应的序号h_M，找到位于该时刻的扫频同步信号电压值V(h_M)；

第三：由序号h_M所对应的扫频同步信号电压值V(h_M)，找到温度反射光强的最大值G(h_M)所对应的光纤光栅温度传感器的Bragg波长；

第四：根据光纤Bragg光栅温度传感器的特性参数，由以上得到的温度反射光强最大值G(h_M)所对应的Bragg波长，计算出当前的海水温度值T_S。

步骤11：由于海水的折射率n_S和温度的变化较为缓慢，可忽略海水的折射率n_S和温度T_S在每个测量时间段内的变化，所述每个测量时间段为两个扫频周期，用该时间段的第一个扫频周期T，即t＝nT至t＝(n+1)T内所得的海水样本的折射率n_S，近似为整个时间段，即t＝nT至t＝(n+2)T内的海水样本的折射率n_S，同理，用该时间段的第二个扫频周期T，即t＝(n+1)T至t＝(n+2)T内所得的海水的温度值T_S，近似为整个时间段，即t＝nT至t＝(n+2)T内的海水的温度值T_S；

步骤12：根据以上得到的当前时间段的海水折射率n_S、温度值T_S及扫频激光光源输出光波的平均波长 求解以下经验方程式：

$n_S=n_0+(n_1+n_2{T}_S+n_3{T_S}^2)S+n_4{T_S}^2+\frac{n_5+n_{6}S+n_7{T}_S}{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}+\frac{n_8}{{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}^2}+\frac{n_9}{{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}^3}$

计算出当前时间段，即t＝nT至t＝(n+2)T之间的海水盐度值S，其中，各系数为，

n₀＝1.31405，n₁＝1.779×10^-4，n₂＝-1.05×10^-6，n₃＝1.6×10^-8，n₄＝-2.02×10^-6，n₅＝15.868，n₆＝0.01155，n₇＝-0.00423，n₈＝-4382，n₉＝1.1455×10⁶；

步骤13：令n＝n+2，对时间段的起始时刻更新，指向下一时间段；

步骤14：重复步骤6-13，如此循环进行，测量在起始时刻t＝0之后的任意时间段，即t＝nT至t＝(n+2)T内的海水盐度值S，实现海水盐度的实时检测，其中，n＝0 12,......。

实施方式2：

实施方式1中的一种海水盐度在线检测方法所用的海水盐度在线检测装置，其包括水上平台、海水折射率检测探头和海水温度检测探头，水上平台通过两根光纤分别与海水折射率检测探头、海水温度检测探头相连，该水上平台输出调频光波，其波长是一个周期性的锯齿波信号，该调频光波分成两路，分别送至海水中的折射率探头和温度探头；该折射率探头为干涉仪结构，该温度探头内置光纤Bragg光栅。

所述水上平台包括扫频激光光源、扫频同步信号源、第一光电探测器、第二光电探测器、A/D转换器、信号处理单元、控制单元、第一1x2光纤耦合器、第一光纤环形器及第二光纤环形器，该海水温度检测探头包括一个光纤Bragg光栅温度传感器，所述的扫频激光光源设有光纤接口和电气接口，该扫频激光光源通过电气接口与扫频同步信号源相连，该扫频激光光源通过光纤接口与第一1x2光纤耦合器的1臂相连，所述的第一1x2光纤耦合器的2臂、3臂分别接第一光纤环形器的1臂、第二光纤环形器的1臂，所述的第一光纤环形器的2臂、3臂分别与折射率检测探头、第一光电探测器相连，所述的第二光纤环形器的2臂、3臂分别与海水温度检测探头、第二光电探测器相连，所述的扫频同步信号源、第一光电探测器及第二光电探测器均与A/D转换器相连，所述的A/D转换器相连与信号处理单元相连，所述的信号处理单元与控制单元相连，所述的控制单元还与扫频同步信号源相连。

该折射率检测探头包括第二1x2光纤耦合器、第一光纤自聚焦透镜、第二光纤自聚焦透镜、海水样本腔、参考介质、第一反射镜、第二反射镜，第二1x2光纤耦合器的1臂与所述的第一光纤环形器的2臂相连，第二1x2光纤耦合器的2臂、3臂分别与第一光纤自聚焦透镜、第二光纤自聚焦透镜相连，所述的海水样本腔位于第一光纤自聚焦透镜与第一反射镜之间，所述的参考介质位于第二光纤自聚焦透镜与第二反射镜之间

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (10)

1.一种海水盐度在线检测方法，其特征在于，包括以下步骤：

步骤1：控制扫频同步信号源，使其输出的扫频同步信号V(t)是扫频周期为T的锯齿波电压信号，其具体形式如下：

$V\left(t\right)=\left\{\begin{array}{cc}{V}_0+\frac{2(V_P-V_0)}{T}(t-nT)& nT\&le;t\&le;nT+\frac{T}{2}\\ V_P+\frac{2(V_P-V_0)}{T}(nT+\frac{T}{2}-t)& nT+\frac{T}{2}\&le;t\&le;nT+T\end{array}\right.$

其中，t为时间变量，n为非负整数，n＝0 1 2,......，V₀为V(t)在每个扫频周期T的起始或终止时刻，即t＝nT或t＝(n+1)T对应的电压值，而V_P为V(t)在每个扫频周期的中间时刻，即对应的电压值，且V₀≥0，V_P＞V₀；

步骤2：在所述的扫频同步信号源的作用下，扫频激光光源所输出的调频光波，其波长λ(t)随扫频同步信号V(t)而线性变化，波长λ(t)也是一个扫频周期为T的锯齿波信号，具体形式如下：

$\mathrm{\&lambda;}\left(t\right)=\left\{\begin{array}{cc}{\mathrm{\&lambda;}}_0+\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}}{T}(t-nT)& nT\&le;t\&le;nT+\frac{T}{2}\\ {\mathrm{\&lambda;}}_0+\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}}{T}(nT+T-t)& nT+\frac{T}{2}\&le;t\&le;(n+1)T\end{array}\right.$

其中，λ₀为波长λ(t)在每个扫频周期T的起始或终止时刻，即t＝nT或t＝(n+1)T对应的波长值，且λ₀为扫频过程中的最小波长，Δλ＞0，Δλ为波长的扫频变化范围，而λ₀+Δλ为λ(t)在每个扫频周期的中间时刻，即对应的波长值，且λ₀+Δλ为扫频过程中的最大波长；

步骤3：扫频激光光源输出的调频光波被光波导器件分成两路光，分别进入放置在被测海水中的折射率检测探头和温度检测探头，该折射率检测探头为干涉仪结构，干涉仪的传感臂主要由被测海水样本与第一反射镜构成，干涉仪的参考臂主要由已知折射率的参考介质与第二反射镜构成，且传感臂与参考臂的长度相等，二者均为l，进入折射率检测探头的光，在干涉仪的作用下，形成传感光和参考光，二者因干涉效应，产生干涉光，该干涉光由折射率检测探头返回；所述的温度检测探头内置光纤Bragg光栅温度传感器，进入温度检测探头的调频光波，在该光纤Bragg光栅温度传感器的作用下，形成反射光及反射光谱，因该反射光谱的Bragg波长与海水温度值相关，将该反射光称为温度反射光；

步骤4：采用光电探测器，测量上述干涉光强度信号I(t)和温度反射光强度信号G(t)，将整个测量时间，分成若干个长度相等的测量时间段，该时间段长度取为2个扫频周期，即为2T，在每个时间段的第一个扫频周期T内进行干涉光强度I(t)的测量，第二个扫频周期T内进行温度反射光强度G(t)的测量；

步骤5：先令n＝0，设置当前时间段的起始时刻为t＝nT；

步骤6：从t＝nT时刻开始，用A/D转换器，对当前时间段第一个扫频周期T内的干涉光强度I(t)信号，以采样周期为T₁，进行离散抽样，得到第一个扫频周期T内，即t＝nT至t＝(n+1)T时间的抽样信号序列I(m)，该序列的长度为L₁，m为序列序号且是非负整数，m＝0 12......L₁-1，将I(m)保存在信号处理单元中；

步骤7：从t＝(n+1)T开始，用A/D转换器，对当前时间段第二个扫频周期T内的扫频同步信号V(t)和温度反射光强度G(t)，以采样周期为T₂，进行离散抽样，得到第二个扫频周期T内，即t＝(n+1)T至t＝(n+2)T时间的扫频同步信号序列V(h)、温度反射光强信号序列G(h)，序列V(h)和G(h)的长度均为L₂，h为序列序号且为非负整数，h＝0 1 2......L₂-1，将V(h)和G(h)保存在信号处理单元中；

步骤8：信号处理单元对当前时间段第一个扫频周期T内的干涉光强度抽样信号序列I(m)，做快速傅立叶变化FFT，得到干涉光强度I(t)在该时间内的频谱分布，由该频谱计算出以上时间内干涉光强度I(t)交流分量I_AC(t)的频率值ω_s；

步骤9：根据干涉光强度交流分量I_AC(t)的频率值ω_s与海水折射率n_S之间的关系，按下式计算出当前时间段第一个扫频周期T内，即t＝nT至t＝(n+1)T内的被测海水样本折射率n_S：

$n_S=n_R+\frac{{\mathrm{\&omega;}}_s{{\mathrm{T\&lambda;}}_0}^2}{8\mathrm{\&pi;}\mathrm{\&Delta;}\mathrm{\&lambda;}l}$

其中，l为折射率检测探头传感臂和参考臂的长度，n_R为该探头参考介质的已知折射率；

步骤10：利用当前时间段第二个扫频周期T内的扫频同步信号序列V(h)和温度反射光强度信号序列G(h)，按照光纤Bragg光栅温度传感解调算法，计算得到第二个扫频周期T内，即t＝(n+1)T至t＝(n+2)T内的被测海水温度值T_S；

步骤11：由于海水的折射率n_S和温度T_S的变化较为缓慢，可忽略海水的折射率n_S和温度T_S在每个测量时间段内的变化，所述每个测量时间段为两个扫频周期，用该时间段的第一个扫频周期T，即t＝nT至t＝(n+1)T内所得的海水样本的折射率n_S，近似为整个时间段，即t＝nT至t＝(n+2)T内的海水样本的折射率n_S，同理，用该时间段的第二个扫频周期T，即t＝(n+1)T至t＝(n+2)T内所得的海水的温度值T_S，近似为整个时间段，即t＝nT至t＝(n+2)T内的海水的温度值T_S；

步骤12：根据以上得到的当前时间段的海水折射率n_S、温度值T_S及扫频激光光源输出光波的平均波长求解以下经验方程式：

$n_S=n_0+(n_1+n_2{T}_S+n_3{T_S}^2)S+n_4{T_S}^2+\frac{n_5+n_{6}S+n_7{T}_S}{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}+\frac{n_8}{{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}^2}+\frac{n_9}{{\stackrel{\&OverBar;}{\mathrm{\&lambda;}}}^3}$

计算出当前时间段，即t＝nT至t＝(n+2)T之间的海水盐度值S，其中，各系数为，

n₀＝1.31405，n₁＝1.779×10^-4，n₂＝-1.05×10^-6，n₃＝1.6×10^-8，

n₄＝-2.02×10^-6，n₅＝15.868，n₆＝0.01155，n₇＝-0.00423，

n₈＝-4382，n₉＝1.1455×10⁶；

步骤13：令n＝n+2，对时间段的起始时刻更新，指向下一时间段；

步骤14：重复步骤6-13，如此循环进行，测量在起始时刻t＝0之后的任意时间段，即t＝nT至t＝(n+2)T内的海水盐度值S，实现海水盐度的实时检测，其中，n＝0 1 2,......。

2.根据权利要求1所述的一种海水盐度在线检测方法，其特征在于：步骤3中所述的温度检测探头内置的光纤Bragg光栅温度传感器，其反射谱Bragg波长随海水温度的变化而移动，且在海水温度的整个变化范围内，该温度检测探头的Bragg波长的移动区间不超过扫频激光光源的扫频波长范围，即该Bragg波长在λ₀至λ₀+Δλ之间左右移动。

3.根据权利要求1所述的一种海水盐度在线检测方法，其特征在于：步骤4中的干涉光强度I(t)：

$I\left(t\right)=\left\{\begin{array}{cc}{I}_S+I_R+2\sqrt{I_S{I}_R}cos\&lsqb;\frac{4\mathrm{\&pi;}l(n_S-n_R)}{{\mathrm{\&lambda;}}_0}(1-\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}(t-nT)}{{\mathrm{T\&lambda;}}_0})\&rsqb;& nT\&le;t\&le;nT+\frac{T}{2}\\ I_S+I_R+2\sqrt{I_S{I}_R}\cos \&lsqb;\frac{4\mathrm{\&pi;}l(n_S-n_R)}{{\mathrm{\&lambda;}}_0}(1-\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}(nT-t)}{{\mathrm{T\&lambda;}}_0})\&rsqb;& nT+\frac{T}{2}\&le;t\&le;(n+1)T\end{array}\right.$

其中，I_S、I_R分别为干涉仪的传感光强和参考光强，如果忽略I_S、I_R在一个扫频周期T内的变化，将一个扫频周期T内的I_S、I_R视为常数，则以上干涉光强度I(t)为直流分量I_DC与交流分量I_AC(t)之和，即：

I(t)＝I_DC+I_AC(t)

其中，干涉光强度I(t)的直流分量I_DC：

I_DC＝I_S+I_R

干涉光强度I(t)的交流分量I_AC(t)为：

$I_{AC}\left(t\right)=\left\{\begin{array}{cc}2\sqrt{I_S{I}_R}cos\&lsqb;\frac{4\mathrm{\&pi;}l(n_S-n_R)}{{\mathrm{\&lambda;}}_0}(1-\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}(t-nT)}{{\mathrm{T\&lambda;}}_0})\&rsqb;& nT\&le;t\&le;nT+\frac{T}{2}\\ 2\sqrt{I_S{I}_R}\cos \&lsqb;\frac{4\mathrm{\&pi;}l(n_S-n_R)}{{\mathrm{\&lambda;}}_0}(1-\frac{2\mathrm{\&Delta;}\mathrm{\&lambda;}(nT-t)}{{\mathrm{T\&lambda;}}_0})\&rsqb;& nT+\frac{T}{2}\&le;t\&le;(n+1)T\end{array}\right.$

上式的λ₀、Δλ、l及n_R为常数，海水折射率n_S一般变化较慢，n_S在一个扫频周期T内的也可视为常数，则该交流分量I_AC(t)为单频信号，其频率值ω_s

${\mathrm{\&omega;}}_s=\frac{8\mathrm{\&pi;}(n_S-n_R)l}{{\mathrm{\&lambda;}}_0}\frac{\mathrm{\&Delta;}\mathrm{\&lambda;}}{{\mathrm{\&lambda;}}_{0}T}.$

4.根据权利要求1所述的一种海水盐度在线检测方法，其特征在于：步骤6中的采样周期T₁应满足采样定理的要求，即满足以下条件：

$T_1<\frac{2\mathrm{\&pi;}}{2{\mathrm{\&omega;}}_s}=\frac{{{\mathrm{\&lambda;}}_0}^{2}T}{8(n_S-n_R)l\mathrm{\&Delta;}\mathrm{\&lambda;}},$

且抽样信号序列I(m)的长度L₁为：

$L_1=\frac{T}{T_1}.$

5.根据权利要求1所述的一种海水盐度在线检测方法，其特征在于：

步骤7中的采样周期T₂为：

6.根据权利要求5所述的一种海水盐度在线检测方法，其特征在于：

步骤7中的扫频同步信号序列V(h)和温度反射光强信号序列G(h)长度均为L₂，且L₂为：

7.根据权利要求1所述的一种海水盐度在线检测方法，其特征在于：步骤10中的光纤Bragg光栅温度传感解调算法，其包括以下步骤：

第一：由温度反射光强度信号序列G(h)，找到温度反射光强的最大值G(h_M)及所对应的序号h_M；

第二：由以上温度反射光强的最大值G(h_M)所对应的序号h_M，找到位于该时刻的扫频同步信号电压值V(h_M)；

第三：由序号h_M所对应的扫频同步信号电压值V(h_M)，找到温度反射光强的最大值G(h_M)所对应的光纤光栅温度传感器的Bragg波长；

第四：根据光纤Bragg光栅温度传感器的特性参数，由以上得到的温度反射光强最大值G(h_M)所对应的Bragg波长，计算出当前的海水温度值T_S。

8.根据权利要求1至7任意一项所述的一种海水盐度在线检测方法所用的海水盐度在线检测装置，其特征在于：其包括水上平台、海水折射率检测探头和海水温度检测探头，水上平台通过两根光纤分别与海水折射率检测探头、海水温度检测探头相连，该水上平台输出调频光波，其波长是一个周期性的锯齿波信号，该调频光波分成两路，分别送至海水中的折射率探头和温度探头；该折射率探头为干涉仪结构，该温度探头内置光纤Bragg光栅。

9.根据权利要求8所述的一种海水盐度在线检测装置，其特征在于：所述水上平台包括扫频激光光源、扫频同步信号源、第一光电探测器、第二光电探测器、A/D转换器、信号处理单元、控制单元、第一1x2光纤耦合器、第一光纤环形器及第二光纤环形器，该海水温度检测探头包括一个光纤Bragg光栅温度传感器，所述的扫频激光光源设有光纤接口和电气接口，该扫频激光光源通过电气接口与扫频同步信号源相连，该扫频激光光源通过光纤接口与第一1x2光纤耦合器的1臂相连，所述的第一1x2光纤耦合器的2臂、3臂分别接第一光纤环形器的1臂、第二光纤环形器的1臂，所述的第一光纤环形器的2臂、3臂分别与折射率检测探头、第一光电探测器相连，所述的第二光纤环形器的2臂、3臂分别与海水温度检测探头、第二光电探测器相连，所述的扫频同步信号源、第一光电探测器及第二光电探测器均与A/D转换器相连，所述的A/D转换器相连与信号处理单元相连，所述的信号处理单元与控制单元相连，所述的控制单元还与扫频同步信号源相连。

10.根据权利要求8所述的一种海水盐度在线检测装置，其特征在于：该折射率检测探头包括第二1x2光纤耦合器、第一光纤自聚焦透镜、第二光纤自聚焦透镜、海水样本腔、参考介质、第一反射镜、第二反射镜，第二1x2光纤耦合器的1臂与所述的第一光纤环形器的2臂相连，第二1x2光纤耦合器的2臂、3臂分别与第一光纤自聚焦透镜、第二光纤自聚焦透镜相连，所述的海水样本腔位于第一光纤自聚焦透镜与第一反射镜之间，所述的参考介质位于第二光纤自聚焦透镜与第二反射镜之间。