我正在模拟钻井液在地层中的渗透半径,可是得到的结果却不是以圆弧形态均匀扩散的(不对称):
但是我看大神做出了这个结果:
这是我的程序:
new
config fluid
set fluid off
;=================================
;建立网络模型
generate zone cshell size 10 10 10 10 ratio 1 1 1 1.2 p0 0 0 0 p1 3 0 0 p2 0 6 0 p3 0 0 3 dimension 0.15,0.15,0.15
;===================================
;显示
plot set rotation 20 0 30
plot set center auto
plot set magnification 1.5
plot add surface blue
plot add axes red
plot show
;===================================
group soil
group dam ran x 0 3 y 0 2 z 0 3 ;设置不透水面
group dam ran x 0 3 y 4 6 z 0 3 ;设置不透水面
hist fltime
hist id=2 gp zd 0.2,3,0.2
hist id=3 gp pp 0.3,3,0.3
hist id=4 gp pp 0.4,3,0.4
hist id=5 gp pp 0.5,3,0.5
hist id=6 gp pp 0.6,3,0.6
hist id=7 gp pp 0.7,3,0.7
hist id=8 gp pp 0.8,3,0.8
hist id=9 gp pp 0.9,3,0.9
hist id=10 gp pp 1,3,1
hist id=11 gp pp 1.2,3,1.2
hist id=12 gp pp 1.4,3,1.4
hist id=13 gp pp 1.6,3,1.6
hist id=14 gp pp 1.8,3,1.8
hist id=15 gp pp 2,3,2
m e
prop bu 4.17e7 sh 10.19e7 ;设置体积模量和剪切模量
ini pp 0 grad 0 0 -10e3 ran y 0 6 ;生成孔隙压力
ini dens 1500 ;定义地层密度
model fl_iso ;使用各项同性渗流模型
prop por 0.02 perm 1e-10;设置孔隙率
ini fden 1025 ften -1e10 ;设置流体密度
ini sat 0.0 ran y 2 4 ;设置饱和度
model fl_null ran gro dam
;ini pp 0 ran gro dam ;设置不透水层为不透水材料模型
fix z ran z -0.01 0.01 ;定义边界条件
fix z ran z 2.99 3.01
fix x ran x -0.01 0.01
fix x ran x 2.99 3.01
fix y ran y -0.01 0.01
fix y ran y 5.99 6.01
set grav 0,10,0 ;设置重力加速度
solve
save elastic1.sav
rest elastic1.sav
set fluid on mech off
ini fmod 2e3 ften 0.0 ran gro soil
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
app pp 30e3 grad 0 0 -10e3 ran y 2 4 x 0
app pp 0 ran y 2 4 x 3
plot con pp out on
step 1000
solve
我估计应该是没有定义钻孔圆弧边界、没有在圆弧面上加载压力,结果导致了钻井液扩散不均匀、不对称!所以请教各路大神支支招!!万分感谢!