kzplanner.c

seismic software,very useful

/* KZMIG planing program */
#include "par.h"

char *sdoc = 
"KZPLANNER - KZMIG Planning to output parameters for kzvelo, rayt3d and kzmig\n"
"\n"
"kzplanner [parameters] 		\n" 
"\n"
"Required parameters:					\n"
"offmax=    maximum offset to migrate	\n"
"feather=   maximum feathering distance of input data \n"
"           (maximum crossline distance from source or receive to midpoint \n"
"minli=     minimum line number of input data \n"
"maxli=     maximum line number of input data \n"
"minti=     minimum trace number of input data \n"
"maxti=     maximum trace number of input data \n"
"minlo=     minimum line number of migration output \n"
"maxlo=     maximum line number of migration output \n"
"minto=     minimum trace number of migration output \n"
"maxto=     maximum trace number of migration output \n"
"Optional Parameters: \n"
"ds=12.5    distance of two inline traces (trace number difference of 1) \n"  
"dl=20      distance of two lines (line number difference of 1) \n"
"apers=5000 inline migration aperture (m or ft) \n"
"aperl=5000 crossline migration aperture (m or ft) \n"
"minzo=0    minimum depth of migration output \n"
"maxzo=8000 maximum depth of migration output \n"
"dzo=10     depth interval of migration output \n"
"ntvel=6    trace number increment of velocity grid \n"
"nlvel=4    line number increment of velocity grid \n"
"dzvel=50   depth interval of velocity grid (m or ft) \n"
"dztt=50    depth interval of travel time grid (m or ft) \n"
"nttt=10    trace number increment of travel time grid  \n"
"nltt=6     line number increment of travel time grid  \n"
"nlray=10   line number increment of ray-tracing source grid \n"
"ntray=16   trace number increment of ray-tracing source grid \n"
"\n"
"Example: \n"
" kzplanner offmax=10000 feather=3280.8 \\ \n"
" minli=5001 maxli=7065 minti=2950 maxti=5500 \\ \n"
" minlo=5899 maxlo=6179 minto=4100 maxto=5050 \\ \n"
" ds=41.01 dl=65.616 apers=16404 aperl=16404 \\ \n"
" minzo=5000 maxzo=28000 dzo=30 ntvel=6 nlvel=4 dzvel=200 \\ \n"
" dztt=200 nttt=10 nltt=20 nlray=10 ntray=16 \n"
"\n"
"\n"
"AUTHOR:		Zhiming Li,       ,	4/27/2000   \n"		    
;

