超聲多普勒測量海底滲漏氣泡流速的實驗研究

劉潤華，鄒大鵬，龍建軍，胡少興

（廣東工業大學　機電工程學院，廣東　廣州 510006）

超聲多普勒測量海底滲漏氣泡流速的實驗研究

劉潤華，鄒大鵬*，龍建軍，胡少興

（廣東工業大學機電工程學院，廣東廣州 510006）

海底滲漏氣泡流速是測量海底冷泉滲漏氣體流量的主要參數。為了測量海底滲漏氣泡流速，利用超聲多普勒測量的原理和氣泡的二次多普勒效應，提出海底滲漏氣泡速度的超聲多普勒測量方法。超聲多普勒測量根據超聲多普勒原理設計了氣泡流量測量裝置，通過多普勒頻移算出海底滲漏氣泡流速，同時用視頻法測出海底滲漏氣泡流速。從多普勒法和視頻法的氣泡速度數據t分布95%的置信水平誤差可知，視頻法的氣泡速度數據的波動大于多普勒法的氣泡速度數據的波動，表明了超聲多普勒測量海底滲漏氣泡流速的有效性和準確性，優于視頻法測量海底滲漏氣泡流速。

海底冷泉滲漏；多普勒；頻移；氣泡流速

海底冷泉滲漏是海洋環境中廣泛分布的自然現象［1］，部分滲漏到海洋水體以及大氣中，天然氣的成分主要為甲烷，甲烷是重要的溫室氣體之一，每年通過海底冷泉天然氣滲漏釋放到海洋水體及大氣中的甲烷的量是巨大的［2-3］；部分滲漏天然氣在適合的環境中將生成天然氣水合物，儲存在水合物中的天然氣有可能成為未來能源。人類有必要對海底冷泉進行觀測，探測出海底冷泉天然氣的滲漏量。

目前測量海底冷泉天然氣的滲漏量應用的方法有：容積收集法［4-6］，利用收集器直接收集海底冷泉天然氣滲漏的氣泡氣體的量，難以用于較深的水域和連續觀測；回波強度檢測氣泡流量法［7-8］，利用氣泡柱反射回波的強度間接檢測海洋中的氣泡的流量，回波信號強度主要取決于被測氣泡柱狀流界面的反射，回波沒有探測到氣泡柱狀流內部信息；攝像裝置觀測滲漏氣泡法［9］，采用高分辨CCD獲取氣泡流動圖像，采用圖像處理軟件計算氣泡的直徑，此方法分辨率高，需要光源耗電大，所以適宜短時觀測；遙測紅外成像法［10］，根據甲烷強烈吸收短波紅外線的特性，采用光譜數據評估釋放到海面的甲烷數量，該方法適用于大范圍的海底冷泉滲漏量定性估計，難以定量測量滲漏量。

本文基于測量前調整氣泡狀態并用透射聲波測量滲漏氣泡流量的方法［11-12］，在開展了超聲測截面含氣率和光學測氣泡流速的實驗研究基礎上，提出應用超聲多普勒測速法測量海底滲漏氣泡流速。

1　超聲多普勒測量原理和方法

1.1超聲多普勒測量的原理

多普勒超聲流量測量的原理是基于物理學中聲波的多普勒效應，當超聲波聲源和反射界面或散射體之間存在相對運動時，被反射或被散射的超聲波信號的頻率將產生變化，該頻率與聲源超聲波頻率之間的差值與相對運動的速度成正比。這一現象是物理學家Christian Johann Doppler于1842年首次發現的，被稱為多普勒效應。接收信號的頻率與聲源的頻率之差稱為多普勒頻移，相應的頻差信號稱為多普勒信號。當氣液兩相中存在著可供反射或散射超聲波的運動粒子（例如氣泡、固體粒子等）時，便可利用多普勒信號確定運動粒子的速度。

1.2超聲測量的二次的多普勒效應

如圖1所示，發射超聲波流量傳感器透穿過氣泡流動的通道是一個矩形通道，是由網狀材料組成，目的是超聲波透射而不改變原來的路徑和運動氣泡能固定在矩形測量通道內，使作用在運動氣泡的超聲波能返回接收超聲波流量傳感器。

圖1　反射式二次多普勒圖解

