Code Matlab
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% Calcul des flux dans le bassin de lagunage %%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
close all
clear all
load Data_Kaedi.mat
% Données
A=1.7E-3; % surface de la lagune en km2
uz=1; % vitesse du vent en m/s
In=40; % Flux entrant en m3/j
Inf=1; % Flux d'infiltration mm/j
Volinit=In*60;
Soleil=[8,9,9,10,10,9,9,9,9,9,8,8];
for i=1:1:12
Tmax(i)=Data_Kaedi(i,1); % Température maximale en °C
Tmin(i)=Data_Kaedi(i,2); % Température minimale en °C
P(i)=Data_Kaedi(i,3); % Précipitation en mm
RH(i)=Data_Kaedi(i,4); % Humidité relative en %
end
% Calcul de la pression de vapeur saturante de l'air en z=O e0
for i=1:1:12
e0max(i)=exp(13.7)*exp(-5120/(Tmax(i)+273.3))*1013*100;
e0min(i)=exp(13.7)*exp(-5120/(Tmin(i)+273.3))*1013*100;
end
% Calcul de la pression de vapeur saturante de l'eau es
for i=1:1:12
esmax(i)=610.8*exp(17.17*Tmax(i)/(273.3+Tmax(i)));
esmin(i)=610.8*exp(17.17*Tmin(i)/(273.3+Tmin(i)));
end
% Calcul (EVP constitue l'évaporation)
for i=1:1:12
EVPmax(i)=1.69E-11*A^(-0.05)*(esmax(i)-RH(i)*e0max(i)/100)*86400000; % en mm/J (2592)
EVPmin(i)=1.69E-11*A^(-0.05)*(esmin(i)-RH(i)*e0min(i)/100)*86400000; % en mm/J (2592)
%multiplié par le nombre de secondes ensoleillées dans une journée !!
end
figure
plot((1:1:12),EVPmax(:),'r')
hold on
plot((1:1:12),EVPmin(:),'b')
title('Evaporation mensuelle à Kaedi')
xlabel('Mois')
ylabel('Evaporation en mm')
legend('Evaporation maximale','Evaporation minimale')
%plot((1:1:12),P(:),'c')
%Calcul de flux pour un jour d'un mois spécifique
for i=1:12
delta(i)=In+P(i)*A*1000/30-EVPmax(i)*A*1000*Soleil(i)/24-Inf*A*1000;
end
figure
plot((1:1:12),delta(:),'g')
title('Différence de flux')
xlabel('Mois')
ylabel('Delta (m3)')
Eff=sum(delta(:))/12 %Valeur du débit constant d'effluent journalier nécessaire pour retrouver un même volume en 1 an
%Calcul du volume final sur 1 an
Vol1=Volinit;
for i=1:12
Vol(i)=Vol1+30*(delta(i)-Eff);
Vol1=Vol(i);
end
figure
plot((1:1:12),Vol(:),'m')
title('Evoluvion du volume')
xlabel('Mois')
ylabel('Volume (m3)')
set(gca,'YLim',[0 3000],'Fontsize',12);
figure
plot((1:1:12),Vol(:)/(A*10^4),'b')
title('Evoluvion de la hauteur d eau dans la lagune')
xlabel('Mois')
ylabel('Hauteur d eau (cm)')
set(gca,'YLim',[0 150],'Fontsize',12);
%Scénario de remplissage commencant début mars
Vol2=31*delta(3)+30*delta(4); %mois de mars et avril
j=0;
while Vol2 < Volinit
Vol2=Vol2+delta(5);
j=j+1; %nombre de jours en mai
end
if j > 31
Vol2=31*delta(3)+30*delta(4)+31*delta(5); %mois de mars, avril et mai
k=0;
while Vol2 < Volinit
Vol2=Vol2+delta(6);
k=k+1; %nombre de jours en juin
end
tr=31+30+31+k;
else
tr=31+30+j;
end