diff --git a/算例1-组分流-不规则非均质计算域/main1.m b/算例1-组分流-不规则非均质计算域/main1.m index e6cd6e4..402a502 100644 --- a/算例1-组分流-不规则非均质计算域/main1.m +++ b/算例1-组分流-不规则非均质计算域/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,1)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx, dy,dz,nx,ny,nz,NTG); diff --git a/算例2-组分流-无效网格/main1.m b/算例2-组分流-无效网格/main1.m index 1759026..ece2201 100644 --- a/算例2-组分流-无效网格/main1.m +++ b/算例2-组分流-无效网格/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,20) ];ny=size(dy,2); dz=[10*ones(1,3)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx,dy,dz,nx,ny,nz,NTG); diff --git a/算例3-单重双重离散裂缝混合模型/main1.m b/算例3-单重双重离散裂缝混合模型/main1.m index e6d8bc2..57a3601 100644 --- a/算例3-单重双重离散裂缝混合模型/main1.m +++ b/算例3-单重双重离散裂缝混合模型/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,1)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx,dy,dz,nx,ny,nz,NTG); diff --git a/算例4-气水两相流-压力导数曲线/main1.m b/算例4-气水两相流-压力导数曲线/main1.m index ca69cef..4c2b4cb 100644 --- a/算例4-气水两相流-压力导数曲线/main1.m +++ b/算例4-气水两相流-压力导数曲线/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,1)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx, dy,dz,nx,ny,nz,NTG); diff --git a/算例5-油水两相流-注焖采/main1.m b/算例5-油水两相流-注焖采/main1.m index 472478f..3e7d182 100644 --- a/算例5-油水两相流-注焖采/main1.m +++ b/算例5-油水两相流-注焖采/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,1)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx,dy,dz,nx,ny,nz,NTG); diff --git a/算例6-油水两相流-直井注压裂水平井采/main1.m b/算例6-油水两相流-直井注压裂水平井采/main1.m index 99df1fa..f4b7afd 100644 --- a/算例6-油水两相流-直井注压裂水平井采/main1.m +++ b/算例6-油水两相流-直井注压裂水平井采/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,2)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx, dy,dz,nx,ny,nz,NTG); diff --git a/算例7-油水两相流-不同裂缝制度不同/main1.m b/算例7-油水两相流-不同裂缝制度不同/main1.m index e13eafe..db459c3 100644 --- a/算例7-油水两相流-不同裂缝制度不同/main1.m +++ b/算例7-油水两相流-不同裂缝制度不同/main1.m @@ -3,8 +3,8 @@ function [ r, Times, OutputRs, Wellpara, trun ] = main1() % 1-- classical EDFM modelflag = 1;% 此处固定为1即可 %% ==================================PART 2: 基质网格定义===================================== -dx=[10*ones(1,100) ];nx=size(dx,2); -dy=[10*ones(1,50) ];ny=size(dy,2); +dx=[50*ones(1,20) ];nx=size(dx,2); +dy=[50*ones(1,10) ];ny=size(dy,2); dz=[10*ones(1,1)];nz=size(dz,2); NTG= 1*ones(nx*ny*nz,1); r = GridProp_pre(dx,dy,dz,nx,ny,nz,NTG);