廖国祥 李艳妮

摘要：為深入研究井喷和管道破损等海底事故后溢油在水下环境中的输移扩散过程，文章利用自制的组合式环形水槽（周长9 . 7m、宽0 . 4 5m、深1m） ，以阿曼原油与消油剂混合物和淡水（含示踪剂）为模拟污染物，初步开展静水和动水环境中海底溢油浮射流的物理模拟实验，并应用基于拉格朗日积分方法的水下溢油浮射流模型进行数值模拟比较分析。研究结果表明：静水环境中，水下溢油浮射流主要沿喷口的垂直中心线向水面输移扩散；非均匀流动水环境中，横流速度越大，浮射流输移轨迹的弯曲程度越明显；数值模拟的浮射流轨迹总体上与实验观测结果符合较好；研究结果可为今后相关物理实验和数值模型的改进研究提供有益参考。

关键词：水下溢油；浮射流；输移扩散；非均匀流；数值模拟

中图分类号： P 7 6； X 5 5 文献标志码： A 文章编号： 1 0 0 5-9 8 5 7（ 2 0 2 0） 1 2-0 0 7 1-0 7

P r e l i m i n a r yE x p e r i m e n t a l S t u d ya n dN u m e r i c a lM o d e l i n g o f t h eU n d e r w a t e rO i l S p i l B u o y a n t J e t s i n S t a t i ca n dF l o w i n gW a t e rE n v i r o n m e n t s

L I AOG u o x i a n g1， 2， L IY a n n i3

（ 1 .N a t i o n a lM a r i n eE n v i r o n m e n t a lM o n i t o r i n gC e n t e r， D a l i a n1 1 6 0 2 3， C h i n a；

2 . S t a t eE n v i r o n m e n t a lP r o t e c t i o nK e yL a b o r a t o r yo fC o a s t a lE c o s y s t e m， D a l i a n1 1 6 0 2 3， C h i n a； 3 . E n v i r o n m e n t a lS c i e n c ea n dE n g i n e e r i n gC o l l e g e， D a l i a nM a r i t i m eU n i v e r s i t y， D a l i a n1 1 6 0 2 6， C h i n a）

A b s t r a c t： I no r d e r t os t u d y t h e t r a n s p o r t a n dd i f f u s i o np r o c e s s e so f o i l s p i l l i nu n d e r w a t e r e n v i r o n -m e n t a f t e rs u b s e aa c c i d e n t s（ e . g . s u b s e ao i l - w e l lb l o w o u t， o i l /g a sp i p e l i n el e a k a g e） ， t h i sp a p e r u s e da na n n u l a r t a n kw i t hd i m e n s i o n so f 9 . 7m×0 . 4 5m×1m（ p e r i m e t e r×w i d t h×d e p t h）a s t h e b a s i ct e s t f a c i l i t y， a n ds e l e c t e dOm a nc r u d eo i l， m i x t u r eo f c r u d eo i l a n dc h e m i c a l d i s p e r s a n t， f r e s h w a t e r（ w i t h i n ka s t r a c e r）a s s i m u l a t e dp o l l u t a n t s， a n dc a r r i e do u t as e r i e so fu n d e r w a t e rb u o y a n t j e tp h y s i c a l e x p e r i m e n t s， a n dc o n d u c t e d t h en u m e r i c a lm o d e l i n ga n dc o m p a r a t i v e a n a l y s i sb yu s i n g aL a g r a n g i a ni n t e g r a lm e t h o db a s e du n d e r w a t e ro i l /g a sb u o y a n t j e tm o d e l . E x p e r i m e n t a lr e s u l t s s h o w e dt h a t t h eb u o y a n t j e t s t r a n s p o r t e da l o n gt h ev e r t i c a l l i n eo f t h en o z z l ea n ds p r e a di nh o r i -z o n t a l d i r e c t i o n s i nt h es t a t i cw a t e re n v i r o n m e n t；w h i l e i nt h en o n - u n i f o r mc r o s sf l o w i n gw a t e r e n v i r o n m e n t s， t h eb u o y a n t j e t t r a j e c t o r i e sw e r eb e n t o v e rb y t h e s t r o n g e r c r o s s f l o w.I t a l s o s h o w e d t h a t t h en u m e r i c a lm o d e l e d t r a j e c t o r i e so fb u o y a n t j e t sg e n e r a l l ya g r e e dw e l lw i t ht h ee x p e r i m e n -t a l o b s e r v e dd a t a . T h er e s u l t so f t h i sp a p e rc a np r o v i d eau s e f u l r e f e r e n c e f o r t h e i m p r o v e m e n to f p h y s i c a l e x p e r i m e n t sa n dn u m e r i c a lm o d e l i n g i nt h e f u t u r e .

