【海上溢油漂移扩散】海上溢油漂移扩散的数值模拟+无人机溢油边缘定位附matlab代码

⛄ 内容介绍
石油是一种最为典型的海洋污染物，随着近年 来海上石油开发和海上运输量的日益增加，给经济 发展带来巨大效益的同时也加大了海洋的油污染 隐患。若发生溢油事故，泄漏的原油就会在水动力 的推动下发生漂移扩散，所以建立科学准确的溢油 模型是对事故进行重演与预测的重要手段，也是溢 油应急反应系统的重要组成部分，模型的建立能够 预测油膜在扩散漂移过程中性质和状态的变化及 最终归宿，能够为应急决策的制定和清除手段的选 择及溢油损害的评估提供定量化的依据。采用油、水分离的 Lagrangian 粒子追踪模型，对渤海海域的溢油漂移扩散问题进行了模拟。并将溢油边界采用无人机进行定位。

⛄ 部分代码

clear all

%% 根据经纬度创建部分海上地图

m_proj('Equidistant Cylindrical','lat',[5 40],'long',[110 155]); %投影方式及绘图范围设定

m_gshhs('lc1','patch',[0.98 0.98 0.98],'EdgeColor',[0.98 0.98 0.98]); %线条及色块颜色设定

m_grid('box','none','linest','none','linewidth',0.1,'tickdir','out','backcolor',[0.7020 0.8588 0.9098]);

set(gcf,'GraphicsSmoothing','off');

%% 创建a，b，c参数，服从正态分布

a=normrnd(0,1/3,[650,1]);

b=normrnd(0,1/3,[650,1]);

c=normrnd(0,1/3,[650,1]);

%% 初始海上溢油位置及浓度

x=135*ones(650,1)*0.00001;

y=20*ones(650,1)*0.00001;

z=1.25*ones(650,1)*0.00001;

nnn=35; %%******************此处修改时间5、10、15、20min*****************************************

%% 扩散模拟

for t=1:nnn

lx=a*(120*20)^0.5;

ly=b*(120*20)^0.5;

lz=c*(120*20)^0.5;

x=(x+12+0.06*y+lx);

r1=x/250;

r1=131+r1;

y=(y+0.02*50*60+ly);

r2=y/250;

r2=r2+15;

z=z+1.2+lz;

r3=abs(z/100);

end

hold on

%m_plot(r1,r2,'*','Markersize',1.5);

XYZ=[r1,r2,r3];

datamin = min(min(XYZ(:, 1)), min(XYZ(:, 2)));

datamax = max(max(XYZ(:, 1)), max(XYZ(:, 2)));

datamin = floor(datamin);

datamax = ceil(datamax);

Length = 600;

Width = 600;

axismin = datamin;

axismax = datamax;

% 绘制密度图

XYZ(:, 3) = XYZ(:, 3) - XYZ(:, 3);

sizeXYZ = size(XYZ);

searchR = 2;

for i = 1 : sizeXYZ(1)

index_i = find(XYZ(:, 1) > XYZ(i, 1) - searchR & XYZ(:, 1) < XYZ(i, 1) + searchR ...

& XYZ(:, 2) > XYZ(i, 2) - searchR & XYZ(:, 2) < XYZ(i, 2) + searchR);

sizeIndexI = size(index_i);

XYZ(i, 3) = sizeIndexI(1);

end

[sortXYZ, sortI] = sort(XYZ(:, 3));

hold on

m_scatter(XYZ(sortI, 1), XYZ(sortI, 2), [5], XYZ(sortI, 3), 'filled');

gc = get(gca);

set(gcf, 'position', [100, 100, Length, Width]);

AAA=[XYZ(sortI,1),XYZ(sortI,2),XYZ(sortI,3)];

bbb=AAA(1:nnn,:)

⛄ 运行结果

5min

15min

25min

35min

⛄ 参考文献

[1] 严志宇, 熊德琪, 殷佩海. 海上溢油风化模型评述[J]. 大连海事大学学报, 2001, 27(4): 36-39. YAN Zhi-yu, XIONG De-qi, YIN Pei-hai. Weathering modeling review of oil spills at sea[J]. Journal of Dalian Maritime University, 2001, 27(4): 36-39.

[2] 李琼. 海上溢油风化预测系统的研制[D]. 大连海事 大学, 大连, 中国, 2002. LI Qiong. Development of oil weathering system[D]. Dalian Maritime University, Dalian, China, 2002.

[3] MILLER, STOLZENBACK K D. A review and evaluation of basic techniques predictiong the behavior of surface oil slicks[R]. Sea Grant Program Report No. MITSG-77-9. Cambridge, United States: Massachusetts Institute of Technology, 1977, 22: 20-30.

[4] WILLIAMS G N, HANN R W. Simulation models for oil spill transport and diffusion[C]. Summer Computer Simulation Conference, La Jolla, California, United States, 1975: 748-752.

[5] 黄娟. 渤海海上溢油漂移扩散数值模拟研究[J]. 海洋 科学, 2014, 38(11): 100-107. HUANG Juan. Numerical simulation of oil spill driftdiffusion in the Bohai sea[J]. Marine Science, 2014, 38(11): 100-107.

[6] 黄娟. 渤海溢油三维漂移数值模拟研究[J]. 海洋科学, 2015, 39(2): 110-117. HUANG Juan. Numerical simulation of oil spill in Bohai[J]. Marine Science, 2015, 39(2): 110-117