2. 中国石化工程建设有限公司, 北京 100101
海洋石油开发及运输伴随着原油泄漏污染问题,给环境和社会带来了极大影响[1,2,3]. 为此,对于污染区域的治理显得非常重要,为防污染区域随海浪作用或海洋风作用飘移,导致污染治理目标区域遗失,可采取 海上浮油跟踪定位浮标定位污染区域[4,5,6],这种装置由卫星定位发射装置和油污吸附承载装置两大部分组成(图 1). 卫星定位发射装置安装在油污吸附承载装置中,为保证发射信号强度,同时降低浮标承受的风载荷,发射装置微露水面;承载装 置的主要作用部分是一种高性能的吸油材料[7,8,9],其主要原材料为羽毛,对各种黏度油污具有强烈的吸附能力,吸油比达到1:18$\sim $1:20,吸附后不易离析,并有良好的耐海水腐蚀能力;次要材料为羽毛固结材料,确保在海洋波浪和海洋风作用下,吸油材料不会分散. 承载装置采用大小圆圈相套外形设计,两圆圈之间利用双层高强度丝网夹持并固定吸油材料,可以确保在海洋风和波浪的双重作用 下总体形态不会发生变化,吸附的原油可以较好地保存在吸附材料中,从而使装置与油污保持性质类似.
在发生海上漏油事件后,可以将海上浮油跟踪定位浮标投放到污染(漏油)区域,承载装置中的吸油材料吸附到海面原油上后,在吸附力作用下会将浮标``锁定''在油污面上,使得浮标与污染区域一起移动,其后根据卫星定位发射装置的定位作用就可以锁定污染区域飘移轨迹,实现溢油跟踪目标,以便后续的石油回收和环境清理.
在海洋风作用下,特别是恶劣海况条件时,若要保证定位浮标与污染区域同步运动,则需要浮标中吸油材料具备足够的承载能力,在不考虑吸油材料的抗分散能力下,吸油材料受力基本条件应该满足能够抵抗海洋风产生的风载荷[10,11,12],即要求吸油材料与海面原油之间的黏附力高于吸油材料承受的风载. 此外,由于定位浮标随船携带,使用保存必须灵活方便,且不能占用太大空间,因此在满足能够抵抗飓风产生风载的条件下,希望浮标的直径越小越好,故要求对吸油材料的直径进行优化设计.
本文拟通过室内实验,首先对吸油材料与原油之间的黏附力随吸油材料直径的变化关系进行研究,并且分析黏附力的组成部分,研究其主要作用力,在此基础上可以对吸油材料的材质进行优选;其次,对吸油材料受到的风载荷随其直径以及风速变化的规律进行研究,以期将变化规律扩展至海洋风达到的风速条件下;最后结合黏附力必须要高于风载荷的原则,在考虑恶劣海况条件下,对吸油材料的直径进行优化设计.
1 吸油材料受力分析吸油材料直接接触的介质包括原油和空气,空气的流通反映到材料上即风载,原油的流动作用于材料即黏附力,风载荷将吸油材料推离其附着区域,而黏附力则起到固定作用. 黏附力为各向同性,水平方向与垂直方向受力等同;风载荷为水平方向的力,将两者均考虑为水平方向的力时,可以不考虑受力方向,因此文中所涉及到的力均为标量.
1.1 风载荷风载荷是风压与受力面积的乘积,油污吸附承载装置为圆盘形,所受风载面积与风向没有关联,受力面积是一定的,风载荷只和风压相关. 风压计算式[13,14]如式(1)
$ P=P_m-P_b=\dfrac 12 \dfrac{\gamma}{g} v_{\rm a}^2 $ | (1) |
1个大气压、15${^\circ}$C和干燥的条件下, $\gamma =0.012 018$ kN/m$^{3}$,重力加速度取值为9.8 m/s$^{2}$,代入式(1)得到风压公式
$P=\dfrac 12 \dfrac \gamma g v^2_{\rm a} =\dfrac{v^2_{\rm a}}{1 630} $ | (2) |
式中,$\gamma $为空气单位体积的重力,$g$为重力加速度,$v_{\rm a}$为风速,$P$的单位为kN/m$^2$.
1.2 黏附力吸油材料与原油之间存在相互作用力,包括分子间作用力、液体表面张力和剪切黏性力,黏附力是这3种力的综合作用[15]. 其中分子间作用力是油分子与材料分子之间的相互作用,液体表面张力主要指吸油材料受到的源于油膜的液体表面张力,剪切黏性力是原油流动导致吸油材料受到的摩擦和压差作用力. 由空气引起的材料表面黏附力主要来源于毛细力,毛细力的大小与水膜厚度、相对湿度以及水膜自由能有关[16,17,18],考虑到实验在常规室内环境下进行,由空气相对湿度引起的黏附力很小,在文中不考虑空气引起的吸油材料表面黏附力.
(1)分子间作用力
分子间作用力分为吸引力和排斥力[19],吸引力分为van der Waals力、Casimir力、排空力、疏水相互作用等,排斥力分为双电层力、溶剂化力、空间位型排斥力等;其中van der Waals力包括取向力、诱导力和色散力[20]. 分子间作用力在宏观层面上无法通过计算公式直接计算得到,本文通过材料受力平衡计算分子间作用力,具体过程见2.2节.
(2)液体表面张力
由于体系的能量越低体系越稳定,故液体表面具有自动收缩的趋势,我们可以把这种趋势看作液体表面分子相互吸引的结果,此吸引力与表面平行,它的大小反映了表面自动收缩的趋势大小,称其为液体表面张力[21,22,23]. 当吸油材料与原油相接触之后,原油液面的原有平衡被打破,原油为了恢复原有状态,将体系能量将至最低,必然会对吸油材料产生一个作用力,即源于液体表面张力的力作用. 源于液体表面张力的力可以用式(3)来计算
$ f=2\sigma L $ | (3) |
式中,$f$为外力,$\sigma $为液体表面张力,$L$为材料周长.
(3)剪切黏性力
在风载荷和海浪的作用下,吸附材料四周以及下部的原油发生流动,对材料产生摩擦作用,并顺着液流方向产生压差作用,将这两种作用产生附加力称为剪切黏性力.
从海洋流体力学角度讲,这个剪切黏性力可以看作是水流对吸油材料的拖曳力[24,25,26,27],拖曳力计算式为
$ F_{\rm D} =\dfrac 12 C_{\rm D} \rho A v^2 $ | (4) |
式中,$C_{\rm D}$为拖曳力系数,圆盘为1.2;$A$为圆盘浸入水体部分垂直于流动方向的投影面积;$\rho $为水密度;$v$为水流速度.
实验中,水流速度是未测量的,可以通过海洋流体力学知识求得. 在风的作用下,水体表面的水流速度为
$ v=\dfrac{T_{\rm A}}{\sqrt{2\rho A_Z \omega \sin \varphi}} $ | (5) |
$T_{\rm A}$可以根据作用力与风速的关系求得,风对水体的作用力与风速有以下关系
$ T_{\rm A}=C\rho_{\rm a} v^2_{\rm a} $ | (6) |
式中,$T_{\rm A}$为风对水的作用力,$C$为系数,$\rho_{\rm a}$为空气密度, $A_{Z}$为垂直涡动动力黏滞系数,$\omega $为地球自转角速度,$\varphi $为实验地点所处纬度.
结合式(4)$\sim $(6)得出剪切黏性力的计算公式
$ F_{\rm D}=\dfrac{C_{\rm D}AC^2 \rho^2_{\rm a}v^4_{\rm a}}{4A_Z\omega \sin \varphi} $ | (7) |
为了得出吸油材料所受风载荷、黏附力与直径之间的关系,将材料制作成直径为8 cm,10 cm,12 cm,16 cm,20 cm, 24 cm,28 cm 的圆盘(图 2),针对不同直径圆盘测量其在不同风速下的风载,并测量其与实验用油之间的黏附力.
风载荷和风压与受力面积相关,受力面积取决于材料圆盘直径和厚度,风压由风速确定.
(1)风速测量
实验中风扇产生的风速分为3个档位:高档、中档和低档,利用皮托管测量风速[28,29,30],如图 3.
由于U形管面板具有一定角度,直接测量的高差 $\Delta H$是倾斜面高差,不能直接用于风速计算,需要将此高差转变为垂直方 向的高差 $\Delta H_{\rm h}$,通过此垂直高差可以计算出风压,再由式(2)得出式(9),即可计算出风速. 高差公式为
$ \Delta H_{\rm h}=\Delta H \sin \phi=\Delta H \sin 5^\circ $ | (8) |
$ \rho g \Delta H_{\rm h} =\dfrac{v^2_{\rm a}}{1.63} $ | (9) |
风速测量实验数据见表 1.
(2)风载测量
由于实验材料所受风载值小,实验中采用杠杆原理放大风载以便测量. 如图 4,给漂浮细长条标上尺寸,在一端用细线系上材料圆盘,另一端使用两个弹簧秤在指定尺寸处,分两边勾住长条. 使长条处于水平状态,以此测定材料所受风载$F$. 风载测量采用量程2 N,最小刻度值为0.01 N的弹簧秤.
$ F=F_2-F_1 $ | (10) |
风载荷实验原始数据如表 2.
表 2中空白处没有值,其一在于圆盘直径小,承受风压面积小;其二在于低档时风速小,产生风压小,使得圆盘受到风载荷小,测量不准确.
图 5是风载与直径关系图,可以看出风载荷与圆盘直径呈线性关系.
在进行风载荷实验过程中,风不仅仅对圆盘有作用力,而且对圆盘周边及下方的水体也有推动作用,所以实际测量得到的风载 中包含有圆盘受到的风载以及流体对圆盘的剪切黏性力,所以上图直线不经过原点. 通过风载荷计算公式将实际风载部分计算出来,计算结果如表 3.
实验中采用十字形硬铁丝固定吸油材料,采用弹簧秤上提的方式,得到提出前后弹簧秤的示数,提出前弹簧秤示数减去提出 后示数为黏附力. 黏附力测量分别采用量程为5 N,2.5 N和2 N,最小刻度值为0.05 N,0.05 N和0.01 N的弹簧秤.
(1)黏附力测量
由于吸油材料吸油能力太强,吸油过多会影响黏附力的测量,所以在实验中限制吸油材料浸入油中不超过5 mm.
表 4中的黏附力就是圆盘提出前测重与提出后测重的差值,图 6是表征圆盘直径与黏附力的关系,其拟合方程为式(11).
$ Y=0.003 4 x^2-0.021 3x +0.023 6 $ | (11) |
图 6说明吸油材料和油之间的黏附力与吸油材料的直径呈二次方程关系,圆盘直径的增大,黏附力会急剧增大.
(2)液体表面张力
查找资料,发现油的液体表面张力一般在40$\sim $46 mN/m之间波动[31],实验原油为辽河油田产出原油,本文计算 液体表面张力取值为42 mN/m.
根据式(3)计算得出源于液体表面张力的力如表 6所示.
(3)剪切黏性力
剪切黏性力是流体流动对圆盘的影响,在风力作用下,流体与圆盘之间有相对运动,产生剪切作用,这个剪切力增大了我们实验中 测量的实际风载. 实际测量风载与理论计算风载的差值就是这个剪切黏性力. 由于实测风载数据并不完整(实验条件的限制),所以对于剪切黏性力的计算,只计算最大风速下的剪切黏性力,剪切黏性力计算结果如表 7.
表 7给出了根据实验结果计算出的剪切黏性力值,而式(7)给出的是理论计算公式,在式(7)中,$C$与$A_{Z}$是未知量常量, 可 以通过实验数据来确定这个未知量.
式(7)中,$\rho_{\rm a}=1.226$ kg/m$^{3}$,$v_{\rm a}$取值为最大风速4.358 m/s,$\omega $取值为7.27$\times $10$^{ - 5}$ rad/s,$\varphi $取值40$^\circ$,$A$取值为圆盘直径$D$与浸入水中厚度的乘积. 圆盘浸入水中的厚度可以根据阿基米德原理求得,其值为6.6$\times $10$^{ - 5}$ m. 最后得到剪切黏性力计算公式为
$ F_{\rm D}=100 \dfrac{c^2}{A_Z} D=K D $ | (12) |
由上式可以发现剪切黏性力最终与圆盘直径$D$呈正比,将表 7中实际剪切黏性力第1行数据代入式(12),可以得到$K =0.118$. 从而
$ F_{\rm D}=0.118 D $ | (13) |
由式(12)可以得到表 7中的理论剪切黏性力.
对比实际与理论剪切黏性力可知,在圆盘直径较大($ > 12$ cm)时,理论计算数值与实际测量计算值十分吻合. 说明式(12)适合于圆盘直径较大时的计算.
(4)分子间作用力
在2.2节中,前文给出了总的黏附力、液体表面张力和剪切黏性力,利用总黏附力减去液体表面张力和剪切黏性力得到分子间作用力,如表 8.
比较表 6$\sim $表 8发现相比于分子间作用力来说,液体表面张力和剪切黏性力太小. 故可以认为在吸油材料与油相互作用中,以分子间作用力为主导影响因素. 根据此结论,可以优选分子间作用力较大的吸油材料.
3 吸油材料圆盘直径优化前文风载荷都是以实验室风扇的风速为基础计算出来的,优化计算中可以扩大风速范围,这里取风速为4 m/s,6 m/s, 8 m/s,10 m/s,12 m/s,14 m/s,16 m/s,18 m/s,20 m/s,22 m/s,24 m/s,26 m/s,28 m/s,30 m/s, 32 m/s,34 m/s. 海洋工程中32.7 m/s以上的风速是14级飓风,在这种情况下,即使发生漏油事件,也不可能去进行定位工作,所以本文最多只计 算到34 m/s的风速.
图 6给出了圆盘直径与黏附力的关系曲线,以式(11)为基础,可以计算圆盘直径从0到60 cm所有的黏附力值.
根据风压计算公式,计算出各个风速下的风压,然后乘以圆盘的受力面积,就可以得到风载荷.图 7给出一定风速下的风载荷和对应圆盘直 径下的剪切黏性力.
从该图可以看出以下几点:
(1)风速越大,风载越大,且同一圆盘直径下,风速差相同时,风载差值变大.
(2)圆盘直径越大,风载呈线性增大,黏附力呈二次关系增大,黏附力增大速率远远高于风载.
(3)黏附力与每一条风载直线都有交点,说明在此处圆盘半径是临界半径,若是圆盘直径大于这个临界直径,则黏附力大于风载,从而吸附装置可以很好的与油污“黏”住. 若圆盘直径小于临界直径,则风载大于黏附力,定位装置会与油污脱离,达不到定位作用.
(4)如果考虑海洋气候的极限情况,即风速达到34 m/s,黏附力曲线与风载曲线的交点的横坐标在38 cm处. 也就是说,在不考虑其他因素影响,只考虑吸油材料本身,当吸油材料的直径超过38 cm时,吸油材料就可以抵抗几乎所有的海洋复杂气候.
4 结 论(1)分析了吸油材料的受力状态,得出吸油材料受到风载荷和黏附力的共同作用,黏附力包括分子间作用力、液体表面张力和剪切黏性力,其中以分子间作用力为主.
(2)得到了风载荷和黏性力的测量方法,并得到了分子间作用力宏观测量方法.
(3)通过实验得出随圆盘直径变化,剪切黏性力和在不同风速下风载荷的变化趋势,并将其扩展于大风速和大圆盘直径下,优化出适合于海况的圆盘直径.
1 | 陈涛. 美国海洋溢油事件的社会学研究. 中国农业大学学报(社 会科学版),2012,02:110-117 (Chen Tao. Sociological study on American oil spill event.China Agricultural University Journal of Social Sciences Edition, 2012, 02: 110-117 (in Chinese)) |
2 | Gill DA, Picou JS, Ritchie LA. The exxon valdez and BP oil spills: A comparison of initial social and psycho-logical impacts.American Behavioral Scientist, 2012, 56(1): 3-23 |
3 | 李京梅,侯怀洲,姚海燕.基于资源等价分析法的海洋溢油生物 资源损害评估.生态学报,2014,34(13):3762-3770 (Li Jinmei, Hou Huaizhou, Yao Haiyan. Marine biological damage assessment of oil spill based on resources equivalency analysis.Acta Ecologica Sinica, 2014, 34(13): 3762-3770 (in Chinese)) |
4 | 杨瑞,顾群,王玉林. 水上溢油应急装备跟踪定位浮标关键技术 研究. 中国海事,2014,04:51-53 (Yang Rui, Gu Qun,Wang Yulin. Research on the key technology of tracking location buoy of the marine oil spill emergency equipment.China Marine, 2014, 04: 51-53 (in Chinese)) |
5 | Hurford N. The use of a satellite-tracked buoy to monitor the movement of an oil slick at sea.Oil and Chemical Pollution, 1990: 6(2):101-112 |
6 | 杨瑞, 顾群. 基于北斗卫星的溢油跟踪浮标. 上海海事大学学报,2014, 03: 23-27 (Yang Rui, Gu Qun. Spilled oil tracking buoys based on Beidou satellites.Journal of Shanghai Maritime University,2014, 03: 23-27 (in Chinese)) |
7 | 彭丽,刘昌见,刘百军. 天然有机纤维吸油材料的研究进展. 化工 进展,2014,02:405-411 (Peng Li, Liu Changjian, Liu Baijun. Research progress of natural organic fibers for oil absorption material.Chemical Industry and Engineering Progress, 2014, 02: 405-411 (in Chinese)) |
8 | Ceylan D, Dogu S, Karacik B, et al. Evaluation of butyl rubber as sorbent material for the removal of oil and polycyclic aromatic hydrocarbons from seawater.Environ Sci Technol,2009, 43(10):3846-3852 |
9 | Wang JT, Zheng Y, Wang AQ. Coated kopak fiber for removal of spilled oil.Marine Pollution Bulletin, 2013, 69: 91-96 |
10 | 岳晓瑞,徐海祥,罗薇. 海洋工程结构物风载荷计算方法比较. 武 汉理工大学学报(交通科学与工程版),2011,03:453-456 (Yue Xiaorui, Xu Haixiang, Luo Wei. Comparison of wind load calculation methods for marine engineering structures.Journal of Wuhan University of Technology, 2011, 03: 453-456 (in Chinese)) |
11 | Haddara MR. Wind loads on marine structures.Marine Structures,1999, (12): 199-209 |
12 | OCIMF. Prediction of wind and current loads on VLCCs. Oil Companies International Marine Forum, 1994 |
13 | 张相庭. 结构风压和风振计算. 上海:同济大学出版社,1985:4-5 (Zhang Xiangting. Calculation of Structural Wind Pressure and Vibration.Shanghai: Tongji University Press,1985: 4-5 (in Chinese)) |
14 | Nyi Nyi Aung, Ye Jihong. Coherence of wind pressure on domes.Journal of Southeast University, 2010, 26(1): 100-106 |
15 | Fay JA. Oil on the Sea. New York: Plenum Press, 1969: 53-63 |
16 | 赵亚溥. 纳米与介观力学. 北京:科学出版社,2014:342-351 (Zhao Yapu. Nano and Mesoscopic Mechanics.Beijing: Science Press, 2014: 342-351 (in Chinese)) |
17 | 卿涛, 邵天敏, 温诗铸. 相对湿度对材料表明粘附力影响的研究. 摩擦学学报, 2006, 26(4): 295-299 (Qing Tao, Shao Tianmin, Wen Shizhu. E ects of relative humidity on surface adhesion.Tribology,2006, 26(4): 295-299 (in Chinese)) |
18 | Patton ST, Eapen KC, Zabinski JS. E ect of adsorbed water and sample aging in air on the 1N lever adhesion force between Si(100) and silicon nitride.Tribology International, 2001, 34: 481-491 |
19 | 赵亚溥. 表面与界面物理力学. 北京:科学出版社,2012:40-77 (Zhao Yapu. Surface and Interface Physics Mechanics.Beijing: Science Press, 2012: 40-77 (in Chinese)) |
20 | 李发胜. 有机化学. 北京: 科学出版社, 2012: 8-9 (Li Fasheng. Organic Chemistry. Beijing: Science Press, 2012: 8-9 (in Chinese)) |
21 | 黄志宇,张太亮,鲁红升. 表面及胶体化学. 北京:石油工业出 版社, 2012:1-4 (Huang Zhiyu, Zhang Tailiang, Lu Hongsheng. Surface and Colloid Chemistry. Beijing: Petroleum Industry Press,2012: 1-4 (in Chinese)) |
22 | 尹东霞, 马沛生, 夏淑倩. 液体表面张力测定方法的研究进展. 科 技通报, 2007, 23(3): 424-429,433 (Yin Dongxia, Ma Peisheng, Xia Shuqian. Progress on methods for measuring surface tension of liquids.Bulletin of Science and Technology, 2007, 23(3): 424-429,433 (in Chinese)) |
23 | Paddy JF, Pitt AR, Pashley RM. Menisci at a free liquid surface: surface tension from the maximum pull on a rod.J Chem Soc Faraday Trans, 1975, 71: 1919-1931 |
24 | 张远君. 流体力学大全. 北京:北京航空航天大学出版社,1991:226-228 (Zhang Yuanjun. Hydrodynamics. Beijing: Beijing University of Aeronautics and Astronautics Press, 1991: 226-228 (in Chinese)) |
25 | 邓邵云. 桩基绕流阻力特性研究现状与展望. 水运工程,2006,1:10-15 (Deng Shaoyun. Current situation and prospect of study on drag force characteristics of piles.Port & Waterway Engineering,2006, 1: 10-15 (in Chinese)) |
26 | 徐远杰, 康欢, 楚锡华. 漂浮物对流体拖曳力系数影响的试验研 究. 中国力学大会论文集, 2011. 1-9 (Xu Yuanjie, Kang Huan, Chu Xihua. Influence of floating-debris on fluid drag force using experimental investigation.In: Chinese Mechanical Assembly, 2011. 1-9 (in Chinese)) |
27 | 陈志清. 河道输沙学导论. 徐氏基金会, 1977, 1: 11-12 (Chen Zhiqing. Introduction to sediment transport in river.Xu Foundation,1977, 1: 11-12 (in Chinese)) |
28 | 王敏,魏根宝,马锦国. 二等皮托管测风误差分析及偏差模拟计 算. 气象科技,2013,05:843-846,851 (Wang Min, Wei Genbao, Ma Jinguo. Measurement error analysis and deviation simulation of second-class pitot tube.Meteorological Science and Technology,2013, 05: 843-846,851 (in Chinese)) |
29 | 李涛, 樊荣. 皮托管在瓦斯抽放煤层气输送计量中的应用. 自动 化仪表, 2015, 01: 95-98 (Li Tao, Fan Rong. Application of pitot tube in transport metering of coalbed methane gas drainage.Process Automation Instrumentation, 2015, 01: 95-98 (in Chinese)) |
30 | 王晓蕾, 沈建华, 张玉存. 皮托管风速测量不确定度的分析与评 定. 计量技术, 2003, 9: 53-54 (Wang Xiaolei, Shenjianhua, Zhang Yucun. Analysis and evaluation of pitot tube wind speed measurement uncertainty.Measurement Technology, 2003, (9): 53-54 (in Chinese)) |
31 | 王家禄. 油藏物理模拟. 北京:石油工业出版社, 2010. 35 (Wang Jialu. Reservoir Physical Simulation. Beijing: Petroleum Industry Press, 2010. 35 (in Chinese)) |
2. Sinopec Engineering Incorporation, Beijing 100101, China