示例#1
0
def simulateddata(sdat, rdat, slo, sd, ss, par):

    # single reflector
    par['kxref'] = 50
    par['lxref'] = par['nximg'] - par['kxref']
    par['kzref'] = par['nzimg'] - par['nzpad'] / par['jzimg'] - 10
    Flow(
        'ref', None, '''
         spike nsp=1 mag=1 
         n1=%(nximg)d d1=%(dximg)g o1=%(oximg)g k1=%(kxref)d l1=%(lxref)d
         n2=%(nzimg)d d2=%(dzimg)g o2=%(ozimg)g k2=%(kzref)d |
         smooth rect1=25 repeat=3 |
         spray axis=2 n=1 o=0 d=1 |
         put label1=x label2=y label3=z 
         ''' % par)
    Plot('ref',
         'window | transp | smooth rect1=3 |' + fdmod.cgrey('pclip=100', par))
    Result('ref', ['ref', ss], 'Overlay')

    # source wavelet
    par['frq'] = 5
    par['kt'] = 50
    Flow(
        'wvl', None, '''
         spike nsp=1 mag=1
         n1=%(nt)d d1=%(dt)g o1=0   k1=%(kt)d |
         ricker1 frequency=%(frq)g |
         scale axis=123 |
         fft1 |
         window squeeze=n n1=%(nw)d min1=%(ow)g
         ''' % par)

    # source data on the surface
    for iexp in range(par['ns']):
        etag = "-e%03d" % iexp

        xsou = par['os'] + (par['fs'] + iexp * par['js']) * par['ds']
        isou = (xsou - par['oximg']) / par['dximg']
        Flow(
            'spk' + etag, 'wvl', '''
             pad beg2=%d n2out=%d |
             put label1=w label2=x label3=y o2=%g d2=%g |
             transp memsize=250 plane=12 |
             transp memsize=250 plane=23 
             ''' % (
                isou,
                par['nximg'],
                par['oximg'],
                par['dximg'],
            ))
    Flow('spk', map(lambda x: 'spk-e%03d' % x, range(par['ns'])),
         'cat space=n axis=4 ${SOURCES[1:%d]}' % par['ns'])

    # datumed source data
    zomig.Cdtone3(sdat, 'spk', sd, par)
    Result(sdat, fdmod.fgrey('', par))

    # datumed receiver data
    spmig.modelPW3(rdat, slo, sdat, 'ref', par)
    Result(rdat, fdmod.fgrey('', par))
示例#2
0
def simulateddata(sdat,rdat,slo,sd,ss,par):

    # single reflector
    par['kxref']=50
    par['lxref']=par['nximg']-par['kxref']
    par['kzref']=par['nzimg']-par['nzpad']/par['jzimg']-10
    Flow('ref',None,
         '''
         spike nsp=1 mag=1 
         n1=%(nximg)d d1=%(dximg)g o1=%(oximg)g k1=%(kxref)d l1=%(lxref)d
         n2=%(nzimg)d d2=%(dzimg)g o2=%(ozimg)g k2=%(kzref)d |
         smooth rect1=25 repeat=3 |
         spray axis=2 n=1 o=0 d=1 |
         put label1=x label2=y label3=z 
         ''' % par )    
    Plot('ref','window | transp | smooth rect1=3 |'
         + fdmod.cgrey('pclip=100',par))
    Result('ref',['ref',ss],'Overlay')

    # source wavelet
    par['frq']=5
    par['kt']=50
    Flow('wvl',None,
         '''
         spike nsp=1 mag=1
         n1=%(nt)d d1=%(dt)g o1=0   k1=%(kt)d |
         ricker1 frequency=%(frq)g |
         scale axis=123 |
         fft1 |
         window squeeze=n n1=%(nw)d min1=%(ow)g
         ''' % par)

    # source data on the surface
    for iexp in range(par['ns']):
        etag = "-e%03d" % iexp

        xsou = par['os'] + (par['fs'] + iexp * par['js']) * par['ds']
        isou = (xsou-par['oximg'])/par['dximg']
        Flow('spk'+etag,'wvl',
             '''
             pad beg2=%d n2out=%d |
             put label1=w label2=x label3=y o2=%g d2=%g |
             transp memsize=250 plane=12 |
             transp memsize=250 plane=23 
             ''' % (isou,
                    par['nximg'],
                    par['oximg'],
                    par['dximg'],) )
    Flow('spk',
         map(lambda x: 'spk-e%03d'  % x,range(par['ns'])),
         'cat space=n axis=4 ${SOURCES[1:%d]}'%par['ns'])

    # datumed source data
    zomig.Cdtone3(sdat,'spk',sd,par)
    Result(sdat,fdmod.fgrey('',par))

    # datumed receiver data
    spmig.modelPW3(rdat,slo,sdat,'ref',par)
    Result(rdat,fdmod.fgrey('',par))
