Pengolahan Sinyal Seismic 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


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

6. 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

7. Hasilnya seperti yang terlihat dalam nomor 5. 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

5. Menampilkan semua stasiun dengan macro berikut. Sehingga kita bisa mengatur kisaran
periode yang akan kita potong (cut range)

SAC> r *II*KAPI*00*LHZ*SAC *II*PALK*00*LHZ*SAC
*II*DGAR*00*LHZ*SAC *II*WRAB*00*LHZ*SAC *II*UOSS*00*LHZ*SAC
SAC> cut 100 1300
SAC> r
SAC> p1

6. Dalam Macro harus disiapkan nilai APV3 Apparent Velocity pangkat 3 dan Distance dalam
meter. Kita harus hitung menggunakan persamaan berikut. Tentukan juga kisaran potongan
sinyal masing-masing stasiun umumnya 100 – 150 detik. Dalam latihan ini digunakan 130 detik
APV3 = (0.16*ı+7.9*1000)^3
Distance = ı* 111.267e+3



===================================================


Bentuk umum pole zeros yang dihasilkan adalah sebagai berikut
* **********************************
* 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




Komentar

Posting Komentar

Postingan populer dari blog ini

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 Bila perintah sac b
Baca selengkapnya

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.     LAI coupling menunjukkan
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