一种基于BOTDR技术的热力耦合解耦方法与流程

本发明涉及一种基于botdr技术的未知环境液体密度计算方法，属于机器分布式光纤检测领域中的工程问题。

背景技术：

光纤传感技术是20世纪70年代伴随着光纤技术和光纤通信技术的发展而兴起的一种新型传感技术。它以光波为传感信号，以光纤为传输介质，感知和探测外界被测信号，在传感方式、传感原理以及信号的探测与处理等方面都与传统的电学传感器有很大差异。光纤本身不带电、体积小、质量轻、易弯曲、抗电磁干扰、抗辐射性好，特别适合在易燃、易爆、空间受严格限制及强电磁干扰等恶劣环境下使用。因此，光纤传感技术一经问世就受到了极大重视，在各个重要领域得到了研究和应用。

光纤传感技术的种类繁多，有多种分类方法。往往同一种被测参量可以用不同类型的传感器测量，而同一原理的传感器又可以测量多种物理量。按照光纤的感知范围分，光纤传感种类可以分为单点式光纤传感器、多点式光纤传感器和全分布式光纤传感器。有些研究对象往往不是一个点或者几个点，而是呈现一定空间分布的场，如温度场、应力场、振动场等，这一类被测对象不仅涉及距离长、范围广，而且呈三维空间连续性分布，此时单点式甚至多点准分布式传感已经无法胜任传感检测，全分布式光纤传感系统应运而生。

热力耦合过程是应力场与温度场两个物理场之间相互影响的过程，即温度对受力变形有影响，同时受力变形对温度变化也有影响。因为热力耦合影响，光纤本身随温度变化，应变和应力对应的关系会发生微小的变化，而本发明就是解耦这种微小变化，使应变和应力关系更加准确。

技术实现要素：

为了克服上述热力耦合现象所导致的现有技术不足，本发明提供了一种基于botdr技术的热力耦合解耦方法，符合温度变化下的botdr技术的“应变-应力”测量。

在全分布式光纤传感系统中，光纤既作为信号传输介质，又是传感单元。即它将整根光纤作为传感单元，传感点是连续分布的，也有人称其为海量传感头，因此该传感方法可以测量光纤沿线任意位置处的信息。根据被测光信号的不同，全分布式光纤传感技术可以分为基于光纤中的瑞利散射、拉曼散射和布里渊散射三种类型；根据信号分析方法，可以分为基于时域和基于频域的全分布式光纤传感技术。截至目前，在应变/温度检测中比较成熟的技术是botdr技术，botdr技术的全称是“brillouinopticaltimedomainreflection”，其中文名称是“布里渊光时域反射”，它是通过布里渊散射的方法检测时域上的光信号，最终达到检测应变/温度的目的。

本发明提供一种基于botdr技术的热力耦合解耦方法，主要包括如下几个步骤：

步骤一、光纤标定；

步骤二、温度和应变关系的分离；

步骤三、应变和应力的关系建立；

步骤四、热力耦合的解耦。

本发明的优点在于：

1.本发明利用botdr技术具有良好的适应性，可以分离光纤温度和应变的关系。

2.本发明在botdr技术的基础上，不仅建立了中心频率和温度/应变之间的关系，同时还建立了应变和应力之间的关系。

3.本发明为了力求应力计算的准确，引入了热力耦合解耦的方法，使应力求解更加准确。

附图说明

图1为本发明中算法系统框图；

图2为本发明中应力/应变曲线模型；

图3为本发明中rbf神经网络拓扑结构图；

图4为本发明中光纤实物图；

图5为本发明中7020设备实物图；

图6为本发明中不同温度下拉力和应变的关系图。

具体实施方式

下面结合附图和实施例对本发明进行详细说明。

本发明提供的一种基于botdr技术的未知环境液体密度计算方法，系统框图如图1所示，具体包括步骤如下：

步骤一、光纤标定；

根据botdr技术的原理，可以很容易得出如下公式：

δf布里渊＝a11·δt+a12·δε

其中δf布里渊是指布里渊散射光中心频率的偏移，δt是指温度的变化，δε是指应变的变化，a11是指温度变化对应的中心频率的变化，a12是指应变变化对应的中心频率的变化。

在实际使用过程中，辅助有温度计、应变表等设备，需要对使用的光纤进行标定，使a11和a12两个系数更加准确。

步骤二、温度和应变关系的分离；

因为在botdr技术中，布里渊中心频率的偏移是伴随着温度和应变同时变化的，如果只需要应变关系的话，需要对温度和应变两个指标进行一下分离。本发明采用同一根光纤中，利用不同的方法进行处理，从而达到最终分离的效果。

在第一步布里渊检测的基础上，利用瑞利的方法进行测量，公式如下：

δf瑞利＝a21·δt+a22·δε

其中δf瑞利是指瑞利散射光频率的偏移，瑞利散光的频率与入射光的频率相等，δt是指温度的变化，δε是指应变的变化，a11是指温度变化对应的频率的变化，a12是指应变变化对应的频率的变化。

根据布里渊和瑞利的检测方法，可以求解二元一次方程组，从而达到温度和应变分离的效果。

步骤三、应变和应力的关系建立；

现行的分布式光纤传感技术，可以测量的应变范围在-3000με至+4000με以内。因此本发明只取0～4000με微应变进行研究。采用应变表和拉力计同时作用，每100个微应变记录一次拉力数据，数据点如图2所示。根据图2的数据点，拟合曲线等到公式如下：

f拉＝a·δε^1.2

步骤四、温度和应变之间的相互耦合作用。

因为温度的变化，导致光纤本身分子结构间作用力会产生微小的变化，因此会导致不同温度下，同样的应变发生的作用力会有变化。本发明称这种现象为热力耦合现象，记做δγ＝f(δt,δε)，这种关系并没有典型的数学公式，因此采用rbf神经网络作为训练模型进行处理。

rbf网络的结构与多层前向网络类似，它是一种三层前向网络。输入层由信号源结点组成；第二层为隐含层，隐单元数视所描述问题的需要而定，隐单元的变换函数是rbf径向基函数，它是对中心点径向对称且衰减的非负非线性函数；第三层为输出层，它对输入模式的作用作出响应。从输人空间到隐含层空间的变换是非线性的，而从隐含层空间到输出层空间变换是线性的。rbf神经网络的拓扑结构如图3所示。

实施例一

步骤一、光纤标定；

选取光纤型号为双芯单模的光纤，如图4所示。采用日本光纳株式社的nbx-7020设备，如图5所示。

通过标定可以知道，基于布里渊散射原理的光纤的“中心频率-应变/温度”系数。如下所示：

a11＝1.07mhz/℃,a12＝0.0497mhz/με

步骤二、温度和应变关系的分离；

通过标定可以知道，基于瑞利散射原理的光纤的“光频率-应变/温度”系数。如下所示：

a21＝-1.379ghz/℃,a22＝-0.1542ghz/με

根据布里渊和瑞利的原理联立组成方程，如下所示：

在光纤末端挂上重物前后，如图3所示，浸入液体前后分别检测一回，在这里温度可以认为是不变的，应变换算公式如下，求解方程组得到，浸入液体前的应变为δε1＝495με，浸入液体后的应变为δε2＝378με。

步骤三、应变和应力的关系建立；

取1米长的光纤在25℃下，固定光纤的一端，在另一端从小到大施加拉力，分别用拉力计和应变表记录数据，得到的数据表格如下：

表1应变和应力关系对应表

对上表的40个点进行描点绘图如图2所示，进行曲线拟合后，得到每个单位长度光纤上应变与应力的曲线方程为：

f拉′＝0.38·δε^1.2

受应变整个范围内，单位光纤的应力得出后，计算出总的受力。

f拉＝f拉′·l

步骤四、温度和应变之间的相互耦合作用。

根据表1数据格式，利用恒温箱环境，分别在20℃，25℃，30℃，35℃，40℃，绘制“应变/应力”曲线如图6所示。根据图上证明，当温度越高时光纤发生同样的应变，所受到的应力越小。经过多次试验统计热力耦合参数如下表所示：

表2