二次的多普勒效應［13］（如圖1）是指：

（1）若聲源靜止、接收點運動，則接收點收到一個波長的時間（周期）變化，因而產生多普勒頻移。運動著的氣泡首先作為接收體即產生了第一次多普勒效應；

（2）若聲源運動、接收點靜止，則由于聲源運動使波在傳播過程中發生波長變化，因而產生多普勒頻移。由于聲波散射（反射），氣泡再作為運動的聲源產生了第二次多普勒效應。

1.3超聲測量的方法

接收信號的頻率f與聲源的頻率f0之差稱為多普勒頻移fd，即fd=f-f0?？梢詫懗墒剑?）［14］：

式中：c為超聲波高頻在水氣介質中傳播速度，可以看成是水介質的聲速，隨著頻率增高，水氣介質的聲速接近于無氣泡時的狀態［15］；v為氣泡流速；θ1為發射超聲波束與氣泡流速方向矢量夾角；θ2為接收超聲波束與氣泡流速方向矢量夾角（其中發收超聲波束是無方向性）。

當θ1與θ2大致相等，由式（1）可得：

由（2）式可知，氣泡速度v與多普勒頻移fd成正比。測出多普勒頻移fd就可以算出氣泡速度v。頻移可以通過短時FFT計算，頻譜分析可以得出超聲換能器接收到頻率成分。

2　超聲多普勒測量海底滲漏氣泡流速裝置的設計

根據超聲波發射和散射定律［16］：當障礙物直徑大于超聲波波長的1/2時，在障礙物表面主要發生反射；當障礙物直徑小于超聲波波長的1/2時，在障礙物表面主要發生散射。反射的回波比散射的回波功率強。因此，超聲多普勒需要設計成反射式的方式。氣泡破碎均勻后的氣泡直徑在3～5 mm，選用中心頻率為973 kHz的超聲換能器，超聲波在水介質傳播的波長為1.5 mm，符合超聲波發射和散射定律的反射條件。

海底滲漏氣泡流速的超聲多普勒測量裝置（圖2）由4個部分組成：模擬海底冷泉滲漏氣泡形成部分、氣泡收集預處理部分、聲波測量部分和視頻驗證部分。

模擬海底冷泉滲漏氣泡形成部分，空壓機（氣源）產生壓縮空氣進入水池中壓縮氣體的流量，壓縮氣體在水中形成氣泡，模擬海底冷泉滲漏氣泡產生過程。

氣泡收集預處理部分，氣泡收集預處理裝置收集產生的氣泡，預先破碎調整氣泡的流動狀態，使其大小基本一致，在水中的分布趨于均勻。

圖2　超聲多普勒測量裝置示意圖

聲波測量部分，聲波測量部分由奧迪威超聲波流量傳感器（AW5Y0980K03L142Z）、TH204B聲波儀、電子電路、PC機組成。通過TH204B聲波儀單次采集激發電子電路產生有限個正弦波作用于奧迪威超聲波流量傳感器，為TH204B聲波儀同步采集接收換能器的數據。

視頻驗證部分主要由工業攝像機、視頻采集卡與PC機組成，用來觀察驗證氣泡預處理的效果以及計算氣泡流速來檢測多普勒頻移法測出的海底滲漏氣泡流速的準確性。

3　超聲多普勒測量海底滲漏氣泡流速的實驗結果

實驗采用有限個（20 000個）正弦波發射頻率取值為奧迪威超聲波流量傳感器的諧振頻率973 kHz作為發射，同步采樣的間隔為0.5 μs（采樣頻率為2 MHz），采樣長度32k字節，θ1取值為60°，水中聲速常溫時1 500 m/s。實驗在無氣泡和有氣泡不同流量的情況下進行對比。

3.1無氣泡時的實驗

因有直達波（聲電泄漏）的存在，實際上接收換能器也會有信號進來，時域波形如圖3，而此時的信號就是發射基頻。對采集到的數據經過短時FFT多次測量讀數可以得到頻率為f0=972.716 4 kHz（圖4）。

圖3　直達波的幅值曲線

圖4　無氣泡時接收數據的頻譜圖

3.2有氣泡時的實驗

在有氣泡時對采集到接收超聲流量換能器的數據（圖5）進行短時FFT的運算可得有氣泡時的頻譜圖（圖6）。有兩個峰對應的頻率是基頻和經過氣泡反射回來的頻率。

