Scribd
Upload
60 views

Coding Matlab Bagus Ansory

This document contains Matlab code for calculating properties of brine and water based on temperature, pressure, and salinity. It defines constants and variables, calculates water velocity and density as a function of temperature and pressure using a polynomial equation, then calculates brine velocity and density by incorporating the effects of salinity into similar polynomial equations. The brine bulk modulus is then calculated from the brine density and velocity.

Uploaded by

Farsa Randitama
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
60 views

Coding Matlab Bagus Ansory

This document contains Matlab code for calculating properties of brine and water based on temperature, pressure, and salinity. It defines constants and variables, calculates water velocity and density as a function of temperature and pressure using a polynomial equation, then calculates brine velocity and density by incorporating the effects of salinity into similar polynomial equations. The brine bulk modulus is then calculated from the brine density and velocity.

Uploaded by

Farsa Randitama
Copyright
© © All Rights Reserved
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

Coding Matlab Bagus Ansory, Susi Susanti, Farsa Randitama

Brine Prop
rho_o = 42;
% Oil gravity (deg API)
GOR = 160.0;
% GOR (L/L)
rho_g = 0.9;
% Gas gravity (API)
T = 150.00;
% Temperature (0 C)
Ppsi = 3200.00; % Pressure (psi)
S = 3800;
% water salinity (ppm)
phi = 0.20;
% porosity (in fraction)
vsh = 0.20;
% Vsh (volume shale in fraction)
isw = 0.40; % initial water saturation (SW)
tsw = 1.00; % target water saturation (in fraction)
ifluid = 1; % initial hydrocarbon is 1(oil), 2(gas)
fluid = 1; % Desired fluid is 1(brine), 2(oil) 3(gas)
vp = 11000.0; % ft/s - from log (initial value)
vs = 6500.0;
% ft/s - from log (initial value)
rho = 2.2;
% gm/c - from log (initial value)
%
% Fixed parameters (e.g., Mavko et al., 1998)
%
k_clay = 20.9;
% Bulk mod (GPa)
k_qtz = 36.6;
rho_clay = 2.58; % gm/cc
rho_qtz = 2.65;
%
% some applied properties
%
div_mill = 1/1000000; % factor used to divide by million
fs2kms = 0.000305; % factor for ft/s to km/s conversion
kms2fs = 3280.84; % factor for km/s to ft/s conversion
v_clay = vsh*0.70;
% Assumption: V_clay = 70% of VSH
v_qtz= 1-v_clay;
% quartz fraction in mineral
ish = 1-isw;
% initial hydrocarbon saturation
tsh = 1-tsw;
% final hydrocarbon saturation
rho_o = 141.5/(rho_o+131.5); % oil gravity in gm/cc (from API)
P = Ppsi*6.894757*0.001; % Press in MPa (from Psi)
S = S*div_mill;
% salinity as weight fraction
vp = vp*fs2kms;
% ft/s to km/s
vs = vs*fs2kms;
% ft/s to km/s
%
% Step 2: water/brine properties (Equations 10 and 11)
%
w(1,1) = 1402.85; w(1,3) = 3.437*10^(-3); % Table 2
w(2,1) = 4.871; w(2,3) = 1.739*10^(-4);
w(3,1) = -0.04783; w(3,3) = -2.135*10^(-6);
w(4,1) = 1.487*10^(-4); w(4,3) = -1.455*10^(-8);
w(5,1) = -2.197*10^(-7); w(5,3) = 5.230*10^(-11);
w(1,2) = 1.524; w(1,4) = -1.197*10^(-5);
w(2,2) = -0.0111; w(2,4) = -1.628*10^(-6);
w(3,2) = 2.747*10^(-4); w(3,4) = 1.237*10^(-8);
w(4,2) = -6.503*10^(-7); w(4,4) = 1.327*10^(-10);
w(5,2) = 7.987*10^(-10); w(5,4) = -4.614*10^(-13);
sum = 0;
for i=1:5
for j=1:4
sum = sum+w(i,j)*T^(i-1)*P^(j-1);
end
end
v_water = sum
v1 = 1170-9.6*T+0.055*T*T-8.5*10^(-5)*T*T*T+2.6*P-0.0029*T*P-0.0476*P*P
v_brine= v_water+S*v1+S^1.5*(780-10*P+0.16*P*P)-1820*S*S
% m/s
r1= 489*P-2*T*P+0.016*T*T*P-1.3*10^(-5)*T*T*T*P-0.333*P*P-0.002*T*P*P
rho_water=1+10^(-6)*(-80*T-3.3*T*T+0.00175*T*T*T+r1)
r2= 300*P-2400*P*S+T*(80+3*T-3300*S-13*P+47*P*S)
rho_brine= rho_water+0.668*S+0.44*S*S+10^(-6)*S*r2
%gm/cc (held const)
k_brine= rho_brine*v_brine*v_brine*div_mill
%GPa(held const)
%

You might also like