Calculate PRFs of a seismic array in batch

This example provide a procedure for batch calculating PRFs of NCISP-III array (Zheng et al., 2007)

Download this example

Download link: ex-batrf.tar.gz

Download and unzip the file to any directory.

wget -O ex-batrf.tar.gz
tar -xzf ex-batrf.tar.gz

The package include Data.ZX with 47 sub folders named as network.station, which include sac files cut out from 300 s after origin time to 1000 s after origin time seismic events. The SACHeader` of station information (``netwk, stnm, stla and stlo) had been written into sac files.

Write a script for batch calculation

Seispy has provided a command as setpar, which will modify values in the .cfg file. Combining with a list of station information with 3 columns (staname, stla, stlo), We can use the setpar and prf command in shell script ( for RF calculation in batch. In this example, please execute in ex-batRF directory for the batch calculation.

while read staname stla stlo
    setpar rf.cfg path datapath ./Data.ZX/$staname
    setpar rf.cfg path rfpath ./RFresult/$staname
    prf rf.cfg
done < sta.lst
  • The sta.lst includes station information.

  • The rf.cfg is the configure file. In addition to the datapath and rfpath, other parameters can be set in advance.

  • If set up a local catalog, the calculation time will be saved.

CCP stacking along a linear seismic array

Download this example

Download link: ex-ccp.tar.gz

Download and unzip the file to any directory.

wget -O ex-ccp.tar.gz
tar -xzf ex-ccp.tar.gz

This package include:

  • A Folder with PRFs (RFresult)

  • A configure file with cfg format (ccp.cfg).

Time-to-depth conversion

specify following parameters in the ccp.cfg for time-to-depth conversion

  • rfpath: Path to PRFs

  • depthdat: Output data structure after time-to-depth conversion

  • stalist: The list of station information with 3 columns (staname, stla, stlo)

Then execute follow command to do the conversion

rf2depth ccp.cfg

After the time-to-depth conversion, the RFdepth.mat will be generated following the path set as depthdat.

CCP stacking

specify following parameters in the ccp.cfg for CCP stacking. See Stack PRFs with Common Conversion points (CCP) method in detail.

  • stackfile: Output file for the CCP stacking

  • stack_sta_list: Stations used for the CCP stacking. In this example, we use the same stations as the stalist.

  • Section [bin] for parameters of stacking bins

  • Section [line] for location of the profile.

  • Section [stack] for Depth samples.

Then execute following command for CCP stacking

ccp_profile ccp.cfg -t

An text file set by stackfile will be generated including 6 columns: bin_lat, bin_lon, bin_distance, depth, amplitude and stack_num.

Plot the stacking image

Seispy dose not provide any functions or script for plotting stacking image. Here we provide a script with GMT6 format for plotting the image. Just run in ex-ccp directory.


CCP stacking along the seismic array


Zheng T, Chen L, Zhao L, et al. Crustal structure across the Yanshan belt at the northern margin of the North China Craton[J]. Physics of the Earth and Planetary Interiors, 2007, 161(1-2): 36-49.