b'\n \n \n
 
\n  
sfmpignfwi (4.0)
index
user/zhiguang/Mmpignfwi.c
\n Acoustic FWI using Gauss-Newton optimization \n

\n \n \n \n \n \n
 
\n Synopsis
       sfmpignfwi Fvel=Fv.rsf Fwavelet=Fw.rsf output=Fdat.rsf output=Finv.rsf Ferr=Ferr.rsf Fmod=Fmod.rsf Fgrad=Fgrad.rsf function=2 verb=n nb=100 coef=0.002 acqui_type=1 ns= ds= s0= sz=3 nr=acpar->nx dr=acpar->dx r0=acpar->x0 rz=3 fhi=0.5/acpar->dt flo=0. frectx=2 frectz=2 onlygrad=n wt1=acpar->t0 wt2=acpar->t0+(acpar->nt-1)*acpar->dt woff1=acpar->r0 woff2=acpar->r0+(acpar->nr-1)*acpar->dr gain=1 waterz=51 grectx=3 grectz=3 drectx=1 drectz=1 nrepeat=1 tangent=0 sigma1=-1 sigma2=-1 v1=0. v2=10. lniter=10 niter= conv_error= nls=20 factor=10 repeat=5 err_type=0
\n

\n \n \n \n \n \n
 
\n Parameters
       \n \n \n
\n  
file Ferr=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fgrad=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fmod=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
file Fvel=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
file Fwavelet=
\tauxiliary input file name
\n
\n \n\n \n \n
\n  
int acqui_type=1
\tif 1, fixed acquisition; if 2, marine acquisition; if 3, symmetric acquisition
\n
\n \n\n \n \n
\n  
float coef=0.002
\tabsorbing boundary coefficient
\n
\n \n\n \n \n
\n  
float conv_error=
\tfinal convergence error
\n
\n \n\n \n \n
\n  
float dr=acpar->dx
\treceiver interval
\n
\n \n\n \n \n
\n  
int drectx=1
\tsmoothing kernel radius in x
\n
\n \n\n \n \n
\n  
int drectz=1
\tsmoothing kernel radius in z
\n
\n \n\n \n \n
\n  
float ds=
\tshot interval
\n
\n \n\n \n \n
\n  
int err_type=0
\tif 0, true misfit function; if 1, both smoothing kernel and original L2 norm misfits
\n
\n \n\n \n \n
\n  
float factor=10
\tstep length increase factor
\n
\n \n\n \n \n
\n  
float fhi=0.5/acpar->dt
\thigh frequency in band, default is Nyquist
\n
\n \n\n \n \n
\n  
float flo=0.
\tlow frequency in band, default is zero
\n
\n \n\n \n \n
\n  
int frectx=2
\tsource smoothing in x
\n
\n \n\n \n \n
\n  
int frectz=2
\tsource smoothing in z
\n
\n \n\n \n \n
\n  
int function=2
\tif 1, forward modeling; if 2, FWI
\n
\n \n\n \n \n
\n  
float gain=1
\tvertical gain power of data residual
\n
\n \n\n \n \n
\n  
int grectx=3
\tgradient smoothing radius in x
\n
\n \n\n \n \n
\n  
int grectz=3
\tgradient smoothing radius in z
\n
\n \n\n \n \n
\n  
int lniter=10
\tCG iteration number
\n
\n \n\n \n \n
\n  
int nb=100
\tboundary width
\n
\n \n\n \n \n
\n  
int niter=
\titeration number
\n
\n \n\n \n \n
\n  
int nls=20
\tline search number
\n
\n \n\n \n \n
\n  
int nr=acpar->nx
\tnumber of receiver
\n
\n \n\n \n \n
\n  
int nrepeat=1
\tsmoothing kernel repeat number
\n
\n \n\n \n \n
\n  
int ns=
\tshot number
\n
\n \n\n \n \n
\n  
bool onlygrad=n [y/n]
\tonly calculate gradident or not
\n
\n \n\n \n \n
\n  
file output=
\tauxiliary output file name
\n
\n \n\n \n \n
\n  
float r0=acpar->x0
\treceiver origin
\n
\n \n\n \n \n
\n  
int repeat=5
\tafter how many iterations the step length goes back to 1
\n
\n \n\n \n \n
\n  
int rz=3
\treceiver depth
\n
\n \n\n \n \n
\n  
float s0=
\tshot origin
\n
\n \n\n \n \n
\n  
float sigma1=-1
\tsmoothing kernel radius moving step in z
\n
\n \n\n \n \n
\n  
float sigma2=-1
\tsmoothing kernel radius moving step in x
\n
\n \n\n \n \n
\n  
int sz=3
\tsource depth
\n
\n \n\n \n \n
\n  
int tangent=0
\tif 1, calculate prediction corrector
\n
\n \n\n \n \n
\n  
float v1=0.
\tlower limit of estimated velocity
\n
\n \n\n \n \n
\n  
float v2=10.
\tupper limit of estimated velocity
\n
\n \n\n \n \n
\n  
bool verb=n [y/n]
\tverbosity flag
\n
\n \n\n \n \n
\n  
int waterz=51
\twater layer depth
\n
\n \n\n \n \n
\n  
float woff1=acpar->r0
\twindow data residual: rmin
\n
\n \n\n \n \n
\n  
float woff2=acpar->r0+(acpar->nr-1)*acpar->dr
\twindow data residual: rmax
\n
\n \n\n \n \n
\n  
float wt1=acpar->t0
\twindow data residual: tmin
\n
\n \n\n \n \n
\n  
float wt2=acpar->t0+(acpar->nt-1)*acpar->dt
\twindow data residual: tmax
\n
\n \n
'