Programming
Terimakasih sebelumnya untuk pak Eddy Hartantyo yang kemarin memberikan clue dalam membuat kurva dispersi..
oke, jadi setelah mengutakatik, kira-kira proses dari data multi-channel refraksi ke phase velocitiy vs frekuensi itu seperti ini :
1. data timeseries tiap receiver (x,t) di fft 1d tiap x menjadi (x,f)
2. lalu data (x,f) di fft2d menjadi (k,f) lalu di shift (di matlab commandnya fftshift) sehingga bagian tengah merupakan frekuensi 0.
3. gunakan rumus v=2*pi*f/k untuk mengubah (k,f) menjadi (v,f)
referensi :AN APPRAISAL OF SURFACE WAVE METHODS FOR SOIL CHARACTERIZATION
By Khiem Tat Tran
See ya,
L
Transformasi Domain F-V Kurva Dispersi Data Refraksi menggunakan Matlab (MASW/SASW)
Hello Digital world.Terimakasih sebelumnya untuk pak Eddy Hartantyo yang kemarin memberikan clue dalam membuat kurva dispersi..
oke, jadi setelah mengutakatik, kira-kira proses dari data multi-channel refraksi ke phase velocitiy vs frekuensi itu seperti ini :
1. data timeseries tiap receiver (x,t) di fft 1d tiap x menjadi (x,f)
2. lalu data (x,f) di fft2d menjadi (k,f) lalu di shift (di matlab commandnya fftshift) sehingga bagian tengah merupakan frekuensi 0.
3. gunakan rumus v=2*pi*f/k untuk mengubah (k,f) menjadi (v,f)
referensi :AN APPRAISAL OF SURFACE WAVE METHODS FOR SOIL CHARACTERIZATION
By Khiem Tat Tran
 |
| Fig 1. Data refraksi yang akan ditransformasikan |
 |
| Fig 2. Transformasi menggunakan software Pickwin |
 |
| Fig 3. Hasil transformasi FK dan FV metode Park et al di Matlab |
See ya,
L
this is best
ReplyDeleteActually, according to Tran's thesis I mentioned above, the best method for dispersion curve for an active MASW data is the cylindrical beamformer method. Though, from my experience using my own data, it doesn't really show much difference between Cylindrical, Park et al, and this method.
ReplyDeletehello , i have one question to ask. i have got figure in f-k domain like figure3 (up) using matlab. but how to transform to f-v domain (like figure 3down) thank you very much
ReplyDeletehello, you can use the formula I mentioned above in the 3rd step.
Deletev=2*pi*f/k where v is stand for phase velocity.
you might need to use quite a lot of nested loop (for .. end) to find these "v" for each frequency and wavenumber.
hello mas leo, saya ingin minta tolong bagaimana persamaan matematika inversi kurva dispersi sehingga diperoleh profil Vs 1D terhadap kedalaman pada metode MASW? terimakasiih sebelumnya mas
ReplyDeleteJAWABAN VERSI PENDEK:
Deleterumus forward pakai dari cara perhitungan di referensi thesis yg saya tuliskan di atas.
rumus inversi pakai Gauss Newton atau sejenis monte carlo / simulated annealing juga bisa.
JAWABAN VERSI PANJANG:
Geofisika adalah bagian dari ilmu yang mempelajari bumi menggunakan kaidah/teori atau prinsip-prinsip fisika (wikipedia).
Teori tiap metode-metode geofisika tersebut umumnya dalam bentuk..
d=f(m), yang artinya kalau misalnya kita masukin parameter fisis di dalam bumi "m" ke dalam rumus teori dasar suatu metode geofisika "f", nanti hasilnya adalah data teoritis "d" yang merupakan respon kalau misalnya kita ngukur di situ (umumnya orang bilang rumus forward engine/ formula/ modelling )
Oke, kita ngukur pakai alat di lapangan nih, kan dapet data respon bumi (data observasi). Gimana caranya tau dalem isi bumi (parameternya)? Kita utak-atik nilai "m" di rumus teori metode tsb sampai nilai "d" hasil rumus (data teoritis) cocok dengan "d" hasil ngukur (data observasi). Kalau cocok berarti ya berarti isi bumi itu "m". Proses ini namanya inversi :) Metodenya ada banyak. I personally recommend baca bukunya Richard Aster Parameter Estimation (Ada latihannya juga).
agak panjang yah? hehe, supaya paham aja :)
permisis bang leo, saya irsyad mahasiswa fisika KBK Geofisika universitas mulawarman, sekarang saya mengambil skripsi dengan judul MASW, saya kesulitan untuk menemukan sofwer untuk membuat kurva dispersi, boleh kah saya minta scrip matlab bang leo, untuk membuat kurva dispersi, mohon bantuannya
ReplyDeleteHi Irsyad, kamu bisa gunakan pickwin. Software tsb bisa dari editing, prosessing, transformasi, hingga inversi.
ReplyDeleteHalo Mas, kenalin saya Mita Dari Geofisika Unhas. Mau nanya, untuk mendapatkan persamaan mathematis cara mendapat kurva dispersi kalkulasi caranya gimana ya Mas? Terima kasih🙏🙂
ReplyDelete