Macro SAC untuk sinyal pergeseran

  • Mengubah percepatan menjadi pergeseran

               r *LHZ
               rtrend
               int
               r more         !  menghilangkan drift “penyok”
               p1



  •  Mengubah sampling sinyal menjadi lebih kecil (down sampling)


            SAC > qdp off
            SAC > bd x
            SAC > synch               # memperjelas sinyal to synchronice the time all components
            SAC > interpolate delta 0.5
            SAC > p1
            SAC > lh
                        DELTA : 5E-1
                        sebelumnya   DELTA : 1E+0

            NyQIST frequency = 1/(2*DELTA) = 1/ (2*0.5) = 1.0
            NyQIST frequency adalah maksimum frequency. DELTA adalah sampling per second
           
            SAC > decimate 2

Downsampling diperlukan untuk mengurangi harga filter sehingga benar-benar mendapatkan harga frekwensi corner yang layak. Nilai downsampling hanya sekitar 2 – 7.

  • Prosedur untuk phase picking gelombang P dan S

              SAC > r coco.sac
              SAC > lh DIST
                        Ubah dalam bentuk degree DEG = DIST / 111.267
                        Lihat dalam Tabel Jeffrey Bullens (Gambar 1) untuk DEGREE 
                        sehingga dapat diketahui kapan waktu kedatangan gel P dan S
                        Misal :     COCO
                                      DEG  = 1.201442e+03 / 111.267 = 10.8
                                    dalam tabel didapatkan waktu gel P = 3.70 dan S = 5.80
                                    
Kemudian observasi sinyal secara manual untuk mendapatkan nilai waktu sebenarnya

            SAC > bd x
            SAC > qdp off
            SAC > r coco.sac
            SAC > synch
            SAC > ppk  ----------> pick gel P dan S
            SAC > lh 

            Bandingkan hasil obs dengan tabel, usahakan deviasi tidak terlalu jauh (+/- 40 dtk)

            SAC > w over   

             menimpa nama file lama

            SAC > ch a -12345 t0 -12345

           menghapus picking dan kembali ke sinyal asli

  • Untuk mengubah dari gambar sac menjadi ps file:

                      bd sgf    # start sac graphics file (sgf)
                       p1
                      sgftops f001.sgf COCO_enr.ps  # convert to sgf to ps file
                      quit from sac
                      rm -f *.sgf   # remove sgf files
                      gs COCO_enr.ps

  • Membuat label judul, absis, dan ordinat

                      title on
                      xlabel "Time @(s@)"
                      ylabel "Velocity @(nm/s@)"
                      * set the xlabel and ylabel
                      ylim all
                      * set the y-axis to be the same
                      xlim o o +8
                set the time axis to be from origin time to 8 sec after the origin time


Gabungkan tiga komponen sekaligus

               rm -f *.sgf
               setbb sacfile "$sta$.*"  #setbb sacfile "$sta$.[zne]"
               setbb psfile "$sta$_3c.eps"  #set three component data

               qdp off
               r COCO.sac
               rmean
               bp p 2 n 4 c 1 10 #4th order bandpass butterworth filter with 1 - 10 Hz

               setbb ds '( concatenate ' &1,kstnm ' ' &1,kzdate ' ' &1,kztime ' ) '
               evaluate to r &1,dist
              title '( concatenate ' %ds% ' ', dist = ' ' %r km ') '
              * get the stn name, reference time, and distance from SAC header
              * and put them into the title
              fileid l ul
              fileid t l kcmpnm
             * display the component name
             * fileid Controls the file id display found on most SAC plots.
             * TYPE LIST hdrlist  :  Define a list of header fields to display in the
             * fileid.
             * LOCATION UL  :  Place file id in upper left hand corner.

Perintah envelope untuk membalikkan sinyal 90 derajat atau menjadikan bernilai positif dengan transformasi Hilbert. 


Gambar 1. Tabel Jeffery Bullens




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