地震動予測方法及び地震動ハザードマップの制作方法
【課題】不整形地盤においても適用可能な地震動予測法、及び、不整形地盤を含む任意の空間における震度ハザードマップの作成方法を提供する。
【解決手段】仮想地震に対して工学的基盤面における地震動波形を予測する。次いで、調査により得られる地盤データにより工学的基盤から地表面までの地震動の伝達関数を計算する。これらの操作をJISによる地域地盤メッシュ毎に行うことによって任意の空間における地震動予測、及び、震度ハザードマップの作成が可能となる。不整形地盤を含む空間における地震動の伝達関数は、メッシュ内の位置と不整形性の程度を表すパラメータとともに、相互に隣接するメッシュの伝達関数を重ね合わせることにより計算する。
【解決手段】仮想地震に対して工学的基盤面における地震動波形を予測する。次いで、調査により得られる地盤データにより工学的基盤から地表面までの地震動の伝達関数を計算する。これらの操作をJISによる地域地盤メッシュ毎に行うことによって任意の空間における地震動予測、及び、震度ハザードマップの作成が可能となる。不整形地盤を含む空間における地震動の伝達関数は、メッシュ内の位置と不整形性の程度を表すパラメータとともに、相互に隣接するメッシュの伝達関数を重ね合わせることにより計算する。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は地震動予測方法及び地震動ハザードマップの制作方法に係り、特に、不整形地盤を備えた地域における地震動伝達関数の補間手法を用いて地震動を容易かつ的確に予測する方法に関する。
【背景技術】
【0002】
一般に、地震動予測は、震源から地表への震動の伝達メカニズムを想定して算出されるが、将来予測される地震による地域の被害予測のための震度ハザードマップを作成する場合には、表層地盤はボーリング調査(点測定)に基づき数100m四方の独立した成層地盤メッシュとして扱われる。この場合、各ボーリング調査による地盤データを用いた一次元の地盤震動解析の分野においては、地震動伝達関数の計算方法として、等価線形化法による地盤震動解析プログラムSHAKEが広く使用されてきた(以下の非特許文献1参照)。また、これを改良してやや軟弱な地盤にまで適用範囲を広げた、周波数依存型等価線形化法FDELも開発されている(以下の非特許文献2参照)。これらは計算に必要な地盤パラメータが比較的容易に得られることから、実用性が高く評価されている。
【0003】
一方、採取された地盤データから地震動パラメータを算出し、この地震動パラメータに基づき、クリギングの手法によって、任意時点又は任意領域において補間された補間地震動パラメータ及びその補間精度を求める方法が知られている(以下の特許文献1参照)。また、以下の特許文献2においても地盤メッシュ間の地震動の補間手法が記載されている。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開2005−156273号公報
【特許文献2】特開2003−287574号公報
【非特許文献】
【0005】
【非特許文献1】P.B.Schnabel, J.Lysmer and H.B.Seed: SHAKE a computer program for earthquake response analysis of holizontally layered sites, EERC, 72-12, 1972.(SHAKE)
【非特許文献2】杉戸真太, 合田尚義, 増田民夫: 周波数特性を考慮した等価ひずみによる地盤の地震応答解析法に関する一考察, 土木学会論文集, No.493/III-27, pp.49-58.1994.6.(FDEL)
【非特許文献3】Lysmer, J., Udaka, T., Tsai, C.-F. and Seed, H.B.: FLUSH a computer program for approximate 3D-analyisis of soil structure interaction problem, Proc. of EERC: 75-30. 1975(FLUSH)
【非特許文献4】Y.Furumoto, M.Sugito, A.Yashima:Frequency-Dependent Equivalent Linearized Technique for FEM Response Analysis of Ground, 12WCEE in Auckland, New Zealand (CD-ROM),2000.2(FDEL−FEM)
【非特許文献5】Sugito,M., Furumoto,Y., and Sugiyama,T., Strong Motion Prediction on Rock Surface by Superposed Evolutionary Spectra, 12th World Conference on Earthquake Engineering, CD-ROM, Auckland, New Zealand, January 2000.(強震動予測法(EMPR))
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、既存の手法で作られる震度ハザードマップ等の分布図は、地盤メッシュ(例えば500m四方)で区切られているのが一般的であり、滑らかなコンターマップとはなっていない。これは、地盤のモデル化が地盤メッシュ内の高々一点を基準にして行われているためで、日本のように起伏の激しい地形(=不整形地盤)では、正確な震度マップの作成方法とはいえない。地表での震度を予測するためには、調査済みの地盤モデルから、基盤から地表までの地震動の伝達関数を計算する必要があるが、既存の手法は、地盤メッシュ内での一次元的な応答計算に留まっている。このため、地盤構造に急激な変動がある不整形地盤を有する地域において予想される波動の反射・屈折は考慮されていない。ところが、不整形地盤では、波動の屈折・反射の影響が無視できないので、隣接するメッシュ領域の伝達関数の影響を受ける。すなわち、従来の地震動予測は隣接する地盤メッシュ間のエネルギー収支の相互作用が考慮されていないため不整形地盤地域の予測には精度的な問題がある。
【0007】
これに対して、上記の不整形地盤を含む領域において正確な地震動伝達関数を得るには、有限要素法により3次元または2次元的モデル化による擬似3次元での解析をする必要がある。しかし、この解析に用いることのできる地盤情報は限られていることと、有限要素解析を行うことに莫大なコストがかかるため、研究目的以外での適用は実用上不可能である。
【0008】
一方、上記特許文献1及び2のような補間手法では、地盤データの採取地点間で地震動パラメータの補間をしたり、地盤メッシュと計測震度との関係に応じて地震動の補間を行ったりしているため、いずれも不整形地盤の分布状況とは無関係に、複数の地震動パラメータを演算的に処理するものに過ぎないことから、予測時の演算処理の効率と予測精度とを高次元で両立することができないという問題点がある。
【0009】
そこで、本発明は上記問題点を解決するものであり、その課題は、不整形地盤による地震動伝達関数の隣接領域間の相互作用を考慮することにより、予測時の処理効率と予測精度を共に向上させることの可能な地震動予測方法を提供し、また、この地震動予測方法を用いて高精度の地震動ハザードマップを効率的に制作する方法を実現することにある。
【課題を解決するための手段】
【0010】
斯かる実情に鑑み、本発明の地震動予測方法は、表層地盤の地質特性に関する地盤データから地震動特性を予測する地震動予測方法において、前記地盤データから当該地盤データの採取地点若しくは該採取地点に対応して設けられた別の地点である特定地点における特定地震動伝達関数を複数の前記特定地点についてそれぞれ算出する特定地震動算出工程と、隣接する2つの前記特定地点間における地盤構成の高低差に対応する不整形性を表す不整形地盤指標を設定する不整形指標設定工程と、前記隣接する2つの特定地点における前記特定地震動伝達関数と、前記隣接する2つの特定地点間に対応する前記不整形地盤指標とに基づいて、前記隣接する2つの特定地点を通過する直線上の複数の予測地点における推定地震動伝達関数若しくは該推定地震動伝達関数から得られる予測地震動特性値を算出する推定地震動算出工程と、を具備し、共通の前記特定地点について前記不整形指標設定工程及び前記推定地震動算出工程を異なる2方向に沿って実行して各方向に対応する前記予測地点についてそれぞれ前記推定地震動伝達関数若しくは前記予測地震動特性値を算出することにより、前記推定地震動伝達関数若しくは前記予測地震動特性値の平面分布を求めることを特徴とする。
【0011】
本発明によれば、隣接する2つの特定地点における特定地震動伝達関数及びこれに対応する不整形地盤指標に基づいて2つの特定地点を通過する直線上の複数の予測地点において推定地震動伝達関数若しくは予測地震動特性値を算出することで、不整形地盤指標によって上記直線に沿った方向の二次元的な影響、すなわち地盤構成の高低差による地震波の屈折や反射の影響を考慮した、地震動に関する地盤特性を簡易に得ることができる。また、或る特定地点について異なる2方向に沿った複数の予測地点(予測地点群)について推定地震動伝達関数若しくは予測地震動特性値を算出することで、それらの平面分布を求めることができるので、地盤構造の平面的な分布に対応する疑似三次元的な地震動に関する地盤特性を得ることができる。なお、隣接する2つの特定地点を通過する直線上の複数の予測地点において推定地震動伝達関数若しくは予測地震動特性値を算出する推定地震動算出工程には、特定地点間の予測地点において補間推定する場合に限らず、特定地点間の外側にある予測地点において外挿推定を行う場合も含む。
【0012】
本発明において、前記不整形地盤指標は、前記隣接する2つの特定地点間の距離と、地表面若しくは工学的基盤の高低差とに対応する直線状の傾斜を備えた仮想傾斜構造を特定する設定値であることが好ましい。これによれば、隣接する2つの特定地点の地盤データのみで地盤の不整形の影響を取り込むことができるとともに、仮想傾斜構造が特定地点間の距離と地表面若しくは工学的基盤の高低差(レベル差)に対応する直線状の傾斜を備えたものとされるため、簡易な傾斜構造であることから、推定地震動伝達関数の算出処理を容易に行うことが可能になる。ここで、上記不整形地盤指標の設定値としては、上記仮想傾斜構造の傾斜方向の長さと傾斜角の他に、同長さと高低差などといった、仮想傾斜構造を特定できる実質的に等価なものであれば種々の設定値を用いることができる。
【0013】
本発明において、前記平面分布は前記予測地震動特性値の平面分布であり、前記予測地震動特性値は、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期またはSI値のいずれかであることが好ましい。予測地震動特性値としては、予測地震動(波形)の特徴を何らかの意味で反映する各種の地震動パラメータを用いることができるが、具体的で一般的なパラメータとしては上記のものが挙げられる。
【0014】
本発明において、所定の入力地震動に基づいて前記推定地震動伝達関数により前記予測地点における予測地震動波形をさらに求めることが好ましい。これによれば、予測地震動波形、例えば数値加速度時刻歴など、を求めることにより、具体的な地震動に基づいて建物の耐震対策などをより確実かつ正確に策定することが可能になる。
【0015】
本発明において、前記特定地点は平面的に配列された複数の地盤メッシュの代表地点であり、前記不整形指標設定工程では、隣接する2つの前記地盤メッシュにおける前記特定地点の間の前記不整形地盤指標が算出され、前記推定地震動算出工程では、前記隣接する2つの地盤メッシュにおける前記特定地点を通過する直線上の、前記特定地点間の距離よりも小さな間隔で配列された複数の前記予測地点における前記推定地震動伝達関数若しくは前記予測地震動特性値が算出されることが好ましい。これによれば、地盤メッシュごとに存在する地盤データを元にその代表地点における特定地震動伝達関数を求めるとともに、隣接する地盤メッシュの代表地点を通過する直線上において代表地点間の距離よりも小さな間隔で配列された複数の予測地点を設定することにより、地盤メッシュ単位よりも細密に、地震動に関する地盤特性の平面分布を得ることができる。
【0016】
本発明の地震動ハザードマップの制作方法は、上記の地震動予測方法を用いて、所定の入力地震動に基づいて前記推定地震動伝達関数により得られる前記予測地震動特性値の平面分布を求め、該予測地震動特性値の平面分布をマップ上に表示することが好ましい。これによれば、予測地震動特性値の平面分布をマップ上に表示することにより、地震動に関する地盤特性の平面分布の態様を地盤メッシュに制約を受けずに細密に示すことが可能になる。
【発明の効果】
【0017】
本発明によれば、不整形地盤による地震動の屈折・反射の影響を反映させることができるとともに、有限要素法のような高いコストを要することなく、簡易に地震動に関する地盤特性の平面分布を得ることができるという優れた効果を奏し得る。特に、ハザードマップやコンターマップなどのマップ上に推定地震動伝達関数若しくは予測地震動特性の平面分布を表示することにより、滑らかな地震動に対する地盤特性の表示態様が実現できる。また、地盤メッシュの地盤データに基づいて推定地震動伝達関数若しくは予測地震動特性の平面分布を求めることにより、地盤メッシュの制約を受けずに、不整形地盤に対応した細密なものとした、地震動に関する地盤特性の平面分布を簡易に得ることができる。
【図面の簡単な説明】
【0018】
【図1】本実施形態において考慮する不整形地盤をタイプごとに模式的に示す地盤モデル構造図であり、(a)は地表面に高低差(傾斜)を有するタイプAの地盤構造に対するもの、(b)は工学的基盤に高低差(傾斜)を有するタイプBの地盤構造に対するもの。
【図2】傾斜基盤を備えた不整形地盤における、(a)ポイントA(図1(b)に示す傾斜上端地点)、(b)ポイントB(同、傾斜中央地点)、(c)ポイントC(同、傾斜下端地点)における地震動伝達関数を示す周波数特性グラフであり、各グラフにおいて、各ポイントにおける、符号1Dを付した線は1次元の重複反射理論による等価線形解析による解析結果、符号2Dを付した線は2次元FEM解析による解析結果、太い2本の灰色線AとCはそれぞれ比較のために示したポイントAとポイントCにおける1次元の重複反射理論による等価線形解析による解析結果である。
【図3】1次元の地震動伝達関数ΩTとΩBを用いた簡易推定手法を説明するための説明図。
【図4】重み係数CTのXLS(=x/Ls)依存性を傾斜角ごとに示すグラフ(a)〜(d)。
【図5】下部に示す二次元有限要素法を適用した分割要素を示す図(分割ブロック構造図)に示すポイントA、ポイントB及びポイントCにおける地震動伝達関数の周波数特性を示すグラフであり、符号1Dを付した線は1次元の重複反射理論による等価線形解析による解析結果、符号2Dを付した線は2次元FEM解析による解析結果、符号SEを付した線は本実施形態に用いた簡易推定法による解析結果である。
【図6】ポイントA、ポイントB及びポイントCにおけるベクトル軌跡を示す図(ナイキスト線図)。
【図7】縦横に隣接する4つの地盤メッシュ1〜4の解析結果としての特定震度と模式化した地盤構造とを対比して4組の隣接する2つの地盤メッシュの対をそれぞれ示す説明図。
【図8】縦横に隣接する4つの地盤メッシュ1〜4の特定震度を濃淡で示すとともに各地盤メッシュの中心点である代表地点1a〜4aを示す説明図。
【図9】4組の隣接する2つの地盤メッシュの対ごとに、2つの代表地点を通過する直線上の予測震度を濃淡と正方形ブロックの高低で示す図。
【図10】東西方向(図示左右方向)に沿って算出した予測震度を同方向に配列された2つの代表地点間の範囲に限定して濃淡で示す左図と、南北方向(図示上下方向)に沿って算出した予測震度を同方向に配列された2つの代表地点間の範囲に限定して濃淡で示す右図とを対比して示す図。
【図11】各地盤ブロックの4つの代表地点で囲まれる正方形の範囲内で東西方向に沿って算出した予測震度と南北方向に沿って算出した予測震度とを重ね合わせた結果を濃淡で表示した図。
【図12】本実施形態の対象地区として選定した範囲を5箇所の区域とともに示すマップ。
【図13】本実施形態による図12に示す対象地区の解析結果である予測震度の平面分布を濃淡で示す図。
【図14】従来の1次元の等価線形化手法による地盤メッシュごとの特定震度の平面分布を濃淡で示す図。
【図15】図13に示す解析結果を、図12に示す三輪の区域について拡大して示す図。
【図16】図14に示す解析結果を、図12に示す三輪の区域について拡大して示す図。
【図17】図13に示す解析結果を、図12に示す富竹の区域について拡大して示す図。
【図18】図14に示す解析結果を、図12に示す富竹の区域について拡大して示す図。
【図19】図13に示す解析結果を、図12に示す安茂里の区域について拡大して示す図。
【図20】図14に示す解析結果を、図12に示す安茂里の区域について拡大して示す図。
【図21】図13に示す解析結果を、図12に示す篠ノ井の区域について拡大して示す図。
【図22】図14に示す解析結果を、図12に示す篠ノ井の区域について拡大して示す図。
【図23】図13に示す解析結果を、図12に示す松代の区域について拡大して示す図。
【図24】図14に示す解析結果を、図12に示す松代の区域について拡大して示す図。
【発明を実施するための形態】
【0019】
次に、添付図面を参照して本発明の実施形態について詳細に説明する。最初に、本実施形態において特定地震動伝達関数や特定地震動特性値の算出に用いられる算出方法の例として、等価線形化手法(SHAKE)及び周波数依存型等価線形化手法(FDEL)について説明する。
【0020】
地盤震動解析法では、重複反射理論を用いることで、複数の地層における地震動の解析を行うことができる。せん断波速度Vsの高い工学的基盤(一般的には300〜700m/s、本実施形態では500m/s以上を基盤層としている。)上に複数の地層(jは1〜nの自然数)が存在する場合には、工学的基盤からSH波の一次元的な上昇伝達について地層境界での反射・屈折を考慮することにより、SH波の変位波形の基礎方程式から、ひずみは以下の数式1により表される。ここで、gjは数式2で表される。なお、γjは地表からj番目の第j層のひずみ、ujは第j層での地震波の振幅、tは時間、Ajは第j層の上昇波の変位振幅、Bjは第j層の下降波の変位振幅、Gjは第j層のせん断剛性、zjは第j層の厚さ、ωは円振動数、ρjは第j層の密度、ξj(=2Gjhj/ω、hjは減衰定数)は第j層の粘性係数である。このとき、(1)地表面でせん断応力が0、(2)地層の境界面でせん断応力は連続、(3)地層の境界面で変位が連続の3つの境界条件を用いると、地表層についてはA1=B1が成立し、これによって、Ajは以下の数式3、Bjは以下の数式4の漸化式によって表される。
【0021】
【数1】
【0022】
【数2】
【0023】
【数3】
【0024】
【数4】
【0025】
この場合、地表の加速度に対する第j層の上面の加速度波の応答倍率である伝達関数Tfjは、以下の数式5で表される。
【0026】
【数5】
【0027】
第j層が工学基盤上の基盤層とすれば、工学基盤への入力波は基盤層の上昇波と考えることができるので、第j層の応答倍率(伝達関数)T*fjは以下の数式6で表される。ここで、a*jは上昇波だけの加速度波である。また、解放基盤波形を求めるときは、応答倍率は上昇波を2倍して以下の数式7とする。
【0028】
【数6】
【0029】
【数7】
【0030】
地表の地震動をフーリエ変換し、周波数領域で上記数式5のTfj又は上記数式6又は7のT*fjを乗じた後にフーリエ逆変換することにより、任意の層の加速度波形を得ることができる。また、同周波数領域で上記数式1のγjを乗じることで、任意の層のひずみ波形を得ることができる。
【0031】
一方、地盤の地震応答解析において用いられてきた等価線形化手法の上記SHAKEでは、以下の数式8に示すように、ひずみの最大値の或る一定の割合を等価平均ひずみγeとし、これより定まるせん断剛性と減衰定数を周波数領域の計算において周波数に依存することなく一律に適用している。ここで、Cは係数、γmaxはひずみの時刻歴波形の最大値である。上記SHAKEでは係数Cを0.65としている。
【0032】
【数8】
【0033】
上記のように、等価線形化手法(SHAKE)でも或る地点の地盤データから一次元的な地震動解析を行うことができる。しかしながら、地震動のレベルが上がりひずみが大きくなると地盤のせん断剛性と減衰定数は非線形に変動し、特に低・中周波数領域では実際の観測記録と比較的よく一致するが、強震動や軟弱地盤の応答解析を行う場合、高周波数領域で計算された結果が観測値を下回るなど、地盤の非線形性の影響による修正を考慮する必要がある。このため、周波数依存型等価線形化手法(FDEL)では、各周波数ごとにひずみ波形に寄与する度合に応じて適切なせん断剛性と減衰定数を与える周波数依存型の等価平均ひずみγf(ω)を以下の数式9のように定義する。ここで、Fγ(ω)はひずみ波形のフーリエスペクトル、Fγmaxはひずみ波形のフーリエスペクトルの最大値である。なお、C=0.65、Fγ(ω)/Fγmax=1.0とすると上記SHAKEと同じになる。係数Cの値はひずみ波形のフーリエスペクトルの卓越周波数周辺の伝達関数に大きく影響し、種々の検討によりC=0.65程度で問題がないことが確認されている。ただし、この係数Cの値は実情に合わせて調整可能である。
【0034】
【数9】
【0035】
ここで、せん断ひずみγとせん断剛性Gとの関係、並びに、せん断ひずみγと減衰定数hとの関係を用いることにより、ひずみに依存する非線形性を考慮した応答を得ることができる。非線形な応力とひずみの関係としては、例えば、拘束圧の影響を考慮してモデル化されたHardin-Drnevichモデルを用いることができる。このモデルは、地層の土質別に種々の拘束圧におけるせん断強度に対する実験に基づいて、各土質ごとに基準化ひずみγγを定義し、ひずみとせん断剛性、ひずみと減衰定数の関係を求めたものである。ここで、数式9に示す周波数依存型の平均ひずみγf(ω)を用いる場合には、以下の数式10によりせん断剛性G(ω)が表され、以下の数式11により減衰定数h(ω)が表され、いずれも周波数の関数になる。ここで、Gmaxは初期せん断剛性、hmaxは最大減衰定数、γγは基準化ひずみ(=τmax/Gmax、τmaxはせん断応力の最大値)である。
【0036】
【数10】
【0037】
【数11】
【0038】
上記の等価線形化手法(SHAKE)や周波数依存型等価線形化手法(FDEL)では、一次元の地盤震動解析により地震動伝達関数を算出することができる。一方、二次元有限要素法(プログラム名:FRUSH)や、これを改良した周波数依存型等価線形化法による二次元FEM地盤震動解析法(プログラム名:FDEL−FEM)も使用されているが、これらは膨大な計算処理を必要とする。そこで、本実施形態では、上述の種々の一次元地盤震動解析により求められた地震動伝達関数に基づいて、不整形地盤の簡易推定法を用いる。まず、この簡易推定法と二次元有限要素法による解析結果とを比較し、簡易推定法の有効性を示す。
【0039】
図1は、簡易推定法と二次元有限要素法の解析結果の比較に用いるモデル地盤として、工学的基盤が平坦で盛土などにより地表面が傾斜している不整形地盤のモデル(a)と、地表面が平坦で工学的基盤が傾斜している不整形地盤のモデル(b)をそれぞれ模式的に示す概略構成図である。ここで、領域1(Region-1)は表層地盤構造であり、領域2(Region-2)は工学的基盤である。二次元FEM地盤震動解析では、上記数式9に示す周波数依存型等価ひずみγf(ω)を導入し、C=0.65とした。また、モデル地盤の形状及び材料物性は、以下の表1及び表2に示すものとした。
【0040】
【表1】
【0041】
【表2】
【0042】
また、地震動伝達関数は地震動の入力強度に依存することから、マグニチュードMと震源距離Rを仮定し、入力レベルが異なるシミュレーション地震動を算出して解析に用いた。入力地震動の諸元を表3に示す。なお、表2において、領域1のせん断波速度(単位m/s)については、二種類(100と300)の数値でそれぞれ計算を行い、解析を行っている。
【0043】
【表3】
【0044】
上記タイプAの不整形地盤における解析結果の一例(Case−3,傾斜角1:5,Vs=300m/sec,入力地震動No.Acc−3)を図2に示す。ポイントAは図1に示す傾斜上端の地点、ポイントCは図1に示す傾斜下端の地点、ポイントBは着目地点である傾斜上の地点である。図中の2Dで示す実線は二次元FEM解析(2D)による伝達関数(地震動の出力スペクトル)、1Dで示す実線は着目地点であるポイントBと同じ深さの成層地盤を仮定して1次元の重複反射理論による等価線形解析(1D)を行った場合の結果である。また、指示線を有しない2本の灰色の太線は傾斜上端のポイントA及び下端のポイントCでの1Dによる伝達関数を示し、それぞれ傾斜上端のグラフをA、傾斜下端のグラフをCとした。
【0045】
傾斜上端の地表面であるポイントAでは、不整形性の影響が比較的少なく、2Dと1Dによる地震動伝達関数は、概ね一致している。傾斜中央のポイントBでは、地盤の不整形性の影響を受け、2Dによる伝達関数の卓越振動数は、1Dによる伝達関数の固有振動数と異なった位置に卓越振動数を持つ。また、傾斜下端のポイントCにおいてもポイントBと同様に不整形性の影響が確認できる。また、上述の結果は、タイプBのCase−1とCase−2についても共通していることが確認された。なお、2Dによる伝達関数は、傾斜上端および下端の1Dによる伝達関数の固有震動数の両方を卓越振動数として含むことから、傾斜部における地表面における伝達関数は、傾斜上端・下端における伝達関数(固有振動数)と密接に関係しているといえる。すなわち、地震動伝達関数は、傾斜部が存在することで傾斜方向に沿った他の地点の地盤構造に影響を受けることがわかる。
【0046】
不整形構造を有する地盤の傾斜部の地震動伝達関数の卓越震動数は、傾斜上端・下端の固有振動数の双方の位置で卓越し、その振幅は傾斜角度や地点に依存していることから、傾斜部上端・下端での1Dによる伝達関数に重みを付け、それらを複素平面上で重ね合わせることで当該地点の伝達関数を補正推定する簡便法を用いることができる。なお、この伝達関数は上述と同様に、SH波の鉛直入射に対する面内震動における入射地点から地表への水平震動成分の増幅率を表す。
【0047】
この簡易推定法では、以下の数式12に示すように、傾斜上端での1Dによる伝達関数ΩTと傾斜下端での1Dによる伝達関数ΩBの重み付き相乗平均として重ね合わせ、図3に示す任意地点(着目地点、傾斜上端からの距離xで規定される。)の伝達関数ΩEを推定する。伝達関数の掛算は複素平面上での回転操作となるため、振幅と位相角を正しく補間できる。一方、足し算を行うと、位相角が相反する場合、位相角と振幅が相殺される場合がある。なお、CTはΩTに対する重み係数、CBはΩBに対する重み係数である。このとき、伝達関数はいずれも複素数で与えられるものとし、重み係数は実数で与えられるものとする。
【0048】
【数12】
【0049】
上記数式12は、x/Ls=+∞のときに重み係数CTが0で右側の自由地盤の伝達関数ΩB(ω)、x/Ls=−∞のときに重み係数CTが1で左側の自由地盤の伝達関数ΩT(ω)とそれぞれ一致し、傾斜面上端からの距離xの変動により滑らかに変化することから、重み係数CTを以下の数式13で表されるロジスティック関数でモデル化する。
【0050】
【数13】
【0051】
傾斜面上端からの距離xと傾斜部の長さLsとの比(x/Ls)と傾斜角θをパラメータとして、数式12におけるΩEと有限要素法により求めたΩ2Dとの振幅の残差が固有振動数で最小となるように回帰分析により決定する。この決定方法の一例としては、例えば、数式13を変形すると数式14となり、この数式14を用いてb0及びb1を得ることができる。その結果を数式15に示す。なお、R′2は自由度調整済みの決定係数である。ただし、本発明では、重み係数の算出方法は上記の数式13乃至15に示す方法に限らない。例えば、上記のように傾斜構造の設定値(θやLs)だけでなく、表層地盤の物性(例えば、せん断波速度、ポアソン比、地盤堆積密度など)、表層地盤の地盤構成(地層数、土質、深さなど)、入力地震動の強度のいずれか一つ若しくは二以上を、陽に変数として表れる形で重み係数の数式13中に組み込んで算出し、数式12を計算してもよい。
【0052】
【数14】
【0053】
【数15】
【0054】
ここで、簡易推定式の構築に際して行った解析条件より、傾斜角の範囲を1:2〜1:5、傾斜部分の水平長さLsが500m程度以下をもって適用範囲とする。重み係数CTの標本値とモデル式を示す曲線(実線)を傾斜角ごとに図4に示す。ここで、図4(a)〜(d)は、表1に示す5つの全てのケースについて、それぞれ同じ傾斜角のものを選んで基礎となる地盤モデルとし、表2の地盤材料と表3の地震入力をそれぞれ適用して得た複数のFEM計算結果から算出した重み係数CTの標本値(薄い網目線)と、上述の方法で係数を決定した重み係数のモデル式(実線)のXLS=x/Lsに対する依存性(分布)を示すものである。図4の各図をみると、モデル式は各傾斜角の分布を表現できていることが確認できる。すなわち、タイプAとタイプBのいずれにおいても傾斜角θごとに上記の共通のモデル式を用いることが可能であることが示されている。
【0055】
簡易推定式を用いた伝達関数と二次元FEM解析による伝達関数の比較を図5に示す。対象地点を傾斜面上端のポイントA、傾斜面中央のポイントB、傾斜面下端のポイントCとする。図5の下図は二次元FEM解析に用いる地盤の分割要素の構成を示す図である。図中のSEで示す実線は簡易推定法により算出した伝達関数、2Dで示す実線は上記2D、1Dで示す実線は上記1Dによる当該地点の伝達関数である。ポイントAにおいて伝達関数の振幅の2Dによる結果と簡易推定法を適用して推定した結果を比較すると、簡易推定法を適用して推定した結果は2Dの伝達関数をほぼ再現できている。ポイントBでは、2Dと1Dの結果が異なることは前節記述した通りであり、推定した伝達関数はポイントAと同様に振幅、位相角の両方について2Dの伝達関数をほぼ再現できている。ポイントCに関しても同様の結果となり、簡易推定式は概ね2Dによる結果を再現できている。
【0056】
図6にはそれぞれのベクトル軌跡(ナイキスト線図)を表す。ベクトル軌跡は、複素平面上に伝達関数を表示したもので、ここでは周波数0Hzから3Hzまでの軌跡を表示した。図5と同様に、簡易推定法(SE)により、有限要素法(2D)の振幅(原点からの距離)をほぼ再現できるのを表すとともに、位相角(実軸と虚軸を基準とした角度)もほぼ同様に推定できていることを示している。
【0057】
ここで、タイプAの不整形地盤については、工学的基盤が同レベルであり、地表面に傾斜が存在するため、地表面の傾斜角を上記数式13のθとすることにより、中間地点xにおける伝達関数ΩEを求めることができる。また、タイプBの不整形地盤については、工学的基盤が傾斜し、地表面が同レベルであるため、工学的基盤の傾斜角を上記数式13のθとすることにより、中間地点xにおける伝達関数ΩEを求めることができる。
【0058】
また、地表面と工学的基盤のいずれにも傾斜がない場合(すなわちθ=0の場合、以下、タイプCという。)には、傾斜による波動の屈折や反射の影響がないと考えられるので、CT=1−x/Ls(0<x<Ls)、CB=1−CT(0<CB,CT<1)を用いて、重み係数CT、CBを線形補間により求めて上記数式12により中間点の波動関数ΩEを求めることができる。
【0059】
なお、本実施形態では、タイプA及びBについては、いずれも数式13によって数式12の重み係数CT及びCBを推定し、タイプCでは重み係数CT及びCBを線形補間して数式12を適用しているが、本発明はこのような態様に限られるものではなく、後述するように地盤タイプの少なくとも一部において不整形地盤指標に基づいて推定が行われさえすればよい。また、本実施形態では、タイプAとBとにおいて、傾斜角θと距離Lsが等しければ同じ方法(数式12及び13)で推定するようにしているので、例えば、地表面が一方側に傾斜角θで傾斜し工学的基盤が平坦なタイプAの地盤モデルは、地表面が平坦で工学的基盤が反対側に傾斜角θで傾斜しているタイプBと実質的に同じ推定処理を行うことになる。したがって、地表面と工学的基盤が共に傾斜し、どちらの傾斜も無視できない場合には、工学的基盤を基準としたときの地表面の高低差によって地表面の傾斜角を求め、工学的基盤が平坦なタイプAの地盤モデルとして計算するか、或いは、地表面を基準としたときの工学的基盤の高低差によって工学的基盤の傾斜角を求め、地表面が平坦なタイプBの地盤モデルとみなして計算すればよい。これによって全ての地盤構造をタイプA〜Cのいずれかの地盤モデルに対応させることができる。ただし、本実施形態の上述の方法に拘わらず、本発明ではタイプAとタイプBで相互に異なる方法(数式)で推定するようにしても構わない。この場合、本実施形態と同様に上述のようなタイプAとタイプBの間で変換可能な関係をもたせてもよいが、例えば地表面の傾斜角と工学的基盤の傾斜角のうちの傾斜角の大きい方の傾斜を残すように、地盤構造を何らかの指標に応じてタイプAとタイプBのいずれかの地盤モデルに分類するようにしてもよい。さらに、上述のように地表面と工学的基盤が共に傾斜し、どちらの傾斜も無視できない場合をタイプDとして別途の方法(数式)で推定するようにしてもよい。
【0060】
本実施形態では、一次元等価線形化手法等により得られた上記地震動伝達関数ΩT及びΩBに相当する特定地震動伝達関数と、上記Ls及びθに相当する不整形地盤指標(上記傾斜部を特定する設定値)とに基づいて、上記の簡易補間推定法を用いて、複数のメッシュ間で上記地震動伝達関数ΩEに相当する推定地震動伝達関数を求め、疑似3次元的な地盤構造の解析を行う。
【0061】
この場合、簡易推定法を適用する2つの特定地点の地盤データの関係に応じて、上記タイプA、タイプB及びタイプCのいずれか一つを選択するとともに、選択されたタイプごとに上記不整形地盤指標を設定し、これらに応じて重み係数CT及びCBを上述のタイプ別にそれぞれに求めることで、上記数式12により推定地震動伝達関数ΩEを求めることができる。なお、本実施形態では、不整形地盤指標は地盤の傾斜(高低差)に関するもののみであるので、傾斜のないタイプCについては実質的に不整形地盤指標が適用されずに推定地震動伝達関数ΩBが求められるが、この場合には、不整形地盤指標のうち傾斜角θ=0の特別な場合において不整形地盤指標とは無関係に算出しているだけであり、傾斜角θを0とみなすことのできないタイプA及びBの場合において不整形地盤指標に基づいた算出が行われうる構成である。このように、全ての地盤において等しく不整形地盤指標を用いる必要はなく、必要な一部の場合にのみ不整形地盤指標を用いるものであっても構わない。
【0062】
次に、この解析方法の原理を説明するために、図7に示すように隣接する4つの地盤メッシュを基にして解析を行う例を示す。図7に示す4つの地盤メッシュ1〜4の地盤構造は、各地盤メッシュ1〜4の代表地点(図示例では各地盤メッシュの中心点)1a〜4aにおける地盤データに対応する地盤構造(工学的基盤上の表層地盤の地層構成)を示す。この場合、解析には、従来法で使用していた地盤メッシュの地盤データをそのまま使用することができる。ここで、工学的基盤上の表層地盤が複数の地層よりなる場合には、上記重複反射理論に基づいて伝達関数を算出する。なお、各地盤メッシュでは、地盤データを採取した採取地点と、上記代表地点(特定地点)とが一致しない場合もあり得る。例えば、地盤データの採取地点が中心点以外にしかない場合、複数の採取地点がある場合であって複数の採取地点の寄与を含めた地盤データを用いる場合などでは、各地盤メッシュの中心点などの代表地点(特定地点)における地盤データは、実際には中心点などの代表地点(特定地点)とは異なる採取地点から取得した地盤データである。この場合、使用する地盤データの採取地点を地盤メッシュの代表地点(特定地点)としてもよく、或いは、採取地点から取得した地盤データを代表地点(特定地点)における地盤データと見做しても構わない。
【0063】
まず、基にする地盤データであるが、図7の下部に示されているような地盤構造が与えられているとし、この地盤構造に対応するせん断波速度、単位体積重量、土質、深さを用いる。ここで、土質とは、上述のHarding-Drnevichモデルで定義される基準化ひずみを与えるもので、具体的には、砂、粘土、礫、岩に対応する土質区分を言う。また、図7の下部の地盤構造の図において、最下部の二重線は工学的基盤との境界を示し、この二重線より上方部分は表層地盤の概略構成(1又は複数の地層の境界と厚み)を示す。この地盤構造では、実際には地表面レベルにはばらつきがあるものの、地表面レベルを基準として工学的基盤の高低差を図示しており、この場合には、上記のタイプBの地盤モデルを想定して計算を行うことになる。これを使い各地盤メッシュ1〜4での上記の一次元地震動解析を行って特定地震動伝達関数Ω1〜Ω4を求め、所定の共通の地震動(地震波スペクトル)を入力として工学的基盤に与えた場合について、図7の上部に示す震度が算出され(濃淡と数値で示してある。)、実際の震度予測図に用いられる平面上に落とすことで図8のようになる。ここで、濃淡と数値で示した各地盤メッシュの震度は、各地盤メッシュ1〜4の代表地点1a〜4aの震度に対応するものであるが、図示上では各地盤メッシュ全域に表している。この震度は、所定の入力地震動に基づいて上記特定地震動伝達関数から得られる特定地震動特性値である。なお、所定の入力地震動に基づいて地震動伝達関数から得られる地震動特性値としては、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期、SI値などが挙げられる。
【0064】
これらの4つの隣接する地盤メッシュ間で、各地盤メッシュ1〜4の代表地点(図示例では中央点)1a〜4aの地盤データに基づく特定地震動伝達関数Ω1〜Ω4と、それぞれの隣り合う2つの地盤メッシュの代表地点1aと2a、2aと3a、3aと4a、4aと1aの間の地表面や工学的基盤の傾斜に応じた上記と同様の不整形地盤指標で、2つの地盤メッシュ1と2、2と3、3と4、4と1の全体を(すなわち、隣接する2つの代表地点の間だけでなく、一方の代表地点に対して他方の代表地点とは反対側にある地点についても)補間・外挿推定し、これらの複数の予測地点について得られた推定地震動伝達関数ΩCに基づいてそれぞれ予測地震動特性値である上記共通の地震動を入力したときの予測震度を求めたものが図9である。上記の算出処理を行った隣接する2つの地盤メッシュの各組では、タイプA及びBの地盤モデルに対応する場合には上記数式12乃至15で示す方法で、タイプCの地盤モデルに対応する場合には上記の線形補間された重み係数による上記数式12により、東西方向(図8の左右方向)又は南北方向(図8の上下方向)をそれぞれ上記のx方向として、1つの地盤メッシュ当たりで等間隔の複数(図示例では10個)のx値(隣接する2つの地盤メッシュでは合計20個)の予測地点について、上記推定地震動伝達関数ΩCをそれぞれ求めた後、上記共通の入力地震動に基づいて各推定地震動伝達関数ΩCによる出力地震動に対応する予測震度をそれぞれのx値ごとに示してある。
【0065】
次に、東西方向・南北方向の2つの結果を地盤メッシュの代表地点(中心点)から隣の地盤メッシュの代表地点(中心点)までの範囲で抽出して示したものが図10である。ここで、図10の左側の図は東西方向の簡易補間推定結果であり、地盤メッシュ1と2、3と4でそれぞれ簡易補間推定を行った結果を東西方向の分布に反映させ、南北方向には一定の値としている。また、図10の右側の図は南北方向の簡易補間推定結果であり、地盤メッシュ1と4、2と3でそれぞれ簡易補間推定を行った結果を南北方向の分布に反映させ、東西方向には一定の値としている。
【0066】
次に、上述の東西方向・南北方向の結果を重ね合わせる。なお、この結果としては、図示のように予測震度などの予測地震動特性値を用いてもよく、予測震度などの予測地震動特性値を求める前の推定地震動伝達関数ΩCを用いてもよい。推定地震動伝達関数ΩCを用いた場合には重ね合わせた後に上記共通の地震動を入力することで予測震度(予測地震動特性値)を求めることができる。この時、東西方向・南北方向では予測震度に大きな差の出る場合がある。このため、地震波の伝達を考え、両方向で得られた2つの結果(2方向で計算した結果)の平均値や加算値を用いるなど、両方向の推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を共に反映させた値とする必要がある。ただし、予測地点の位置に応じて、両方向の算出値の重み付け(寄与)を適宜に設定して統合値(合成値)に反映させても構わない。
【0067】
また、図10に示した処理では、東西方向の簡易補間推定では、東西方向の位置(x値)が同じであれば同一の地盤メッシュ中の南北方向に沿った全ての予測地点(図面上では帯状の範囲)において同じ推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、また、南北方向の位置(x値)が同じであれば同一の地盤メッシュ中の東西方向に沿った全ての予測地点(図面上では帯状の範囲)において同じ推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、これらを平面的に重ね合わせることで合成している。しかし、東西方向と南北方向の簡易補間推定において、それぞれ対応する2つの代表地点を通過する直線上にある予測地点のみが算出された推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、両方向の直線上の予測地点が重ならないように設定して、これらの予測地点ごとに対応する予測地震動特性値がそのまま平面分布を構成するようにしてもよい。
【0068】
なお、地震動伝達関数から震度などの地震動特性値を算出するためには、入力地震動を決める必要があるが、本実施例の場合には、特定地震動特性値と予測地震動特性値のいずれについても、マグニチュード7.3相当の直下型地震を想定し、予測地区(長野市内)が震源からの距離20〜30km以内の震源域にあると仮定し、予測地区全体にわたり同一の基盤入力があったとして、同一の入力地震動を用いて予測震度を求めている。ただし、仮定された震源位置からの距離に応じた地震動の減衰を考慮し、上記特許文献5に示す強震動予測方法(EMPR)を用いて、地盤メッシュごとに、或いは、複数の地盤メッシュを含むより大きな区域ごとに工学的基盤へ入力される地震動波形を計算し、地盤メッシュ若しくは上記区域ごとに入力地震動を設定してもよい。本明細書では、上記のいずれの場合でも、共通の入力地震動に基づいて予測地震動特性値を求めているものとする。
【0069】
以上の過程から、図11の4つの地盤メッシュの中心1a〜4aを結んだ四角形の中の予測震度の表示(平面分布の作成)がなされる。この方法を対象地域内で隣接する地盤メッシュ全てで行い繋げていくことにより、図12に示す予測地区全域における予測震度を図13に示すように表示することができる。図13に示す簡易補間推定を行った結果は、隣接する各代表地点間、或いは、1つの地盤メッシュの東西方向及び南北方向の一辺当たり10個の予測地点を設定し、その結果を同一辺あたりで東西方向・南北方向それぞれに10段階で表示している。そのため、図14に示す従来法(1D)による特定地震動特性値のみを表示した場合に比べて1地盤メッシュあたり10×10で100倍の細かさで表示することができる。
【0070】
上記結果の詳細な比較・検討を行うため、図12に示す各区域の推定結果を示す。これらの各区域は以下の表4に示す5つの区域である。また、各区域における本実施形態による簡易補間推定結果及び従来法の表示を、図15と図16(三輪)、図17と図18(富竹)、図19と図20(安茂里)、図21と図22(篠ノ井)、図23と図24(松代)にそれぞれ対比して示す。なお、これらの各区域は一辺が3kmの正方形状で、各区域内を縦横それぞれ6つに分割してなる図中のメッシュサイズは500m×500mである。
【0071】
【表4】
【0072】
まず、図13の長野市全体での予測震度の表示結果(平面分布)を図14の従来手法と比較して述べる。本実施形態の方法により得られた地盤メッシュの境界における予測震度の表示は急激な震度の差が小さくなっていることが見て取れる。これは、予測震度が大きい場所から小さい場所にかけて段階的に変化している様子から、簡易補間推定法による正しい結果が表れているものと思われる。段階的に変化している部分では、従来法に比べ表示されている予測震度が大きくなっている場所がある。
【0073】
この長野市全体における結果を詳しく調べるため、次に比較対象地においての従来法との比較結果を見ていく。従来法では、地盤メッシュ内に山が入っている地点において山における地盤データが基になっているため、予測震度が小さく表示されている(例、安茂里2−5、松代4−4、※それぞれ左から数えたメッシュ数−下から数えたメッシュ数を示す。以下同様。)。それにより、例えば松代の例では上の地盤メッシュとの予測震度の差が大きくなり、境界や山の麓での予測震度が不明確であった。しかし、簡易補間推定法により予測震度の差は小さくなり、実際に起こりうる震度に近づけることが出来た。さらに、地盤メッシュに関係なく震度表示できるため、地盤メッシュ内での震度の変化が表され、山の形に合った震度表示がなされている。地盤メッシュの境界付近では地盤メッシュの移り変わりに伴い段階的に変化している様子が表されている(例、三輪5−2と5−3の境界、篠ノ井4−2と4−3の境界)。これは地盤の不整形性を考慮した結果であり、エネルギー収支の関係性が影響していることが示されている。
【0074】
どの地点においても隣り合う地盤メッシュ間でおきるエネルギーの収支の状態が表れており、地盤メッシュ境界での予測震度の大きな差がなくなった。これにより、地盤メッシュ境界付近の予測震度が小さかった地域でも安全側の予測結果を表示できる。各図において、対象地域の外縁で震度表示がなされていないが、これは簡易補間推定を行っているためである。簡易補間推定は隣り合う地盤メッシュのデータを基にメッシュ間の震度表示を行うため、隣り合う地盤メッシュのデータが存在しないと外縁の震度表示はされない。対象地域内の震度表示をするには対象地域に接するメッシュデータも必要となる。
【0075】
地盤の不整形性を考慮することにより、今まで予測震度の小さかった地域においても危険にさらされる可能性のあることが確認できた。また、従来法では地盤メッシュの境界で大きく震度が変化していたが、地盤メッシュに関係なく予測震度の表示が出来るようになったため、段階的に変化している状態が表されている。
【0076】
以上より、地盤の不整形性を考慮した長野市の詳細な震度予測図が作成できたと言える。今回、予測震度の大きくなる場所が多数見られ、このような場所に主要交通網や盛土斜面がある場合、震災時の予測被害は大きくなる。従来法による震度予測図では想定されなかった被害も想定できるようになったことから、地震対策を進める上で必要な情報として貢献できるものとなった。本実施形態で行った簡易補間推定法による解析は、1次元解析を行っていた今までの地盤データを使用することができるため、作成したプログラムを用いることで自治体レベルの広域での地震動予測を行える。本実施形態では入力地震動を対象地域内の工学的基盤面において同じ値を入力しているが、これを信濃川断層等の被害が予測される地震でのEMPR(地震動予測モデル)の結果を用いることで、実際に地震が発生したときの震度表示をすることが可能である。
【0077】
上記実施例では、推定地震動伝達関数ΩCから所定の入力地震動に基づいて予測地震動特性値(予測震度)を求めているが、さらに詳細な地震動予測を行うために、所定の入力地震動に基づいて推定地震動伝達関数ΩCから予測地点の地表面における予測地震動波形を求めるようにしてもよい。この地震動波形を求めることにより、より詳細な地震動の態様が判明するので、具体的な防災対策を策定する場合に有効である。この場合には、上記入力地震動として予想される震源域、規模、震動特性に応じたものを用いることが好ましい。
【0078】
尚、本発明の照明装置及び電気光学装置は、上述の図示例にのみ限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。例えば、上記実施例では、地盤メッシュをベースとして、隣接する地盤メッシュ間で簡易補間推定を実施しているが、例えば地盤データの採取地点間に上記方法を適用して複数の予測地点における算出値の平面分布を求めてもよいなど、本発明は地盤メッシュとは無関係に、隣接する特定地点間の簡易補間推定を行う場合を広く包含する。また、上記実施例では、隣接する2つの特定地点の特定地震動伝達関数に基づいて予測地点の推定地震動伝達関数を求め、この推定地震動伝達関数から予測地震動特性値である予測震度を導出する場合について説明したが、本発明は、地盤傾斜に関連付けられた不整形地盤指標に基づくものであれば、特定地震動特性値から予測地震動特性値を直接導出する場合をも包含する。さらに、上記実施形態では地表面と工学的基盤面の高低差に関する不整形地盤指標を取り扱っているが、表層地盤内に複数の地層が存在する場合、地表面と工学的基盤面以外の他の地層面の傾斜について考慮してもよい。
【符号の説明】
【0079】
ΩT…傾斜上端の地盤データに対応する自由端の地震動伝達関数
ΩB…傾斜下端の地盤データに対応する自由端の地震動伝達関数
ΩE…推定地震動伝達関数
Ls…傾斜部分の水平距離
x…傾斜上端からの水平距離
A…傾斜上端の地点
B…傾斜上の地点
C…傾斜下端の地点
θ…傾斜角
Ω1〜Ω4…特定地震動伝達関数
ΩC…推定地震動伝達関数
【技術分野】
【0001】
本発明は地震動予測方法及び地震動ハザードマップの制作方法に係り、特に、不整形地盤を備えた地域における地震動伝達関数の補間手法を用いて地震動を容易かつ的確に予測する方法に関する。
【背景技術】
【0002】
一般に、地震動予測は、震源から地表への震動の伝達メカニズムを想定して算出されるが、将来予測される地震による地域の被害予測のための震度ハザードマップを作成する場合には、表層地盤はボーリング調査(点測定)に基づき数100m四方の独立した成層地盤メッシュとして扱われる。この場合、各ボーリング調査による地盤データを用いた一次元の地盤震動解析の分野においては、地震動伝達関数の計算方法として、等価線形化法による地盤震動解析プログラムSHAKEが広く使用されてきた(以下の非特許文献1参照)。また、これを改良してやや軟弱な地盤にまで適用範囲を広げた、周波数依存型等価線形化法FDELも開発されている(以下の非特許文献2参照)。これらは計算に必要な地盤パラメータが比較的容易に得られることから、実用性が高く評価されている。
【0003】
一方、採取された地盤データから地震動パラメータを算出し、この地震動パラメータに基づき、クリギングの手法によって、任意時点又は任意領域において補間された補間地震動パラメータ及びその補間精度を求める方法が知られている(以下の特許文献1参照)。また、以下の特許文献2においても地盤メッシュ間の地震動の補間手法が記載されている。
【先行技術文献】
【特許文献】
【0004】
【特許文献1】特開2005−156273号公報
【特許文献2】特開2003−287574号公報
【非特許文献】
【0005】
【非特許文献1】P.B.Schnabel, J.Lysmer and H.B.Seed: SHAKE a computer program for earthquake response analysis of holizontally layered sites, EERC, 72-12, 1972.(SHAKE)
【非特許文献2】杉戸真太, 合田尚義, 増田民夫: 周波数特性を考慮した等価ひずみによる地盤の地震応答解析法に関する一考察, 土木学会論文集, No.493/III-27, pp.49-58.1994.6.(FDEL)
【非特許文献3】Lysmer, J., Udaka, T., Tsai, C.-F. and Seed, H.B.: FLUSH a computer program for approximate 3D-analyisis of soil structure interaction problem, Proc. of EERC: 75-30. 1975(FLUSH)
【非特許文献4】Y.Furumoto, M.Sugito, A.Yashima:Frequency-Dependent Equivalent Linearized Technique for FEM Response Analysis of Ground, 12WCEE in Auckland, New Zealand (CD-ROM),2000.2(FDEL−FEM)
【非特許文献5】Sugito,M., Furumoto,Y., and Sugiyama,T., Strong Motion Prediction on Rock Surface by Superposed Evolutionary Spectra, 12th World Conference on Earthquake Engineering, CD-ROM, Auckland, New Zealand, January 2000.(強震動予測法(EMPR))
【発明の概要】
【発明が解決しようとする課題】
【0006】
しかしながら、既存の手法で作られる震度ハザードマップ等の分布図は、地盤メッシュ(例えば500m四方)で区切られているのが一般的であり、滑らかなコンターマップとはなっていない。これは、地盤のモデル化が地盤メッシュ内の高々一点を基準にして行われているためで、日本のように起伏の激しい地形(=不整形地盤)では、正確な震度マップの作成方法とはいえない。地表での震度を予測するためには、調査済みの地盤モデルから、基盤から地表までの地震動の伝達関数を計算する必要があるが、既存の手法は、地盤メッシュ内での一次元的な応答計算に留まっている。このため、地盤構造に急激な変動がある不整形地盤を有する地域において予想される波動の反射・屈折は考慮されていない。ところが、不整形地盤では、波動の屈折・反射の影響が無視できないので、隣接するメッシュ領域の伝達関数の影響を受ける。すなわち、従来の地震動予測は隣接する地盤メッシュ間のエネルギー収支の相互作用が考慮されていないため不整形地盤地域の予測には精度的な問題がある。
【0007】
これに対して、上記の不整形地盤を含む領域において正確な地震動伝達関数を得るには、有限要素法により3次元または2次元的モデル化による擬似3次元での解析をする必要がある。しかし、この解析に用いることのできる地盤情報は限られていることと、有限要素解析を行うことに莫大なコストがかかるため、研究目的以外での適用は実用上不可能である。
【0008】
一方、上記特許文献1及び2のような補間手法では、地盤データの採取地点間で地震動パラメータの補間をしたり、地盤メッシュと計測震度との関係に応じて地震動の補間を行ったりしているため、いずれも不整形地盤の分布状況とは無関係に、複数の地震動パラメータを演算的に処理するものに過ぎないことから、予測時の演算処理の効率と予測精度とを高次元で両立することができないという問題点がある。
【0009】
そこで、本発明は上記問題点を解決するものであり、その課題は、不整形地盤による地震動伝達関数の隣接領域間の相互作用を考慮することにより、予測時の処理効率と予測精度を共に向上させることの可能な地震動予測方法を提供し、また、この地震動予測方法を用いて高精度の地震動ハザードマップを効率的に制作する方法を実現することにある。
【課題を解決するための手段】
【0010】
斯かる実情に鑑み、本発明の地震動予測方法は、表層地盤の地質特性に関する地盤データから地震動特性を予測する地震動予測方法において、前記地盤データから当該地盤データの採取地点若しくは該採取地点に対応して設けられた別の地点である特定地点における特定地震動伝達関数を複数の前記特定地点についてそれぞれ算出する特定地震動算出工程と、隣接する2つの前記特定地点間における地盤構成の高低差に対応する不整形性を表す不整形地盤指標を設定する不整形指標設定工程と、前記隣接する2つの特定地点における前記特定地震動伝達関数と、前記隣接する2つの特定地点間に対応する前記不整形地盤指標とに基づいて、前記隣接する2つの特定地点を通過する直線上の複数の予測地点における推定地震動伝達関数若しくは該推定地震動伝達関数から得られる予測地震動特性値を算出する推定地震動算出工程と、を具備し、共通の前記特定地点について前記不整形指標設定工程及び前記推定地震動算出工程を異なる2方向に沿って実行して各方向に対応する前記予測地点についてそれぞれ前記推定地震動伝達関数若しくは前記予測地震動特性値を算出することにより、前記推定地震動伝達関数若しくは前記予測地震動特性値の平面分布を求めることを特徴とする。
【0011】
本発明によれば、隣接する2つの特定地点における特定地震動伝達関数及びこれに対応する不整形地盤指標に基づいて2つの特定地点を通過する直線上の複数の予測地点において推定地震動伝達関数若しくは予測地震動特性値を算出することで、不整形地盤指標によって上記直線に沿った方向の二次元的な影響、すなわち地盤構成の高低差による地震波の屈折や反射の影響を考慮した、地震動に関する地盤特性を簡易に得ることができる。また、或る特定地点について異なる2方向に沿った複数の予測地点(予測地点群)について推定地震動伝達関数若しくは予測地震動特性値を算出することで、それらの平面分布を求めることができるので、地盤構造の平面的な分布に対応する疑似三次元的な地震動に関する地盤特性を得ることができる。なお、隣接する2つの特定地点を通過する直線上の複数の予測地点において推定地震動伝達関数若しくは予測地震動特性値を算出する推定地震動算出工程には、特定地点間の予測地点において補間推定する場合に限らず、特定地点間の外側にある予測地点において外挿推定を行う場合も含む。
【0012】
本発明において、前記不整形地盤指標は、前記隣接する2つの特定地点間の距離と、地表面若しくは工学的基盤の高低差とに対応する直線状の傾斜を備えた仮想傾斜構造を特定する設定値であることが好ましい。これによれば、隣接する2つの特定地点の地盤データのみで地盤の不整形の影響を取り込むことができるとともに、仮想傾斜構造が特定地点間の距離と地表面若しくは工学的基盤の高低差(レベル差)に対応する直線状の傾斜を備えたものとされるため、簡易な傾斜構造であることから、推定地震動伝達関数の算出処理を容易に行うことが可能になる。ここで、上記不整形地盤指標の設定値としては、上記仮想傾斜構造の傾斜方向の長さと傾斜角の他に、同長さと高低差などといった、仮想傾斜構造を特定できる実質的に等価なものであれば種々の設定値を用いることができる。
【0013】
本発明において、前記平面分布は前記予測地震動特性値の平面分布であり、前記予測地震動特性値は、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期またはSI値のいずれかであることが好ましい。予測地震動特性値としては、予測地震動(波形)の特徴を何らかの意味で反映する各種の地震動パラメータを用いることができるが、具体的で一般的なパラメータとしては上記のものが挙げられる。
【0014】
本発明において、所定の入力地震動に基づいて前記推定地震動伝達関数により前記予測地点における予測地震動波形をさらに求めることが好ましい。これによれば、予測地震動波形、例えば数値加速度時刻歴など、を求めることにより、具体的な地震動に基づいて建物の耐震対策などをより確実かつ正確に策定することが可能になる。
【0015】
本発明において、前記特定地点は平面的に配列された複数の地盤メッシュの代表地点であり、前記不整形指標設定工程では、隣接する2つの前記地盤メッシュにおける前記特定地点の間の前記不整形地盤指標が算出され、前記推定地震動算出工程では、前記隣接する2つの地盤メッシュにおける前記特定地点を通過する直線上の、前記特定地点間の距離よりも小さな間隔で配列された複数の前記予測地点における前記推定地震動伝達関数若しくは前記予測地震動特性値が算出されることが好ましい。これによれば、地盤メッシュごとに存在する地盤データを元にその代表地点における特定地震動伝達関数を求めるとともに、隣接する地盤メッシュの代表地点を通過する直線上において代表地点間の距離よりも小さな間隔で配列された複数の予測地点を設定することにより、地盤メッシュ単位よりも細密に、地震動に関する地盤特性の平面分布を得ることができる。
【0016】
本発明の地震動ハザードマップの制作方法は、上記の地震動予測方法を用いて、所定の入力地震動に基づいて前記推定地震動伝達関数により得られる前記予測地震動特性値の平面分布を求め、該予測地震動特性値の平面分布をマップ上に表示することが好ましい。これによれば、予測地震動特性値の平面分布をマップ上に表示することにより、地震動に関する地盤特性の平面分布の態様を地盤メッシュに制約を受けずに細密に示すことが可能になる。
【発明の効果】
【0017】
本発明によれば、不整形地盤による地震動の屈折・反射の影響を反映させることができるとともに、有限要素法のような高いコストを要することなく、簡易に地震動に関する地盤特性の平面分布を得ることができるという優れた効果を奏し得る。特に、ハザードマップやコンターマップなどのマップ上に推定地震動伝達関数若しくは予測地震動特性の平面分布を表示することにより、滑らかな地震動に対する地盤特性の表示態様が実現できる。また、地盤メッシュの地盤データに基づいて推定地震動伝達関数若しくは予測地震動特性の平面分布を求めることにより、地盤メッシュの制約を受けずに、不整形地盤に対応した細密なものとした、地震動に関する地盤特性の平面分布を簡易に得ることができる。
【図面の簡単な説明】
【0018】
【図1】本実施形態において考慮する不整形地盤をタイプごとに模式的に示す地盤モデル構造図であり、(a)は地表面に高低差(傾斜)を有するタイプAの地盤構造に対するもの、(b)は工学的基盤に高低差(傾斜)を有するタイプBの地盤構造に対するもの。
【図2】傾斜基盤を備えた不整形地盤における、(a)ポイントA(図1(b)に示す傾斜上端地点)、(b)ポイントB(同、傾斜中央地点)、(c)ポイントC(同、傾斜下端地点)における地震動伝達関数を示す周波数特性グラフであり、各グラフにおいて、各ポイントにおける、符号1Dを付した線は1次元の重複反射理論による等価線形解析による解析結果、符号2Dを付した線は2次元FEM解析による解析結果、太い2本の灰色線AとCはそれぞれ比較のために示したポイントAとポイントCにおける1次元の重複反射理論による等価線形解析による解析結果である。
【図3】1次元の地震動伝達関数ΩTとΩBを用いた簡易推定手法を説明するための説明図。
【図4】重み係数CTのXLS(=x/Ls)依存性を傾斜角ごとに示すグラフ(a)〜(d)。
【図5】下部に示す二次元有限要素法を適用した分割要素を示す図(分割ブロック構造図)に示すポイントA、ポイントB及びポイントCにおける地震動伝達関数の周波数特性を示すグラフであり、符号1Dを付した線は1次元の重複反射理論による等価線形解析による解析結果、符号2Dを付した線は2次元FEM解析による解析結果、符号SEを付した線は本実施形態に用いた簡易推定法による解析結果である。
【図6】ポイントA、ポイントB及びポイントCにおけるベクトル軌跡を示す図(ナイキスト線図)。
【図7】縦横に隣接する4つの地盤メッシュ1〜4の解析結果としての特定震度と模式化した地盤構造とを対比して4組の隣接する2つの地盤メッシュの対をそれぞれ示す説明図。
【図8】縦横に隣接する4つの地盤メッシュ1〜4の特定震度を濃淡で示すとともに各地盤メッシュの中心点である代表地点1a〜4aを示す説明図。
【図9】4組の隣接する2つの地盤メッシュの対ごとに、2つの代表地点を通過する直線上の予測震度を濃淡と正方形ブロックの高低で示す図。
【図10】東西方向(図示左右方向)に沿って算出した予測震度を同方向に配列された2つの代表地点間の範囲に限定して濃淡で示す左図と、南北方向(図示上下方向)に沿って算出した予測震度を同方向に配列された2つの代表地点間の範囲に限定して濃淡で示す右図とを対比して示す図。
【図11】各地盤ブロックの4つの代表地点で囲まれる正方形の範囲内で東西方向に沿って算出した予測震度と南北方向に沿って算出した予測震度とを重ね合わせた結果を濃淡で表示した図。
【図12】本実施形態の対象地区として選定した範囲を5箇所の区域とともに示すマップ。
【図13】本実施形態による図12に示す対象地区の解析結果である予測震度の平面分布を濃淡で示す図。
【図14】従来の1次元の等価線形化手法による地盤メッシュごとの特定震度の平面分布を濃淡で示す図。
【図15】図13に示す解析結果を、図12に示す三輪の区域について拡大して示す図。
【図16】図14に示す解析結果を、図12に示す三輪の区域について拡大して示す図。
【図17】図13に示す解析結果を、図12に示す富竹の区域について拡大して示す図。
【図18】図14に示す解析結果を、図12に示す富竹の区域について拡大して示す図。
【図19】図13に示す解析結果を、図12に示す安茂里の区域について拡大して示す図。
【図20】図14に示す解析結果を、図12に示す安茂里の区域について拡大して示す図。
【図21】図13に示す解析結果を、図12に示す篠ノ井の区域について拡大して示す図。
【図22】図14に示す解析結果を、図12に示す篠ノ井の区域について拡大して示す図。
【図23】図13に示す解析結果を、図12に示す松代の区域について拡大して示す図。
【図24】図14に示す解析結果を、図12に示す松代の区域について拡大して示す図。
【発明を実施するための形態】
【0019】
次に、添付図面を参照して本発明の実施形態について詳細に説明する。最初に、本実施形態において特定地震動伝達関数や特定地震動特性値の算出に用いられる算出方法の例として、等価線形化手法(SHAKE)及び周波数依存型等価線形化手法(FDEL)について説明する。
【0020】
地盤震動解析法では、重複反射理論を用いることで、複数の地層における地震動の解析を行うことができる。せん断波速度Vsの高い工学的基盤(一般的には300〜700m/s、本実施形態では500m/s以上を基盤層としている。)上に複数の地層(jは1〜nの自然数)が存在する場合には、工学的基盤からSH波の一次元的な上昇伝達について地層境界での反射・屈折を考慮することにより、SH波の変位波形の基礎方程式から、ひずみは以下の数式1により表される。ここで、gjは数式2で表される。なお、γjは地表からj番目の第j層のひずみ、ujは第j層での地震波の振幅、tは時間、Ajは第j層の上昇波の変位振幅、Bjは第j層の下降波の変位振幅、Gjは第j層のせん断剛性、zjは第j層の厚さ、ωは円振動数、ρjは第j層の密度、ξj(=2Gjhj/ω、hjは減衰定数)は第j層の粘性係数である。このとき、(1)地表面でせん断応力が0、(2)地層の境界面でせん断応力は連続、(3)地層の境界面で変位が連続の3つの境界条件を用いると、地表層についてはA1=B1が成立し、これによって、Ajは以下の数式3、Bjは以下の数式4の漸化式によって表される。
【0021】
【数1】
【0022】
【数2】
【0023】
【数3】
【0024】
【数4】
【0025】
この場合、地表の加速度に対する第j層の上面の加速度波の応答倍率である伝達関数Tfjは、以下の数式5で表される。
【0026】
【数5】
【0027】
第j層が工学基盤上の基盤層とすれば、工学基盤への入力波は基盤層の上昇波と考えることができるので、第j層の応答倍率(伝達関数)T*fjは以下の数式6で表される。ここで、a*jは上昇波だけの加速度波である。また、解放基盤波形を求めるときは、応答倍率は上昇波を2倍して以下の数式7とする。
【0028】
【数6】
【0029】
【数7】
【0030】
地表の地震動をフーリエ変換し、周波数領域で上記数式5のTfj又は上記数式6又は7のT*fjを乗じた後にフーリエ逆変換することにより、任意の層の加速度波形を得ることができる。また、同周波数領域で上記数式1のγjを乗じることで、任意の層のひずみ波形を得ることができる。
【0031】
一方、地盤の地震応答解析において用いられてきた等価線形化手法の上記SHAKEでは、以下の数式8に示すように、ひずみの最大値の或る一定の割合を等価平均ひずみγeとし、これより定まるせん断剛性と減衰定数を周波数領域の計算において周波数に依存することなく一律に適用している。ここで、Cは係数、γmaxはひずみの時刻歴波形の最大値である。上記SHAKEでは係数Cを0.65としている。
【0032】
【数8】
【0033】
上記のように、等価線形化手法(SHAKE)でも或る地点の地盤データから一次元的な地震動解析を行うことができる。しかしながら、地震動のレベルが上がりひずみが大きくなると地盤のせん断剛性と減衰定数は非線形に変動し、特に低・中周波数領域では実際の観測記録と比較的よく一致するが、強震動や軟弱地盤の応答解析を行う場合、高周波数領域で計算された結果が観測値を下回るなど、地盤の非線形性の影響による修正を考慮する必要がある。このため、周波数依存型等価線形化手法(FDEL)では、各周波数ごとにひずみ波形に寄与する度合に応じて適切なせん断剛性と減衰定数を与える周波数依存型の等価平均ひずみγf(ω)を以下の数式9のように定義する。ここで、Fγ(ω)はひずみ波形のフーリエスペクトル、Fγmaxはひずみ波形のフーリエスペクトルの最大値である。なお、C=0.65、Fγ(ω)/Fγmax=1.0とすると上記SHAKEと同じになる。係数Cの値はひずみ波形のフーリエスペクトルの卓越周波数周辺の伝達関数に大きく影響し、種々の検討によりC=0.65程度で問題がないことが確認されている。ただし、この係数Cの値は実情に合わせて調整可能である。
【0034】
【数9】
【0035】
ここで、せん断ひずみγとせん断剛性Gとの関係、並びに、せん断ひずみγと減衰定数hとの関係を用いることにより、ひずみに依存する非線形性を考慮した応答を得ることができる。非線形な応力とひずみの関係としては、例えば、拘束圧の影響を考慮してモデル化されたHardin-Drnevichモデルを用いることができる。このモデルは、地層の土質別に種々の拘束圧におけるせん断強度に対する実験に基づいて、各土質ごとに基準化ひずみγγを定義し、ひずみとせん断剛性、ひずみと減衰定数の関係を求めたものである。ここで、数式9に示す周波数依存型の平均ひずみγf(ω)を用いる場合には、以下の数式10によりせん断剛性G(ω)が表され、以下の数式11により減衰定数h(ω)が表され、いずれも周波数の関数になる。ここで、Gmaxは初期せん断剛性、hmaxは最大減衰定数、γγは基準化ひずみ(=τmax/Gmax、τmaxはせん断応力の最大値)である。
【0036】
【数10】
【0037】
【数11】
【0038】
上記の等価線形化手法(SHAKE)や周波数依存型等価線形化手法(FDEL)では、一次元の地盤震動解析により地震動伝達関数を算出することができる。一方、二次元有限要素法(プログラム名:FRUSH)や、これを改良した周波数依存型等価線形化法による二次元FEM地盤震動解析法(プログラム名:FDEL−FEM)も使用されているが、これらは膨大な計算処理を必要とする。そこで、本実施形態では、上述の種々の一次元地盤震動解析により求められた地震動伝達関数に基づいて、不整形地盤の簡易推定法を用いる。まず、この簡易推定法と二次元有限要素法による解析結果とを比較し、簡易推定法の有効性を示す。
【0039】
図1は、簡易推定法と二次元有限要素法の解析結果の比較に用いるモデル地盤として、工学的基盤が平坦で盛土などにより地表面が傾斜している不整形地盤のモデル(a)と、地表面が平坦で工学的基盤が傾斜している不整形地盤のモデル(b)をそれぞれ模式的に示す概略構成図である。ここで、領域1(Region-1)は表層地盤構造であり、領域2(Region-2)は工学的基盤である。二次元FEM地盤震動解析では、上記数式9に示す周波数依存型等価ひずみγf(ω)を導入し、C=0.65とした。また、モデル地盤の形状及び材料物性は、以下の表1及び表2に示すものとした。
【0040】
【表1】
【0041】
【表2】
【0042】
また、地震動伝達関数は地震動の入力強度に依存することから、マグニチュードMと震源距離Rを仮定し、入力レベルが異なるシミュレーション地震動を算出して解析に用いた。入力地震動の諸元を表3に示す。なお、表2において、領域1のせん断波速度(単位m/s)については、二種類(100と300)の数値でそれぞれ計算を行い、解析を行っている。
【0043】
【表3】
【0044】
上記タイプAの不整形地盤における解析結果の一例(Case−3,傾斜角1:5,Vs=300m/sec,入力地震動No.Acc−3)を図2に示す。ポイントAは図1に示す傾斜上端の地点、ポイントCは図1に示す傾斜下端の地点、ポイントBは着目地点である傾斜上の地点である。図中の2Dで示す実線は二次元FEM解析(2D)による伝達関数(地震動の出力スペクトル)、1Dで示す実線は着目地点であるポイントBと同じ深さの成層地盤を仮定して1次元の重複反射理論による等価線形解析(1D)を行った場合の結果である。また、指示線を有しない2本の灰色の太線は傾斜上端のポイントA及び下端のポイントCでの1Dによる伝達関数を示し、それぞれ傾斜上端のグラフをA、傾斜下端のグラフをCとした。
【0045】
傾斜上端の地表面であるポイントAでは、不整形性の影響が比較的少なく、2Dと1Dによる地震動伝達関数は、概ね一致している。傾斜中央のポイントBでは、地盤の不整形性の影響を受け、2Dによる伝達関数の卓越振動数は、1Dによる伝達関数の固有振動数と異なった位置に卓越振動数を持つ。また、傾斜下端のポイントCにおいてもポイントBと同様に不整形性の影響が確認できる。また、上述の結果は、タイプBのCase−1とCase−2についても共通していることが確認された。なお、2Dによる伝達関数は、傾斜上端および下端の1Dによる伝達関数の固有震動数の両方を卓越振動数として含むことから、傾斜部における地表面における伝達関数は、傾斜上端・下端における伝達関数(固有振動数)と密接に関係しているといえる。すなわち、地震動伝達関数は、傾斜部が存在することで傾斜方向に沿った他の地点の地盤構造に影響を受けることがわかる。
【0046】
不整形構造を有する地盤の傾斜部の地震動伝達関数の卓越震動数は、傾斜上端・下端の固有振動数の双方の位置で卓越し、その振幅は傾斜角度や地点に依存していることから、傾斜部上端・下端での1Dによる伝達関数に重みを付け、それらを複素平面上で重ね合わせることで当該地点の伝達関数を補正推定する簡便法を用いることができる。なお、この伝達関数は上述と同様に、SH波の鉛直入射に対する面内震動における入射地点から地表への水平震動成分の増幅率を表す。
【0047】
この簡易推定法では、以下の数式12に示すように、傾斜上端での1Dによる伝達関数ΩTと傾斜下端での1Dによる伝達関数ΩBの重み付き相乗平均として重ね合わせ、図3に示す任意地点(着目地点、傾斜上端からの距離xで規定される。)の伝達関数ΩEを推定する。伝達関数の掛算は複素平面上での回転操作となるため、振幅と位相角を正しく補間できる。一方、足し算を行うと、位相角が相反する場合、位相角と振幅が相殺される場合がある。なお、CTはΩTに対する重み係数、CBはΩBに対する重み係数である。このとき、伝達関数はいずれも複素数で与えられるものとし、重み係数は実数で与えられるものとする。
【0048】
【数12】
【0049】
上記数式12は、x/Ls=+∞のときに重み係数CTが0で右側の自由地盤の伝達関数ΩB(ω)、x/Ls=−∞のときに重み係数CTが1で左側の自由地盤の伝達関数ΩT(ω)とそれぞれ一致し、傾斜面上端からの距離xの変動により滑らかに変化することから、重み係数CTを以下の数式13で表されるロジスティック関数でモデル化する。
【0050】
【数13】
【0051】
傾斜面上端からの距離xと傾斜部の長さLsとの比(x/Ls)と傾斜角θをパラメータとして、数式12におけるΩEと有限要素法により求めたΩ2Dとの振幅の残差が固有振動数で最小となるように回帰分析により決定する。この決定方法の一例としては、例えば、数式13を変形すると数式14となり、この数式14を用いてb0及びb1を得ることができる。その結果を数式15に示す。なお、R′2は自由度調整済みの決定係数である。ただし、本発明では、重み係数の算出方法は上記の数式13乃至15に示す方法に限らない。例えば、上記のように傾斜構造の設定値(θやLs)だけでなく、表層地盤の物性(例えば、せん断波速度、ポアソン比、地盤堆積密度など)、表層地盤の地盤構成(地層数、土質、深さなど)、入力地震動の強度のいずれか一つ若しくは二以上を、陽に変数として表れる形で重み係数の数式13中に組み込んで算出し、数式12を計算してもよい。
【0052】
【数14】
【0053】
【数15】
【0054】
ここで、簡易推定式の構築に際して行った解析条件より、傾斜角の範囲を1:2〜1:5、傾斜部分の水平長さLsが500m程度以下をもって適用範囲とする。重み係数CTの標本値とモデル式を示す曲線(実線)を傾斜角ごとに図4に示す。ここで、図4(a)〜(d)は、表1に示す5つの全てのケースについて、それぞれ同じ傾斜角のものを選んで基礎となる地盤モデルとし、表2の地盤材料と表3の地震入力をそれぞれ適用して得た複数のFEM計算結果から算出した重み係数CTの標本値(薄い網目線)と、上述の方法で係数を決定した重み係数のモデル式(実線)のXLS=x/Lsに対する依存性(分布)を示すものである。図4の各図をみると、モデル式は各傾斜角の分布を表現できていることが確認できる。すなわち、タイプAとタイプBのいずれにおいても傾斜角θごとに上記の共通のモデル式を用いることが可能であることが示されている。
【0055】
簡易推定式を用いた伝達関数と二次元FEM解析による伝達関数の比較を図5に示す。対象地点を傾斜面上端のポイントA、傾斜面中央のポイントB、傾斜面下端のポイントCとする。図5の下図は二次元FEM解析に用いる地盤の分割要素の構成を示す図である。図中のSEで示す実線は簡易推定法により算出した伝達関数、2Dで示す実線は上記2D、1Dで示す実線は上記1Dによる当該地点の伝達関数である。ポイントAにおいて伝達関数の振幅の2Dによる結果と簡易推定法を適用して推定した結果を比較すると、簡易推定法を適用して推定した結果は2Dの伝達関数をほぼ再現できている。ポイントBでは、2Dと1Dの結果が異なることは前節記述した通りであり、推定した伝達関数はポイントAと同様に振幅、位相角の両方について2Dの伝達関数をほぼ再現できている。ポイントCに関しても同様の結果となり、簡易推定式は概ね2Dによる結果を再現できている。
【0056】
図6にはそれぞれのベクトル軌跡(ナイキスト線図)を表す。ベクトル軌跡は、複素平面上に伝達関数を表示したもので、ここでは周波数0Hzから3Hzまでの軌跡を表示した。図5と同様に、簡易推定法(SE)により、有限要素法(2D)の振幅(原点からの距離)をほぼ再現できるのを表すとともに、位相角(実軸と虚軸を基準とした角度)もほぼ同様に推定できていることを示している。
【0057】
ここで、タイプAの不整形地盤については、工学的基盤が同レベルであり、地表面に傾斜が存在するため、地表面の傾斜角を上記数式13のθとすることにより、中間地点xにおける伝達関数ΩEを求めることができる。また、タイプBの不整形地盤については、工学的基盤が傾斜し、地表面が同レベルであるため、工学的基盤の傾斜角を上記数式13のθとすることにより、中間地点xにおける伝達関数ΩEを求めることができる。
【0058】
また、地表面と工学的基盤のいずれにも傾斜がない場合(すなわちθ=0の場合、以下、タイプCという。)には、傾斜による波動の屈折や反射の影響がないと考えられるので、CT=1−x/Ls(0<x<Ls)、CB=1−CT(0<CB,CT<1)を用いて、重み係数CT、CBを線形補間により求めて上記数式12により中間点の波動関数ΩEを求めることができる。
【0059】
なお、本実施形態では、タイプA及びBについては、いずれも数式13によって数式12の重み係数CT及びCBを推定し、タイプCでは重み係数CT及びCBを線形補間して数式12を適用しているが、本発明はこのような態様に限られるものではなく、後述するように地盤タイプの少なくとも一部において不整形地盤指標に基づいて推定が行われさえすればよい。また、本実施形態では、タイプAとBとにおいて、傾斜角θと距離Lsが等しければ同じ方法(数式12及び13)で推定するようにしているので、例えば、地表面が一方側に傾斜角θで傾斜し工学的基盤が平坦なタイプAの地盤モデルは、地表面が平坦で工学的基盤が反対側に傾斜角θで傾斜しているタイプBと実質的に同じ推定処理を行うことになる。したがって、地表面と工学的基盤が共に傾斜し、どちらの傾斜も無視できない場合には、工学的基盤を基準としたときの地表面の高低差によって地表面の傾斜角を求め、工学的基盤が平坦なタイプAの地盤モデルとして計算するか、或いは、地表面を基準としたときの工学的基盤の高低差によって工学的基盤の傾斜角を求め、地表面が平坦なタイプBの地盤モデルとみなして計算すればよい。これによって全ての地盤構造をタイプA〜Cのいずれかの地盤モデルに対応させることができる。ただし、本実施形態の上述の方法に拘わらず、本発明ではタイプAとタイプBで相互に異なる方法(数式)で推定するようにしても構わない。この場合、本実施形態と同様に上述のようなタイプAとタイプBの間で変換可能な関係をもたせてもよいが、例えば地表面の傾斜角と工学的基盤の傾斜角のうちの傾斜角の大きい方の傾斜を残すように、地盤構造を何らかの指標に応じてタイプAとタイプBのいずれかの地盤モデルに分類するようにしてもよい。さらに、上述のように地表面と工学的基盤が共に傾斜し、どちらの傾斜も無視できない場合をタイプDとして別途の方法(数式)で推定するようにしてもよい。
【0060】
本実施形態では、一次元等価線形化手法等により得られた上記地震動伝達関数ΩT及びΩBに相当する特定地震動伝達関数と、上記Ls及びθに相当する不整形地盤指標(上記傾斜部を特定する設定値)とに基づいて、上記の簡易補間推定法を用いて、複数のメッシュ間で上記地震動伝達関数ΩEに相当する推定地震動伝達関数を求め、疑似3次元的な地盤構造の解析を行う。
【0061】
この場合、簡易推定法を適用する2つの特定地点の地盤データの関係に応じて、上記タイプA、タイプB及びタイプCのいずれか一つを選択するとともに、選択されたタイプごとに上記不整形地盤指標を設定し、これらに応じて重み係数CT及びCBを上述のタイプ別にそれぞれに求めることで、上記数式12により推定地震動伝達関数ΩEを求めることができる。なお、本実施形態では、不整形地盤指標は地盤の傾斜(高低差)に関するもののみであるので、傾斜のないタイプCについては実質的に不整形地盤指標が適用されずに推定地震動伝達関数ΩBが求められるが、この場合には、不整形地盤指標のうち傾斜角θ=0の特別な場合において不整形地盤指標とは無関係に算出しているだけであり、傾斜角θを0とみなすことのできないタイプA及びBの場合において不整形地盤指標に基づいた算出が行われうる構成である。このように、全ての地盤において等しく不整形地盤指標を用いる必要はなく、必要な一部の場合にのみ不整形地盤指標を用いるものであっても構わない。
【0062】
次に、この解析方法の原理を説明するために、図7に示すように隣接する4つの地盤メッシュを基にして解析を行う例を示す。図7に示す4つの地盤メッシュ1〜4の地盤構造は、各地盤メッシュ1〜4の代表地点(図示例では各地盤メッシュの中心点)1a〜4aにおける地盤データに対応する地盤構造(工学的基盤上の表層地盤の地層構成)を示す。この場合、解析には、従来法で使用していた地盤メッシュの地盤データをそのまま使用することができる。ここで、工学的基盤上の表層地盤が複数の地層よりなる場合には、上記重複反射理論に基づいて伝達関数を算出する。なお、各地盤メッシュでは、地盤データを採取した採取地点と、上記代表地点(特定地点)とが一致しない場合もあり得る。例えば、地盤データの採取地点が中心点以外にしかない場合、複数の採取地点がある場合であって複数の採取地点の寄与を含めた地盤データを用いる場合などでは、各地盤メッシュの中心点などの代表地点(特定地点)における地盤データは、実際には中心点などの代表地点(特定地点)とは異なる採取地点から取得した地盤データである。この場合、使用する地盤データの採取地点を地盤メッシュの代表地点(特定地点)としてもよく、或いは、採取地点から取得した地盤データを代表地点(特定地点)における地盤データと見做しても構わない。
【0063】
まず、基にする地盤データであるが、図7の下部に示されているような地盤構造が与えられているとし、この地盤構造に対応するせん断波速度、単位体積重量、土質、深さを用いる。ここで、土質とは、上述のHarding-Drnevichモデルで定義される基準化ひずみを与えるもので、具体的には、砂、粘土、礫、岩に対応する土質区分を言う。また、図7の下部の地盤構造の図において、最下部の二重線は工学的基盤との境界を示し、この二重線より上方部分は表層地盤の概略構成(1又は複数の地層の境界と厚み)を示す。この地盤構造では、実際には地表面レベルにはばらつきがあるものの、地表面レベルを基準として工学的基盤の高低差を図示しており、この場合には、上記のタイプBの地盤モデルを想定して計算を行うことになる。これを使い各地盤メッシュ1〜4での上記の一次元地震動解析を行って特定地震動伝達関数Ω1〜Ω4を求め、所定の共通の地震動(地震波スペクトル)を入力として工学的基盤に与えた場合について、図7の上部に示す震度が算出され(濃淡と数値で示してある。)、実際の震度予測図に用いられる平面上に落とすことで図8のようになる。ここで、濃淡と数値で示した各地盤メッシュの震度は、各地盤メッシュ1〜4の代表地点1a〜4aの震度に対応するものであるが、図示上では各地盤メッシュ全域に表している。この震度は、所定の入力地震動に基づいて上記特定地震動伝達関数から得られる特定地震動特性値である。なお、所定の入力地震動に基づいて地震動伝達関数から得られる地震動特性値としては、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期、SI値などが挙げられる。
【0064】
これらの4つの隣接する地盤メッシュ間で、各地盤メッシュ1〜4の代表地点(図示例では中央点)1a〜4aの地盤データに基づく特定地震動伝達関数Ω1〜Ω4と、それぞれの隣り合う2つの地盤メッシュの代表地点1aと2a、2aと3a、3aと4a、4aと1aの間の地表面や工学的基盤の傾斜に応じた上記と同様の不整形地盤指標で、2つの地盤メッシュ1と2、2と3、3と4、4と1の全体を(すなわち、隣接する2つの代表地点の間だけでなく、一方の代表地点に対して他方の代表地点とは反対側にある地点についても)補間・外挿推定し、これらの複数の予測地点について得られた推定地震動伝達関数ΩCに基づいてそれぞれ予測地震動特性値である上記共通の地震動を入力したときの予測震度を求めたものが図9である。上記の算出処理を行った隣接する2つの地盤メッシュの各組では、タイプA及びBの地盤モデルに対応する場合には上記数式12乃至15で示す方法で、タイプCの地盤モデルに対応する場合には上記の線形補間された重み係数による上記数式12により、東西方向(図8の左右方向)又は南北方向(図8の上下方向)をそれぞれ上記のx方向として、1つの地盤メッシュ当たりで等間隔の複数(図示例では10個)のx値(隣接する2つの地盤メッシュでは合計20個)の予測地点について、上記推定地震動伝達関数ΩCをそれぞれ求めた後、上記共通の入力地震動に基づいて各推定地震動伝達関数ΩCによる出力地震動に対応する予測震度をそれぞれのx値ごとに示してある。
【0065】
次に、東西方向・南北方向の2つの結果を地盤メッシュの代表地点(中心点)から隣の地盤メッシュの代表地点(中心点)までの範囲で抽出して示したものが図10である。ここで、図10の左側の図は東西方向の簡易補間推定結果であり、地盤メッシュ1と2、3と4でそれぞれ簡易補間推定を行った結果を東西方向の分布に反映させ、南北方向には一定の値としている。また、図10の右側の図は南北方向の簡易補間推定結果であり、地盤メッシュ1と4、2と3でそれぞれ簡易補間推定を行った結果を南北方向の分布に反映させ、東西方向には一定の値としている。
【0066】
次に、上述の東西方向・南北方向の結果を重ね合わせる。なお、この結果としては、図示のように予測震度などの予測地震動特性値を用いてもよく、予測震度などの予測地震動特性値を求める前の推定地震動伝達関数ΩCを用いてもよい。推定地震動伝達関数ΩCを用いた場合には重ね合わせた後に上記共通の地震動を入力することで予測震度(予測地震動特性値)を求めることができる。この時、東西方向・南北方向では予測震度に大きな差の出る場合がある。このため、地震波の伝達を考え、両方向で得られた2つの結果(2方向で計算した結果)の平均値や加算値を用いるなど、両方向の推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を共に反映させた値とする必要がある。ただし、予測地点の位置に応じて、両方向の算出値の重み付け(寄与)を適宜に設定して統合値(合成値)に反映させても構わない。
【0067】
また、図10に示した処理では、東西方向の簡易補間推定では、東西方向の位置(x値)が同じであれば同一の地盤メッシュ中の南北方向に沿った全ての予測地点(図面上では帯状の範囲)において同じ推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、また、南北方向の位置(x値)が同じであれば同一の地盤メッシュ中の東西方向に沿った全ての予測地点(図面上では帯状の範囲)において同じ推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、これらを平面的に重ね合わせることで合成している。しかし、東西方向と南北方向の簡易補間推定において、それぞれ対応する2つの代表地点を通過する直線上にある予測地点のみが算出された推定地震動伝達関数ΩC若しくは予測地震動特性値(予測震度)を有するものとし、両方向の直線上の予測地点が重ならないように設定して、これらの予測地点ごとに対応する予測地震動特性値がそのまま平面分布を構成するようにしてもよい。
【0068】
なお、地震動伝達関数から震度などの地震動特性値を算出するためには、入力地震動を決める必要があるが、本実施例の場合には、特定地震動特性値と予測地震動特性値のいずれについても、マグニチュード7.3相当の直下型地震を想定し、予測地区(長野市内)が震源からの距離20〜30km以内の震源域にあると仮定し、予測地区全体にわたり同一の基盤入力があったとして、同一の入力地震動を用いて予測震度を求めている。ただし、仮定された震源位置からの距離に応じた地震動の減衰を考慮し、上記特許文献5に示す強震動予測方法(EMPR)を用いて、地盤メッシュごとに、或いは、複数の地盤メッシュを含むより大きな区域ごとに工学的基盤へ入力される地震動波形を計算し、地盤メッシュ若しくは上記区域ごとに入力地震動を設定してもよい。本明細書では、上記のいずれの場合でも、共通の入力地震動に基づいて予測地震動特性値を求めているものとする。
【0069】
以上の過程から、図11の4つの地盤メッシュの中心1a〜4aを結んだ四角形の中の予測震度の表示(平面分布の作成)がなされる。この方法を対象地域内で隣接する地盤メッシュ全てで行い繋げていくことにより、図12に示す予測地区全域における予測震度を図13に示すように表示することができる。図13に示す簡易補間推定を行った結果は、隣接する各代表地点間、或いは、1つの地盤メッシュの東西方向及び南北方向の一辺当たり10個の予測地点を設定し、その結果を同一辺あたりで東西方向・南北方向それぞれに10段階で表示している。そのため、図14に示す従来法(1D)による特定地震動特性値のみを表示した場合に比べて1地盤メッシュあたり10×10で100倍の細かさで表示することができる。
【0070】
上記結果の詳細な比較・検討を行うため、図12に示す各区域の推定結果を示す。これらの各区域は以下の表4に示す5つの区域である。また、各区域における本実施形態による簡易補間推定結果及び従来法の表示を、図15と図16(三輪)、図17と図18(富竹)、図19と図20(安茂里)、図21と図22(篠ノ井)、図23と図24(松代)にそれぞれ対比して示す。なお、これらの各区域は一辺が3kmの正方形状で、各区域内を縦横それぞれ6つに分割してなる図中のメッシュサイズは500m×500mである。
【0071】
【表4】
【0072】
まず、図13の長野市全体での予測震度の表示結果(平面分布)を図14の従来手法と比較して述べる。本実施形態の方法により得られた地盤メッシュの境界における予測震度の表示は急激な震度の差が小さくなっていることが見て取れる。これは、予測震度が大きい場所から小さい場所にかけて段階的に変化している様子から、簡易補間推定法による正しい結果が表れているものと思われる。段階的に変化している部分では、従来法に比べ表示されている予測震度が大きくなっている場所がある。
【0073】
この長野市全体における結果を詳しく調べるため、次に比較対象地においての従来法との比較結果を見ていく。従来法では、地盤メッシュ内に山が入っている地点において山における地盤データが基になっているため、予測震度が小さく表示されている(例、安茂里2−5、松代4−4、※それぞれ左から数えたメッシュ数−下から数えたメッシュ数を示す。以下同様。)。それにより、例えば松代の例では上の地盤メッシュとの予測震度の差が大きくなり、境界や山の麓での予測震度が不明確であった。しかし、簡易補間推定法により予測震度の差は小さくなり、実際に起こりうる震度に近づけることが出来た。さらに、地盤メッシュに関係なく震度表示できるため、地盤メッシュ内での震度の変化が表され、山の形に合った震度表示がなされている。地盤メッシュの境界付近では地盤メッシュの移り変わりに伴い段階的に変化している様子が表されている(例、三輪5−2と5−3の境界、篠ノ井4−2と4−3の境界)。これは地盤の不整形性を考慮した結果であり、エネルギー収支の関係性が影響していることが示されている。
【0074】
どの地点においても隣り合う地盤メッシュ間でおきるエネルギーの収支の状態が表れており、地盤メッシュ境界での予測震度の大きな差がなくなった。これにより、地盤メッシュ境界付近の予測震度が小さかった地域でも安全側の予測結果を表示できる。各図において、対象地域の外縁で震度表示がなされていないが、これは簡易補間推定を行っているためである。簡易補間推定は隣り合う地盤メッシュのデータを基にメッシュ間の震度表示を行うため、隣り合う地盤メッシュのデータが存在しないと外縁の震度表示はされない。対象地域内の震度表示をするには対象地域に接するメッシュデータも必要となる。
【0075】
地盤の不整形性を考慮することにより、今まで予測震度の小さかった地域においても危険にさらされる可能性のあることが確認できた。また、従来法では地盤メッシュの境界で大きく震度が変化していたが、地盤メッシュに関係なく予測震度の表示が出来るようになったため、段階的に変化している状態が表されている。
【0076】
以上より、地盤の不整形性を考慮した長野市の詳細な震度予測図が作成できたと言える。今回、予測震度の大きくなる場所が多数見られ、このような場所に主要交通網や盛土斜面がある場合、震災時の予測被害は大きくなる。従来法による震度予測図では想定されなかった被害も想定できるようになったことから、地震対策を進める上で必要な情報として貢献できるものとなった。本実施形態で行った簡易補間推定法による解析は、1次元解析を行っていた今までの地盤データを使用することができるため、作成したプログラムを用いることで自治体レベルの広域での地震動予測を行える。本実施形態では入力地震動を対象地域内の工学的基盤面において同じ値を入力しているが、これを信濃川断層等の被害が予測される地震でのEMPR(地震動予測モデル)の結果を用いることで、実際に地震が発生したときの震度表示をすることが可能である。
【0077】
上記実施例では、推定地震動伝達関数ΩCから所定の入力地震動に基づいて予測地震動特性値(予測震度)を求めているが、さらに詳細な地震動予測を行うために、所定の入力地震動に基づいて推定地震動伝達関数ΩCから予測地点の地表面における予測地震動波形を求めるようにしてもよい。この地震動波形を求めることにより、より詳細な地震動の態様が判明するので、具体的な防災対策を策定する場合に有効である。この場合には、上記入力地震動として予想される震源域、規模、震動特性に応じたものを用いることが好ましい。
【0078】
尚、本発明の照明装置及び電気光学装置は、上述の図示例にのみ限定されるものではなく、本発明の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。例えば、上記実施例では、地盤メッシュをベースとして、隣接する地盤メッシュ間で簡易補間推定を実施しているが、例えば地盤データの採取地点間に上記方法を適用して複数の予測地点における算出値の平面分布を求めてもよいなど、本発明は地盤メッシュとは無関係に、隣接する特定地点間の簡易補間推定を行う場合を広く包含する。また、上記実施例では、隣接する2つの特定地点の特定地震動伝達関数に基づいて予測地点の推定地震動伝達関数を求め、この推定地震動伝達関数から予測地震動特性値である予測震度を導出する場合について説明したが、本発明は、地盤傾斜に関連付けられた不整形地盤指標に基づくものであれば、特定地震動特性値から予測地震動特性値を直接導出する場合をも包含する。さらに、上記実施形態では地表面と工学的基盤面の高低差に関する不整形地盤指標を取り扱っているが、表層地盤内に複数の地層が存在する場合、地表面と工学的基盤面以外の他の地層面の傾斜について考慮してもよい。
【符号の説明】
【0079】
ΩT…傾斜上端の地盤データに対応する自由端の地震動伝達関数
ΩB…傾斜下端の地盤データに対応する自由端の地震動伝達関数
ΩE…推定地震動伝達関数
Ls…傾斜部分の水平距離
x…傾斜上端からの水平距離
A…傾斜上端の地点
B…傾斜上の地点
C…傾斜下端の地点
θ…傾斜角
Ω1〜Ω4…特定地震動伝達関数
ΩC…推定地震動伝達関数
【特許請求の範囲】
【請求項1】
表層地盤の地質特性に関する地盤データから地震動特性を予測する地震動予測方法において、
前記地盤データから当該地盤データの採取地点若しくは該採取地点に対応して設けられた別の地点である特定地点における特定地震動伝達関数若しくは該特定地震動伝達関数から得られる特定地震動特性値を複数の前記特定地点についてそれぞれ算出する特定地震動算出工程と、
隣接する2つの前記特定地点間における地盤構成の高低差に対応する不整形性を表す不整形地盤指標を設定する不整形指標設定工程と、
前記隣接する2つの特定地点における前記特定地震動伝達関数と、前記隣接する2つの特定地点間に対応する前記不整形地盤指標とに基づいて、前記隣接する2つの特定地点を通過する直線上の複数の予測地点における推定地震動伝達関数若しくは該推定地震動伝達関数から得られる予測地震動特性値を算出する推定地震動算出工程と、
を具備し、
共通の前記特定地点について前記不整形指標設定工程及び前記推定地震動算出工程を異なる2方向に沿って実行して各方向に対応する前記予測地点についてそれぞれ前記推定地震動伝達関数若しくは前記予測地震動特性値を算出することにより、前記推定地震動伝達関数若しくは前記予測地震動特性値の平面分布を求めることを特徴とする地震動予測方法。
【請求項2】
前記不整形地盤指標は、前記隣接する2つの特定地点間の距離と、地表面若しくは工学的基盤の高低差とに対応する直線状の傾斜を備えた仮想傾斜構造を特定する設定値であることを特徴とする請求項1に記載の地震動予測方法。
【請求項3】
前記平面分布は前記予測地震動特性値の平面分布であり、前記予測地震動特性値は、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期またはSI値のいずれかであることを特徴とする請求項1または2に記載の地震動予測方法。
【請求項4】
所定の入力地震動に基づいて前記推定地震動伝達関数により前記予測地点における推定地震動波形をさらに求めることを特徴とする請求項1ないし3のいずれか一項に記載の地震予測方法。
【請求項5】
前記特定地点は平面的に配列された複数の地盤メッシュの代表地点であり、
前記不整形指標設定工程では、隣接する2つの前記地盤メッシュにおける前記特定地点の間の前記不整形地盤指標が算出され、
前記推定地震動算出工程では、前記隣接する2つの地盤メッシュにおける前記特定地点を通過する直線上の、前記特定地点間の距離よりも小さな間隔で配列された複数の前記予測地点における前記推定地震動伝達関数若しくは前記予測地震動特性値が算出されることを特徴とする請求項1ないし3のいずれか一項に記載の地震動予測方法。
【請求項6】
請求項5に記載の地震動予測方法を用いて、所定の入力地震動に基づいて前記推定地震動伝達関数により得られる前記予測地震動特性値の平面分布を求め、該予測地震動特性値の平面分布をマップ上に表示することを特徴とする地震動ハザードマップの制作方法。
【請求項1】
表層地盤の地質特性に関する地盤データから地震動特性を予測する地震動予測方法において、
前記地盤データから当該地盤データの採取地点若しくは該採取地点に対応して設けられた別の地点である特定地点における特定地震動伝達関数若しくは該特定地震動伝達関数から得られる特定地震動特性値を複数の前記特定地点についてそれぞれ算出する特定地震動算出工程と、
隣接する2つの前記特定地点間における地盤構成の高低差に対応する不整形性を表す不整形地盤指標を設定する不整形指標設定工程と、
前記隣接する2つの特定地点における前記特定地震動伝達関数と、前記隣接する2つの特定地点間に対応する前記不整形地盤指標とに基づいて、前記隣接する2つの特定地点を通過する直線上の複数の予測地点における推定地震動伝達関数若しくは該推定地震動伝達関数から得られる予測地震動特性値を算出する推定地震動算出工程と、
を具備し、
共通の前記特定地点について前記不整形指標設定工程及び前記推定地震動算出工程を異なる2方向に沿って実行して各方向に対応する前記予測地点についてそれぞれ前記推定地震動伝達関数若しくは前記予測地震動特性値を算出することにより、前記推定地震動伝達関数若しくは前記予測地震動特性値の平面分布を求めることを特徴とする地震動予測方法。
【請求項2】
前記不整形地盤指標は、前記隣接する2つの特定地点間の距離と、地表面若しくは工学的基盤の高低差とに対応する直線状の傾斜を備えた仮想傾斜構造を特定する設定値であることを特徴とする請求項1に記載の地震動予測方法。
【請求項3】
前記平面分布は前記予測地震動特性値の平面分布であり、前記予測地震動特性値は、地震動の最大加速度、最大速度、最大変位、震度、実効加速度、卓越周期またはSI値のいずれかであることを特徴とする請求項1または2に記載の地震動予測方法。
【請求項4】
所定の入力地震動に基づいて前記推定地震動伝達関数により前記予測地点における推定地震動波形をさらに求めることを特徴とする請求項1ないし3のいずれか一項に記載の地震予測方法。
【請求項5】
前記特定地点は平面的に配列された複数の地盤メッシュの代表地点であり、
前記不整形指標設定工程では、隣接する2つの前記地盤メッシュにおける前記特定地点の間の前記不整形地盤指標が算出され、
前記推定地震動算出工程では、前記隣接する2つの地盤メッシュにおける前記特定地点を通過する直線上の、前記特定地点間の距離よりも小さな間隔で配列された複数の前記予測地点における前記推定地震動伝達関数若しくは前記予測地震動特性値が算出されることを特徴とする請求項1ないし3のいずれか一項に記載の地震動予測方法。
【請求項6】
請求項5に記載の地震動予測方法を用いて、所定の入力地震動に基づいて前記推定地震動伝達関数により得られる前記予測地震動特性値の平面分布を求め、該予測地震動特性値の平面分布をマップ上に表示することを特徴とする地震動ハザードマップの制作方法。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【公開番号】特開2012−251815(P2012−251815A)
【公開日】平成24年12月20日(2012.12.20)
【国際特許分類】
【出願番号】特願2011−123262(P2011−123262)
【出願日】平成23年6月1日(2011.6.1)
【出願人】(504237050)独立行政法人国立高等専門学校機構 (656)
【公開日】平成24年12月20日(2012.12.20)
【国際特許分類】
【出願日】平成23年6月1日(2011.6.1)
【出願人】(504237050)独立行政法人国立高等専門学校機構 (656)
[ Back to top ]