Pengolahan Sinyal Seismic SAC

PENDAHULUAN SAC 

1.    Install SAC terlebih dahulu dengan mengetikkan :

% echo $SACAUX
                                          /usr/local/geophysics/sac/sac/aux
                                         % echo $PATH

2.    Tambahkan dalam lingkungan (environment) bashrc

${HOME}/.bashrc 
Atau
${HOME} vi bashrc
                       
                       # added for Seismic Analysis Code (SAC)

                        export SACAUX=/usr/local/geophysics/sac/sac/aux
export PATH=/usr/local/geophysics/bin:/usr/local/geophysics/sac/
local_bin:/usr/local/geophysics/sac/sac/bin:$PATH
setenv SACHOME /usr/local/geophysics/sac/sac
                        setenv SACAUX $SACHOME/aux

                        set path=($SACHOME/bin $path)

3.    Buka file package dari SAC

% tar zxvf /usr/local/geophysics/sac/Sac_Tutorial_2009.tar.gz
% cd SAC_Tutorial/SAC_Tutorial_WF
% wget http://geophysics.eas.gatech.edu/people/zpeng/Teaching/
Sac_Tutorial_2009.tar.gz
% tar zxvf Sac_Tutorial_2009.tar.gz
% cd SAC_Tutorial/SAC_Tutorial_WF

4.    Setelah sukses instalasi, ketikkan ls untuk menampilkan file SAC komponen LHZ jaringan II

Perintah ls :
 sugeng@sugeng-laptop:~/DATE$ ls *II*00*LHZ*SAC

2010.298.14.43.17.0695_2010.298.14.56.41.0695.II.KAPI.00.LHZ.R.SAC
2010.298.14.43.48.0695_2010.298.14.57.28.0695.II.PALK.00.LHZ.R.SAC

5. Ketikkan ini untuk menampilkan file Pole Zero pilihan komponen LHZ jaringan II

ls *PZs_II*LHZ_00* -l

SAC_PZs_II_AAK_LHZ_00_2010.179.19.30.00.0000_99999.9999.24.60.60.99999
SAC_PZs_II_ABPO_LHZ_00_2010.272.00.00.00.0000_99999.9999.24.60.60.99999

Bentuk umum pole zeros yang dihasilkan :
* **********************************
* Network   (KNETWK): II
* Station   (KSTNM) : KAPI
* Component (KCMPNM): LHZ
* Location  (KHOLE) : 00
* Start             : 1999,037,00:00:00
* End               : N/A
* Input Units       : METERS
* Output Units      : COUNTS
* Latitude          :  -5.014200
* Longitude         :  119.751700
* Elevation         :   300.0
* Depth             : 300.0
* Dip               : -90.0
* Azimuth           :   0.0
* Instrument Comment: S/N #110
* Channel Flag      : CG
* **********************************
ZEROS 3
POLES 4
-0.0048  0.0000
-0.0739  0.0000
-23.0559  22.8027
-23.0559  -22.8027
CONSTANT 4.171975e+12

  • Kita akan merubah susunan dan nama file pole zero SAC_PZs_II_KAPI_BHZ_00_2005.001.00.00. 

  • Buat macro program dan simpan dengan nama file seperti look.m. Sebelumnya cek dulu file pole zeros karena tidak semua jumlah barisnya sama

sed -n '25p' SAC_PZs_II_KAPI_LHZ_00*  > look1
sed -n '19,24p' SAC_PZs_II_KAPI_LHZ_00* > look2
cat look1 look2 > look3
cp look3 PZSKAPI

  • Hasilnya seperti yang terlihat dibawah ini. Lakukan untuk seluruh stasiun yang akan diolah

CONSTANT 4.171975e+12
ZEROS 3
POLES 4
-0.0048  0.0000
-0.0739  0.0000
-23.0559  22.8027
-23.0559  -22.8027

Komentar

Posting Komentar

Postingan populer dari blog ini

Metode Magnet Sebagai Prekursor Gempabumi

Gambar
Suaidi Ahadi, ST, MT 1 & Dr. Sugeng Pribadi 2   1 Peneliti BMKG & Kandidat Doktor Magnet Bumi ITB 2 Doktor Bidang Seismologi dan Tsunami BMKG Atmosfer mempunyai sifat dinamis dalam hal pemantauan pergerakan bumi seperti pada pengamatan iklim global. Ketika terjadi gempabumi, gaya stress batuan memancarkan gelombang elektromagnetik dan sejumlah ion listrik saat melintasi lapisan bumi atau dikenal dengan LAI ( Lithosphere Atmosphere Ionosphere ) coupling seperti pada Gambar 1. Perubahan densitas elektron di atmosfer dijadikan dasar prekursor gempabumi menggunakan metode elektromagnetik (Kamogawa, 2006). Sebelum terjadinya gempabumi, batuan di kedalaman sumber mengalami peningkatan gaya-gaya stress sampai akhirnya pecah sehingga menimbulkan pergeseran pada bidang litosfer serta getaran di permukaan (Johnston, 2002). Efek gempabumi juga mengakibatkan perubahan ion kemagnetan bumi atau dikenal dengan piezomagnetisme.  Gambar 1.     LA...
Baca selengkapnya

Install SAC

Siapkan sistem berbasis Linux sebagai media aplikasi SAC (misal; Ubuntu) Unduh software SAC. Sesuaikan dengan prosesor komputer anda: SAC 64 bit atau 32 bit http://www.iris.edu/dms/nodes/dmc/software/downloads/sac Masuk ke dalam root $ sudo su password:   Simpan dalam directory  $ cd /usr/local/ Buka extract file $ tar -xvf  sac-101.5c-linux_x86_32.tar.gz Buka environment setting dan atur PATH SAC (bash shell). $ vi .bashrc Perhatian agar tidak merubah baris dan kolom dalam bashrc .  Untuk mengetahui nomor baris , ketik $ s e n u Penambahan PATH selalu diletakkan di akhir baris export SACHOME=/usr/local/sac export PATH=${PATH}:${SACHOME}/bin export SACAUX=${SACHOME}/aux export SAC_PPK_LARGE_CROSSHAIRS=1 Simpan dan compile ba shrc $ chmod 755 .bashrc $ ./bashrc  atau $ .. bash rc Perintah SAC dapat dibuka di semua folder. $ sac SAC > Untuk mengetahui dimana letak file sac $ which sac ...
Baca selengkapnya

Download data IRIS Wilber orfeus II

Gambar
Buka http://www.orfeus-eu.org/ Tentukan event yang akan kita ambil berdasarkan triwulan (Gambar 1) Klik list all event Klik precise event yang kita tentukan (Gambar 2) Klik Select All, kemudian Proceed (Gambar 3). Artinya kita pilih seluruh jaringan seismik global  Atur pilihan berdasarkan Station, Network & Distance. Untuk pengolahan W phase, gunakan pilihan Distance mulai jarak 10 derajat. Untuk pilihan channel pilih ALL untuk keperluan W phase. Isilah personal information dan pilih format data SEED tar file. Kemudian klik Process Request (Gambar 4).  Cek notifikasi email anda. Bila unduhan selesai IRIS akan memberitahu anda. Buka web unduhan dan Save file as SEED.tar (WinRAR) (Gambar 5). Berikut ini adalah Gambar 1-5.
Baca selengkapnya