JP2006266940A

JP2006266940A - 入力地震動の推定方法

Info

Publication number: JP2006266940A
Application number: JP2005086989A
Authority: JP
Inventors: Shingo Nishimura; 新吾西村
Original assignee: Sekisui Chemical Co Ltd
Current assignee: Sekisui Chemical Co Ltd
Priority date: 2005-03-24
Filing date: 2005-03-24
Publication date: 2006-10-05
Anticipated expiration: 2025-03-24
Also published as: JP4339806B2

Abstract

【課題】簡易かつ高精度に入力地震動を推定することのできる入力地震動の推定方法を提供すること。
【解決手段】公開された基盤面における地震波形データを使用して、複数の解析方法に基づいて地表面または表層地盤の適宜の深度における地震動波形を作成し、作成された地震動波形の中から適宜の地震動波形を選択することによって入力地震動が推定される方法である。ここで、複数の解析方法には、表層地盤の材料特性を線形として扱う線形解析法と、表層地盤の歪み依存性が考慮された等価線形解析法と、表層地盤の応力とひずみの関係を表したヒステリシスに基づく非線形解析法とが少なくとも含まれている。公開された基盤面における地震波形データが所定規模以上の地震波形データである場合には、例えば、非線形解析法に基づいて作成された地震動波形が入力地震動であると推定できる。
【選択図】図７

Description

本発明は、基盤面上に堆積する表層地盤上または表層地盤内に建設される構造物の地震応答解析モデルに入力される入力地震動の推定方法に係り、特に、簡易かつ高精度に入力地震動を推定することのできる入力地震動の推定方法に関するものである。

地上や地中に建設される構造物の設計に際しては、特に地震時に十分な耐力（靭性なども含む）が確保できるように構造物の設計がなされるが、近時においては、近い将来にその発生が危惧されている東海地震や東南海地震などを想定した構造解析の必要性が叫ばれている。

構造物の地震応答解析に際しては、該構造物が建設される地域の地盤をどのようにモデル化するかによって構造物に入力される地震動波形が多様に変化する。基盤面（工学的基盤）に入射された地震動は、表層地盤の地盤特性によっては増幅されることがあり、増幅された地震動が地上構造物に載荷される。構造解析の際に使用される入力地震動が相違することによって応答解析結果も多様に相違することとなり、それに応じて実際の構造物の構造や使用される構造部材の仕様も相違することとなる。したがって、構造物の解析モデルに入力する地震動をどのように推定するのか、さらにはそれに付随する要素として地盤特性をどのように評価するのかといった問題は地震応答解析にとって極めて重要なものである。

耐震設計における入力地震動の推定方法には種々の方法がある。その１つは理論的方法と呼ばれる方法であり、この方法は、断層の破壊現象を時間的／空間的に考慮して地震動を推定する方法である。他の方法として小地震の重ね合わせによる半経験的な地震動評価法によるものであり、この方法は、地震波の伝播経路の影響や地盤の増幅特性を同じ伝播経路を持つ小地震記録を用いて表現し、それらを大地震の破壊過程に従って重ね合わせることにより地震動を評価するものである。さらにその他の方法としては、距離減衰式に断層の広がりの効果を考慮した方法や、観測された強震記録に基づく方法などがある。

特許文献１には、模擬地震動作成方法やその装置に関する発明が開示されている。かかる方法は、構造物の建設対象地域近傍の歴史地震や活断層を特定し、解析パラメータや断層パラメータ、破断パラメータ、基盤模擬地震動パラメータに基づいて基盤模擬地震動を作成し、作成された基盤模擬地震動と建設地域の地盤の地盤物性値に基づいて構造物に入力される模擬地震動を作成するものである。

特開平１０−３２０４５０号公報

上記する従来の地震動の推定方法や特許文献１に開示の模擬地震動の作成方法では、断層パラメータをはじめとする多数のパラメータの設定が必要となり、パラメータの設定に多大な時間を要することに加えて、解析する技術者ごとに地震応答解析結果が相違する可能性が高くなる。また、上記する観測された強震記録に基づく方法によれば、極めて簡便に地震動波形を推定することが可能となる一方で、強震記録には建設地域の地盤特性が考慮されていないことから地震応答解析結果の精度は極めて低いものにならざるを得ない。

本発明の入力地震動の推定方法は、上記する問題に鑑みてなされたものであり、汎用性
があって簡便かつ高精度な入力地震動の推定方法を提供することを目的とする。また、解析をおこなう技術者の技量などに左右されることなく、一般的な地震応答解析にかかる知識を有する技術者であればその解析結果にほとんど差異が生じ得ない入力地震動の推定方法を提供することを目的とする。