K e yw o r d s： U n d e r w a t e r o i l s p i l l， B u o y a n t j e t， T r a n s p o r t a n dd i f f u s i o n， N o n - u n i f o r mc r o s s f l o w， N u -m e r i c a lm o d e l i n g

0 引言

自2 0 1 0年墨西哥湾“深水地平线”深海溢油事故以来，水下溢油问题受到世界各海洋石油开发国家的高度重视和关注。近年来，我国不仅加强对渤海和南海的浅海油气资源开采，而且积极推进南海深水油气资源的勘探开发，发生水下溢油事故的风险日益增高。在此背景下，通过开展物理实验和数值模拟研究，深入认识水下溢油污染物在水下环境中的输移扩散过程，可为事故前的风险评估、事故时的应急决策和事故后的损害评估等提供重要的科学和技术支持。

国外学者自2 0世纪7 0年代起就开展了静水条件下的水下溢油浮射流实验，用于支持水下溢油模型的构建和验证[ 1-3]。挪威学者分别于1 9 9 6年和2 0 0 0年先后开展1 0 7m和8 4 4m水下大型溢油现场实验[ 4-5]。由于水下溢油现场实验的成本高昂，且不易观测污染物输移扩散的全过程，现有研究仍以实验室模拟实验研究为主。S o c o l o f s k y等[ 6]开展均匀流环境中的多相流实验，研究油气混合物的输移和分离特性；M a s u t a n i等[ 7]在高为1 . 3 0m、长和宽均为0 . 5 5m的长方体水槽中开展原油样品喷射实验，根据射流形态确立油品破碎成油滴的边界； B r a n d v i k等[ 8]利用直径为3m、高为6m的大型垂直水槽开展系列实验，研究泄漏条件、油品特性和消油剂对油滴粒径分布的影响； Z h a o等[ 9]在美国O h m s e t t大型波浪槽中，利用直径为2 5 . 4mm、泄漏量为6 . 3L/s的水下横管开展大型溢油实验，分析水下溢油的动力特性和油滴粒径分布。

近年来，我国学者积极开展水下溢油实验研究，并取得显著进展。王晶[ 1 0]和肖杰[ 1 1]开展水下单一孔口形成的气泡/油滴实验，研究海底管线微孔泄漏的油滴尺寸和上升速度；臧晓刚等[ 1 2]在2 . 0m×1 . 0m×1 . 8m的水箱中进行水下管道微孔溢油实验，研究在不同泄漏量时溢油油滴上升至水面的时间以及泄漏产生的浮羽射流的形态；钱国栋等[ 1 3]在1m×1m×2m的垂直水槽中开展“蓬莱1 9 - 3”水下溢油模拟实验，并评估水下溢油的喷口直径和喷射速度对油滴破碎形成的油滴粒径分布的影响。

与此同时，真实海洋环境通常是非均匀的分层水流环境，而现有实验室实验研究以静水和均匀流环境为主，难以深入研究更接近真实海洋環境条件的水下溢油污染物输移扩散过程。因此，本研究基于自制的组合式环形水槽，初步开展静水和动水（非均匀横流）环境中水下溢油浮射流的物理模拟实验和数值模拟分析，以期为水下溢油事故污染物输移扩散研究提供有益参考。

1 水下浮射流实验

1 . 1 材料与方法

1 . 1 . 1 物理实验装置

水下溢油浮射流物理模拟实验装置为1台自制的具有初步造波和造流功能的组合式环形溢油水槽（图1）。

水槽周长9 . 7m、宽0 . 4 5m、深1m，其中槽壁材质为玻璃，易于观测实验过程。水槽由2个直段和2个弯段组成，选取水槽的一侧直段作为水下浮射流物理模拟实验区域，并利用1台空气压缩机和1个油罐制造水下溢油浮射流。须指出的是，物理模拟实验重点关注水下溢油浮射流的输移扩散过程，因此以水槽底部喷口中心为原点建立三维空间坐标系，其中X轴和Y轴分别为水平面的E向和N向（具体表示为污染物输移距离，单位为m） ， Z轴为垂直方向（具体表示为实验水深，单位为m）。

1 . 1 . 2 水流环境模拟

