13 lines
340 B
Matlab
13 lines
340 B
Matlab
function well1 = handle_well1(well1,r)
|
|
% 需要找到射孔段具体所在的基质网格编号
|
|
n1=size(well1,1);
|
|
for i=1:n1
|
|
perf=well1{i,3};
|
|
nperf=well1{i,2};
|
|
perf(:,3)=(r.nz+1)*ones(nperf,1)-perf(:,3);
|
|
A=repmat([0 1 1],nperf,1);
|
|
B=(perf-A)*[1;r.nx;r.nx*r.ny];
|
|
B=B'; %转换为行向量
|
|
well1{i,3}=B;
|
|
end
|
|
end |