示例#3
0
def originaldata(sdat,rdat,sd,par):

    Flow(  'wav',None,'spike nsp=1 mag=1 k1=1 n1=%(nt)d d1=%(dt)g o1=%(ot)g' % par)
    Result('wav','window n1=500 |' + fdmod.waveplot('',par))

    sigsbee.getdata('data',par)
    sigsbee.makeshots('sall','data',par)
        
    Flow('swin',
         'sall',
         '''
         window n3=%(ns)d f3=%(fs)d j3=%(js)d |
         bandpass flo=2 fhi=10
         ''' % par)
    
    for file in (['swin','sall']):
        for iexp in range(par['ns']):
            etag = "-e%03d" % iexp

            Plot(file+etag,file,
                 'window n3=1 f3=%d | put o2=%g |' % (iexp,par['jos']+iexp*par['jds'])
                 + fdmod.dgrey('pclip=99',par))
            Result(file,map(lambda x: file+'-e%03d' % x,range(par['ns'])),'Movie')


    encode.shot2grid('sdata-t','rdata-t','wav','swin',par)
    Flow('S-sou-t','sdata-t','window | transp' %par)
    Flow('S-rec-t','rdata-t','window | transp' %par)

    encode.time2freq('S-sou-t','S-sou-w',par)
    encode.time2freq('S-rec-t','S-rec-w',par)
    
    # datum through water
    spmig.datum3('S-dfs',
                 'S-dfr',
                 sd,
                 'S-sou-w',
                 'S-rec-w',par) 

    # window datumed data (all shots)
    Flow(sdat,'S-dfs','window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
    Flow(rdat,'S-dfr','window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)

    Result(sdat,fdmod.fgrey('pclip=99',par))
    Result(rdat,fdmod.fgrey('pclip=99',par))
示例#4
0
def originaldata(sdat, rdat, sd, par):

    Flow('wav', None,
         'spike nsp=1 mag=1 k1=1 n1=%(nt)d d1=%(dt)g o1=%(ot)g' % par)
    Result('wav', 'window n1=500 |' + fdmod.waveplot('', par))

    sigsbee.getdata('data', par)
    sigsbee.makeshots('sall', 'data', par)

    Flow(
        'swin', 'sall', '''
         window n3=%(ns)d f3=%(fs)d j3=%(js)d |
         bandpass flo=2 fhi=10
         ''' % par)

    for file in (['swin', 'sall']):
        for iexp in range(par['ns']):
            etag = "-e%03d" % iexp

            Plot(
                file + etag, file, 'window n3=1 f3=%d | put o2=%g |' %
                (iexp, par['jos'] + iexp * par['jds']) +
                fdmod.dgrey('pclip=99', par))
            Result(file, map(lambda x: file + '-e%03d' % x, range(par['ns'])),
                   'Movie')

    encode.shot2grid('sdata-t', 'rdata-t', 'wav', 'swin', par)
    Flow('S-sou-t', 'sdata-t', 'window | transp' % par)
    Flow('S-rec-t', 'rdata-t', 'window | transp' % par)

    encode.time2freq('S-sou-t', 'S-sou-w', par)
    encode.time2freq('S-rec-t', 'S-rec-w', par)

    # datum through water
    spmig.datum3('S-dfs', 'S-dfr', sd, 'S-sou-w', 'S-rec-w', par)

    # window datumed data (all shots)
    Flow(sdat, 'S-dfs', 'window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)
    Flow(rdat, 'S-dfr', 'window squeeze=n min1=%(oximg)g n1=%(nximg)g' % par)

    Result(sdat, fdmod.fgrey('pclip=99', par))
    Result(rdat, fdmod.fgrey('pclip=99', par))