sconstruct

国外免费地震资料处理软件包

from rsfproj import *

 ####################################### 
 # Sigsbee 2A velocity model construct # 
 #######################################

PRFX = 'sigsbee2a_'
SUFX = '.sgy' 

for c in ('migration_velocity','stratigraphy','reflection_coefficients'):
    if(c=='migration_velocity'):
        v='vmig2A'
        o=3.055
        d=.0143
        s=.0003048
        l='Migration\ Velocity'
        a='y'
        u='\(km\/s\)'
    if(c=='stratigraphy'):
        v='vstr2A'
        o=3.048
        d=.00762
        s=.0003048
        l='Stratigraphic\ Velocity'
        a='y'
        u='\(km\/s\)'
    if(c=='reflection_coefficients'):
        v='reflectionCoefficients'
        o=3.048
        d=.00762
        s=1
        l='Reflection\ Coefficients'
        a='n'
        u=''
    h = c + 'head'
    t = PRFX + c + SUFX
    Fetch(t,'sigsbee')

    Flow([v,h],t,
         '''
         segyread tape=$SOURCE tfile=${TARGETS[1]} |
         put
         o1=0  d1=.00762   label1=Depth    unit1=km
         o2=%f d2=%f       label2=Distance unit2=km
         | scale rscale=%f
         ''' %(o,d,s) ,stdin=0)

    Result(v,v,'''grey color=i scalebar=y allpos=%s screenratio=.375 screenht=3 
                  wheretitle=t title=%s labelsz=4 titlesz=6 barlabel=%s''' % (a,l,u))

End()