前記目的を達成すべく、請求項１に記載の発明による入力地震動の推定方法は、基盤面上に堆積する表層地盤上の地表面または表層地盤内に建設される構造物の地震応答解析モデルに入力される入力地震動の推定方法であって、前記入力地震動の推定方法は、公開された基盤面における地震波形データを使用して、複数の解析方法に基づいて地表面または表層地盤の適宜の深度における地震動波形を作成し、作成された地震動波形の中から適宜の地震動波形を選択することによって入力地震動が推定される方法であり、該複数の解析方法には、表層地盤の材料特性を線形として扱う線形解析法と、表層地盤の歪み依存性が考慮された等価線形解析法と、表層地盤の応力とひずみの関係を表したヒステリシスに基づく非線形解析法とが少なくとも含まれていることを特徴とする。

ここで、公開された基盤面における地震動波形とは、特に阪神淡路大震災以降に各種の学会や協会、各種の研究機関（大学等）にて作成され、公開されている適宜の地震動波形を意味する。かかる公開された地震動波形を使用することで、断層パラメータや断層地点から建設予定地点までの距離減衰パラメータ等の多様なパラメータを設定する必要がなくなり、極めて簡便に基盤面における地震動波形を特定することが可能となる。尤も、建設予定地域に近い断層等を想定して作成された公開地震動波形を使用するのが望ましい。

この公開された基盤面における地震動波形を使用して、基盤面上に堆積する表層地盤の地盤特性を勘案しながら地中の任意の地点または地表における地震動波形（構造物の解析モデルに入力する入力地震動波形）を解析する。

ここで、基盤面における地震動波形から地中の任意の地点または地表における地震動波形を解析する方法としては、表層地盤の地盤モデルや表層地盤の材料特性（材料特性モデル）に応じて大きく３つの方法がある。

その１つは、多層構造または単層構造の表層地盤を１次元の成層構造モデルとし、表層地盤の材料特性はその歪み依存性を考慮しない線形であると仮定して表層地盤をモデル化し、かかる表層地盤モデルの下方（基盤面）から公開地震動波形を入力して地震動波形を解析する方法である。本発明においては、この方法を線形解析法としている。

他の１つは、多層構造または単層構造の表層地盤を１次元の成層構造モデルとするとともに、各層の境界面における地震動の屈折と反射を考慮した重複反射理論を適用した地盤モデルを設定し、表層地盤の材料特性は地盤の応力と歪みの関係を部分的に線形として扱いながら擬似的に非線形性を考慮した等価線形モデルとして解析する方法である。基盤面における公開地震動波形は、表層地盤内で地盤の歪み依存特性に応じて増幅され、地中の任意の地点または地表における地震動波形が解析される。

さらに他の１つは、表層地盤を減衰特性を有するいくつかの質点とばねに置き換えた地盤モデルとし、表層地盤の材料特性は表層地盤の応力とひずみの関係を表したヒステリシスに基づく非線形特性を備えた材料特性として解析する方法である。かかる応力とひずみの関係を表したヒステリシスモデルとしては、公知のバイリニアモデルやランバーグ・オスグット（Ｒａｍｂｅｒｇ−Ｏｓｇｏｏｄ）モデル、ハーディン・ドルネビッチ（Ｈａｒｄｉｎ−Ｄｒｎｅｖｉｃｈ）モデルなどを使用することができる。

本発明においては、少なくとも上記する解析方法にて入力地震動の解析をおこない、例えば経済性や建設される構造物の重要度などに応じて適宜の入力地震動を地震応答解析に使用する入力地震動として推定するものである。

また、請求項２に記載の発明による入力地震動の推定方法は、前記公開された基盤面における地震波形データが所定規模以上の地震波形データである場合には、前記非線形解析法に基づいて作成された地震動波形が入力地震動であると推定されることを特徴とする。

大規模な地震力が作用する場合においては、地表面付近の歪みレベルも大きくなり、場合によっては１０^−３を超えることもある。かかる場合には土のダイレタンシーを考慮する必要が生じ、歪みレベルごとのヒステリシスループを等価のせん断弾性係数と減衰定数に置き換える等価線形解析法が適用できない。尤も、かかる場合に線形解析法が適用できないことは言うまでもない。