真实海洋环境通常是非均匀流环境： 1 9 9 6年挪威“D e e p S p i l l”现场实验的结果表明水下3m、2 0m、5 0m和8 0m等不同深度呈不均匀状态，且随时间变化而动态变化，流速范围为0 . 0 3～0 . 1 0m/s[ 4]；我国学者对南海海流垂直结构的调查研究结果也展示了南海水下环境非均匀流的特点[ 1 4-1 6]。因此，本研究结合组合式环形水槽的实际条件，通过调节水泵的功率和位置，制造2种不同平均速度的流动环境。在实验前，待水槽内的水流循环超过1 5m i n且速度稳定后，利用超声波水流测定仪测定实验区域的水流速度，并记录1m i n的时均速度值、最小值和最大值（表1和图2） ，用于后续的数值模拟分析。

1 . 1 . 3 实验条件和方案设计

参考国内外实验研究经验，结合本研究的物理实验装置实际情况，对模拟污染物、喷口直径、初始喷射速度和水流环境等实验条件参数进行分析与设计。

（ 1）模拟污染物。国内外多采用密度较小和黏度较低的原油开展浮射流实验： B r a n d v i k等[ 8]采用密度为8 3 9 . 3k g/m3的O s e b e r g原油作为实验油品；臧晓刚等[ 1 2]采用密度为8 1 5k g/m3的0号柴油作为实验油品。2 0 1 0年墨西哥湾“深水地平线”溢油事故发生后，挪威和美国学者为研究水下消油剂处理对溢油输移的影响，将原油和消油剂预混合后开展污染物浮射流实验研究[ 8-9]。考虑到油品的获取难度，本研究选取大连港运输较多且密度较小和黏度较低的阿曼原油（ 2 0℃时的密度为8 7 2k g/m3）作为主要实验油品，并以阿曼原油与消油剂的混合物（剂油比为1∶2 0）作为水下浮射流实验的模拟污染物。为开展对比研究，选取淡水（添加蓝色或红色墨水作为示踪剂）作为模拟污染物。

（ 2）喷口直径。B r a n d v i k等[ 8]开展的溢油浮射流系列实验的喷口直径为0 . 5～3 . 0mm，并在原油与消油剂混合物浮射流实验中选取1 . 5mm的喷口直径；臧晓刚等[ 1 2]选取的实验喷口直径为2～ 3mm；钱国栋等[ 1 3]选取的实验喷口直径为1～ 3mm。经综合比较，本研究选取1 . 5mm和2 . 0mm作为实验喷口直径。

（ 3）初始喷射速度和水流环境。主要通过前述的组合式实验装置以及水泵和气泵等设备来调节和设定。

基于上述實验条件分析，制定水下溢油浮射流物理模拟实验方案（表2）。

1 . 1 . 4 实验步骤和测定方法

实验开始：向水槽内注入淡水（温度为室温） ，使水深达到0 . 6 5m（喷口位于水槽底部高5c m处，即喷口与水面的距离为0 . 6m） ；将配置好的模拟污染物（阿曼原油和淡水）添加至体积为2L的“油罐”中；开启压缩气泵，通过调节管路的压力达到实验所需的压强，待压强稳定后打开待测溢油管道的开关，同时利用高清摄像机（分辨率为19 2 0×10 8 0像素，录像频率为3 0帧/s）摄录浮射流输移扩散的过程；污染物的流量和喷射速度根据模拟污染物的体积和持续泄漏时间计算确定。

实验结束：将实验过程视频的高清截图导入I m a g e - P r oP l u s软件，经过空间校正，标记浮射流左边界和右边界的主要拐点，并计算确定每个拐点的水平和垂直坐标；综合所有拐点的坐标，得到浮射流实验污染物输移轨迹和扩散范围的数据。

1 . 2 结果与讨论

静水和动水环境中原油（阿曼原油以及阿曼原油与消油剂混合物）和淡水（含示踪剂）浮射流输移扩散的实验结果如图3和图4所示。需要说明的是，污染物的输移扩散范围是根据处于相对稳定状态且边界容易识别的浮射流实验图像处理分析后确定的。

由图3和图4可以看出：①静水环境中，污染物浮射流基本沿着与喷口垂直的中心线向上运动；随着周围水体不断被卷吸进入浮射流，浮射流的垂直上升速度快速下降，同时其水平扩散范围逐渐变大。②动水环境中，由于初始喷射速度大于横向流速，泄漏初期的浮射流轨迹基本垂直向上；随着射流的浮升，周围水体被卷吸进入浮射流，浮射流中心轨迹逐渐弯曲；对比2种动水环境（ Uc=3 . 4 8c m/s和Uc= 6 . 9 0c m/s）中的淡水浮射流轨迹，横向水流速度越大，浮射流中心轨迹开始弯曲的高度越低，弯曲程度越明显。

