* フーリエスペクトル解析 [#heb33753]

（フーリエ）スペクトル解析は，ランダムデータから有用なシグナルを抽出するために，地震記録の分析においても頻繁に使われる基本的な解析手法である．通常の地震観測記録は，揺れの振幅の時間変化を記録するいわゆる時系列データであるが，地盤や建築構造物はある特定の周期（固有周期）の地震動を増幅させる特徴があることから，時系列波形データを周期ごとの振幅，すなわちスペクトルに分解することにより，地盤の増幅特性や構造物の応答特性を明らかにすることができる．また地盤や建物に対する入力としての地震動は，震源断層のせん断すべり運動によって発生し，地表に伝播したものであるが，こうした入力地震動に対するスペクトル解析によって，震源の破壊過程や伝播経路における減衰特性の解明が期待される．

** 1. フーリエスペクトルと位相スペクトル [#j210608f]
** 1. フーリエスペクトル [#j210608f]
*** 1.1　フーリエ級数 [#z324e51b]

時間#tex(t)の関数#tex(x(t))が区間#tex(\left[-\frac{T}{2},\frac{T}{2}\right])を基本とする任意の周期#tex(T)の関数であるとき，次式のように三角関数#tex(\cos \frac{2n\pai t}{T})と#tex(\sin \frac{2n\pai t}{T})の線形和で級数展開することができる．
時間$t$の関数$x(t)$が区間$\left[-T/2,T/2\right]$を基本とする任意の周期$T$の関数であるとき，次式のように三角関数$\cos \frac{2n\pi t}{T}$と$\sin \frac{2n\pi t}{T}$の線形和で級数展開することができる．

** 2. 震源情報の関連付け [#n7eacf9a]
#tex(x(t)\sim \frac{a_0}{2}+\sum_{n=1}^\infty \left( a_n\cos\frac{2n\pi t}{T}+b_n\sin\frac{2n\pi t}{T} \right) )

取得された強震観測記録は，各種の解析に用いるために地震の震源情報と関連付けられる必要がある．ここでの震源情報とは，地震の発生日時，震源位置（緯度，経度，深さ），マグニチュード，などであるが，これらの震源情報には地震発生直後に気象庁から発表される速報値と，大学やその他公的機関のデータを併せて再推定され，地震発生の翌日ごろに出される暫定値（一元化震源とも呼ばれる），および地震発生からおよそ1年後に発刊される気象庁地震・火山月報（カタログ編）に記載された確定値があり，最終的には確定値が公式の震源情報として用いられる．このため，データベースにおける震源情報についてもこれに対応して更新していくことが重要である．なお，2012年6月現在，速報値は[[気象庁HP:http://www.jma.go.jp/jp/quake/]]より，暫定値（一元化震源要素）は防災科学技術研究所[[Hi-netのHP:http://www.hinet.bosai.go.jp/exclusive_sites/]]より取得可能である((一元化震源要素の取得には防災科学技術研究所へのユーザ登録が必要))．また，地震・火山月報（カタログ編）は，[[一般財団法人気象業務支援センター:http://www.jmbsc.or.jp/index.html]]より購入することができる．
このとき，右辺の級数をフーリエ級数またはフーリエ展開，その係数$a_n, b_n$をフーリエ係数と呼ぶ．フーリエ係数は次式により求められる．

強震観測データと震源情報の関連付けには，実際はかなりしっかりした確認作業が必要である．基本的にはデータの記録開始時刻と地震の発震時刻との時間差と，震源距離から求められる見かけの地震波伝播速度がP波速度の場合は5〜7 km/s前後，S波速度で3〜4 km/s前後であれば両者は関連付けられる可能性が高いが，そのほかにもマグニチュードから推測される振幅レベルの確認なども必要となる．また大地震後の余震や群発地震など，同等規模の地震が続発した場合にはより慎重な作業が要求される．いずれにしても，作業の完全な自動化は困難であり，観測者の目視等による選別が必要となる．
#tex(a_n=\frac{2}{T} \int_{-T/2}^{T/2}x(t)\cos\frac{2n\pi t}{T}dt \\ b_n=\frac{2}{T} \int_{-T/2}^{T/2}x(t)\sin\frac{2n\pi t}{T}dt)

** 3. 複数の観測点データのグルーピング [#i5aa7650]
なおこのとき，関数$\cos\frac{2n\pi t}{T}$と$\sin\frac{2n\pi t}{T}$は直交関数列であり，次のような特徴を持つ．

比較的距離の近い複数の観測点で観測をしている場合には，ペーストアップ図を作成することで関連付けした震源情報の確からしさを検証することができる．特に2011年東北地方太平洋沖地震(Mw 9.0)のような大規模地震における遠方の観測点記録の場合，近傍の小地震の記録としてグルーピングされてしまう可能性がある．さらにノイズの除去も重要な作業である．特に地表に設置されている観測点の場合，交通振動や波浪によるノイズ，さらには空振ノイズ（火山噴火，航空機や砲撃などに伴う空気振動で生じるノイズ）などで強震計がトリガーされる場合がある．これらのノイズ波形は卓越周期が通常の地震波形と大きく異なる，P波とS波が分かれていないなどの点から区別することが可能であるが，震源距離が非常に短い地震記録でもP波とS波の到達時間差がほとんどなく，判別が難しいケースがあることには留意する必要がある．
#tex(\begin{align} \int_{-T/2}^{T/2}\cos\frac{2m\pi t}{T}\cos\frac{2n\pi t}{T}dt =\int_{-T/2}^{T/2}\sin\frac{2m\pi t}{T}\sin\frac{2n\pi t}{T}dt &=0 (m \neq n) \\ &=T/2 (m=n) \end{align})

** 4. 観測点情報 [#w7524582]
さらにm=nのときの積分値を1に正規化すると（この場合は係数$\sqrt{2/T}$を乗じる）正規直交関数列となる．

震源情報と同様に，強震記録が取得された観測点に関する情報も，データ解析を行う際には重要となる．地盤観測の場合には，地震計の設置位置（地表or地中，地中の場合は設置深度）や地盤種別，可能であれば地盤の密度や弾性波速度などの検層データがあると望ましい．建物の観測では，設置建物の構造種別や階数，詳細な設置場所などの情報が重要である．両者に共通する情報としては設置地点の緯度経度や高度，強震計の周波数応答特性，観測期間などがある．また，長期にわたって観測が続けられている場合，担当者の変更などによって観測期間中の強震計設置位置の移動や欠測期間，観測点の追加，削減などの履歴が失われる危険性があるため、これらのログを確実に保存しておくことが重要である．
三角関数を指数関数に置き換えることにより，次のように複素フーリエ級数が導かれる．
#tex(x(t)\sim \sum_{-\infty}^\infty c_ne^{i2\pi nt/T}\\ c_n=\frac{1}{T} \int_{-T/2}^{T/2}x(t)e^{-i2\pi nt/T})

** 5. データベースの作成と公開 [#sc79e4f2]
*** 1.2 フーリエ積分 [#af303a72]
複素フーリエ級数の周期区間を無限大とする極限を考えたとき，次の複素フーリエ積分が成り立つ．また，関数$x(t)$の無限積分が収束し（絶対積分可能であり），全区間で連続，かつ区分的に滑らかであるならば，フーリエ積分はもとの関数$x(t)$に一致する（フーリエの積分定理）．

取得された強震データは，観測者が直接活用・分析する場合を除き，第三者に使いやすい形式で整理し，保存しておくことが望ましい．データのフォーマットについても，可能であれば現在一般的に広く使われている防災科学技術研究所のK−NETフォーマットなどに統一しておくと，他機関の観測データと併せて用いる際に便利であり，また防災科学技術研究所が公開しているユーティリティプログラムを適用することが可能となる．データベースには震源情報と観測点情報が必須であり，個々の震源，および観測点情報にリンクした観測データのテーブルを作成しておくことが有用である．震源情報，観測点情報として必要な要素は2,4で挙げた．さらに震源情報と観測点情報のそれぞれから観測データを呼び出すための検索機能を付加するとよい．現在は，[[防災科学技術研究所:http://www.seis.bosai.go.jp/]]や[[港湾空港技術研究所:http://www.mlit.go.jp/kowan/kyosin/eq.htm]]など，ウェブ上で多数の強震観測データベースが公開されているので，これらを参考にデータベースHPを作成するのも有用である．データの公開先については，観測機関内部限定の場合から，広く一般にインターネット経由で公開する場合まで，機関ごとに適切な形態をとることができる．
#tex(x(t)=\int_{-\infty}^{\infty}X(\omega)e^{i\omega t}d\omega \\ X(\omega)=\frac{1}{2\pi}\int_{-\infty}^{\infty}x(t)e^{-i\omega t}dt)

このとき，$\omega=2\pi f$は角周波数，$f$は周波数で周期$T$の逆数である．

*** 2.3 フーリエ振幅スペクトルと位相スペクトル，および位相差分スペクトル [#xf32fba9]
複素フーリエ係数$c_n$やフーリエ成分$X(\omega)$は複素数であり，実数$a_n, b_n$，または$a(\omega), b(\omega)$を用いて次式のように表すことができる．

#tex(c_n=a_n+ib_n=\sqrt{a_n^2+b_n^2} e^{i\theta_n}\\ \theta_n=\tan^{-1}(b_n/a_n)) 
#tex(X(\omega)=a(\omega)+ib(\omega)=\sqrt{a(\omega)^2+b(\omega)^2} e^{i\theta(\omega)}\\ \theta(\omega)=\tan^{-1}(b(\omega)/a(\omega))) 

このときの絶対値$|X(\omega)|=\sqrt{a(\omega)^2+b(\omega)^2}$はフーリエ振幅スペクトルを構成し，角周波数$\omega$の波の振幅を表す．また偏角（位相角）$\theta(\omega)$は位相スペクトルを成す．$x(t)$が実数であるとき，$|X(\omega)|=|X(-\omega)|$，$\theta(\omega)=-\theta(-\omega)$が成り立ち，$x(t)$は次式のような実数関数の積分として表現される．

#tex(x(t)=2\int_0^\infty|X(\omega)|\cos(\omega t+\theta(\omega))d\omega)

上式は，任意の時系列データ$x(t)$が，位相差$\theta(\omega)$で重なり合う振幅$|X(\omega)|$，周期$T=2\pi/\omega$の余弦関数の群に分解されることを示している．

位相スペクトルに対して，隣り合う周波数の位相角の差（位相差分）の頻度分布を位相差分スペクトルと呼ぶ．位相差分は負の値，すなわち0から$-2\pi$の間で定義される．地震波形記録のような非定常波の位相差分スペクトルは，元となった時系列データの包絡形状によく似ることが指摘されている．

** 3 自己相関関数とパワースペクトル [#s372e05d]
時間の関数$x(t)$に対し，自己相関関数$C(\tau)$は次式に示すような$\tau$時間離れた関数の積の平均値として定義される．

#tex(C(\tau)=E[x(t)x(t+\tau)])

定常確率過程の場合，上式は時間平均で置き換えられ，次のように時間ずれ$\tau$のみの関数となる．

#tex(\begin{align} C(\tau)&=\overline{x(t)x(t+\tau)}\\ &=\lim_{T \to \infty} \frac{1}{T} \int_{-T/2}^{T/2}x(t)x(t+\tau)dt \end{align})


** 4 相互相関関数とクロススペクトル，コヒーレンス [#o864d274]

** 5 フーリエスペクトル比と伝達関数 [#v46e570f]

参考文献 大崎順彦：新・地震動のスペクトル解析入門, 鹿島出版会　日野幹雄：スペクトル解析, 朝倉書店　秋山宏：建築物の耐震極限設計, 東京大学出版会
トップ   新規 一覧 単語検索 最終更新   ヘルプ   最終更新のRSS