したがって、例えば、基盤面の地震動波形が所定規模以上の地震である場合には、非線形解析法にて解析された地震動波形を入力地震動とするものである。なお、ここでいう所定規模とは、例えば地震動の加速度波形における最大値などを意味しており、所定規模を２５０ｇａｌや３００ｇａｌなど適宜の値に設定することができる。

なお、比較的大きな地震力が作用する場合であって等価線形解析の適用が可能な場合には、非線形解析法ではなく等価線形解析法を適用するのが好ましい。非線形解析法では過度な減衰を示すことがあるため、地表面における地震動波形が現実のものよりも小さくなる可能性があるためである。

また、請求項３に記載の発明による入力地震動の推定方法は、基盤面上に堆積する表層地盤上の地表面または表層地盤内に建設される構造物の地震応答解析モデルに入力される入力地震動の推定方法であって、前記入力地震動の推定方法は、公開された基盤面における地震波形データを使用して、複数の解析方法に基づいて地表面または表層地盤の適宜の深度における地震動波形を作成し、作成された地震動波形の大小を比較し、適宜の地震動波形を選択することによって入力地震動が推定される方法であり、該複数の解析方法には、表層地盤の材料特性を線形として扱う線形解析法と、表層地盤の歪み依存性が考慮された等価線形解析法とが少なくとも含まれていることを特徴とする。

例えば、線形解析法によって地表の地震動波形を解析する一方で、等価線形解析法によっても別途、地表の地震動波形を解析する。双方の地震動波形の規模を比較し、過大な設計をしないという設計思想に基づく場合には小さい方の地震動波形を地震応答解析に使用し、安全側の設計をおこなう場合には大きい方の地震動波形を地震応答解析に使用する。

なお、本発明においては、その他の解析方法として、既述する非線形解析法をも含むこともでき、この場合には３つの解析法に基づいて解析されたそれぞれの地震動波形の大小比較をおこない、設計思想に合った地震動波形を選定することができる。

また、請求項４に記載の発明による入力地震動の推定方法は、前記公開された基盤面における地震波形データと、選択された適宜の地震動波形とを比較し、基盤面における地震波形データから選択された適宜の地震動波形への増幅率が１を超える場合には選択された適宜の地震動波形が入力地震動と推定され、増幅率が１以下の場合には基盤面における地震波形データが入力地震動と推定されることを特徴とする。

本発明は、表層地盤における増幅率を考慮して入力地震動を決定（推定）するものである。例えば、等価線形解析法によって解析された地震動波形と線形解析法によって解析さ
れた地震動波形を比較し、過大な設計をしないとする設計思想に基づいて小さい規模の地震動波形を選定する。選定された地震動波形と、この地震動波形のもとになっている基盤面の公開地震動波形とを比較し、公開地震動波形の方が大きい規模の場合には公開地震動波形を選定し、公開地震動波形の方が小さい規模の場合には解析後の地震動波形を選定するものである。選定された地震動波形は、構造物の地震応答解析モデルに直接入力されることとなる。

かかる方法によれば、解析方法によって異なる規模の地震動波形が作成された場合でも過大な地震動波形は選定せず、その一方で表層地盤の地盤特性（ここでは地盤の増幅率）が勘案された地震動波形を入力地震動とすることが可能となる。尤も、地盤の増幅率が１よりも小さな場合には、地盤の場所ごとの不確定性などを勘案して基盤面における公開地震動波形を選定するものである。

さらに、請求項５に記載の発明による入力地震動の推定方法は、前記公開された地震波形データが、内閣府中央防災会議の東海地震に関する専門調査会及び東海地震対策専門調査会において作成された「想定東海地震に係る地震動データ」であることを特徴とする。

上記する内閣府中央防災会議の「東海地震に関する専門調査会」及び「東海地震対策専門調査会」において作成された想定東海地震に係る地震動データは、基盤面（工学的基盤）での想定東海地震波形と、基盤面から地表までの表層地盤モデルとからなるデータのことである。かかる公開データは、埼玉県、東京都、神奈川県、山梨県、静岡県、長野県、岐阜県、愛知県、三重県の１都８県に対して１ｋｍメッシュ単位で作成されている。

この公開データを積極的に使用することで、関東〜近畿にかけての広範囲にわたる地域で建設される構造物の設計に際し、東海地震を対象とした地震応答解析を簡易に実施することが可能となる。

