1.一种海底管道偏移定位方法，其特征在于，包括以下步骤：S1、将组合式智能清管器组安装在海底管道的入口处，组合式智能清管器组上的声源发射装置通过自适应声波频率选择方法，根据海底管中的不同介质对其所发生的声波发射频率进行自适应调节；

S2、AUV跟随组合式智能清管器组移动，实时接收组合式智能清管器组发出的信号，通过加权融合算法确定AUV与组合式智能清管器组之间的信号传播距离 ；

S3、通过卡尔曼滤波算法更新AUV自身的位置估计 ，从而确定组合式智能清管器所在的海底管道的位置；

步骤S1中频率优化算法的具体实现过程如下所述：S1.1、环境数据的采集与实时反馈；

S1.2、根据监测到的环境数据，利用频率优化算法选择最佳频率，调节声源发射装置的声波发射频率；

S1.3、采用迭代的方式实时动态地找到最优频率；

步骤S1.2中，确定最佳声波发射频率范围fopt∈[300Hz,7kHz]；

声波在管道介质中的传播速度Vsignal的计算公式为：，

其中， 表示实时采集的介质温度值；  表示实时采集的介质压力值；

表示介质修正项，其取值为：

，

声波强度随距离d的衰减公式：，

其中，d表示实际测量点与声源之间的距离，单位为m； 表示初始信号强度，即声源处信号强度，单位为dB； 表示参考距离，用于归一化的参考标准；n(f)表示衰减因子，其公式为：n(f)=af+b，

当介质为石油时，a=0.02,b=1.1；当介质为天然气时，a=0.05,b=1.6；

α(f)表示介质对声波信号的吸收系数，单位为db/m，其计算公式为：；

定义综合性能指标J(f)，最大化信噪比 并最小化衰减：,

其中，dmax表示管道最大探测距离；λ表示权重系数，取值0.5；

信噪比SNR（f）的计算公式为： ，

其中， 表示环境噪声功率，单位为dB；

步骤S1.3中，通过梯度上升法实时求解最优频率：，

采用迭代的方式实时动态地求解最优频率，，

其中， 表示第k+1次迭代后的频率更新值； 表示当前第k次迭代时的频率值；为表示频率调整步长； 表示综合性能指标 在频率 处的梯度；

设定停止条件为：

；

其中， 表示在第k次迭代对应频率 时，计算出的综合性能指标； 表示下一次迭代后的频率 对应的性能指标值。

2.根据权利要求1所述的海底管道偏移定位方法，其特征在于，步骤S2中，组合式智能清管器组发射的声波信号在传播过程中因多路径效应产生回波信号，根据回波信号的时间延迟 ，计算基于传播时间延迟估算的AUV与组合式智能清管器组之间的距离 ：,

根据接收信号强度反推传播距离 的公式为：，

其中， 表示流速对信号衰减的修正系数；Vflow表示当前介质流速，单位为m/s；

定义时间延迟权重因子为 ：

，

定义信号强度权重因子为 ：

,

采用加权融合方法计算最终的信号传播距离 ：,

其中，Imax表示环境中接收到的最大信号强度。

3.根据权利要求1所述的海底管道偏移定位方法，其特征在于，步骤S3中，在确定组合式智能清管器以及所在的海底管道的位置时，以AUV自身的位置估计 为基准，通过步骤S2得到的dtotal与方向向量 ，计算清管器在海底三维空间中的位置(xp,yp,zp)：。

4.根据权利要求1所述的海底管道偏移定位方法，其特征在于，步骤S3中，通过卡尔曼滤波算法更新AUV自身的位置估计 的具体过程如下所述：根据AUV上的传感器获取当前AUV的状态向量 ，其中x,y,z表示AUV三维空间中的位置坐标， 表示速度分量；

状态预测与协方差更新公式为：，

，

其中， 表示预测的状态； 表示控制输入； 表示预测的协方差矩阵；A表示系统状态转移矩阵；B表示控制输入矩阵；uk表示控制输入量；Pk表示状态估计协方差矩阵；Q表示过程噪声协方差矩阵；

根据AUV的运动状态和AUV的环境变化，动态调整过程噪声矩阵Q：，

其中，表示调整系数；VfIow表示当前介质流速；Vmax表示预设的最大流速值；Q0表示静态参考过程噪声协方差矩阵；

卡尔曼增益Kk+1的计算公式为：，

更新状态估计的计算公式为：

，

其中，ZK+1表示测量值，由USBL系统、深度计、声学反射定位系统提供；H表示观测矩阵，R表示观测噪声协方差矩阵；

在更新后的状态向量中提取前三个分量，可获得AUV当前的实时三维位置估计：，

其中， 为更新状态向量的第一个元素，表示AUV当前时刻的x坐标；

为更新状态向量的第二个元素，表示AUV当前时刻的y坐标； 为更新状态向量的第三个元素，表示AUV当前时刻的z坐标。

5.一种自适应声源频率调节模块，其特征在于，该模块包括：控制模块，用于根据声源发射频率配置DDS信号生成模块的相位寄存器；与其他模块通信，进行信号控制和实时反馈；

DDS信号生成模块，用于调节频率以满足对于海底管道不同介质的适合发射频率,生成

300Hz到7kHz之间可调的频率信号；通过SPI接口设置DDS的频率寄存器，生成目标频率；通过相位累加器用于控制输出信号的频率；通过更新相位寄存器，准确控制频率；选择使用一个低频信号作为基准信号，利用DDS的频率调制功能生成其他频率；

信号滤波与放大模块，用于对DDS信号生成模块生成的信号进行滤波和增益调整，去除高次谐波并增强信号幅度；去除不必要的高频谐波；放大信号到所需的幅度；

功率放大模块，用于提升信号功率，使其能够驱动声学换能器；

输出模块，用于驱动声波信号的发射，采用声学换能器，通过声学换能器，将电信号转化为声波；

环境适应性模块；用于确保海底使用的可靠性；

自适应声源频率调节模块能够根据环境实时变化，精确选择300Hz至7kHz之间的最佳频率，声波在管道介质中的传播速度Vsignal的计算公式为：，

其中， 表示实时采集的介质温度值；  表示实时采集的介质压力值；

表示介质修正项，其取值为：

，

声波强度随距离d的衰减公式：，

其中，d表示实际测量点与声源之间的距离，单位为m； 表示初始信号强度，即声源处信号强度，单位为dB； 表示参考距离，用于归一化的参考标准；n(f)表示衰减因子，其公式为：n(f)=af+b，

当介质为石油时，a=0.02,b=1.1；当介质为天然气时，a=0.05,b=1.6；

α(f)表示介质对声波信号的吸收系数，单位为db/m，其计算公式为：；

定义综合性能指标J(f)，最大化信噪比 并最小化衰减：,

其中，dmax表示管道最大探测距离；λ表示权重系数，取值0.5；

信噪比SNR（f）的计算公式为：，

其中， 表示环境噪声功率，单位为dB；

通过梯度上升法实时求解最优频率：，

采用迭代的方式实时动态地求解最优频率，，

其中， 表示第k+1次迭代后的频率更新值； 表示当前第k次迭代时的频率值；为表示频率调整步长； 表示综合性能指标 在频率 处的梯度；

设定停止条件为：

；

其中， 表示在第k次迭代对应频率 时，计算出的综合性能指标； 表示下一次迭代后的频率 对应的性能指标值。