如何提取地震数据里的同相轴？

乐观

初学者，还不太明白。
不知道有没有程序。

meil32

地震资料处理流程与方法介绍，另一个压缩包是画地震信号同相轴或者其他剖面时候的必备代码。

meil32

fortran编写简单的形成地震同相轴，用的是雷克子波

1.       program main
   
2.         parameter (lt=100, nmx=60, dt=0.002)
   
3.         real  trace(lt), txdata(lt, nmx)
   
4.   
   
5.       f0=25.0
   
6.          
   
7.         open(1,file='synthetic_100_60.dat',access='direct',recl=4*lt)  
   
8.         do ix=1, nmx
   
9.            do it=1,lt
   
10.              trace(it)=0.0
    
11.            enddo
    
12.            do it=1,lt
    
13.               if(it.eq.30)then
    
14.                 call RICKER_WAVELET(trace(it), f0, lt-it, dt)
    
15.             elseif(it.eq.3*ix)then
    
16.                 call RICKER_WAVELET(trace(it), f0, lt-it, dt)
    
17.               elseif(it.eq.-3*ix+60)then
    
18.                 call RICKER_WAVELET(trace(it), f0, lt-it, dt)
    
19.               endif
    
20.             txdata(it,ix)=trace(it)
    
21.            enddo
    
22.         enddo
    
23.       do ix=1,nmx
    
24.           write(1,rec=ix)(txdata(it,ix),it=1,lt)
    
25.         enddo
    
26.       close(1)
    
27.         end
    
28. ******************************************************
    
29. c    this is the program to make a ricker wavelet with a main frequence
    
30.       SUBROUTINE RICKER_WAVELET(TRACE, F0, LT, DT)
    
31.       PARAMETER (PI=3.14159265359)
    
32.       REAL FO
    
33.       real TRACE(LT)
    
34.       DO IT=1,LT
    
35.         W=PI*F0*IT*DT
    
36.         TRACE(IT)=0.0
    
37.         TRACE(IT)=(1.0-2*W*W)*EXP(-W*W)
    
38.       END DO
    
39.    
    
40.       RETURN
    
41.    
    
42.       END

复制代码

乐观

感激不尽！！