具体实施方式

下面结合附图和具体实施例，进一步阐明本发明，应理解这些实施例仅用于说明本发明而不用于限制本发明的范围，在阅读了本发明之后，本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。

本发明提供一种用于抑制钢轨漏磁检测提离干扰的滤波方法，包括如下步骤：

S1：如图1和图2所示，在被测钢轨轨头表面，磁轭的正下方处沿检测方向间隔排布两列磁敏传感器，分别用于测量x和z方向的磁场，两列各N个，两列磁敏传感器间隔的距离为l；

S2：对每个磁敏传感器进行采样，每个磁敏传感器的采样点数为M；

S3：对采样数据前后各补B个数据点，从实际采样数据的第一个点开始向后推进计算，取该点前后各B个采样点；

S4：以同一行分别检测x方向和z方向磁场的两个磁敏传感器为一组，共有N组磁敏传感器，对每组采样数据进行作差处理，找到差值的绝对值最大的点；

S5：将步骤S4获取的N组绝对值最大的点进行比较，找到其中绝对值最小的点所对应的一组通道；

S6：将步骤S5获取的通道所对应的采样点作为参考信号，将原信号和参考信号作差，得到滤波后的结果。

基于上述方案，这里结合现有方法，对本发明滤波方法进行如下分析：

由于伤损的漏磁场比较微弱，分布范围较小，单个磁敏传感器的检测范围有限。为了使检测范围覆盖整个钢轨轨头表面，现有方法是沿着y方向排布一排磁敏传感器，这样的话，磁敏传感器的输出不仅受伤损信号的影响，还会收到提离干扰的影响。由于伤损和提离的变化都会使磁敏传感器采样数据变化，容易造成误判，且提离带来的干扰难以通过一般的滤波方式消除。

由于磁敏传感器距离较近，提离变化引起的数据变化是类似的。钢轨轨头表面伤损一般分布在轨距角处或顶面中部，一般不会出现贯穿整个表面的伤损。即使这种伤损出现了，一般也不会垂直于x方向。所以一排磁敏传感器在通过伤损部位时，有些磁敏传感器能测到伤损的漏磁场，其输出数据会因为伤损发生变化；有些磁敏传感器不能测到伤损的漏磁场，其数据不因为伤损发生变化。检测到漏磁场的磁敏传感器的数据受到伤损漏磁场和提离变化的共同影响，而没有检测到漏磁场的磁敏传感器的数据只受到提离变化的影响。

由于提离变化和伤损漏磁场对磁敏传感器的共同影响有可能是正向叠加的，也有可能是负向抵消，所以对一排磁敏传感器的检测数据，单纯的取最小值作为参考信号是不可取的。根据x方向的磁敏传感器伤损信号和z方向磁敏传感器伤损信号的特征来考虑，由于x方向磁敏传感器信号是近似单峰的，且在伤损的正上方有最大值，而z方向磁敏传感器信号是双峰的。x方向的伤损信号的峰值和z方向信号伤损的负峰值是反向的。在伤损处，若一个传感器测量到x方向的漏磁信号的正峰值，另一个传感器测量到z方向的漏磁信号的负峰值，两个传感器输出信号求差之后绝对值将比单路信号的绝对值大。

在无伤损处，x方向和z方向的提离干扰信号是同向的。当提离变小时，x方向和z方向的漏磁信号都变大；当提离变大时，x方向和z方向的漏磁信号都变小。在无伤损处，若一个传感器测量到x方向漏磁信号，另一个传感器测量到z方向的漏磁信号，两个传感器输出信号求差之后绝对值将比单路信号的绝对值小。

显然，如果以分别测量x方向和z方向漏磁信号的两个传感器为一对，在y方向依次布设若干对传感器，如图1和图2所示，只要合理设置一对传感器的间距，其输出信号求差之后绝对值越大，没有测量到伤损的可能性越小；绝对值越小，没有测量到伤损的可能性越大。如前文所述，每次采样时，一排磁敏传感器中总有若干个没有测到伤损信号，可以认为输出信号求差之后绝对值最小的一对传感器没有测到伤损，以其输出作为不包含伤损信号，只包含提离干扰信号的参考信号，其它各对传感器的输出信号减去参考信号作为滤波后的结果。

为了保证找到合适的参考信号，一对传感器在伤损处输出信号求差之后绝对值应该尽可能大。如图5所示，若要求测出的最小伤损宽度为2a，深度为b，伤损上方某点P(x，z)的漏磁场强度为H(x，z)，其x方向分量为H_x(x，z)，z方向分量为H_z(x，z)，可分别由式(1)、(2)求得，伤损的变化趋势如图4所示。

σ_ms为伤损侧面的面磁荷密度，可由下式计算：

式中，μ为材料的磁导率，H为外加的磁场强度。

当x方向传感器测到正峰值时，若z方向传感器测到负峰值，一对传感器在伤损处输出信号求差之后绝对值最大，设为Q。

根据式(1)，若提离z＝z₀不变，当x＝0时H_x(x，z₀)有最大值。根据图4，H_z(x，z₀)的最小值出现在x的正半轴，对式(2)求导并令H′_z(x，z₀)＝0，求解得到其大于零的解l。显然，当两个传感器的间距为l时，当x方向传感器测到正峰值时，若z方向传感器测到负峰值。