值得注意的是，在Uc=3 . 4 8c m/s的动水环境中，实验观测到的原油浮射流扩散范围比原油与消油剂混合物浮射流的大。这是因为原油自喷口喷射进入水体后破碎成为粒径更大的油滴，而粒径较大的油滴更易于脱离浮射流并浮升和扩散。此外，实验观测到的原油与消油剂混合物浮射流分散在水体中的油滴粒径明显更小，但受观测设备的限制，本研究未能对实验中油滴粒径的分布特征进行有效的定量分析。

2 数值模拟

2 . 1 水下溢油浮射流模型及其参数设置

本研究仅给出数值模型的质量守恒方程和动量守恒方程，详细介绍可参见文献[ 1 7]和[ 1 8]。

2 . 2 数值模拟结果与实验观测数据的比较

数值模拟结果与实验观测数据的比较如图6和图7所示。

由图6和图7可以看出：无论是原油还是淡水，静水环境中的污染物浮射流输移轨迹和扩散范围的数值模拟结果与实验观测数据的吻合程度均较高；与静水环境相比，在平均速度为3 . 4 8c m/s和6 . 9 0c m/s的非均匀流动水环境中，污染物浮射流在初始射流动量占主导作用阶段（自喷口至其以上的0 . 1m处）的输移轨迹与扩散范围基本吻合；随着非均匀流动水体被卷吸进入浮射流，尽管输移轨迹的数值模拟结果与实验观测数据总体吻合，但扩散范围存在一些差异。考虑到数值模型输入的环境流速是垂向5层实测流速的平均值而不是实时测定的分层流速等因素，本研究所采用的水下溢油浮射流模型总体上能够较好地模拟预测污染物浮射流的输移扩散趋势。

但须指出的是，由于实验过程中未能获取油滴粒径分布数据，更加精细化的溢油浮射扩散过程以及流动环境中大粒径油滴脱离浮射流后输移扩散过程尚未得到充分考虑。这也意味着现有水下溢油浮射流数值模型仍须结合更加精细的油滴粒径分布实验数据进一步验证和改进，以提高模拟预测结果的准确度。

3 结语

本研究尝试通过物理模拟实验和数值模拟分析相结合的方式，初步探析不同流动环境中水下溢油浮射流的输移扩散特征。开展的主要工作包括：应用自制的组合式环形水槽，以原油（阿曼原油与消油剂混合物）和淡水（含示踪剂）作为模拟污染物，开展静水和动水环境中水下浮射流的物理模拟实验；应用基于拉格朗日积分方法的数值模型，开展数值模拟并结合物理实验结果进行比较和分析。

由于实验装置和观测设备的限制，浮射流实验过程中非均匀流的实时流场以及油滴粒径的动态分布等参数未能有效获得，导致模型数值模拟与实验观测的污染物浮射流扩散范围存在一定的差异。但总体来说，本研究初步获得的物理模拟实验和数值模拟分析结果可为今后水下溢油物理实验和数值模型的改进研究提供有益参考。

参考文献

[ 1] F ANN E L？PT K， S J？E N K. H y d r o d y n a m i c so fu n d e r w a t e r b l o w o u t s[M]. T r o n d h e i m： S h i pR e s e a r c hI n s t i t u t eo fN o r w a y， M a r i n eT e c h n o l o g yC e n t r e， 1 9 8 0.

[ 2] Y A P APD， Z H E N GL. S i m u l a t i o no fo i l s p i l s f r o mu n d e rw a t e r a c c i d e n t sⅠ： m o d e l d e v e l o p m e n t[ J]. J o u r n a lo fH y d r a u l i cR e -s e a r c h， 1 9 9 7， 3 5（ 5） ： 6 7 3-6 8 7.

[ 3] Z H E N GL， Y A P APD.S i m u l a t i o no f o i l s p i l s f r o mu n d e r w a t e r a c -c i d e n t sⅡ： m o d e l v e r i f i c a t i o n[ J].J o u r n a l o fH y d r a u l i cR e s e a r c h， 1 9 9 8， 3 6（ 1） ： 6 7 3-6 8 8.

