19/05/14

== What we did ==

前回からの問題点のおさらい
 * geophoneの生信号が小さい~600cnts (BSのgeophoneは ~1000cnts)
 * calibration後のgeophoneの信号が少しずつ発散していく

=== calibration後のgeophoneの信号が少しずつ発散していく問題 ===
谷岡さんと奥富さん曰く，Type-Aでも同じような問題が過去にあったそうだ。

19/05/13 奥富さんのメールから抜粋

''これはFoton上でzpk filterを作ると零点／極の周波数が
RTモデルのsampling周波数を基準に記述されるという仕様によるものらしく、
あまりsampling周波数から離れた位置に零点／極を置くと
sampling: ~ 10^3^ Hzの基準に対して < 10^-3^ Hzの零点／極による信号処理は
周波数領域で6桁以上の隔たりとなってしまうため、
計算機上ではあまりよろしくないことが起きるっぽいです''

これの解決法は，零点/極の位置を高周波側にずらす。RT modelのsampling rateを下げること。
(Type-Aは，RT modelのsamplingを16kHz→2kHzにして，極を3mHzくらいにしたそうだ)

とりあえず極の位置をずらすことから始める。
現状のgeophoneのToVel filterは，こんな感じだった

{{attachment:ToVelfilter.png||width=500}}

DC cutoffの極の位置が，0.1Hzである。(この時点で，Type-Aのうまく行った場合より隔たりが１桁ほど低い気もするけど)

V1のみ，DC cutoffの極の位置を0.3Hzに変更して，ゲインを高周波側で１になるように調整。変更していないV2,V3と同時に測定をしてみて，発散するかを見る。

(測定開始時の信号の様子をスクショするのと，Amplitudeをメモしておけば，よかった。。)

 * １時間後の信号の様子

{{attachment:signal_Test_1hour.png||width=800}}

左下がV1,左上がV2,右下がV3,右上がcaribration前のV1の生信号。V2,V3は今まで通り形が変になっている。一方，V1は今の所，そのようには見えない。

 * ２時間後の信号の様子

{{attachment:signal_test_2hour.png||width=800}}

V2,V3は1時間前に比べて，Amplitudeが大きくなっているのがわかる。一方，V1はその傾向が見られない。また形も変になっていない。

 * ３時間後の信号の様子

{{attachment:signal_Test_3hour.png||width=800}}

V2,V3はAmplitudeが大きくなっている。また，V1の形も変になり始めた。うーん。

そもそも，cutoffを0.3Hzにすると，0.2Hz付近のmicro seismic noiseの見かけが小さくなるんじゃ。。。

{{attachment:190514spec_V1_seis.png||width=800}}

cutoffを入れているので，数百mHzの領域が地震計の値に比べて，小さくなっている。これでは意味ない気がする。

=== まとめと考察 ===
 * DC cutoffのzpkを高周波側にしてみたが，時間が経つと結局発散していく。
 * 数百mHzの領域を変えると，そのあたりの地面振動が測定できなくなる。
 * そもそもType-Aの場合よりも，DC cutoffが高い周波数である。
 * zpkをこれ以上変更するよりも，RT modelのsampling rateを変更する方が，隔たり？という意味では大きく減らせる。