横流环境中水下溢油输移扩散的实验研究

文章基于自制的组合式环形水槽（周长9.7 m、宽0.45 m、深1.0 m）,分别以阿曼原油及其消油剂混合物和淡水（含示踪剂）为模拟污染物,开展水下溢油的物理模拟实验,以浮射流输移轨迹、污染物扩散范围和油滴粒径分布为考察指标,研究横流环境和消油剂的使用对水下溢油输移扩散的影响。实验结果表明:在横流环境中,浮射流输移轨迹开始弯曲的高度随着流速的增加而降低;与淡水浮射流主要在水中输移扩散的情况不同,当污染物为原油时,大粒径油滴脱离浮射流主体并上浮至水面,导致扩散范围更大;消油剂的添加会使原油浮射流内部油滴的体积中值粒径变小,油滴粒径分布曲线向小尺寸方向偏移。实验结果可为后续的物理模拟实验和数值模拟研究提供参考。     

南海深水井喷溢油输移扩散的数值模拟研究

文章以假设发生在南海北部深水油气勘探开发海域的1279m深水井喷溢油事故为例,采用深海溢油输移扩散模型,以2020年冬季（1月）和夏季（7月）三维流场和海面风场为主要环境动力,数值模拟研究未实施和实施海底消油剂处理情景下的溢油在海洋环境中的输移扩散全过程。模拟结果显示:不同深度的海流对溢油输移轨迹有着显著的影响;深水环...     

“深水地平线”事故深海溢油输移扩散的数值模拟

为深入认识深海溢油输运过程和提高深海溢油事故的应急响应能力,文章以2010年墨西哥湾“深水地平线”事故为例,采用深海溢油输移扩散模型,以三维流场和海面风场为主要环境动力,数值模拟溢油深海泄漏后的浮射扩散、水体中输移扩散以及在风场和流场共同作用下在海面上输移扩散的全过程,同时模拟实施海底消油剂喷注处理措施后溢油输移轨迹和扩散范围的变化。研究结果表明：数值模拟结果与相关报道的悬浮油带实际观测结果以及美国国家海洋与大气管理局的海水异常遥感监测结果总体相符,可为更加全面和精细的深海溢油输移扩散数值模拟研究奠定良好的基础。     

深海环境中溢油输移扩散的初步数值模拟

通过对溢油在深海环境中的输移过程及行为特点的分析,初步建立基于拉格朗日积分法的深海溢油模型.该模型除了能够模拟油气混合物在真实深海环境中的共同输移与分离输移扩散过程,还考虑了石油溶解、气体溶解、天然气水合物形成与分解等行为变化对溢油运动轨迹的影响.应用该模型初步数值模拟了一次实际深海溢油试验,结果表明溢油在水下的空间分...     

基于VOF方法的海底管道溢油扩散数值模拟研究

在海底输油管道运行过程中,管道渗漏、穿孔及破碎都会导致原油泄漏。对溢油运动的轨迹及其扩散范围作出预报可为溢油事故的处理提供及时、准确的信息,指导应急处理的正确实施。基于工程实际需求,采用有限体积法,结合k-ε紊流模型,建立了海流作用下海底输油管道溢油扩散数值模型。采用VOF方法(volume of fluid method)追踪多相流界面。首先,将数值模拟结果与Fan的实验值及Zheng和Yapa的数值结果进行了对比,验证数值模型的可靠性;其次,研究了不同原油溢出速度与环境水深对不同时刻溢油轨迹、到达海面时间、横向漂移距离与海面扩散范围的影响。研究表明:随原油溢出速度增大,溢油到达海面时间逐渐减小,溢油横向漂移距离与海面扩散范围则逐渐增大;随环境水深增大,溢油到达海面时间逐渐增大,且其变化接近线性分布。     

水下溢油数值模拟研究

基于Lagrange积分法和Lagrange粒子追踪法建立了一个水下溢油数值模型。该溢油模型由两个子模型组成:羽流动力模型和对流扩散模型,其中羽流动力模型用以模拟溢油的喷发阶段和浮力羽流阶段;对流扩散模型用以模拟溢油的对流扩散阶段。通过数值实验,结合实验室水槽实验和水下溢油现场实验的观测资料进行模型验证。实验结果表明,模拟结果与观测资料一致性较好,从而验证了本文溢油模型的合理性和准确性;羽流动力模型为对流扩散模型提供源,海流、海水的垂向密度结构和油滴的直径分布是影响溢油在对流扩散阶段运动和分布的主要因素。     

平直斜坡上破波带内质量输移流对污染物输移影响的数值研究