在实际工程应用中根据要求测出的最小伤损宽度和深度确定两个传感器的间距l。不同尺寸的伤损漏磁场的分布如图6所示。间距l是根据最小伤损来计算的，而对于尺寸更大的伤损，x方向的正峰值更大，且当x方向传感器测到正峰值时，间距l的z方向传感器虽然没有测到负峰值，但是其信号比最小伤损的负峰值更小，也就是说一对传感器输出信号求差之后绝对值大于Q，数值更大，不会对参考信号的选择造成干扰。

从图4可见，伤损的漏磁场有一定的范围，为了准确的识别和重构伤损，希望尽可能完整的测量到漏磁场。但是在距离伤损较远的位置，漏磁场比较弱，难以从各种干扰中区分出来。以x方向漏磁场信号绝对值达到x方向峰值绝对值5％为伤损漏磁场的边界条件。根据式(1)，若提离z＝z₀不变，当x＝0时H_x(x，z₀)有最大值，即x方向的漏磁场在伤损中心正上方最大，最大值为H_x(0，z₀)。离伤损中心越远，H_x(x，z₀)越小。当H_x(x，z₀)减少到时，可以认为没有伤损x方向的漏磁场分布，即认为x方向的漏磁场分布在以伤损为中心，2L范围之内。根据式(1)，最小伤损x的边界L由下式计算：

若巡检速度为v米/秒，采样速度为s点/秒，则x方向传感器的第j₀个采样点如果在伤损的漏磁场边界之内，则在个采样点内必然包括x方向漏磁场的正峰值。由于x方向传感器在正峰值处，一对传感器输出信号求差之后绝对值最大，最有利于找出参考信号，所以在对j₀采样点进行滤波时，在[j₀-B，j₀+B]，个采样点内，找出一对传感器输出信号求差之后绝对值最大值MAXi(i＝1、2、……、N)，再找出MAXi最小的一对传感器，以其输出作为参考信号。

基于上述分析，本实施例对上述滤波方法进行详细说明，提供一种用于抑制钢轨漏磁检测提离干扰的滤波方法，如图3所示，其包括如下步骤：

步骤1：以同一行分别检测x方向和z方向磁场的两个磁敏传感器S_x[i]、S_z[i](i＝1、2、……、N)为一组，共有N组磁敏传感器；对每个磁敏传感器进行采样，每个磁敏传感器的采样点数为M，S_x[i]、S_z[i]的采样结果为数组：S_x[i，j]、S_z[i，j](j＝1、2、……、M)；

两列磁敏传感器的间隔距离l是根据伤损最小尺寸、巡检速度和采样速度得到，具体是该间隔距离为通过方程H′_z(x，z₀)＝0大于零的解：

其中，要求测出的最小伤损宽度为2a，深度为b，传感器没有发生振动时的提离为z₀，G_ms为伤损侧面的面磁荷密度。

步骤2：计算数组R_x[i，k]、R_z[i，k](k＝1、2、……、M+2B)：

巡检速度为v米/秒，采样速度为s点/秒，L由下式计算：

如图5所示，z₀为钢轨探伤车静止时，磁敏传感器与钢轨表面的提离，a为1/2的应检最小损伤宽度，b为应检最小损伤深度，所以要求测出的最小伤损宽度为2a，深度为b，传感器没有发生振动时的提离为z₀。

设置循环变量q，初值为B+1。

若检测探头的提离z₀＝1mm，a＝1mm，b＝1mm，则该伤损x和z方向的磁场分布如图4所示。显然，若提离z＝z₀不变，当x＝0时H_x(x，z)有最大值，即x方向的漏磁场在伤损中心正上方最大，最大值为H_x(0，z₀)。离伤损中心越远，H_x(x，z)越小。当H_x(x，z)减少到时，可以认为没有伤损x方向的漏磁场分布，即认为x方向的漏磁场分布在2L之内。

步骤3：求

|R_x[i，q-B]-R_z[i，q-B]|、|R_x[i，q-B+1]-R_z[i，q-B+1]|……、|R_x[i，q+B]-R_z[i，q+B]|中的最大值，记为MAX_i。

步骤4：找出MAX₁、MAX₂、……、MAX_N中的最小值，记录最小值的下标i₀；

步骤5：S_x[i，q-B]＝R_x[i，q]-R_x[i₀，q]、s_z[i，q-B]＝R_z[i，q]-R_z[i₀，q]，i＝1、2、……、N；

步骤6：q＝q+1，若q≤M+B，转步骤3，否则执行步骤7；

步骤7：以S_x[i，j]、S_z[i，j]，(i＝1、2、……、N，j＝1、2、……、M)为滤波后的结果。

本实施例中将上述滤波方法进行实例应用，具体如下：

设某钢轨顶面伤损漏磁检测系统在垂直于列车行进方向上并排布设了16路x方向和z方向磁敏传感器，要求检出的最小伤损为a＝1mm、b＝1mm、z₀＝1mm，根据上式利用MATLAB工具可以求出L＝8mm，若v＝1m/s，采样速率f＝10khz，则采样点

将16路磁敏传感器组输出的信号经过采样后保存，在采样数据的前后各补齐80个采样点，保证从第一个点到最后一个点，都能对其前后各80个数据进行计算；分别求每个传感器组的每个实际采样点前后80个x方向和z方向采样点的差值绝对值的最大值Maxi(i＝1，2，3……16)，将该最大值数组的最小值所对应的通道的采样点作为参考信号，将构成的参考信号作为噪声通道和原始数据进行自适应滤波，得出滤波后的数据。