圖5　有氣泡時的幅值曲線

圖6　有氣泡時接收數據的頻譜圖

通過有無氣泡對比實驗（如圖4和圖6）對比驗證了超聲多普勒頻移的存在，多普勒頻移值fd就是圖6兩峰對應頻率差值。

3.3不同氣流量下的頻移值及海底滲漏氣泡速度

為了算出具體的多普勒頻移值，本實驗在不同氣體流量情況下，如在20 sccm，40 sccm，60 sccm，80 sccm，100 sccm，120 sccm，140 sccm，160 sccm，180 sccm，200 sccm（其中sccm為mL/min）時，進行5次測量求出多普勒頻移的平均值，根據（2）式求得海底滲漏氣泡的平均速度。在不同流量下具體對應的多普勒頻移值和氣泡速度值如表1所示。

表1　不同流量下海底滲漏氣泡的速度

測得每隔20 sccm流量下具體對應的多普勒頻移值、多普勒法海底滲漏氣泡速度值數，繪制出氣流量-多普勒頻移圖及擬合曲線，如圖7所示。

圖7　氣體流量—多普勒頻移圖

對比多普勒法氣泡速度數據和視頻法的海底滲漏氣泡速度數據，兩種方法的測量結果接近，趨勢相同，不過視頻法在大流量下因為紊流的因素，引起速度計算減小。兩種方法的海底滲漏氣泡速度數據t分布95%的置信水平誤差如圖8。

圖8　氣體流量—氣泡速度t分布95%的置信水平誤差棒圖

4　結語

本實驗超聲多普勒測量海底滲漏氣泡多普勒頻移在0.30～0.55 kHz之間，氣泡平均速度0.28～0.47 m/s左右，與胡柳等［12］文中用視頻法測量氣泡流速結果一致。本實驗用視頻法的結果和多普勒法進行了對比，多普勒法的數據和視頻法的數據t分布置信水平為95%的預測區間交叉重合度較大，從多普勒法和視頻法的氣泡速度數據t分布95%的置信水平誤差可知，視頻法的氣泡速度數據的波動大于多普勒法的氣泡速度數據的波動，表明了超聲多普勒測量海底滲漏氣泡流速方法的精度更高。

基于多普勒頻移法測量海底滲漏氣泡速度和視頻法測量海底滲漏氣泡速度存在一些差異。原因有：（1）氣泡上升的不只是單純的垂直上升，而是在上升過程中存在著微小的左右擺動［17］，以致超聲波聲束與氣泡流速方向夾角的難以準確測量，只能以一個估算值代替此夾角；（2）多普勒法測的是海底滲漏氣泡的瞬時速度（即是垂直上升方向和水平方向的合成速度）；視頻法測出的數據是投影速度（即垂直上升的速度），忽略了水平方向的運動；（3）短時FFT頻率分辨率有待提高，頻率分辨率（Δf=fs/N，其中fs為采樣頻率，N為數據點數）為0.061 kHz，對應海底滲漏氣泡速度分辨率為0.054 m/s。

以上研究表明，多普勒頻移法測量海底滲漏氣泡速度優于視頻法測量海底滲漏氣泡流速，可以為測量海底冷泉滲漏流量提供一種新的方法。

［1］Judd AG.The Global Importance and Context ofMethane Escape fromthe Seabed［J］.Geo-Marine Letters，2003，23:147-154.

［2］Amouroux D，Roberts G，Rapsomanikis S，et al.Biogenic Gas Emission to the Atmosphere from Near-Shore of the North-Western Black Sea［J］.Estuarine，Coastal&ShelfScience，2002，54（3）:575-587.

［3］劉勇健，李彰明，張麗娟，等.未來新能源可燃冰的成因與環境巖土問題分析［J］.廣東工業大學學報，2010，27（3）:83-87.

［4］Cable J E，Burnett W C，Chanton J P，et al.Field Evaluation of Seepage Meters in the Coastal Marine Environment［J］.Estuarine，Coastal and ShelfScience，1997，45（3）:367-375.