破波带内外都有质量输移流存在，其对破波带内污染物输移有怎样的影响，需要进一步深入研究。本文基于实验以及考虑质量输移流的对流扩散数学模型研究了平直斜坡上破波带内质量输移流对污染物输移影响。数学模型包括波浪模型、近岸流模型以及对流扩散模型。首先建立了破波带内污染物输移数学模型，其中波浪场基于波能守恒方程来计算，波导流场基于Longuet-Higgins提出的辐射应力模拟，污染物对流扩散方程中考虑了质量输移流的影响，并利用算例验证该数学模型。其次简要介绍了平直斜坡上破波带内污染物输移实验，并分析了污染物输移特性。连续投放污染物会形成污染带，本文分析了两种波况下不同时刻污染带与岸线夹角的变化，以及污染物在垂直岸线和沿岸线方向的输移速度，结果表明对两种波况来说在初始10-40s污染团向岸线方向输移速度分别约为0.05m/s、0.017m/s，之后速度分别减小为0.001m/s、0.011m/s。数值模拟结果与实验结果比较表明：考虑质量输移流的模拟结果与实验更为吻合。因而，通过实验以及数模研究表明破波带内质量输移流对破波带内污染物在垂直岸线方向的输移有重要影响，而对沿岸方向的输移则影响较小。     

渤海溢油三维漂移数值模拟研究

国家海洋局北海预报中心于2011年开发了渤海三维溢油模型,该模型在分析国内外溢油模型现状的基础上,借鉴当今流行的数值模拟方法,使用油粒子模型与油膜扩展模型相结合的方式,用拉格朗日方法追踪每个带有一定油量的油粒子的轨迹,针对每一个油粒子则使用油膜扩展理论计算其油膜扩展过程。该模型可实现对溢油油污上升及水平输运过程、海表面油污浓度的预报,通过三组理想试验和2012年的海上溢油实验数据,对模型的各项功能、稳定性及模型精度进行了对比验证,结果较好,模型可实现对渤海海域海底或水下发生溢油的数值模拟。该模型解决了以往二维溢油模型在模拟钻井平台及海底输油管道泄漏等溢油事故方面的不足,可更好地为溢油灾害对海洋环境影响的估计提供有效参考信息。     

水下输油管微孔溢油的实验研究

目前大部分溢油研究都是针对大孔径泄漏进行的。文中对水下管道微孔(直径小于5 mm的孔)泄漏溢油进行了实验研究,获得了不同尺寸微孔泄漏与出流系数的关系、不同泄漏量时溢油油滴上升至水面的时间以及泄漏产生的浮羽射流的形态。实验结果表明,出流系数会随着泄漏孔孔径的增大而变大,而泄漏量的增加则会明显缩短溢油上浮至水面的时间,增大羽射流宽度。研究对水下微孔泄漏的溢油扩散提供了实验参考,为今后微孔溢油的相关研究打下基础。     

水下管道微孔溢油扩散特性及油滴尺寸分布

对水下油管微孔泄漏进行了有关数值模拟及实验研究。首先采用计算流体动力学(CFD)方法对微孔泄漏后的油气两相混合物的运动特性进行了数值模拟,模拟结果表明随着含气率的增大,溢油在海表面横向扩散范围逐渐增大,且其变化接近线性分布;泄漏速度对溢油羽射流形态影响较大,随着泄漏速度的增大,水下溢油羽射流的截面宽度逐渐增大,油滴上升速度较快,且油滴直径不断减小;其次,对泄漏后的油滴尺寸进行了实验测试,实验结果表明泄漏量对油滴尺寸大小有较大的影响,随着泄漏量的增大,油滴获得了更大的出流速度,形成更多更小的油滴,当泄漏量增大到一定程度时油滴逐渐趋于雾化。     

Hydrodynamics of oil jets without and with dispersant: Experimental and numerical characterization

In this paper, we present the analysis of an underwater horizontal oil jet experimental measurement and Computational Fluid Dynamics (CFD) using the Reynolds Averaged Navier Stokes (RANS) equations. Two oil subsurface releases were conducted: one with crude oil and another with crude oil premixed with dispersant at the dispersant to oil ratio (DOR) of 1:20. The jet profile was captured by a camera at moderate resolution, and the instantaneous velocity was measured by a Vectrino Profiler. The velocity components, turbulence kinetic energy, and turbulence dissipation rate from the experiment agreed well with those from the CFD simulation using the k-epsilon turbulence model. The spread angle of the jet was found to be around 21° and 24° from the experiment measurement, for oil without dispersant and oil with dispersant, respectively. The latter is close to the angle of miscible jets at 23°. The jet profile of oil with dispersant had a smaller buoyancy than that without dispersant, which is probably due to the large water entrainment for the oil with dispersant jet. The cross sections of the jet for both cases gradually became flattened with distance, as the plume turned upward.     

海上油气开发水下井喷溢油的化学分散技术研究

为有效应对海上油气开发水下井喷溢油事故,降低海底溢油对海洋环境的危害,文章通过模拟试验研究消油剂的类型、使用量和喷注位置对溢油分散效果的影响,并提出深水水下消油剂使用技术体系。研究结果表明：GM-2消油剂对原油呈现良好的分散性能,且分散效果随剂油比的增大而增强;水下喷注消油剂的使用量为溢油量的10%时即可获得良好的分散效果;在一定范围内,消油剂喷口与水下井口的水平距离对溢油分散效果的影响极小,而垂直高度对溢油分散效果的影响较大。     

