[发明专利]一种基于环状交织阵列的悬浮物动态监测方法与装置

权利要求书

1.一种基于环状交织阵列的悬浮物动态监测方法，其特征在于，包括以下步骤：

S1、发送测量信号；

测量装置相对于地面处于静止状态，在含有悬浮物的水中竖直向下发送测量信号s(t)：

式中，A是信号的幅度，W是信号的单边带宽，f₀是信号的中心频率，其功率谱图在正半轴上为幅度为0.5A、带宽为2W、中心频率为f₀的矩形脉冲：

f₀-Wff₀+W；

每隔周期T发送一次，即发送频率为F＝1/T，发送时长为T₀，T₀T；

s(t)频率范围应包含在所用发射换能器和接收换能器的频率范围中，且发射换能器正下方有一个反射体，以保证接收换能器能充分接收到反射信号；

S2、计算信号到达反射体时的频率；

选取地面作为参考系，以超声发射探头S的位置为坐标原点，垂直地面向下为z方向建立空间直角坐标系，则S的坐标为(0,0,0)，反射体G的坐标为(0,0,z)，相应的多个接收换能器均和超声发射探头S在同一平面上；

发射换能器S向z轴正方向，即方向发射测量信号，发送频率为F；z方向的单位向量可表示为设水流速度为其中v_wx，v_wy，v_wz分别代表水流速度在坐标轴x、y、z三个方向上的分量，则水流速度在S→G方向的速度为：

对于水中的悬浮物，由于悬浮物体积和质量都较小，在水中的流速主要与水流速度有关，因此估计悬浮物的沉降速度或流动速度可转换为计算悬浮水体的流速；假设声波在水下的传播速度为c，由于信号发射点S静止，根据多普勒效应可得，信号到达反射体时的频率为：

S3、计算信号到达接收探头时的频率；

L表示测量信号从发射换能器S出发传至反射体G最后传至接收换能器R_i的距离，为测量信号从反射体G传至接收换能器R_i的传播方向，α_i、β_i、γ_i分别为向量与坐标轴x、y、z的夹角，可得到向量的单位向量为则水流速度在G→R_i方向上的速度为：

声波由反射体G至接收点R_i的过程中，受到水流的作用，因此测量信号会受到多普勒效应；对于测量信号s(t)，在发送端每隔周期T发送一次，信号在水中传播的过程中由于多普勒效应产生多普勒频移，因此接收端R_i最终接收到的反射信号r_i(t)的周期为接收频率由于接收端没有运动，因此信号到达第i个接收换能器R_i时的频率为：

整理，得到测量信号s(t)从发送到接收的过程中，接收频率与发送频率之间的关系为：

其中，△_i为多普勒因子；

S4、需要M个接收换能器列出如上式所示的M个的等式，M≥3，即可求出上述v_wx、v_wy和v_wz三个未知数，由此求得悬浮水体的三维流速，即为悬浮物的速度；一个发射换能器和M个接收换能器组成一个收发阵列；

S5、估计悬浮物的浓度和线度分布；

利用信号功率谱对悬浮物浓度及线度分布做估计过程如下：

S5-1、采样；

由于测量信号s(t)在水中受到衰减，同时受多普勒效应的影响，因此对第i个接收换能器R_i接收到的信号r_i(t)，i＝1,2,..,M，采样后得到的信号为：

其中，B_i、△_i分别为接收信号的幅度和多普勒因子，w_i(n)为信号噪声，采样频率f_s＝1/T_s，其中f_s2f(f₀+W)，△_i可由步骤S3求出；

S5-2、对接收信号r_i(n)的自相关函数R_i(n)作2N点离散傅里叶变换，得：

S5-3、对该收发阵列的信号进行分析，发射换能器S的系统函数为H₁(f)，接收换能器的系统函数均为H₂(f)，同一个收发阵列中的接收换能器参数相同，且系统函数H₁(f)和H₂(f)可通过查看相应的换能器参数获得；

测量信号s(t)发射出去时先经过系统H₁(f)，同时，接收到的信号r_i(t)也先经过系统H₂(f)，因此在水中传播前和在水中传播过后的测量信号功率谱分别为H₁(f)S(f)和S_Ri(k)/H₂(f)；

S5-4、由步骤S5-3，测量信号在接收换能器R_i相应的f_sk/N频点上的离散信号功率谱为测量信号在发射换能器S相应的f_sk/N频点上的离散信号功率谱为由此得到衰减参数Q_k：

表示测量信号经过水体前后的衰减参数；

S5-5、对于M个接收换能器接收的信号，对信号相应区间离散功率谱进行叠加，根据衰减参数Q_k和传输距离L反演出悬浮物的浓度，Q_k的下标k＝0,1,…,N-1，而测量信号的每个频点f_sk/N则对应了四分之一波长线度的悬浮物的浓度。