以上の説明から理解できるように、本発明の入力地震動の推定方法によれば、公開された基盤面における地震動波形を複数の解析方法を介して地震動波形を解析し、設計思想などに応じて適宜の地震動波形を地震応答解析用の入力地震動とすることができるため、簡便で高精度な入力地震動を推定でき、さらには選定自由度の高い入力地震動の推定をおこなうことができる。また、本発明の入力地震動の推定方法によれば、公開されたデータとして想定東海地震に係る地震動データを使用することで、関東〜近畿にかけての広範囲にわたる地域で建設される構造物に対して、東海地震を対象とした地震応答解析を簡易におこなうことができる。さらに、本発明の入力地震動の推定方法によれば、地震応答解析をおこなう技術者の技量に左右されることなく、ほぼ同じ解析結果を得ることができる。

以下、図面を参照して本発明の実施の形態を説明する。図１は、基盤面における公開地震動波形データの一実施形態を示した図を、図２は、公開された任意地点の表層地盤の物性値を示した図をそれぞれ示している。図３は、等価線形解析法を説明したフローを、図４は、等価線形解析法において使用される歪み依存曲線をそれぞれ示している。図５は、非線形解析法における表層地盤のモデル化を説明した模式図である。図６は、非線形解析法における表層地盤の応力とひずみの関係を表したヒステリシスを示した図であり、（ａ）は、バイリニアモデルを、（ｂ）は、ランバーグ・オスグットモデルを、（ｃ）は、ハーディン・ドルネビッチモデルをそれぞれ示している。図７は、本発明の入力地震動の推定方法の一実施形態を示したフローを、図８は、本発明の入力地震動の推定方法の他の実施形態を示したフローをそれぞれ示している。図９ａは、は、基盤面における公開地震動の時刻歴波形を示した図であり、図９ｂは、図８の入力地震動の推定方法によって解析さ
れた地震動の時刻歴波形を示した図である。

図１は、内閣府中央防災会議の東海地震に関する専門調査会及び東海地震対策専門調査会において作成された「想定東海地震に係る地震動データ」の一部を示したものである。図中のデータは地震動の加速度（ｇａｌ）を示しており、１／１００秒刻みのサンプルとなっている。データ数（ＮＤＡＴＡ）が１６３８４あるため、かかるデータに基づいて作成される地震動の時刻歴波形は１６３．８３秒までの波形が作成される。

図２は、公開地震動波形とともに公開された任意地域の表層地盤の物性値データを示したものである。ここで、動的変形特性番号について粘性土（図中の番号１）や砂質土（図中の番号２）などを指定すると、粘性土や砂質土に応じた歪み依存特性などが勘案されることとなる。一方、線形を設定するのは（図中の番号が４，５）、土質が比較的硬質であって地震時に地盤の歪みに応じて地震動の地盤内増幅がなされないと判断されるためである。

次に、図３，４に基づいて、等価線形解析法について説明する。表層地盤は、粘性土や砂質土など物性の相違する複数の土層から構成されるのが一般的である。そこで、等価線形解析法で解析をおこなう際には、予め各土層の歪み依存曲線が求められ、ステップ１（Ｓ１）にて初期物性値の設定がおこなわれる。ステップ２（Ｓ２）では、各土層における剛性の評価がなされ、ステップ３（Ｓ３）において周波数応答計算がおこなわれる。かかる周波数応答計算は複素応答解析と同一の手順でおこなわれるものである。周波数応答計算がおこなわれた後、ステップ４（Ｓ４）で有効歪みｅγが計算され、ステップ５（Ｓ５）では、計算に使用された物性値（各土層のせん断剛性と減衰定数）とステップ４で計算された物性値とが適合するか否かの判定がおこなわれる。双方が適合する場合にはそこで計算は終了となるが、適合しない場合にはステップ６（Ｓ６）で歪み依存曲線により計算された歪みに適合するように物性値の設定が再度なされ、ステップ２に戻る。かかる手順を繰返しながら解の収束を図るものである。図４はその一例である（Ｘ線はせん断剛性：Ｇを示す曲線、Ｙ線は減衰定数：ｈを示す曲線）が、最初に設定した物性値（Ｇ０，ｈ０）と計算された歪み：ｅγ１とが適合しないため、次に物性値をＧ１，ｈ１に設定し直して再度計算する。計算した結果の歪み：ｅγ２はやはり適合しないため、以後同様に繰返し計算をおこない、設定された物性値に適合する歪みとなるまで繰返し計算を実施する。