Underwater spreading and surface drifting of oil spilled from a submarine pipeline under the combined action of wave and current

Tremendous economic loss and environmental damages are caused by oil-spilling accidents in sea. Accurate prediction of the underwater spreading and surface drifting of oil spills is important for the emergency response. In the present study, numerical investigation on the underwater spread and surface drift of oil spilled from a submarine pipeline under the combined action of wave and current was carried out to examine the effects of physical ocean environment, leaking flux and spilled oil density and viscosity. Reynolds-Averaged-Navier-Stokes (RANS) equations, realizable k-ε turbulence model and volume of fluid (VOF) model are employed to describe the multiphase flow, and velocity-boundary wave-making technique combined with the sponge layer damping absorber technique realizes the numerical wave flume. Oil spill experiments were conducted to validate the numerical model. The calculation results indicate that compared with the environmental conditions of still water, only current and only wave, a larger scope of underwater spreading and relatively slower rising rate and relatively faster drifting rate of oil droplets are observed under the combined action of wave and current. The leaking flux affects the floating time and dispersion concentration, while the ocean environment affects the horizontal migration and surface drifting. Under the specific conditions of present work, oil density has obvious effect on the underwater spread but limited effect on the surface drifting, while oil viscosity has little effect on both the two processes.     

Surface manifestations of internal waves investigated by a subsurface buoyant Jet: 1. The mechanism of internal-wave generation

In a large test reservoir with artificial temperature stratification at the Institute of Applied Physics, Russian Academy of Sciences, we have performed a major laboratory simulation of the nonstationary dynamics of buoyant turbulent jets generated by wastewater flows from underwater collector diffusers. The interaction of buoyant jets with the pycnocline leads to an active generation of internal waves. An analysis of the dependence of wave amplitude on the control parameter proportional to the rate of liquid flow from the collector diffuser has indicated that this dependence is adequately described by a function that is characteristic for the presence in the Hopf bifurcation system, which occurs for a soft actuation mode of self-oscillations of the globally instable mode. To check the conditions for the actuation of the globally instable mode, we have performed an auxiliary experiment in a small reservoir with a salt stratification formulated similar to the experiment in the big reservoir. Using the particle image velocimetry (PIV) method, we have measured the velocity field in the buoyant jet and constructed the profiles of transverse velocity in several sections. When the jet approaches the pycnocline, a counterflow is generated at the edges. A stability analysis for the resulting profiles of flow velocities performed by the method of normal modes has revealed that, for the jet portions with counterflow, the condition of absolute instability by the Briggs criterion for axisymmetric jet oscillations is satisfied, which testifies to the fact that the globally instable mode is actuated. The estimates for oscillation frequencies of the globally instable mode are well consistent quantitatively with the measured spectrum of jet oscillations.     

Numerical study on the flow characteristics of tidal jets induced by a tidal-jet-generator

The flow characteristics of tidal jets induced by a Tidal-Jet Generator (TJG) are investigated using a finite-difference numerical scheme, named Navier–Stokes (NS)–Marker and Cell (MAC)-TIDE, based on the fully 3D NS equations. The TJG is an enclosed rectangular breakwater, which has vertical opening and a large enclosed volume inside. During both phases of tide, strong and uni-directional jets can be obtained locally from the inlet of the TJG, due to the water level difference between the inner and outer sides of TJG.The computed results are extensively compared with three other independently developed numerical models; 3D-ADI, DVM, and CIP-CSF. These models are based on quasi-3D, 2D depth-averaged, and fully 3D NS equations, respectively. It is seen that the present fully 3D numerical model NS–MAC-TIDE can predict the maximum intensity of inlet velocity with higher accuracy than the other numerical models when compared with the empirical function proposed from the experiments. The numerical simulations based on NS–MAC-TIDE can reproduce successfully the processes of generation, development, and dissipation of tidal jets. The effects of gap opening on the main characteristics of the tidal jet flow are assessed. Through numerical assessment, it is also clearly demonstrated that the residual time of a pollutant distributed around the front of the TJG can be decreased by significant amount due to the locally induced tidal jet. The TJG can thus utilize tidal energy for water purification in local marine environment by providing a flushing mechanism.     

The interaction between multiple high pressure combustion gas jets and water in a water-filled vessel

A static experimental system based on the background of the underwater gun firing process is designed to simulate the expansion of conical-distributed multiple gas jets in a water-filled vessel. The gas is generated from the gunpowder combustion and injected into water through five orifices on a mock bullet. A high-speed digital camera is taken in experiment to record such expansion and the Euler-Euler multiphase model is utilized in simulation to describe the gas-liquid flow. Results show that the contraction of Taylor cavity in liquid and the gas-liquid entrainment both affect the shock structure in gas. The expansion zone near the orifice shrinks and the gas speed downs when the central gas jet has a contraction, but lateral jets have no contractions due to the obstruction from the wall. The contraction has a greater impact on the shock structure than the gas-liquid entrainment and occurs earlier under a higher injection pressure.