main(int argc, char **argv)
{

float ds, dl, apers, aperl, offmax, feather;
int minli, maxli, minti, maxti, minlo, maxlo, minto, maxto;
float minzo,maxzo,dzo;
int ntvel,nlvel; 
float dzvel, dztt;
int nttt, nltt, ntray, nlray;

float tmp;
int dstrace, dsline;
int smins,smaxs,lmins,lmaxs;
int sminv,smaxv,lminv,lmaxv;
int smint,smaxt,lmint,lmaxt;
int nzo,nss,nls,nst,nlt,nsv,nlv,nzv,smin;
int timin,timax,limin,limax,nzt;
float zmint;

/* initialization */
initargs(argc,argv);
askdoc(1);

/* get required parameters */
if (!getparfloat("offmax", &offmax)) err(" offmax missing \n");
if (!getparfloat("feather", &feather)) err(" feather missing \n");
if (!getparint("minli", &minli)) err(" minli missing \n");
if (!getparint("maxli", &maxli)) err(" maxli missing \n");
if (!getparint("minti", &minti)) err(" minti missing \n");
if (!getparint("maxti", &maxti)) err(" maxti missing \n");
if (!getparint("minlo", &minlo)) err(" minlo missing \n");
if (!getparint("maxlo", &maxlo)) err(" maxlo missing \n");
if (!getparint("minto", &minto)) err(" minto missing \n");
if (!getparint("maxto", &maxto)) err(" maxto missing \n");

/* optional parameters */
if (!getparfloat("ds", &ds)) ds = 12.5;
if (!getparfloat("dl", &dl)) dl = 20.;
if (!getparfloat("apers", &apers)) apers = 5000.;
if (!getparfloat("aperl", &aperl)) aperl = 5000.;
if (!getparfloat("minzo", &minzo)) minzo = 0.;
if (!getparfloat("maxzo", &maxzo)) maxzo = 0.;
if (!getparfloat("dzo", &dzo)) dzo = 10.;
if (!getparint("ntvel", &ntvel)) ntvel = 6;
if (!getparint("nlvel", &nlvel)) nlvel = 4;
if (!getparfloat("dzvel", &dzvel)) dzvel = 50.;
if (!getparfloat("dztt", &dztt)) dztt = 50.;
if (!getparint("nttt", &nttt)) nttt = 10;
if (!getparint("nltt", &nltt)) nltt = 6;
if (!getparint("nlray", &nlray)) nlray = 10;
if (!getparint("ntray", &ntray)) ntray = 16;

tmp = (apers+offmax/2.)/ds + 0.5;
dstrace = tmp;
tmp = (aperl+feather)/dl + 0.5;
dsline = tmp;

/* limits of input source and receive locations */
tmp = offmax/2./ds+0.5; 
timin = minti - (int)tmp;
timax = maxti + (int)tmp;
tmp = feather/dl+0.5; 
limin = minli - (int)tmp;
limax = maxli + (int)tmp;


/* output travel time grid limit */
tmp = (float)(maxto-minto)/nttt + 1.9999;
nst = tmp;
smint=minto;
smaxt=minto+(nst-1)*nttt;
tmp = (float)(maxlo-minlo)/nltt + 1.9999;
nlt = tmp;
lmint=minlo;
lmaxt=minlo+(nlt-1)*nltt;
zmint = minzo;
tmp = (maxzo-minzo)/dztt + 1.9999;
nzt = tmp; 


/* raytracing source limit */
smins = smint - dstrace;
if(smins<timin) smins = timin;
smaxs = smaxt + dstrace;
if(smaxs>timax) smaxs = timax;
lmins = lmint - dsline;
if(lmins<limin) lmins = limin;
lmaxs = lmaxt + dsline;
if(lmaxs>limax) lmaxs = limax;

tmp = (float)(smaxs-smins)/ntray + 1.9999;
nss = (int) tmp; 
tmp = (float)(lmaxs-lmins)/nlray + 1.9999;
nls = (int) tmp; 

smaxs = smins + (nss-1)*ntray;
lmaxs = lmins + (nls-1)*nlray;


/* velocity grid limit */
tmp = (float)(smaxs-smins)/ntvel + 1.999;
nsv = (int) tmp;
tmp = (float)(lmaxs-lmins)/nlvel + 1.999;
nlv = (int) tmp;
tmp = maxzo/dzvel + 1.9999;
nzv = (int) tmp;


/* ouput migration limit */
tmp = (maxzo-minzo)/dzo + 1.9999;
nzo = (int) tmp;


fprintf(stderr," ===  KZPLANNER OUTPUT === \n");
fprintf(stderr," =========================  \n");
fprintf(stderr,"  \n");

fprintf(stderr," ===  KZVELO SETUP === \n");
fprintf(stderr," number of depth samples of velocity output=%d \n", nzv);
fprintf(stderr," starting depth of velocity output =0 \n");
fprintf(stderr," depth sampling interval of velocity output =%g \n",dzvel);
fprintf(stderr," number of traces per line of velocity output=%d \n", nsv+2);
fprintf(stderr," starting trace number of velocity output=%d \n", smins-ntvel);
fprintf(stderr," trace number increment of velocity output=%d \n", ntvel);
fprintf(stderr," number of lines of velocity output=%d \n", nlv+2);
fprintf(stderr," starting line number of velocity output=%d \n", lmins-nlvel);
fprintf(stderr," line number increment of velocity output=%d \n", nlvel);
fprintf(stderr," \n");


fprintf(stderr," ===  RAYT3D SETUP === \n");
fprintf(stderr," first trace number to output travel time=%d \n", smint);
fprintf(stderr," trace number increment to output travel time=%d \n", nttt);
fprintf(stderr," number of traces per line to output travel time=%d \n", nst);
fprintf(stderr," first line number to output travel time=%d \n", lmint);
fprintf(stderr," line number increment to output travel time=%d \n", nltt);
fprintf(stderr," number of lines to output travel time=%d \n", nlt);
fprintf(stderr," minimum depth to output travel time=%g \n", minzo);
fprintf(stderr," depth interval to output travel time=%g \n", dztt);
fprintf(stderr," number of depth samples to output travel time=%d \n", nzt);
fprintf(stderr," trace number to first raytracing surface source=%d \n", smins);
fprintf(stderr," trace number increment of raytracing surface source=%d \n", ntray);
fprintf(stderr," number of traces per line of raytracing surface source=%d \n", nss);
fprintf(stderr," line number to first raytracing surface source=%d \n", lmins);
fprintf(stderr," line number increment of raytracing surface source=%d \n", nlray);
fprintf(stderr," number of lines of raytracing surface source=%d \n", nls);
fprintf(stderr," inline raytracing distance from source=%g \n", dstrace*ds);
fprintf(stderr," crossline raytracing distance from source=%g \n", dsline*dl);
fprintf(stderr," \n");

fprintf(stderr," ===  KZMIG SETUP === \n");
fprintf(stderr," depth of first output sample=%g \n", minzo);
fprintf(stderr," depth interval to output =%g \n", dzo);
fprintf(stderr," number of depth samples to output=%d \n", nzo);
fprintf(stderr," inline aperture of migration=%g \n", apers);
fprintf(stderr," crossline aperture of migration=%g \n", aperl);


return 0;

}