; tit FLAC Modelling tt11 @ 1800m ; same as tt8 but more steps to try and get more tensile failure and props same as nerys' config extra 5 gw grid 100 180 gen 0.00e0,0.00e0 0.00e0,1.00e3 1.00e3,1.00e3 1.00e3,0.00e0 i=1,101 j=1,181 mark i=1 mark j=1 mark i = 101 mark j = 181 ; *************************************** model mohr water bulk 2e9 dens 1000 prop poro=0.3 ; to be elaborated below ;; 1 mech step = 1 fluid step; modify if necessary!!! set nmech 1 ngw 1 flow on *************************************** ; tab 1 0 , 904.11 1000 , 874.82 ; ; tab 2 518.63 , 886.53 0 , 0 382.17 , 0 tab 2 924.11 , 876.77 ; ; tab 3 647.29 , 884.58 119 , 0 255.46 , 0 tab 3 752.56 , 798.68 803.24 , 880.68 ; ; gen tab 1 gen tab 2 gen tab 3 ; ; ; Properties for each sub-regions: ; Shale pro dens= 2500 bulk=2.81e10 shear=1.69e10 reg= 41 , 164 pro coh=5.0e8 tens=5.0e8 reg= 41 , 164 pro fric=30 dil= 3 reg= 41 , 164 pro perm=1.00e-17 porosity=0.2 reg= 41 , 164 ; conglomerate pro dens= 2720 bulk=4.09e10 shear=2.45e10 reg= 15 , 108 pro coh=5.0e8 tens=5.0e8 reg= 15 , 108 pro fric=30 dil= 3 reg= 15 , 108 pro perm=1.00e-16 porosity=0.2 reg= 15 , 108 ; conglomerate pro dens= 2720 bulk=4.09e10 shear=2.45e10 reg= 84 , 64 pro coh=5.0e8 tens=5.0e8 reg= 84 , 64 pro fric=30 dil= 3 reg= 84 , 64 pro perm=1.00e-16 porosity=0.2 reg= 84 , 64 ; Intrusive Suite pro dens= 2700 bulk=4.67e10 shear=2.80e10 reg= 53 , 98 pro coh=2.0e7 tens=2e6 reg= 53 , 98 pro fric=30 dil= 2 reg= 53 , 98 pro perm=1.00e-13 porosity=0.3 reg= 53 , 98 ; Alteration Envelope pro dens= 2720 bulk=4.09e10 shear=2.45e10 reg= 57 , 71 pro coh=2.0e7 tens=2e6 reg= 57 , 71 pro fric=30 dil= 2 reg= 57 , 71 pro perm=1.00e-13 porosity=0.3 reg= 57 , 71 ; Alteration Envelope pro dens= 2720 bulk=4.09e10 shear=2.45e10 reg= 40 , 107 pro coh=2.0e7 tens=2e6 reg= 40 , 107 pro fric=30 dil= 2 reg= 40 , 107 pro perm=1.00e-13 porosity=0.3 reg= 40 , 107 apply syy = -4.41e7 j= 181 ; ***** def inistr_phg loop i (1,igp-1) loop j (1,jgp-1) jno = jgp - 1 If j = jno Then hh1 = (abs(y(i,j)-y(i,j+1))+abs(y(i+1,j)-y(i+1,j+1)))*0.25 sss1 = density(i, j) * 9.8 * hh1 sss2 = sss1 Else jno1 = j jno2 = j + 1 hh1 = (abs(y(i,j)-y(i,j+1))+abs(y(i+1,j)-y(i+1,j+1)))*0.25 sss1 = density(i, j) * 9.8 * hh1 sss2 = sss1 loop jj (jno2,jgp-1) aaa = abs(y(i, jj) - y(i, jj + 1)) bbb = abs(y(i + 1, jj) - y(i + 1, jj + 1)) ccc = 0.5 * (aaa + bbb) sss3 = density(i, jj) * 9.8 * ccc sss2 = sss2 + sss3 end_loop end_if syy(i, j) = -sss2-4.41e7 sxx(i, j) = -sss2-4.41e7 szz(i, j) = -sss2-4.41e7 end_loop end_loop end inistr_phg * ; * def ini_pp loop i (1,igp) loop j (1,jgp) if j = jgp then sss2 = 4.41e7 else jno1 = j sss2 = 4.41e7 loop jj (jno1,jgp-1) aaa = abs(y(i, jj) - y(i, jj + 1)) if x(i, jj + 1) > x(i, jj) then ddden = density(i - 1, jj) end_if if x(i, jj + 1) < x(i, jj) then ddden = density(i, jj) end_if if x(i, jj + 1) = x(i, jj) then if i = 1 then ii1 = 1 ii2 = 1 end_if if i > 1 then if i < igp then ii1 = i ii2 = i - 1 end_if end_if if i = igp then ii1 = igp - 1 ii2 = igp - 1 end_if ddden = 0.5 * (density(ii1, jj) + density(ii2, jj)) end_if sss3 = ddden * 9.8 * aaa sss2 = sss2 + sss3 end_loop end_if gpp(i, j) = sss2 end_loop end_loop end ini_pp * ; fix pp j=181 ini sat 1 fix sat ; set gravity=9.8 ; ;return set large ini xvel 4e-4 var= -8e-4 ,0 i=1, 101 j=1, 181 ;ini yvel 0.002 var=-0.004 ,0 i=1, 101 j=1, 181 fix x i=1 fix x i=101 fix y j=1 ;fix y j=181 * fix i=1 101 fix j=1 181 * ;apply dis= 1e-7 i= 13 24 j= 1 apply dis= 1e-7 i= 2 38 j= 1 ;set nmech 1 ngw 1 flow on set nmech 1 ngw 10 flow on save tt11.sav; pre compression & pre influx ; step 500 save tt11a.sav; 1 % deformation step 500 save tt11b.sav; 2 % deformation step 500 save tt11c.sav; 3 % deformation step 500 save tt11d.sav; 4 % deformation step 500 save tt11e.sav; 5 % deformation step 500 save tt11f.sav; 6 % deformation step 500 save tt11g.sav; 7 % deformation step 500 save tt11h.sav; 8 % deformation step 500 save tt11i.sav; 9 % deformation step 500 save tt11j.sav 10 % deformation return stop