config fluid
plot show
;set fluid off
gen zon brick p0 0 0 0 p1 300 0 0 p2 0 200 0 p3 0 0 18 size 30 20 18 group 石灰岩
gen zon brick p0 0 0 18 p1 300 0 18 p2 0 200 18 p3 0 0 29 size 30 20 11 group 砂岩
gen zon brick p0 0 0 29 p1 300 0 29 p2 0 200 29 p3 0 0 66 size 60 40 37 group 泥岩
gen zon brick p0 0 0 66 p1 300 0 66 p2 0 200 66 p3 0 0 88 size 60 40 22 group 砂岩
gen zon brick p0 0 0 88 p1 300 0 88 p2 0 200 88 p3 0 0 92 size 60 40 4 group 泥岩
gen zon brick p0 0 0 92 p1 300 0 92 p2 0 200 92 p3 0 0 98 size 60 40 6 group 煤层
gen zon brick p0 0 0 98 p1 300 0 98 p2 0 200 98 p3 0 0 105 size 60 40 7 group 泥岩
gen zon brick p0 0 0 105 p1 300 0 105 p2 0 200 105 p3 0 0 141 size 60 40 36 group 砂岩
gen zon brick p0 0 0 141 p1 300 0 141 p2 0 200 141 p3 0 0 159 size 30 20 18 group 泥岩
gen zon brick p0 0 0 159 p1 300 0 159 p2 0 200 159 p3 0 0 180 size 30 20 21 group 砂岩
group 陷落柱 ran x 180 220 y 80 120 z 0 180
model fl_iso
prop perm 1.5e-12 poro 0.2 ran group 砂岩
prop perm 4.5e-12 poro 0.15 ran group 泥岩
prop perm 0.1e-10 poro 0.26 ran group 煤层
prop perm 0.36e-10 poro 0.25 ran group 石灰岩
ini fden 1000 fmod 2e9 ften -1e10
ini pp 3e6 grad 0 0 -1.7e4 ；加孔压
fix pp ran z -0.1 0.1
m e
prop bu 3e10 sh 1.8e10 ran group 砂岩
prop bu 2.7e10 sh 1.5e10 ran group 泥岩
prop bu 7e9 sh 2e9 ran group 煤层
prop bu 3e10 sh 2e10 ran group 石灰岩
prop bu 1.8e9 sh 0.94e9 ran group 陷落柱
ini density 2700 ran group 砂岩
ini density 2300 ran group 泥岩
ini density 1400 ran group 煤层
ini density 2000 ran group 陷落柱
ini density 2600 ran group 石灰岩
fix x ran x -0.1 0.1
fix x ran x -299.9 299.9
fix y ran y -0.1 0.1
fix y ran y 199.9 200.1
fix z ran z 179.9.1 180.1
ini szz -8.1e6 ran z -0.1 180.1
ini sxx -12.05e6 ran z -0.1 180.1
ini syy -12.05e6 ran z -0.1 180.1
apply nsterss -4e6 ran z 179.9 180.1
set grav 0 0 -10
solve
save 原始.sav
