Q =0.0025375;%('输入体积流率(m^3/s)');
L =609;%('水平井筒长度(m)');
qw=Q/L;%延水平井筒方向单位长度上的注入流量(m^3/s*m);
Cw = 4.2*10^6;%('注入流体定压比热容 (J/m^3*K)');
Cr =2.35*10^6;%('储集层岩石定压比热容 (J/m^3*K)');
pw =1.08;%('注入流体密度(g/m^3)');
pr =2.15;%('储集层岩石密度(g/m^3)');
O =0.2;%('孔隙度(小数)');
KT =2.7778;%('储集层平均导热系数(J/m*K*s)');
T1=100;
Te=100;
re =3;%('油藏径向半径(m)');
rw =0.06858;%('水平井井筒半径(m)');
Aw=pi*rw^2;%水平井筒横截面积(m^2)
x =10;%('水平井筒水平网格数');
Xi=L/x;%(水平井筒水平步长(m);
z=2;%(水平井筒加相邻地层在j方向的网格数)
Zi=1;
ti =5;%('时间增量min');
NR=10;%('径向网格数');
a=(re/rw)^(1/(NR-1));%('径向网格转化为线性网格系数')
A0=(Cw*Aw*pw*Xi*Zi+Cw*pw*qw*Xi*Zi*ti-Cw*pw*qw*Zi*ti+(2*pi*KT*a*Xi*ti)/log(a))/(Cw*Aw*pw*Xi*Zi);%差分方程T(i,j)系数
B0=(Cw*pw*qw*Zi*ti)/(Cw*Aw*pw*Xi*Zi); %差分方程T(i+1,j)系数
C0=((2*pi*KT*a*Xi*ti)/log(a))/(Cw*Aw*pw*Xi*Zi);%差分方程T(i,j+1)系数
%稀疏矩阵的创建
V0=ones([(x+1)*2,1])*A0;%稀疏矩阵主对角线上的系数,+1是在趾部地层多加一个网格;
V1=ones([(x+1)*2-1,1])*B0;%稀疏矩阵泛1对角线系数,+1是在 趾部地层多加一个网格;
V2=ones([(x+1)*2-10,1])*C0;%稀疏矩阵泛10对角线系数,+1是在趾部地层多加一个网格;
V=diag(V0)+diag(V1,1)+diag(V2,10);%线性方程组稀疏矩阵
jj=1;
T=rand((x+1)*2,1);
D=rand((x+1)*2,1);%D:线性方程组的右端项
D=D((x+1)*2,1);
for ii=1:10;
D(ii,jj)=T1;
end
for ii=11;
D(ii,jj)=Te;
end
for ii=11:21;
D(ii,jj)=(T1+Te)/2;
end
for ii=22;
D(ii,jj)=Te;
end
for i=1
T(i)=V\D;
for i=2:42
T(i)=V\T(i-1);
end
end
运行程序说是
In an assignment A(I) = B, the number of elements in B and I must be the same.
Error in haole (line 48)
T(i)=V\D;
我这里的V,和D都是列向量,维数也是相同的,不明白为什么会出错?