CN111610520A - 一种基于声学信号的气泡释放通量量化方法 - Google Patents

一种基于声学信号的气泡释放通量量化方法 Download PDF

Info

Publication number
CN111610520A
CN111610520A CN202010571291.5A CN202010571291A CN111610520A CN 111610520 A CN111610520 A CN 111610520A CN 202010571291 A CN202010571291 A CN 202010571291A CN 111610520 A CN111610520 A CN 111610520A
Authority
CN
China
Prior art keywords
bubbles
volume
bubble
acoustic
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202010571291.5A
Other languages
English (en)
Other versions
CN111610520B (zh
Inventor
武智
李捷
李新辉
朱书礼
杨计平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Pearl River Fisheries Research Institute CAFS
Original Assignee
Pearl River Fisheries Research Institute CAFS
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pearl River Fisheries Research Institute CAFS filed Critical Pearl River Fisheries Research Institute CAFS
Priority to CN202010571291.5A priority Critical patent/CN111610520B/zh
Publication of CN111610520A publication Critical patent/CN111610520A/zh
Application granted granted Critical
Publication of CN111610520B publication Critical patent/CN111610520B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88—Radar or analogous systems specially adapted for specific applications
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52077—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开了一种基于声学信号的气泡释放通量量化方法，包括以下步骤，步骤一，混响信号去除；步骤二，气泡声学特性获取；步骤三，气泡固有属性计算；步骤四，气泡释放通量计算；其中，步骤四的计算公式为：
Figure DDA0002549698800000011
公式中，ρ(mL/m3)为水中的气泡浓度，v(m/s)为气泡的上升速度，P为采样点实测大气压，P0为标准大气压(1013.25hPa)。所述步骤一包括以下步骤：第一步，观察回波图，确定气泡分布范围，通过设定阈值，将该范围外的噪声排除掉，实现首次去噪；第二步，手动定义无效数据区域，进行二次去噪，去除行船噪声和鱼类信号。本发明采用声学技术对气泡进行信号采集，能够在大水域实现连续监测和全面捕捉气泡排放信息，其效率高、覆盖面积大。

Description

一种基于声学信号的气泡释放通量量化方法
技术领域
本发明涉及声学探测技术领域，具体为一种基于声学信号的气泡释放通量量化方法。
背景技术
水体中的甲烷由有机物厌氧降解产生，并且主要以气泡形式逐渐上升，之后进入大气。我国淡水河流、湖泊、水库众多，植被丰富，大量有机质随水流进入水体；此外，人类活动导致水体富营养化，盛长的藻类不能通过食物链输出，沉积于水底。水底有机质分解，会排放大量的温室气体，导致全球气温不断升高，温室效应越来越显著，因此，对淡水水体气泡排放进行研究并量化，对于研究碳循环和碳管理等具有重要意义。
由于气泡释放过程中在宏观和微观尺度上的高度异质性和间歇性，监测并量化水中的气泡是非常困难的。传统的气泡释放统计一般采用气体收集、光学视频或者图片拍摄，如公开号为CN109598739A公开的一种水体中渗漏气泡数目的统计方法及处理终端就公开了通过光学视频拍摄进行气泡数目统计的方法，该方法将拍摄的光学视频的每一帧进行气泡个数统计，进而得到气泡释放通量，然而上述方法仅限于浅水及小范围水域研究，无法在大水面实现连续监测和全面捕捉气泡排放信息，并且其计算的方式繁琐，效率低下，目前仍无方法可以对大水面气泡释放进行量化。
发明内容
本发明针对上述问题，提供一种基于声学信号的气泡释放通量量化方法。
为实现以上目的，本发明通过以下技术方案予以实现：
一种基于声学信号的气泡释放通量量化方法，其特征在于：包括以下步骤：
S1：混响信号去除；
S2：气泡声学特性获取；
S3：气泡固有属性计算；
S4：气泡释放通量计算；
其中，步骤S4的计算公式为：
Figure BDA0002549698780000021
公式(5)中，ρ(mL/m3)为水中的气泡浓度，v(m/s)为气泡的上升速度，P为采样点实测大气压，P0为标准大气压(1013.25hPa)。
进一步地，所述步骤S1包括以下步骤：
B1：观察回波图，确定气泡分布范围，通过设定阈值，将该范围外的噪声排除掉，实现首次去噪；
B2：手动定义无效数据区域，进行二次去噪，去除行船噪声和鱼类信号。
进一步地，所述步骤S2中气泡声学特性包括气泡平均目标强度TS(dB)、体积反向散射强度Sv(dB)以及体积反向散射系数sv(m2/m3)；所述步骤S2包括以下步骤：
C1:通过单体目标信号检测，求得气泡平均目标强度TS(dB)；
C2:通过回波积分获取体积反向散射强度Sv(dB)；
C3:根据体积散射强度Sv(dB)计算体积散射系数sv(m2/m3)；
其中，步骤C3的计算公式为
Figure BDA0002549698780000022
进一步地，所述步骤S3中的气泡固有属性包括单体气泡体积V(mL)、气泡上升速度v(m/s)以及气泡浓度ρ(mL/m3)，所述步骤S3包括以下步骤：
D1：水体中气泡体积范围在0.005-20mL之间，对于该区间的气泡，根据气泡平均目标强度TS(dB)计算单体气泡体积V(mL)；
D2：体积为0.005mL-0.035mL之间的气泡随着半径的增大而上升速度减小，0.035-0.22mL体积的气泡上升速度基本上是恒定的，大于0.5-1mL的气泡上升速度随着半径增大而增大，根据气泡体积计算气泡上升速度v(m/s)；
D3：根据体积散射系数sv(m2/m3)计算气泡浓度ρ(mL/m3)；
其中，D1中的计算公式为：V＝995600e0.3092×TS (1)；
D2中的计算公式为：v＝-(0.00219±0.00093)(log V)4-(0.000737±0.00091)(log V)3+(0.04413±0.005)(log V)2+(0.0662+0.0034)log V+(0.2663+0.005)(2)；
D3中的计算公式为：ρ≈13500×sv (3)。
本发明的有益效果：
本发明提供的基于声学信号的气泡释放通量量化方法采用声学技术对气泡进行信号采集，能够在大水域实现连续监测和全面捕捉气泡排放信息，其效率高、覆盖面积大，可以对气泡的空间异质性进行无干扰评估，该方法可为定量研究水库温室气体排放，为减少水力工程和水电站开发等过程中温室气体的排放提供参考，也可为调查大水面温室气体释放提供技术手段。
附图说明
图1为本发明基于声学信号的气泡释放通量计算方法的流程框图；
图2为本发明混响信号去除前的声学数据图；
图3为本发明混响信号去除后的声学数据图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。
以2019年8月1日长洲水利枢纽库区监测数据为例。本次利用声学仪器进行气泡信号采集的时间长度约9分钟，数据量为95M。
请参考图1，一种基于声学信号的气泡释放通量量化方法，包括以下步骤：
S1：混响信号去除；
S2：气泡声学特性获取；
S3：气泡固有属性计算；
S4：气泡释放通量计算。
其中，步骤S1包括以下步骤：
A1：设置去噪阈值进行首次去噪，去噪阈值为为(-72，-36)dB；
A2：手动定义无效数据区域，进行二次去噪，去除行船噪声和鱼类信号。
请参考图2和图3，首先将.raw的声学原始数据导入Echoview软件中，通过观察回波图，发现气泡目标强度分布范围在-72～-36dB之间，因此设定输出阈值为(-72，-36)dB，将大于-36dB和小于-72dB的信号排除在外；另外通过手动定义无效数据区域，将行船噪声、鱼类信号排除掉。
其中，步骤S2中气泡声学特性包括气泡平均目标强度TS(dB)、体积反向散射强度Sv(dB)以及体积反向散射系数sv(m2/m3)，步骤S2包括以下步骤：
C1:通过单体目标信号检测，求得气泡平均目标强度TS(dB)，气泡的平均目标强度为-54.59dB。
C2:通过回波积分获取体积反向散射强度Sv(dB)，气泡体积反向散射强度的大小为-57.05dB。
C3:根据体积散射强度Sv(dB)计算体积散射系数sv(m2/m3)，步骤C3的计算公式为：
Figure BDA0002549698780000051
计算出的体积散射系数大小为1.97×10-6。
其中，步骤S3中的气泡固有属性包括单体气泡体积V(mL)、气泡上升速度v(m/s)以及气泡浓度ρ(mL/m3)。
步骤S3包括以下步骤：
D1：根据气泡平均目标强度TS(dB)计算单体气泡体积V(mL)，文献资料表明，水体中甲烷气泡单体体积大多数位于0.005-20mL之间。甲烷气泡单体体积的计算公式为：V＝995600e0.3092×TS (1)，
通过公式(1)计算得出气泡单体体积为0.057ml。
D2：根据气泡单体体积计算气泡上升速度v(m/s)，气泡在上升过程中，速度并非恒定的，而是根据体积不同上升速度也有差异。体积为0.005mL-0.035mL之间的气泡随着半径的增大而上升速度减小，体积为0.035-0.22mL体积的气泡上升速度基本上是恒定的，体积大于0.5-1mL的气泡上升速度随着半径增大而增大，气泡上升速度v(m/s)的计算公式为：
v＝-(0.00219±0.00093)(log V)4-(0.000737±0.00091)(log V)3
+(0.04413±0.005)(log V)2+(0.0662+0.0034)log V
+(0.2663+0.005)
(2)
通过公式(2)计算可得，气泡平均上升速度为：0.240-0.256m/s。
D3：根据体积散射系数sv(dB)计算气泡浓度，计算公式为：
ρ≈13500×sv (3)
通过公式(3)计算可得，气泡浓度为：0.027mL/m3。
其中，步骤S4中气泡释放通量计算公式为：
Figure BDA0002549698780000061
公式(5)中P为采样点实测大气压，P0为标准大气压(1013.25hPa)，根据气象网可以查得，调查期间采样点大气压为1012hPa，根据公式(5)计算得出，长洲水利枢纽库区示例数据每天释放甲烷气泡，即气泡释放通量大小为24.96-26.63mmol/m2。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (4)

1.一种基于声学信号的气泡释放通量量化方法，其特征在于：包括以下步骤：
S1：混响信号去除；
S2：气泡声学特性获取；
S3：气泡固有属性计算；
S4：气泡释放通量计算；
其中，步骤S4的计算公式为：
Figure FDA0002549698770000011
公式(5)中，ρ(mL/m3)为水中的气泡浓度，v(m/s)为气泡的上升速度，P为采样点实测大气压，P0为标准大气压(1013.25hPa)。
2.根据权利要求1所述的基于声学信号的气泡释放通量量化方法，其特征在于：所述步骤S1包括以下步骤：
B1：观察回波图，确定气泡分布范围，通过设定阈值，将该范围外的噪声排除掉，实现首次去噪；
B2：手动定义无效数据区域，进行二次去噪，去除行船噪声和鱼类信号。
3.根据权利要求1所述的基于声学信号的气泡释放通量量化方法，其特征在于：所述步骤S2中气泡声学特性包括气泡平均目标强度TS(dB)、体积反向散射强度Sv(dB)以及体积反向散射系数sv(m2/m3)；所述步骤S2包括以下步骤：
C1:通过单体目标信号检测，求得气泡平均目标强度TS(dB)；
C2:通过回波积分获取体积反向散射强度Sv(dB)；
C3:根据体积散射强度Sv(dB)计算体积散射系数sv(m2/m3)；
其中，步骤C3的计算公式为
Figure FDA0002549698770000012
4.根据权利要求3所述的基于声学信号的气泡释放通量量化方法，其特征在于：所述步骤S3中的气泡固有属性包括单体气泡体积V(mL)、气泡上升速度v(m/s)以及气泡浓度ρ(mL/m3)，所述步骤S3包括以下步骤：
D1：水体中气泡体积范围在0.005-20mL之间，对于该区间的气泡，根据气泡平均目标强度TS(dB)计算单体气泡体积V(mL)；
D2：体积为0.005mL-0.035mL之间的气泡随着半径的增大而上升速度减小，0.035-0.22mL体积的气泡上升速度基本上是恒定的，大于0.5-1mL的气泡上升速度随着半径增大而增大，根据气泡体积计算气泡上升速度v(m/s)；
D3：根据体积散射系数sv(m2/m3)计算气泡浓度ρ(mL/m3)；
其中，D1中的计算公式为：V＝995600e0.3092×TS (1)；
D2中的计算公式为：
v＝
-(0.00219±0.00093)(logV)4-(0.000737±0.00091)(logV)3+(0.04413±0.005)(logV)2+(0.0662+0.0034)logV+(0.2663+0.005) (2)；
D3中的计算公式为：ρ≈13500×sv (3)。
CN202010571291.5A 2020-06-22 2020-06-22 一种基于声学信号的气泡释放通量量化方法 Active CN111610520B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010571291.5A CN111610520B (zh) 2020-06-22 2020-06-22 一种基于声学信号的气泡释放通量量化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010571291.5A CN111610520B (zh) 2020-06-22 2020-06-22 一种基于声学信号的气泡释放通量量化方法

Publications (2)

Publication Number Publication Date
CN111610520A true CN111610520A (zh) 2020-09-01
CN111610520B CN111610520B (zh) 2021-08-17

Family

ID=72200831

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010571291.5A Active CN111610520B (zh) 2020-06-22 2020-06-22 一种基于声学信号的气泡释放通量量化方法

Country Status (1)

Country Link
CN (1) CN111610520B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114042391A (zh) * 2021-10-13 2022-02-15 三峡大学 一种水下气泡制造装置及方法
CN114275114A (zh) * 2022-01-04 2022-04-05 江南造船(集团)有限责任公司 一种船底气泡的检测方法及气泡层监测仪

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007000580A1 (en) * 2005-06-25 2007-01-04 University Of Southampton Contrast enhancement between linear and nonlinear scatterers
WO2014106039A2 (en) * 2012-12-28 2014-07-03 Exxonmobil Research And Engineering Company Methods for processing and interpreting signals from static and acoustic probes in fluidized bed reactor systems
CN108590636A (zh) * 2018-06-22 2018-09-28 浙江大学 海底天然气水合物气泡泄露监测装置
CN108931777A (zh) * 2018-04-09 2018-12-04 青岛海洋地质研究所 海底冷泉输出气体运移通量的测量方法
CN109598739A (zh) * 2018-12-14 2019-04-09 广州海洋地质调查局 一种水体中渗漏气泡数目的统计方法及处理终端
CN110531363A (zh) * 2019-09-04 2019-12-03 中国水产科学研究院珠江水产研究所 一种淡水水体气泡的声学信号识别方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2007000580A1 (en) * 2005-06-25 2007-01-04 University Of Southampton Contrast enhancement between linear and nonlinear scatterers
WO2014106039A2 (en) * 2012-12-28 2014-07-03 Exxonmobil Research And Engineering Company Methods for processing and interpreting signals from static and acoustic probes in fluidized bed reactor systems
CN108931777A (zh) * 2018-04-09 2018-12-04 青岛海洋地质研究所 海底冷泉输出气体运移通量的测量方法
CN108590636A (zh) * 2018-06-22 2018-09-28 浙江大学 海底天然气水合物气泡泄露监测装置
CN109598739A (zh) * 2018-12-14 2019-04-09 广州海洋地质调查局 一种水体中渗漏气泡数目的统计方法及处理终端
CN110531363A (zh) * 2019-09-04 2019-12-03 中国水产科学研究院珠江水产研究所 一种淡水水体气泡的声学信号识别方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
J. SCHNEIDER VON DEIMLING ET AL.: "Acoustic imaging of natural gas seepage in the North Sea:Sensing bubbles controlled by variable currents", 《LIMNOLOGY AND OCEANOGRAPHY-METHODS》 *
梅赛等: "甲烷羽状流水体声学探测及气体运移通量测算", 《海洋地质前沿》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114042391A (zh) * 2021-10-13 2022-02-15 三峡大学 一种水下气泡制造装置及方法
CN114275114A (zh) * 2022-01-04 2022-04-05 江南造船(集团)有限责任公司 一种船底气泡的检测方法及气泡层监测仪
CN114275114B (zh) * 2022-01-04 2023-03-14 江南造船(集团)有限责任公司 一种船底气泡的检测方法及气泡层监测仪

Also Published As

Publication number Publication date
CN111610520B (zh) 2021-08-17

Similar Documents

Publication Publication Date Title
CN111610520B (zh) 一种基于声学信号的气泡释放通量量化方法
Crumpton et al. Effects of differential growth and mortality in the seasonal succession of phytoplankton populations in Lawrence Lake, Michigan
Gao et al. New insight into global blue carbon estimation under human activity in land-sea interaction area: A case study of China
Zhang et al. Nitrous oxide in the Changjiang (Yangtze River) Estuary and its adjacent marine area: Riverine input, sediment release and atmospheric fluxes
Wyatt et al. Particulate nutrient fluxes over a fringing coral reef: Source‐sink dynamics inferred from carbon to nitrogen ratios and stable isotopes
Sun et al. Temporal and spatial variability of methane emissions in a northern temperate marsh
WO2015149408A1 (zh) 一种河湖沉积物中重金属污染的生态风险确定方法
Long et al. Ebullition of oxygen from seagrasses under supersaturated conditions
Tang et al. Rapid urbanization effects on partial pressure and emission of CO2 in three rivers with different urban intensities
CN113673737B (zh) 基于卫星遥感影像的藻型湖泊水体溶解二氧化碳估算方法
CN109959619B (zh) 一种测定贝类养殖蓝碳在沉积物中的累积速率的方法
Hoffmann et al. A simple calculation algorithm to separate high-resolution CH 4 flux measurements into ebullition-and diffusion-derived components
CN115062071B (zh) 一种用于河道的水资源污染扩散分析方法及系统
CN111750283A (zh) 基于深度学习的强背景噪声环境下的气体管道泄漏识别方法
CN102453743A (zh) 一种制备水源水水质急性生物毒性检测试剂的方法及其检测方法
CN117077542B (zh) 基于仿生红树林模型拖曳力系数的红树林消浪评估方法
CN110263876B (zh) 一种雨声信号降噪处理方法和系统
Shao et al. Version 2 of the global oceanic diazotroph database
CN111204842A (zh) 一种通过神经网络实现超滤膜污染评价的方法、装置、系统
CN207081564U (zh) 着生生物采样器
Müller et al. Fate of peat-derived carbon and associated CO 2 and CO emissions from two Southeast Asian estuaries.
CN101033617A (zh) 一种基础结构测试信号的动态检测方法
CN104237318B (zh) 一种利用电子鼻检测根茎类植物根部病虫害的方法
CN108760890A (zh) 一种使用信号能谱系数监测脆性材料损伤演化过程的方法
CN105002273B (zh) 用于实时荧光定量pcr方法监测链状亚历山大藻细胞爆发性生长的引物组及监测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant