clear clc close all size = 200; % %%% 555 T=[logspace(-8,log10(1e-03),size)]; % Dimensionless time T=[logspace(-8,log10(5e-4),2000)]'; q = 100/(60*60); lambda = 608.29875*1000; % Dimensionless filtration coefficient lambda=10; beta = 110875; % Formation damage coeffcient alpha = 0.0047; %U = linspace(0,0.025,size); U = linspace(0,0.1,size); [J, c_cumulative,C_cumulative] = SY2Ver2_C_J_file(q,lambda,beta,alpha,T); figure(1) plot(T,J,'r','LineWidth',3); xlabel('T','FontSize',22,'fontname','Times New Roman'); ylabel('J','fontsize',22,'fontname','Times New Roman'); % figure(2) % plot(T,c_cumulative,'b','LineWidth',3); % xlabel('T','FontSize',22,'fontname','Times New Roman'); % ylabel('Cacc','fontsize',22,'fontname','Times New Roman'); % figure(1) % plot(T_D,J_D,'r','LineWidth',3); % hold on % plot(T_D,J_B,'b','LineWidth',3); % hold on % plot(T_D,J_Overall,'k','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('J','fontsize',20,'fontname','Trebuchet'); % legend('Detrital particles','Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast'); % ylim([1 25]); % % figure(2) % plot(T_D,J_D,'r','LineWidth',3); % hold on % plot(T_D,J_B,'b','LineWidth',3); % hold on % plot(T_D,J_Overall,'k','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('J','fontsize',20,'fontname','Trebuchet'); % legend('1/3 Detrital particles','2/3 Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast'); % ylim([1 25]); % % % figure(3) % plot(T_D,c_D,'r','LineWidth',3); % hold on % plot(T_D,c_B,'b','LineWidth',3); % hold on % plot(T_D,c_Overall,'k','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('C_{acc}','fontsize',20,'fontname','Trebuchet'); % legend('Detrital particles','Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast'); % % figure(4) % plot(T_D,c_D,'r','LineWidth',3); % hold on % plot(T_D,c_B,'b','LineWidth',3); % hold on % plot(T_D,c_Overall,'k','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('C_{acc}','fontsize',20,'fontname','Trebuchet'); % legend('1/3 Detrital particles','2/3 Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast'); % % % % %% % figure(5) % plot(T_S1_DLVO,J_S1_DLVO,'r','LineWidth',3); % hold on % plot(T_S1_DLVO,J_S1_Breakage,'b','LineWidth',3); % hold on % plot(T_S1_DLVO,J_S1_Overall,'k--','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('J','fontsize',20,'fontname','Trebuchet'); % legend('Detrital particles','Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast'); % % % figure(6) % plot(T_S1_DLVO,c_S1_DLVO,'r','LineWidth',3); % hold on % plot(T_S1_DLVO,c_S1_Breakage,'b','LineWidth',3); % hold on % plot(T_S1_DLVO,c_S1_Overall,'k--','LineWidth',3); % set(gca,'fontname','Trebuchet','FontSize',12); % xlabel('T, PV','FontSize',20,'fontname','Trebuchet'); % ylabel('C_{acc}','fontsize',20,'fontname','Trebuchet'); % legend('Detrital particles','Authigenic particles','Total particles','fontsize',12,'fontname','Trebuchet','Location','Southeast');