[ 4] R Y EH， B RAN D V I KPJ， S T R OMT. S u b s u r f a c eb l o w o u t s： r e -s u l t s f r o mf i e l de x p e r i m e n t s[ J]. S p i l lS c i e n c e & T e c h n o l o g y B u l l e t i n， 1 9 9 7， 4（ 4） ： 2 3 9-2 5 6.

[ 5] J OHAN S E N？. D e v e l o p m e n ta n dv e r i f i c a t i o no fd e e p - w a t e r b l o w o u t m o d e l s[ J]. M a r i n e P o l l u t i o n B u l l e t i n，2 0 0 3，4 7： 3 6 0-3 6 8.

[ 6] S O C O L O F S KYSA， A D AM SEE. M u l t i - p h a s ep l u m e s i nu n i -f o r m a n d s t r a t i f i e d c r o s s f l o w[J]. J o u r n a l o f H y d r a u l i c R e s e a r c h， 2 0 0 2， 4 0（ 6） ： 6 6 1-6 7 2.

[ 7] MA S UT AN ISM， A D AM SEE. E x p e r i m e n t a l s t u d yo fm u l t i -p h a s ep l u m e sw i t ha p p l i c a t i o nt od e e po c e a no i l s p i l s[ R]. H a -w a i ：H a w a i N a t u r a l E n e r g y I n s t i t u t e，U n i v e r s i t y o f H a w a i ， 2 0 0 1.

[ 8] B R AN D V I K PJ， J OHAN S E N？， L E I R V I K F， e ta l . D r o p l e t b r e a k u p i ns u b s u r f a c e o i l r e l e a s e s - P a r t 1： e x p e r i m e n t a l s t u d yo f d r o p l e tb r e a k u pa n de f f e c t i v e n e s so fd i s p e r s a n ti n j e c t i o n[ J]. M a r i n eP o l l u t i o nB u l e t i n， 2 0 1 3， 7 3： 3 1 9-3 2 6.

[ 9] Z HAOL， S HA F F E RF，R O B I N S ON B， e ta l . U n d e r w a t e ro i l j e t： h y d r o d y n a m i c sa n dd r o p l e ts i z ed i s t r i b u t i o n[ J]. C h e m i c a l E n g i n e e r i n gJ o u r n a l， 2 0 1 6， 2 9 9： 2 9 2-3 0 3.

[ 1 0] 王晶.海底管線溢油在水体中的运移扩散过程研究[ D].大连：大连理工大学， 2 0 0 6.

[ 1 1] 肖杰.海底管线微孔泄漏油滴尺寸与上升速度研究[ D].大连：大连理工大学， 2 0 0 7.

[ 1 2] 臧晓刚，张军，陈怀民，等.水下输油管微孔溢油的实验研究[ J].海洋技术学报， 2 0 1 6， 3 5（ 2） ： 1 0 4-1 0 8.

[ 1 3] 钱国栋，赵宇鹏，于顺.蓬莱1 9 - 3原油水下溢油的模拟实验[ J].海洋环境科学， 2 0 1 6， 3 5（ 6） ： 9 0 1-9 0 7.

[ 1 4] 邱章，方文东.南海北部春季海流的垂向变化[ J].热带海洋学报， 1 9 9 9， 1 8（ 4） ： 3 2-3 9.

[ 1 5] 杜岩，王东晓，陈荣裕，等.南海西边界A D C P观测海流的垂直结构[ J].海洋工程， 2 0 0 4， 2 2（ 2） ： 3 1-3 8.

[ 1 6] 王雪竹，李培良，张婷婷.南海中部深水海盆的潮流垂直结构[ J].海洋学报， 2 0 0 9， 3 1（ 1） ： 2 0-2 7.

[ 1 7] 廖国祥，高振会.水下溢油事故污染物输移扩散的数值模拟研究[ J].海洋环境科学， 2 0 1 1， 3 0（ 4） ： 5 7 8-5 8 2.

[ 1 8] 廖国祥，高振会，熊德琪.水下油气溢漏事故污染物输移预测模型[ J].大连海事大学学报， 2 0 1 0（ 4） ： 1 1 5-1 2 0.

收稿日期： 2 0 2 0 - 0 3 - 1 3；修订日期： 2 0 2 0 - 1 1 - 2 7

基金项目：国家自然科学基金项目（ 4 2 0 7 6 2 1 5、4 1 3 0 6 0 9 9） ；国家重点研发计划项目（ 2 0 1 8 Y F D 0 9 0 0 6 0 6）.

作者简介：廖国祥，副研究员，博士，研究方向为海洋溢油数值模拟与风险评估、滨海湿地与海洋自然保护地监测评估