Scribd
Upload
27 views3 pages

%core #1 %normalised SW %normalised Kro %normalised KRW %plotting The Data For This Core

The document contains code for plotting core flood data from 3 rock cores. It normalizes and fits saturation-relative permeability data from each core using Fourier series. The normalized saturation and relative permeability curves are plotted for each core individually, with different markers to distinguish them. It then calculates and plots the average relative permeability curves across the 3 cores. Finally, it displays the normalized saturation values alongside the average relative permeability values.

Uploaded by

MohamedElAbbasy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as RTF, PDF, TXT or read online on Scribd
27 views3 pages

%core #1 %normalised SW %normalised Kro %normalised KRW %plotting The Data For This Core

The document contains code for plotting core flood data from 3 rock cores. It normalizes and fits saturation-relative permeability data from each core using Fourier series. The normalized saturation and relative permeability curves are plotted for each core individually, with different markers to distinguish them. It then calculates and plots the average relative permeability curves across the 3 cores. Finally, it displays the normalized saturation values alongside the average relative permeability values.

Uploaded by

MohamedElAbbasy
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as RTF, PDF, TXT or read online on Scribd
You are on page 1/ 3

%Core #1

%Normalised Sw
sw1=[0 0.125 0.375 0.625 0.875 1]';
%Normalised Kro
kro1=[1 0.887 0.655 0.414 0.154 0]';
%Normalised Krw
krw1=[0 0.045 0.23 0.495 0.818 1]';
%plotting the data for this core
c1=fit(sw1,kro1,'fourier2');
plot(c1,sw1,kro1)
hold on
plot(sw1,kro1,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','k',...
'MarkerSize',5)

axis([0 1 0 1])

hold on
c2=fit(sw1,krw1,'fourier2');
plot(c2,sw1,krw1)
plot(sw1,krw1,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','k',...
'MarkerSize',5)

%-----------------------------------------%
%Core 2

%Normalised Sw
sw2=[0 0.238 0.476 0.714 0.833 1]';
%Normalised Kro
kro2=[1 0.741 0.491 0.252 0.139 0]';
%Normalised Krw
krw2=[0 0.154 0.382 0.646 0.788 1]';
%plotting data
c3=fit(sw2,kro2,'fourier2');
plot(c3,sw2,kro2)
plot(sw2,kro2,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','m',...
'MarkerSize',5)

c4=fit(sw2,krw2,'fourier2');
plot(c4,sw2,krw2)

plot(sw2,krw2,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','m',...
'MarkerSize',5)

legend off

%---------------------------------------------%
%Core 3

%Normaised Sw
sw3=[0 0.111 0.222 0.444 0.667 0.889 1]';
%Normalised Kro
kro3=[1 0.872 0.839 0.663 0.463 0.215 0]';
%Normalised Krw
krw3=[0 0.023 0.077 0.251 0.503 0.817 1]';

%plotting data
c5=fit(sw3,kro3,'fourier2');
plot(c5,sw3,kro3)
plot(sw3,kro3,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','g',...
'MarkerSize',5)
%
%plot(sw3,kro3,'o')
c6=fit(sw3,krw3,'fourier2');
plot(c6,sw3,krw3);
%plot(sw3,krw3,'o','Color','b')
plot(sw3,krw3,'o',...
'MarkerEdgeColor','k',...
'MarkerFaceColor','g',...
'MarkerSize',5)
legend off

text(0.5,0.95,' Core #1','FontSize',10,'Color','k')


text(0.5,0.90,' Core #2','FontSize',10,'Color','m')
text(0.5,0.85,' Core #3','FontSize',10,'Color','g')
figure
sw=linspace(0,1);
plot(sw,(c1(sw)+c3(sw)+c5(sw))/3) %for oil
hold all
plot(sw,(c2(sw)+c4(sw)+c6(sw))/3) %for water

axis([0 1 0 1])

% Displaying the normalised sw and average Kro , Krw


sw=0:0.1:1;
disp('Normalised Sw |avg.Kro | avg.Krw ')
p=[sw',(c1(sw)+c3(sw)+c5(sw))/3,(c2(sw)+c4(sw)+c6(sw))/3];
disp(p)

You might also like