view release on metacpan or  search on metacpan

lib/App/SeismicUnixGui/sunix/inversion/suinvzco3d.pm  view on Meta::CPAN

 vfile=                  file containing velocity array v[nz]		

 nz=                    number of z samples (1st dimension) in velocity

 nxm=			number of midpoints of input traces		

 ny=			number of lines 				



 Optional Parameters:							

 dt= or from header (dt) 	time sampling interval of input data	

 offs= or from header (offset) 	source-receiver offset	 	

 dxm= or from header (d2) 	sampling interval of midpoints 		

 fxm=0                  first midpoint in input trace			

 nxd=5			skipped number of midpoints (see note)		

 dx=50.0                x sampling interval of velocity		

 fx=0.0                 first x sample of velocity			

 dz=50.0                z sampling interval of velocity		

 nxb=nx/2		band centered at midpoints (see note)		

 fxo=0.0                x-coordinate of first output trace 		

 dxo=15.0		horizontal spacing of output trace 		

 nxo=101                number of output traces 			",	

 fyo=0.0		y-coordinate of first output trace		

 dyo=15.0		y-coordinate spacing of output trace		

 nyo=101		number of output traces in y-direction		

 fzo=0.0                z-coordinate of first point in output trace 	

 dzo=15.0               vertical spacing of output trace 		

 nzo=101                number of points in output trace		",	

 fmax=0.25/dt		Maximum frequency set for operator antialiasing 

 ang=180		Maximum dip angle allowed in the image		

 verbose=1              =1 to print some useful information		



 Notes:									



 This algorithm is based on formula (50) in Geophysics Vol. 51, 	

 1552-1558, by Cohen, J., Hagin, F., and Bleistein, N.			



 Traveltime and amplitude are calculated by ray tracing.		

 Interpolation is used to calculate traveltime and amplitude.		", 

 For each midpoint, traveltime and amplitude are calculated in the 	

 horizontal range of (xm-nxb*dx, xm+nxb*dx). Velocity is changed by 	

 linear interpolation in two upper trianglar corners whose width is 	

 nxc*dx and height is nzc*dz.						",	



 Eikonal equation will fail to solve if there is a polar turned ray.	

 In this case, the program shows the related geometric information. 	

 

 Offsets are signed - may be positive or negative. 			", 



=head2 User's notes (Juan Lorenzo)
untested

=cut


=head2 CHANGES and their DATES

=cut

use Moose;
our $VERSION = '0.0.1';


=head2 Import packages

=cut

use aliased 'App::SeismicUnixGui::misc::L_SU_global_constants';

use App::SeismicUnixGui::misc::SeismicUnix qw($go $in $off $on $out $ps $to $suffix_ascii $suffix_bin $suffix_ps $suffix_segy $suffix_su);
use aliased 'App::SeismicUnixGui::configs::big_streams::Project_config';


=head2 instantiation of packages

=cut

my $get					= L_SU_global_constants->new();
my $Project				= Project_config->new();
my $DATA_SEISMIC_SU		= $Project->DATA_SEISMIC_SU();