# History file 'sal-inv2.h' to fit the SEG/EAGE 3-D Salt Model by a model # with interfaces. # Input data for the salt interface are taken from the velocity field. # Optionally, the data from the GOCAD file 'salt.ts', splitted using # 'sal-salt.h' into three input files, may be used. # Dimension of array RAM in include file 'ram.inc' should be # MRAM=15000000 or greater for this history file. # ===================================================================== # Input files required chk.pl: "data/sal/" "sal-top.dat" chk.pl: "data/sal/" "sal-bot.dat" chk.pl: "data/sal/" "sal-m2.dat" chk.pl: "model/" "sob22.dat" chk.pl: "forms/" "inv.cal" chk.pl: "forms/" "addsob.cal" chk.pl: "forms/" "sqrt.cal" chk.pl: "data/sal/" "sal-icb.cal" # Output of 'sal-ts.h' chk.pl: "" "top.pts" chk.pl: "" "hrz_1.pts" chk.pl: "" "gp2.pts" chk.pl: "" "bottom.pts" # Output of 'sal-vel.h' chk.pl: "" "s-vc810r.out" chk.pl: "" "s-sc810r.out" # Initial model and its consistency check # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MODEL='sal-m2.dat' MODIN='sal-m2.dat' N1=20 N2=20 N3=10 LFREE=1 MODLOG='modchk.out' #modchk: # Form of the files with matrices # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ FORMM='unformatted' MATFORM='unformatted' # Inversion of interfaces # ~~~~~~~~~~~~~~~~~~~~~~~ ICLASS=1 # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsob.out' SOBW00=1 invsoft: MODSOB= SOBW00= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' MPAR=0 PTS='top.pts' INDFUN=1 ERRMUL=2.000 invpts: # 4 points M2IN='m2.out' # PTS='hrz_1.pts' INDFUN=2 ERRMUL=24.52 invpts: # 601 points PTS='gp2.pts' INDFUN=3 ERRMUL=42.19 invpts: # 1780 points PTS='bottom.pts' INDFUN=4 ERRMUL=2.000 invpts: # 4 points # Using the data from the velocity field PTS='sal-top.dat' INDFUN=5 ERRMUL=80.16 invpts: # 6425 points PTS='sal-bot.dat' INDFUN=6 ERRMUL=79.60 invpts: # 6336 points M2IN=' ' PTS=' ' INDFUN= ERRMUL= # 15150 total # Using the data from the GOCAD surfaces # PTS='salt-t1.pts' INDFUN=5 ERRMUL=28.16 invpts: # 793 points # PTS='salt-t2.pts' INDFUN=5 ERRMUL=14.04 invpts: # 197 points # PTS='salt-b1.pts' INDFUN=6 ERRMUL=24.02 invpts: # 577 points # M2IN=' ' PTS=' ' INDFUN= ERRMUL= # 3956 total # Matrix operations MATIN1='dm1.out' MATOUT='dm2.out' MATFUN='inv' matfun: MATIN1='dm2.out' MATIN2='gm1.out' MATOUT='dm2gm1.out' SYMMETRY= MATT1= MATT2=1 matmul: MATIN1='gm1.out' MATIN2='dm2gm1.out' MATOUT='sm1.out' SYMMETRY='sym' MATT1= MATT2= matmul: MATIN1='sm1.out' MATIN2='modsob.out' MATOUT='sm2.out' COEF1= COEF2=169000000 matlin: MATIN1='sm2.out' MATOUT='sm3.out' matinv: MATIN1='dm2.out' MATIN2='gm2.out' MATOUT='gm3.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm1.out' MATIN2='gm3.out' MATOUT='gm4.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='sm3.out' MATIN2=' ' MATOUT='sm3g.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='sm3g.out' MATIN2='gm4.out' MATOUT='gm5.out' SYMMETRY=' ' MATT1= MATT2= matmul: # Updating the model M1='m1.out' MODNEW='gm5.out' MODOUT='sal-m2.tmp' modmod: # Indices of complex blocks relevant to the velocities # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ N1=68 N2=68 N3=20 D1=200 D2=200 D3=200 O1=0 O2=0 O3=40 N1=46 N2=46 N3=14 D1=300 D2=300 D3=300 O1=0 O2=0 O3=100 N1=34 N2=34 N3=20 D1=400 D2=400 D3=200 O1=200 O2=200 O3=40 N1=23 N2=23 N3=21 D1=600 D2=600 D3=200 O1=200 O2=200 O3=100 N1=17 N2=17 N3=20 D1=800 D2=800 D3=200 O1=400 O2=400 O3=40 N1=11 N2=11 N3=20 D1=1340 D2=1340 D3=200 O1=40 O2=40 O3=40 N1=17 N2=17 N3=39 D1=800 D2=800 D3=100 O1=400 O2=400 O3=40 ISRF=3 VEL='sal-s3.out' grid: ISRF=4 VEL='sal-s4.out' grid: ISRF= CAL='sal-icb.cal' GRD2='sal-s3.out' GRD3='sal-s4.out' GRD4='sal-icb.out' SALTV1=1000 SALTV2=1501 SALTV5=4449 SALTD3=152.4 SALTD4=1 GRD1='s-va200.out' GRD1='s-vc200.out' GRD1='s-va200r.out' GRD1='s-vc200r.out' GRD1='s-va300.out' GRD1='s-va420.out' GRD1='s-va620.out' GRD1='s-va820.out' GRD1='s-va134.out' GRD1='s-vc810r.out' grdcal: # Inversion of velocities # ~~~~~~~~~~~~~~~~~~~~~~~ MODEL='sal-m2.tmp' MODIN='sal-m2.tmp' ICLASS=2 NEGPAR=1 # Calculating matrices for inversion M1='m1.out' M2='m2.out' SOBOLEV='sob22.dat' MODSOB='modsob.out' SOBW01=1 invsoft: MODSOB= SOBW01= GM1='gm1.out' GM2='gm2.out' GM3=' ' DM1='dm1.out' GRDICB='sal-icb.out' MPAR=1 POWERM=-1 GRD='s-sa200.out' ERRMUL=304.11 #92480 points GRD='s-sc200.out' ERRMUL=304.11 #92480 points GRD='s-sa200r.out' ERRMUL=304.11 #92480 points GRD='s-sc200r.out' ERRMUL=304.11 #92480 points GRD='s-sa300.out' ERRMUL=172.12 #29624 points GRD='s-sa420.out' ERRMUL=152.05 #23120 points GRD='s-sa620.out' ERRMUL=105.40 #11109 points GRD='s-sa820.out' ERRMUL=76.026 #5780 points GRD='s-sa134.out' ERRMUL=50.408 #2541 points GRD='s-sc810r.out' ERRMUL=106.16 #11271 points invpts: GRD= GRDICB= MPAR= POWERM= ERRMUL= # Matrix operations MATIN1='dm1.out' MATOUT='dm2.out' MATFUN='inv' matfun: MATIN1='dm2.out' MATIN2='gm1.out' MATOUT='dm2gm1.out' SYMMETRY= MATT1= MATT2=1 matmul: MATIN1='gm1.out' MATIN2='dm2gm1.out' MATOUT='sm1.out' SYMMETRY='sym' MATT1= MATT2= matmul: MATIN1='sm1.out' MATIN2='modsob.out' MATOUT='sm2.out' COEF1= COEF2=25000000 matlin: MATIN1='sm2.out' MATOUT='sm3.out' matinv: MATIN1='dm2.out' MATIN2='gm2.out' MATOUT='gm3.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm1.out' MATIN2='gm3.out' MATOUT='gm4.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='sm3.out' MATIN2='gm4.out' MATOUT='gm5.out' SYMMETRY=' ' MATT1= MATT2= matmul: # Updating the model M1='m1.out' MODNEW='gm5.out' MODOUT='sal-m2.out' modmod: # Sobolev norm of the slowness in the updated model # ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ MATIN1='modsob.out' MATIN2='gm5.out' MATOUT='gm6.out' SYMMETRY=' ' MATT1= MATT2= matmul: MATIN1='gm5.out' MATIN2='gm6.out' MATOUT='gm7.out' SYMMETRY='diag' MATT1=1 MATT2=0 matmul: MATIN1='gm7.out' MATOUT='sal-sns.out' SPARSE=-1 MATFUN='sqrt' MATFORM='formatted' matfun: # List of output files # ~~~~~~~~~~~~~~~~~~~~ # 'sal-m2.tmp' ... Inverted model after inversion of interfaces. # 'sal-m2.out' ... Inverted model after inversions of interfaces and # of velocities. # 'sal-sns.ou@'... Sobolev norm of the slowness in the model.