［5］Sassen Roger，Milkov Alexei V，Ozgul Ercin，et al.Gas Venting and Subsurface Charge in the Green Canyon Area，Gulf of Mexico Continental Slope:Evidence ofa Deep Bacterial Methane Source［J］.Organic Geochemistry，2003，34（10）:1455-1464.

［6］Washburn Libe，Johnson Cyril，Gotschalk Chris C，et al.A Gas-Capture Buoy for Measuring Bubbling Gas Flux in Oceans and Lakes［J］.Journal ofAtmospheric and Oceanic Technology，2001，18（8）:1411-1420.

［7］Greinert J，Nützel B.Hydroacoustic Experiments to Establish a Method for the Determination of Bubble Fluxes at Cold Seeps［J］.Geo-Marine Letters，2004，24（2）:75-85.

［8］J Greinert，BNutzel，余平，陶軍（譯）.冷泉中甲烷氣泡流量檢測方法的水聲試驗［J］.海洋地質，2004（3）:53-66.

［9］LeiferI，DeLeeuwG，CohenLH.OpticalMeasurementofBubbles:SystemDesignandApplication［J］.JournalofAtmosphericandOceanic Technology，2003，20（9）:1317-1332.

［10］Roberts D A，Bradley E S，Cheung R，et al.Mapping Methane Emissions from a Marine Geological Seep Source Using Imaging Spectrometry［J］.Remote SensingofEnvironment，2010，114（3）:592-606.

［11］龍建軍，黃為，鄒大鵬，等，海底天然氣滲漏流量聲學測量方法及初步實驗研究［J］.熱帶海洋學報，2012（5）:100-105.

［12］胡柳.冷泉滲漏聲波測量裝置主體研制與氣泡—水聲學特性的實驗研究［D］.廣州：廣東工業大學，2014.

［13］王威琪，邵謙明.血流Doppler效應的起源問題［J］.復旦學報：自然科學版，1992，03:241-246.

［14］王威琪.超聲多普勒血流測量技術［J］.醫療器械，1983，02:24-33.

［15］張憶，孫輝.含氣泡水的聲速及聲吸收系數［C］//中國聲學學會.中國聲學學會2009年青年學術會議［CYCA’09］論文集，中國聲學學會，2009.

［16］徐修萍，楊明.基于超聲多普勒人工心臟血栓檢測系統的研究［J］.電子測量與儀器學報，2010，04:396-401.

［17］宋禹林，譚思超，付學寬.晃蕩對氣泡上升運動影響的數值研究［J］.核動力工程，2014，S1:71-74.

Experimental Study on Ultrasonic Doppler Measurement of the Velocity of Seafloor Leakage Bubbles

LIU Run-hua，ZOU Da-peng，LONG Jian-jun，HU Shao-xing
School of Electromechanical Engineering，Guangdong University of Technology，Guangzhou 510006，Guangdong Province，China

The velocity of seafloor leakage bubbles is an important parameter to measure the gas flow from seafloor cold spring leakage.In order to measure the velocity of seafloor leakage bubbles，an ultrasonic Doppler measurement method is proposed in this paper using the principle of ultrasonic Doppler measurement and secondary Doppler effect.In accordance with the principles of ultrasonic Doppler measurement，a device is designed for measuring bubble flow，so as to calculate the velocity of seafloor leakage bubbles through Doppler frequency shift.Meanwhile，the velocity of leakage bubbles is also measured by the video method.From the confidence level error of 95%t distribution of bubble velocity data form the Doppler and video methods，it can be concluded that the fluctuation of bubble velocity data from the video method is larger than that from the Doppler method，proving the effectiveness and accuracy of ultrasonic Doppler measurement of the velocity of seafloor leakage bubbles，which is better than the method of video measuring.

seafloor cold spring leakage；Doppler；frequency shift；bubble velocity

O429

A

1003-2029（2016）02-0046-05

10.3969/j.issn.1003-2029.2016.02.008

2015-03-10

三亞深?？茖W與工程研究所知識創新工程領域前沿項目資助（SIDSSE-201208）；廣東省科技計劃資助項目（2014A040402008）

劉潤華（1987-），男，碩士研究生，主要從事海洋聲學檢測技術研究。E-mail:294863052@qq.com

鄒大鵬（1977-），男，副教授，主要從事工程裝備和聲學檢測技術研究。E-mail:anthonyzou@126.com