等価線形解析法においては、土層モデルは、１次元の成層地盤モデルが適用され、各土層ごとに既述する収束計算がおこなわれることで地盤の歪み依存性が考慮された任意地点（地表面や適宜の深度）における地震動波形が解析される。

一方、図５，６には、非線形解析法における地盤モデルと、地盤の材料特性を表すヒステリシスがそれぞれ示されている。図５に示す地盤モデルは、各土層（Ｇ１〜Ｇ４、Ｇ０は基盤）を減衰特性を有する幾つかの質点（Ｍ１，Ｍ２，Ｍ３，Ｍ４）とばね（Ｋ１，Ｋ２，Ｋ３，Ｋ４）に置き換えた、いわゆるばね−質点系モデルである。ここで、地盤は集中質点に置き換えられ、その間はせん断ばねにて結合される。地震時に地盤はせん断波（Ｓ波）が卓越することから、地盤を１次元のせん断振動系にモデル化するのが簡易でかつ効果的なモデル化であると言える。

このばね−質点系モデルにおけるせん断ばねには、以下に示す３つの材料特性を表すヒステリシスを適用することができる。

その１つは、図６ａに示すバイリニアモデルである。このモデルの骨格曲線は、原点から降伏ひずみ：γ_ｆまでは勾配Ｇ_０の直線とし、γ_ｙ以上では勾配Ｇ_ｆの直線となる。

他の１つは、図６ｂに示すランバーグ・オスグットモデルである。このモデルの骨格曲線（図中のＸ線）は、降伏歪み：γｙ、せん断応力：τおよび定数：αを用いて式１で表すことができる。なお、修正ランバーグ・オスグットモデルを使用することもできる。

[式１]
Ｇ_０γ＝τ＋α｜τ｜^γ／（Ｇ_０γ_ｙ）^γ−１

さらに図６ｃに示すモデルは、ハーディン・ドルネビッチモデルであり、微小歪みレベルでのせん断弾性係数：Ｇ_０、γ_ｒ（＝τ_ｆ／Ｇ_０）で定義される基準歪みを用いてその骨格曲線（図中のＹ線）は、式２で表すことができる。

[式２]
τ＝Ｇ_０γ／（１＋γ／γ_ｒ）

基盤面における地震動波形が比較的大きな場合であって、表層地盤内の土質が比較的軟質である場合には、地盤の歪みレベルが大きくなるため、線形解析は勿論のこと、等価線形解析法による方法でも地盤の歪みを十分に反映させることができない。そこで、かかるケースの場合においては、上記するばね−質点系モデルを適用し、ばねには上記する材料特性を反映したヒステリシスを導入して非線形解析法により地震動波形の解析をおこなうことが望ましい。

図７は、入力地震動の推定方法の一実施形態のフローを示したものである。まず、ステップ１において、基盤面における公開地震動波形の加速度の閾値（フロー中のＸｇａｌ）を設定する。このＸの値は２５０や３００など、例えば、表層地盤の地盤特性を勘案して、等価線形解析法と非線形解析法のどちらを選定するのが解析結果として好ましいかを技術者が判断して設定することができる。閾値であるＸｇａｌをステップ１（Ｓ１）で設定し、適用する公開地震動波形の加速度がその値以上の場合にはステップ２（Ｓ２）のルート、すなわち非線形解析法を採用してその解析された地震動波形を入力地震動として推定する（ステップ３（Ｓ３））。一方、公開地震動波形の加速度が閾値未満の場合には、ステップ４（Ｓ４）のルートを選定し、等価線形解析法にて解析された地震動波形を入力地震動として推定する。

また、図８には、入力地震動の推定方法の他の実施形態のフローを示している。ここでは、公開地震動波形を使用し、一方では線形解析法を実施して地震動波形を解析し（ステップ１１（Ｓ１１））、他方では等価線形解析法を実施して地震動波形を解析する（ステップ１２（Ｓ１２））。

双方の解析結果を比較し、例えば、等価線形解析法による結果（地震動波形の加速度）が線形解析法による結果よりも大きいか否かを判断する（ステップ１３（Ｓ１３））。図示する実施例では、過大な設計をしないという設計思想に基づいて小さな地震動波形を採用する場合を示しており、ステップ１４，１５では加速度の小さくなる地震動波形が選定される。

