1.一种水下输气管道泄漏源动态参数反演方法，其特征在于，所述方法包括：根据气泡速度点集和气泡半径点集，得到气泡粒子速度曲线和气泡速度半径曲线，并由速度微分方程(1)和半径微分方程(2)绘制速度和半径参数的理论曲线，对绘制的两个理论曲线分别进行参数拟合；

其中，vb为气泡粒子速度，vL为波浪速度，R为气泡半径，ρ1为常数，P为破浪场中广义压力，t为扩散时间，τ为时间，Cd为阻力系数， Re为雷诺系数，z是流体运动粘滞系数；

其中，ρg为气泡密度，R为气泡半径，N为普适气体常量与气体摩尔质量的比值，tg为水温，σ为液体表面张力系数，P为破浪场中广义压力，通过曲线匹配的方式对气泡受力情况进行推导，得到压强分布规律方程(3)；

kz

P(x,y,τ)＝‑gρ1+Aρ1kce cosk(x‑cτ)+Patm       (3)其中，c为波速，k为波数，H为波高，g为重力加速度，ρ1为常数，Patm为水面处标准大气压；

根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到水下输气管道泄漏源参数，所述水下输气管道泄漏源参数包括泄漏的气体达到水平面所需要的时间t和利用气体扩散水平迁移距离L中的一种及其组合。

2.根据权利要求1所述的方法，其特征在于，所述根据气泡速度点集和气泡半径点集，得到气泡粒子速度曲线和气泡速度半径曲线，包括：根据对应点集数量分别对所述气泡速度点集和气泡半径点集进行过滤，得到气泡粒子速度曲线和气泡速度半径曲线。

3.根据权利要求1所述的方法，其特征在于，对绘制的两个理论曲线分别采用最小二乘法进行参数拟合。

4.根据权利要求1所述的方法，其特征在于，根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到泄漏的气体达到水平面所需要的时间t和管道内压p、水流速度u、泄漏孔径d之间对关系如下式(4)所示：其中，v表示气体泄漏速率，H表示水深，p’表示环境水压。

5.根据权利要求1所述的方法，其特征在于，根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到利用气体扩散水平迁移距离L、泄漏孔径d、水流速度u和管道内压p的关系如下式(5)所示：其中，D表示水下管道直径，p’表示环境水压。

6.一种水下输气管道泄漏源动态参数反演装置，其特征在于，所述装置包括处理器，所述处理器被配置为：根据气泡速度点集和气泡半径点集，得到气泡粒子速度曲线和气泡速度半径曲线，并由速度微分方程(1)和半径微分方程(2)绘制速度和半径参数的理论曲线，对绘制的两个理论曲线分别进行参数拟合；

其中，vb为气泡粒子速度，vL为波浪速度，R为气泡半径，ρ1为常数，P为破浪场中广义压力，t为扩散时间，τ为时间，Cd为阻力系数， Re为雷诺系数，z是流体运动粘滞系数；；

其中，ρg为气泡密度，R为气泡半径，N为普适气体常量与气体摩尔质量的比值，tg为水温，σ为液体表面张力系数，P为破浪场中广义压力，通过曲线匹配的方式对气泡受力情况进行推导，得到压强分布规律方程(3)；

kz

P(x,y,τ)＝‑gρ1+Aρ1kce cosk(x‑cτ)+Patm                (3)其中，c为波速，k为波数，H为波高，g为重力加速度，ρ1为常数，Patm为水面处标准大气压；

根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到水下输气管道泄漏源参数，所述水下输气管道泄漏源参数包括泄漏的气体达到水平面所需要的时间t和利用气体扩散水平迁移距离L中的一种及其组合。

7.根据权利要求6所述的装置，其特征在于，所述处理器被进一步配置为根据对应点集数量分别对所述气泡速度点集和气泡半径点集进行过滤，得到气泡粒子速度曲线和气泡速度半径曲线。

8.根据权利要求6所述的装置，其特征在于，所述处理器被进一步配置为根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到泄漏的气体达到水平面所需要的时间t和管道内压p、水流速度u、泄漏孔径d之间对关系如下式(4)所示：其中，v表示气体泄漏速率，H表示水深，p’表示环境水压。

9.根据权利要求6所述的装置，其特征在于，所述处理器被进一步配置为根据压强分布规律方程(3)确定压力场的初始值，不断迭代修正，得到利用气体扩散水平迁移距离L、泄漏孔径d、水流速度u和管道内压p的关系如下式(5)所示：其中，D表示水下管道直径。

10.一种存储有指令的非暂时性计算机可读存储介质，当所述指令由处理器执行时，执行根据权利要求1至5中任一项所述的方法。