一种海底输气管道微泄漏声学检测定位方法

1.本发明涉及一种海洋领域的水声信号处理方法，具体涉及一种海底输气管道微泄漏声学检测定位方法。


背景技术：

2.海底管道油气运输具有运输成本低、稳定安全的特点，是海洋油气资源长距离运输经济可靠的方式。长输管道是我国油气资源的主要输送方式，是保证国家能源和战略安全的有效手段。海底管道往往会穿越或临近公共航道、船舶作业区、临时锚地等通航环境复杂的水域。须定期或长期在线对水下结构进行检测，及时发现并处理微泄漏等故障隐患，保障水下生产系统的安全运行。
3.分布式光纤监测方法、序惯概率法以及次声波法等现有的管道泄漏检测方法只能检测管道存在明显泄漏的情况，并且难以做到实时泄漏监测定位，更无法实现长期对管道结构的输气管道微泄漏的检测定位。


技术实现要素：

4.本发明的目的在于克服现有技术难以实现对管道结构的输气管道微泄漏实时检测定位的问题，提出了一种海底输气管道微泄漏声学检测定位方法。
5.本发明提出一种海底输气管道微泄漏声学检测定位方法，该方法是一种基于高频高分辨率前视成像声纳的检测定位方法，该方法首先采用高分辨率算法获取气体微泄漏波束域数据，接着根据区域分块原则估计自适应门限阈值进行微泄漏动态检测，结合泄漏气泡的区域分布特征剔除杂点干扰，最后对数据进行微泄漏特征约束并大致判断气泡密集程度，检测出气泡微泄漏并对微泄漏位置进行定位。
6.本发明提出了一种海底输气管道微泄漏声学检测定位方法，该方法基于高频高分辨率前视成像声纳对海底管道微泄露位置进行检测定位，首先，采用高分辨率算法获取气体微泄漏波束域数据，然后，根据区域分块原则估计自适应门限阈值进行微泄漏动态检测，结合泄漏气泡的区域分布特征剔除杂点干扰，最后，对数据进行微泄漏特征约束并大致判断气泡密集程度，检测出气泡微泄漏并对微泄漏位置进行定位。
7.作为上述技术方案的改进之一，所述方法具体包括以下步骤：
8.1)首先，利用水下无人机器人搭载高频前视成像声纳，获取海底管道微泄漏声纳探测回波数据，利用常规波束形成方法形成多波束域数据，并采用高分辨率算法处理波束域数据以获取高清晰的泄漏气泡数据；
9.2)对于步骤1)获得的泄露气泡数据进行分块处理，分块的大小主要依据微泄漏气泡分布区域面积获取；求取每块区域的波束域数据幅值的均值及泄漏气泡分布密集区域，所述泄露气泡分布密集区域为分块区域幅值均值的最大值；并将分块区域幅值均值的最大值乘以常数以获得自适应门限阈值；
10.3)根据步骤2)获取的自适应门限阈值，进行疑似微泄漏气泡点的检测，并根据检
测点的分散程度剔除杂点干扰；
11.4)对于步骤3)中剔除杂点干扰的疑似微泄漏气泡点，进一步采用区域特征约束，根据前视声纳对微泄漏探测的特征，将疑似微泄漏气泡点约束在矩形区域内，并结合区域内气泡分布的密集程度，完成海底输气管道的微泄漏检测；
12.统计约束矩形区域内的疑似微泄漏气泡检测点的数据(这个数据是从上述步骤筛选出来的数据，可以理解为一定区域内疑似微泄露气泡检测点的个数或分布程度)mc，若mc大于设定的区域气泡密集程度门限值，则此区域为管道微泄漏区域，并得出微泄漏矩形框的中心点位置(θ0,r0)，其中，θ0为方位值，r0为距离值；
13.5)根据步骤4)获得的微泄漏矩形区域的方位值和距离值(θ0,r0)，并结合水下无人机器人的经纬度值，在地球坐标系下实现对于海底输气管道泄露的定位。
14.作为上述技术方案的改进之一，所述步骤2)中求取每块区域的波束域数据幅值均值的表达式为：
[0015][0016]
其中，bf(i,j)是第(i,j)点的波束域幅值，n
l1
表示每块区域方位向长度为l1对应的波束数目，n
l2
表示每块区域距离向长度距离向长度为l2对应的采样点数目；
[0017]
同时，获得自适应门限阈值的表达式为：
[0018]
c＝ηmax(p
i,j
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0019]
其中，c为自适应门限阈值；η为常数，是经验值，由统计得来的，用于调整门限阈值。
[0020]
作为上述技术方案的改进之一，所述步骤3)中根据检测点的分散程度剔除杂点干扰的表达式为：
[0021][0022]
其中，r(j)是疑似微泄露气泡检测点的距离向值，是超过步骤2)自适应门限阈值检测点的距离向均值，μ为比例系数，是经验值，由统计得来的，l2为每块区域距离向长度。
[0023]
作为上述技术方案的改进之一，所述步骤4)中将疑似微泄漏气泡点约束在矩形区域内的表达式为：
[0024][0025]
其中，θ(i)是疑似微泄漏气泡检测点的方位向值，是步骤3)剔除杂点干扰的疑似微泄漏气泡检测点的方位向均值，是步骤3)剔除杂点干扰的疑似微泄漏气泡检测点的距离向均值；l1为每块区域方位向长度。
[0026]
作为上述技术方案的改进之一，所述步骤5)中在地球坐标系下实现海底输气管道泄露定位的表达式为：
[0027]
long2＝long1 r0*sinθ0/[γ*cos(lat1)*2π/360]
ꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0028]
lat2＝lat1 r0*cosθ0/(γ*2π/360)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0029]
其中，*是相乘运算，(long1,lat1)为水下无人机器人的经纬度，(long2,lat2)为
海底管道微泄漏位置的经纬度，γ是地球的平均半径。
[0030]
本发明方法能够实现海底管道微泄漏的检测和定位，检测正确率较高，且计算量较低，较易工程实现。
附图说明
[0031]
图1为本发明所述方法流程图；
[0032]
图2为管道微泄漏声图；
[0033]
图3为采用本发明所述方法的管道微泄漏检测定位结果示意图。
具体实施方式
[0034]
以下结合实施例进一步说明本发明所提供的技术方案。
[0035]
本发明涉及一种海底管道微泄漏声学检测定位方法，是一种基于高频高分辨率前视成像声纳的检测定位方法，该方法首先采用高分辨率算法获取气体微泄漏波束域数据，接着根据区域分块原则估计自适应门限阈值进行微泄漏动态检测，结合泄漏气泡的区域分布特征剔除杂点干扰，最后，对数据进行微泄漏特征约束并大致判断气泡密集程度，检测出气泡微泄漏并对微泄漏位置定位。
[0036]
如图1所示，为本发明所述方法流程图，具体包括以下步骤：采用高分辨率算法获取到气体微泄漏波束域数据后，对其进行常规预成多波束形成处理后再进行高分辨率波束形成处理；对于处理的波束数据进行分块处理并求取每块区域的波束域数据幅值的均值，并根据分块均值的最大值进行自适应门限阈值估计；然后再对疑似泄露检测点进行杂点干扰剔除；对于剔除干扰点的疑似泄露检测点进行泄露区域特征约束，最后进行泄露区域定位。
[0037]
具体地，为了实现上述目的，本发明提供了一种海底管道微泄漏声学检测定位方法，所述方法包含：
[0038]
1)首先，利用水下无人机器人搭载高频前视成像声纳，获取海底管道微泄漏声纳探测回波数据，利用常规波束形成方法形成多波束域数据，波束数目为m，并采用高分辨率算法处理波束域数据以获取高清晰的泄漏气泡数据。
[0039]
2)对高分辨率算法处理之后的波束域数据进行分块处理，分块的大小主要依据微泄漏气泡分布区域面积获取，设每块区域方位向长度为l1，方位向分块数目为h1，距离向长度为l2，距离向分块数目为h2。求取每块区域的波束域数据幅值的均值，泄漏气泡分布密集区域，即分块均值最大值，并将分块均值最大值乘以常数获得自适应阈值。
[0040][0041]
c＝ηmax(p
i,j
)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0042]
其中，bf(i,j)是第(i,j)点的波束域幅值，n
l1
表示每块区域方位向长度为l1对应的波束数目，n
l2
表示每块区域距离向长度距离向长度为l2对应的采样点数目，η为一常数，用于调整门限阈值。
[0043]
3)根据步骤2)获取的自适应门限阈值，进行疑似微泄漏气泡的检测，并根据检测点的分散程度剔除杂点干扰。
[0044][0045]
其中，是步骤2)超过自适应门限阈值检测点的距离向均值，μ为比例系数。
[0046]
4)对步骤3)中剔除杂点干扰的疑似微泄漏气泡点，进一步采用区域特征约束，根据前视声纳对微泄漏探测的特征，将疑似微泄漏气泡点约束在矩形区域内，并结合区域内气泡分布的密集程度，完成海底输气管道的微泄漏检测。
[0047][0048]
其中，θ(i)是疑似微泄漏气泡检测点的方位向值，是步骤3)剔除杂点干扰的疑似微泄漏气泡检测点的方位向均值，r(j)是疑似微泄漏气泡检测点的距离向值，是步骤3)剔除杂点干扰的疑似微泄漏气泡检测点的距离向均值。
[0049]
统计满足公式(4)约束的疑似微泄漏气泡检测点的数据为mc，若mc大于设定的区域气泡密集程度门限值，则认为此区域为管道微泄漏区域，并得出微泄漏矩形框的中心点位置(θ0,r0)。
[0050]
5)根据步骤4)获得的微泄漏矩形区域的方位值和距离值(θ0,r0)，并结合水下无人潜器上的经纬度值，在地球坐标系下实现海底输气管道的定位。
[0051]
long2＝long1 r0*sinθ0/[γ*cos(lat1)*2π/360]
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0052]
lat2＝lat1 r0*cosθ0/(γ*2π/360)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(6)
[0053]
其中，(long1,lat1)为载体的经纬度，(long2,lat2)为海底管道微泄漏位置的经纬度，γ是地球的平均半径。
[0054]
在本实施例中，系统参数为：采用频率为900khz的高频前视声呐探测管微泄漏，进行海底输气管道消声水池模拟试验，前视声呐布放于水下约3米，输气管道布放于水下约4米，管道的直径为15.24cm，管道的泄漏孔直径为0.5mm，管道内压力为3mpa，只处理一帧数据，处理结果如下。
[0055]
如图2所示，为管道微泄漏声图；
[0056]
如图3所示，为采用本发明所述方法的管道微泄漏检测定位结果示意图(图中，*为检测出的疑似微泄漏气泡，矩形框为泄漏定位区域)。由图2和图3可知，本发明方法检测的微泄露定位区域能够准确覆盖管道微泄露区域。
[0057]
从上述对本发明的具体描述可以看出，本发明方法能够实现海底管道微泄漏的检测和定位，检测正确率较高，且计算量较低，较易工程实现。
[0058]
最后所应说明的是，以上实施例仅用以说明本发明的技术方案而非限制。尽管参照实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，对本发明的技术方案进行修改或者等同替换，都不脱离本发明技术方案的精神和范围，其均应涵盖在本发明的权利要求范围当中。