70 lines
2.7 KiB
Matlab
70 lines
2.7 KiB
Matlab
function r = grid_discretization_DP_model(modelflag, r, f, frac_information, fellip, nf)
|
|
% 数据读取
|
|
nx = r.nx; ny = r.ny; nz = r.nz; NTG = r.NTG; cell_mid_coords = r.cell_mid_coords;
|
|
% 裂缝层网格渗渗透率
|
|
kx = 5*1e-3 * ones(nx*ny*nz, 1); % 达西，D
|
|
% load('kx.mat');
|
|
ky = 5*1e-3 * ones(nx*ny*nz, 1);%
|
|
kz = 5*1e-3 * ones(nx*ny*nz, 1);%
|
|
pori = 0.01 * ones(nx*ny*nz, 1);%
|
|
prpor = 20; % 参考压力，MPa
|
|
cpor = 1.0e-5;% 基质压缩系数，MPa
|
|
NTG= 1 * ones(nx*ny*nz, 1);
|
|
% 基质层网格渗透率
|
|
kx_matrixLayer = 0.01*1e-3 * ones(nx*ny*nz, 1); % 达西，D
|
|
ky_matrixLayer = 0.01*1e-3 * ones(nx*ny*nz, 1);%
|
|
kz_matrixLayer = 0.01*1e-3 * ones(nx*ny*nz, 1);%
|
|
pori_matrixLayer = 0.1 * ones(nx*ny*nz, 1);%
|
|
% prpor = 20; % 参考压力，MPa
|
|
% cpor = 1.0e-5;% 基质压缩系数，MPa
|
|
% 形状因子
|
|
sigma = 0.08*ones(nx*ny*nz, 1);
|
|
for k = 1:nz
|
|
for j= 1:ny
|
|
for i= 1:floor(nx/2)
|
|
grid_numbering = (k-1)*(nx*ny)+(j-1)*nx+i;
|
|
sigma(grid_numbering,1) = 0.08;
|
|
end
|
|
end
|
|
end
|
|
% 有效网格
|
|
valid_grids = ones(nx*ny*nz, 1);
|
|
% SRV区域定义
|
|
% % % base_x=23:78; nx_SRV=length(base_x);
|
|
% % % base_y=15:36; ny_SRV=length(base_y);
|
|
% % % base_z=1:1;nz_SRV=length(base_z);
|
|
% % % SRV=[];
|
|
% % % for k=1:nz_SRV
|
|
% % % for j=1:ny_SRV
|
|
% % % for i=1:nx_SRV
|
|
% % % SRV=[SRV;(base_z(k)-1)*nx*ny+(base_y(j)-1)*nx+base_x(i)];
|
|
% % % end
|
|
% % % end
|
|
% % % end
|
|
% % % % SRV渗透率赋值
|
|
% % % kx(SRV)=0.001*1e-3 * ones(length(SRV), 1);
|
|
% % % ky(SRV)=0.001*1e-3 * ones(length(SRV), 1);
|
|
% % % kz(SRV)=0.001*1e-3 * ones(length(SRV), 1);
|
|
% 裂缝孔渗
|
|
Kf = 10000*1e-3*ones(1,nf);% 裂缝渗透率
|
|
Wf =1e-2*ones(1,nf);% 裂缝开度
|
|
Porf = 0.30*ones(1,nf);% 裂缝孔隙度
|
|
prporf = 20;% 裂缝系统参考压力
|
|
cporf = 1.0e-5;% 裂缝系统压缩系数
|
|
% 应力敏感系数
|
|
stress_factor_fracture = 0.000; % 1/MPa 0.001
|
|
stress_factor_matrix = 0.000; % 1/MPa
|
|
stress_factor_ref_pressure = 20; % 一般取为原始地层压力
|
|
%
|
|
Rpt = 2;cf = 1;ca = 1;
|
|
%% 结合基质、裂缝的几何及物性信息开展前处理
|
|
% r = GridProp_new(modelflag, nx, ny, nz, dx, dy, dz, kx, ky, kz, f, fellip, Kf, Wf, pori, prpor, cpor, Porf, prporf, cporf, cf, ca,co, NTG);
|
|
r = GridProp_DP(modelflag, r, kx, ky, kz, kx_matrixLayer, ky_matrixLayer, kz_matrixLayer, f, frac_information, fellip, Kf, Wf, pori, pori_matrixLayer, prpor, cpor, Porf, prporf, cporf, cf, ca,NTG,sigma,valid_grids);
|
|
rock_density = 2700; % 岩石密度，kg/m^3
|
|
r.sigma = sigma;
|
|
r.rock_density = rock_density;
|
|
r.valid_grids = valid_grids;
|
|
r.stress_factor_fracture = stress_factor_fracture;
|
|
r.stress_factor_matrix = stress_factor_matrix;
|
|
r.stress_factor_ref_pressure = stress_factor_ref_pressure;
|
|
end |