次に、表層地盤の増幅率を勘案するため、増幅率の値が１．０よりも大きいかそれ未満かの判定をおこなう（ステップ１６（Ｓ１６））。この表層地盤の増幅率は、基盤面から地表まで地震動が伝播する過程での地震動の増幅を表すものであり、表層地盤の地盤特性を反映した重要な要素である。そこで、ステップ１６での判定は、基盤面における公開地震動波形の加速度とステップ１４またはステップ１５で解析されて選定された地震動波形の加速度とを比較しておこなわれる。例えば、図９ａに示す地震動の時刻歴波形が基盤面における公開地震動波形であり、選定された地震動波形が図９ｂに示す地震動の時刻歴波
形である。図９ａによれば、加速度の最大値（絶対値の最大）は−５２５ｇａｌであり、図９ｂによれば、加速度の最大値は１１５９ｇａｌである。したがって、かかる実施例では、表層地盤内を地震動が伝播する過程で地表面においては約２倍に増幅された地震動が形成されることとなり、増幅率は約２ということになる。したがって、この場合には、図８のフローにおいて、ステップ１７（Ｓ１７）が選択されて、解析後の地震動波形が入力地震動として推定されることとなる。一方、増幅率が１未満の場合には、地盤の場所ごとの不確定性などを勘案して、公開地震動波形を入力地震動として推定する安全側の推定がなされることとなる（ステップ１８（Ｓ１８））。

以上、本発明の実施の形態を図面を用いて詳述してきたが、具体的な構成はこの実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲における設計変更等があっても、それらは本発明に含まれるものである。

基盤面における公開地震動波形データの一実施形態を示した図。公開された任意地点の表層地盤の物性値を示した図。等価線形解析法を説明したフロー。等価線形解析法において使用される歪み依存曲線。非線形解析法における表層地盤のモデル化を説明した模式図。非線形解析法における表層地盤の応力とひずみの関係を表したヒステリシスを示した図であり、（ａ）は、バイリニアモデルを示した図。（ｂ）は、ランバーグ・オスグットモデルを示した図。（ｃ）は、ハーディン・ドルネビッチモデルを示した図。本発明の入力地震動の推定方法の一実施形態を示したフロー。本発明の入力地震動の推定方法の他の実施形態を示したフロー。（ａ）は、基盤面における公開地震動の時刻歴波形を示した図。（ｂ）は、図８の入力地震動の推定方法によって解析された地震動の時刻歴波形を示した図。

Claims

基盤面上に堆積する表層地盤上の地表面または表層地盤内に建設される構造物の地震応答解析モデルに入力される入力地震動の推定方法であって、
前記入力地震動の推定方法は、公開された基盤面における地震波形データを使用して、複数の解析方法に基づいて地表面または表層地盤の適宜の深度における地震動波形を作成し、作成された地震動波形の中から適宜の地震動波形を選択することによって入力地震動が推定される方法であり、該複数の解析方法には、表層地盤の材料特性を線形として扱う線形解析法と、表層地盤の歪み依存性が考慮された等価線形解析法と、表層地盤の応力とひずみの関係を表したヒステリシスに基づく非線形解析法とが少なくとも含まれていることを特徴とする入力地震動の推定方法。
前記公開された基盤面における地震波形データが所定規模以上の地震波形データである場合には、前記非線形解析法に基づいて作成された地震動波形が入力地震動であると推定されることを特徴とする請求項１に記載の入力地震動の推定方法。
基盤面上に堆積する表層地盤上の地表面または表層地盤内に建設される構造物の地震応答解析モデルに入力される入力地震動の推定方法であって、
前記入力地震動の推定方法は、公開された基盤面における地震波形データを使用して、複数の解析方法に基づいて地表面または表層地盤の適宜の深度における地震動波形を作成し、作成された地震動波形の大小を比較し、適宜の地震動波形を選択することによって入力地震動が推定される方法であり、該複数の解析方法には、表層地盤の材料特性を線形として扱う線形解析法と、表層地盤の歪み依存性が考慮された等価線形解析法とが少なくとも含まれていることを特徴とする入力地震動の推定方法。
前記公開された基盤面における地震波形データと、選択された適宜の地震動波形とを比較し、基盤面における地震波形データから選択された適宜の地震動波形への増幅率が１を超える場合には選択された適宜の地震動波形が入力地震動と推定され、増幅率が１以下の場合には基盤面における地震波形データが入力地震動と推定されることを特徴とする請求項３に記載の入力地震動の推定方法。
前記公開された地震波形データが、内閣府中央防災会議の東海地震に関する専門調査会及び東海地震対策専門調査会において作成された「想定東海地震に係る地震動データ」であることを特徴とする請求項１〜４のいずれかに記載の入力地震動の推定方法。