地震防災のための地盤情報と今後の展開
Future Initiatives of Underground Structure Information for Earthquake Disaster Prevention.
先名 重樹＊
Shigeki Senna

つくば事業部第四技術部では、独立行政法人防災科学技術研究所（以下、「防災科研」と呼ぶ）研究
業務支援等により、地盤情報（地質情報）データベースおよび地盤構造モデル情報作成、地震動予測
地図作成作業における地盤情報の取り扱い手法等を検討してきた。本論文では、上記技術を利用した、将来の地盤情報サービスの展開や方向性についてそれぞれ述べる。
I worked continuously in the disaster prevention section for more than the past ten years as a
researcher（now, visiting researcher）of NIED. The contents of the measure are the directivity of
deployment of the foundation information in a underground structure information（geologicalinformation）database, structural model information creation, and seismic hazard map creation
work, and a future underground structure information service.

で公開している。J-SHISからは、全国250mメッシュサ

1．はじめにイズで、工学的基盤における最大速度と、その最大速度

我が国は、1995年兵庫県南部地震や2011年東北地方太を今後30年あるいは50年の間に超える確率、表層地盤の平洋沖地震など、過去において地震による甚大な被害を微地形区分と増幅率などのデータを取得することができ経験してきた。そして、将来においては南海トラフの巨る。これらのデータと住宅の建物情報を用いることによ

大地震や首都直下地震など、極めて甚大な被害を生じる

り、その場所においての地盤のゆれやすさや大きな揺れ

地震の発生が予想されている。そこで、国は2000年に建に見舞われる確率を知ることができ、そのとき予想され築基準法の改正を行い、建築物の基礎の構造は、国交省

る建物被害リスク（以下、地震リスク）を評価すること

告示1347号1） に示す、「地盤の長期許容応力度に応じてができる。暮らしている地域の地震・地盤環境の理解に

基礎の種類を選ぶ」ことや、常時微動探査による地盤の役立つとともに、将来的な耐震診断の促進に資するもの

卓越周期から地盤種別の推定を行い増幅特性として準用と考える。しかしながら、より詳細で耐震診断および設することが可能（国交省告示1113号 ）となっている等、計に使用できるような地盤情報に関する公開は不十分な

地盤調査が必要となってきている。一方で、地震被害を状況である。本論文では、近年検討されている地盤情報

10年で半減させるという地震防災戦略を2005年３月に策に関わる研究・開発や地下構造データベース、地盤構造

定し、さらに2020年までに耐震化率を95％にするというモデル作成・微動システムの研究および今後の新たな展

目標が新成長戦略として閣議決定され、住宅をはじめ地開等を交えて報告する。

2）

盤を考慮した建物の耐震化は喫緊の課題となっている。

2．これまでの地盤情報の現状

しかしながら、地盤情報に基づく建物の耐震化をはじめ、耐震診断手法は十分に進んでいないのが現状であ

地下構造の情報がきわめて重要であることは前述にもる。地盤や揺れやすさの情報は、市民にとって関心の高述べたが、本章では、これまでの地盤情報に関する取りい事柄であり、自宅の耐震性等について考える機会を捉組みの例として、防災科研による強震動予測のための地

えて地震に関する知識を提供すれば、効果的に知識の普盤構造モデル作成と、地下構造データベースの構築につ及と地震防災に対する意識の啓発につながると考える。

いて述べる。

防災科研では、国の地震調査研究推進本部が作成・公表
している全国地震動予測地図のデータを地震ハザードス
テーション（J-SHIS：http://www.j-shis.bosai.go.jp）3）＊つくば事業部 第四技術部
1

MSS技報・Vol.24

2.1

強震動予測のための地盤構造モデルの作成

ル」で、主に地質情報から構築された防災科研による（全国一次地下構造モデル）「全国初期モデル」6） であり、次は「0.5次地下構造モデ

工学的基盤以深の高精度な三次元地下構造モデル（深ル」で、2005年版「全国を概観した地震動予測地図」に

部地盤構造モデル）として、各種探査や地震観測の結果

おいてシナリオ地震予測が行われた地域などである程度を総合した全国１次地下構造モデル作成の取り組みがあ

高精度化されたモデル、３段階目の「１次地下構造モデ

る。モデル構築に際しては、Koketsu et al.（2009） にル」は、防災科研が平成19年度に構築した九州地域全域

より提案された、わが国における地域規模の三次元地下のモデル、大都市大震災軽減化特別プロジェクト7）や重構造（速度構造）の標準的なモデル化手法（図１）に従点的調査観測の一環として首都圏・近畿圏や宮城県域なって、屈折法／反射法地震探査、重力探査、表層地質、どで構築された地盤モデルを取り込んで作成されたもの

ボーリング、微動探査等のさまざまなデータを利用し、である。なお、この全国１次地下構造モデルは、現在、面的調査データや自然地震観測記録等により調整し、さ長周期地震動予測地図作成に使用されている。一方、工らに地震動の観測記録とシミュレーション結果を比較す

学的基盤以浅の浅部地盤構造モデルは広域的には作成さることにより調整されている。この手法を首都圏の地下れておらず、推本における「警固断層帯（南東部）の強構造に適用して構築された参照速度構造モデル は、主震動評価」で実施された強震動計算8）において、福岡市に長周期地震動予測地図作成のために用いられ、手法のの一部地域のみで評価が実施されている状態でしかな妥当性が確認されている。なお、全国地下構造モデルの

い。

4）

5）

構築の段階は３段階あり、１つめは「0次地下構造モデ

図１

強震動評価のための深部地盤構造モデル作成の流れ

2

MSS技報・Vol.24

2.2

地下構造データベースの構築

3．現状の地震動予測のための地盤モデルの作成について

日本では、地震災害をはじめ各種自然災害によるリス
クが高く、その対策を進めるために、基礎資料となる地

強震動予測を行うために重要な情報としては、図３に下構造に関するデータの利用に対するニーズが高くなっ示すように、震源特性・伝播経路特性・地盤増幅（サイている。

ト）特性の３つの特性が重要であるが、中でも地盤増幅

過去には膨大な地下構造調査が実施されてきている。

特性については、防災上重要な周期特性や増幅特性（長しかしながら、これらの調査は、様々な目的に基づいて

周期地震動および地上における震度予測等）を評価する

様々な機関により実施されているため、限定された目的ために最も重要な情報である。現状公開されている地震の利用に留まり、データは各府省・地方公共団体・関係

動予測地図（震源を特定した地震動予測地図）では、図機関等に散在している。これらのデータは他には十分利

３に示す工学的基盤まで（深部地盤構造）をモデル化

用されないまま死蔵の状態となっているものが多く、中し、その上面までの詳細法による地震波形計算10）を行っには散逸の危機にあるものもある。データの散逸を防ている。一方、浅部地盤構造については、地震動計算用

ぎ、誰でもが利用可能なデータベースを構築することはの速度構造を作成することなく、約250mメッシュ単位の大変重要な課題となっている。

微地形区分11） に基づいたAVS30による地盤増幅率12） で、こうした背景の下、平成18年７月より5カ年間に渡り、工学的基盤の最大速度に増幅率をかけることで、地上の科学技術振興調整費重要課題解決型研究「統合化地下構最大速度（震度）としている。しかしながら、地盤・耐造データベースの構築」 が開始され、防災科研をはじ

震工学の分野での利用を考えると、地表の地震動におけめ６つの研究機関が参加して、地震防災など社会の安る周波数・増幅特性を求める必要がある。以上のことか

全・安心と持続的発展を実現するため、国土の地下構造

ら、国においても地表までの地盤構造モデルを作成するに関する情報を「国民共有の公的財産」と位置づけ、収必要性が出てきている。推本においては、地盤構造モデ

集・整備し、表層地盤から深部に至る系統的・網羅的なル高度化への取り組みが平成23年度から開始し、0.1秒

統合化地下構造データベースを構築することを目指しから10秒程度の広帯域の地震動特性を評価できるようなた。検討では、データの標準化や地質と物性特性など異地盤モデルの構築を念頭に、これまで別々にモデル化をなるデータ間の統合化を進め、各機関に分散したデータ

実施してきた浅部地盤モデル13）と深部地盤モデル5）を統

ベースの相互利用・公開を進めるための分散管理型ネッ

合し、地震・常時微動等の観測記録を再現できるような

トワークシステムを開発してきた。さらに、地震防災にモデルの作成を進めている（図４）。

9）資する地下構造モデルの高度化やハザードマップ・リス

ちなみに、既往の研究では地表まで作成された３次元

ク評価への利活用について研究し、データベースの内容の地盤モデルにおいての地震動の検証が十分になされての高度化や、その有効性を示すため、利活用に関する検いる状況にはない。上記の理由から、防災科研では地盤討を具体的な事例に基づいて実施している（統合化のイ

構造モデルの高度化の研究として、まずは、国・自治体

メージ：図２）。

等のボーリングデータの収集し、初期の地質構造モデル
を地表から工学的基盤までを作成し、既往の深部地盤モ
デルを結合して「初期地盤モデル」を作成している。次に、このモデルを初期値として、K-NET、KiK-net、気象庁および自治体等の地震観測点における地震記録と、多数の常時微動探査の記録より、地盤のS波速度構造・増幅特性（スペクトル増幅率）等を求め、地盤モデルの高度化を行っている。また、面的な補間方法の検討も行い、最終的には約250mメッシュ単位の地盤構造モデル（通称：浅部・深部統合地盤モデル）14）,15）を作成している（図５）
なお、地盤モデル作成結果の検証については、周期２秒よりも短周期側については、１次元重複反射法による
図２

地震観測記録とサイト増幅特性等との検証を行い、周期統合化のイメージ

２秒よりも長周期側については、有限差分法16）による地
震観測記録との結果の検証を行い、既往の地盤モデルに

3

MSS技報・Vol.24

図３

地盤構造のＳ波速度と各構造モデルの関係概念図

図３（補足）

図４

各地盤モデル作成に関する対象周期との関係概念図

浅部・深部統合地盤モデル作成概念図

4

MSS技報・Vol.24

図５

浅部・深部統合地盤モデル作成結果（南関東地域暫定結果）

対しての高精度化の検証を行っている。この検討によにコストがかかる。一方で、近年、直接的な地震観測以

り、日本の広帯域地震動予測用のスタンダードな地盤モ

外で地震時の揺れの周期特性や地盤の速度構造を良くとデルを構築・公開していくことで、将来的に本検討の地

らえることのできる常時微動観測が行われてきている。

盤モデルを基本として新たな地盤モデルの研究や民間企

常時微動観測は、地震動予測のための地盤構造モデル作業が利用できる地盤情報としての利活用を考えており、成に非常に有用な物理探査手法である。その観測の手軽

本研究の最終的な目標としている。

さにより、ここ10年で各大学・研究期間では非常に多くの研究がなされているたとえば17）。しかしながら、研究その

4．地盤モデル作成のための微動観測システムの構築にものはバラバラの機材・考え方により各大学の小さな研ついて究素材としてとじられており、データの共有等について前節に示すように、地震動予測地図の高度化のためにも、ボーリングデータ以上に集約の難しい情報となっては、高精度な地盤モデルの作成は必要不可欠であるが、しまっている。ボーリングデータはその物性の性質（堅そのためのボーリング調査や反射法地震波探査は、莫大

さ柔らかさ）や層序に関しての情報は抽出できるもの

図６

常時微動システムの概要

5

MSS技報・Vol.24

の、地震動にとって重要な一部の物性値構造は推定で与観測を実施したことがない等、大部分の研究機関は観測えるしかない。動的な特性を得るには常時微動やPS検を行いたい希望はあっても、二の足を踏み研究活動が出層・表面波探査のような観測が必要であるが、特にPS

来ないところが多い。全国的にくまなく常時微動データ

検層や表面波探査についてはコストがかかる上に観測者の集約をはかっていくには、このような地域の研究者にの能力や解析結果の解釈も難しく多くの観測は望めな

機材貸与および技術協力を行いデータ集約の仕組みに賛い。

同していただき、観測者の裾野の拡大を行って「ビック常時微動観測は現在、先名・他（2006,2008）のデータ」収集の足がかりにしていくことが、新たな精度微動観測システムの開発により、素人でも数時間のレクの高い全国ハザードマップ（ゆれやすさマップ）を作成

チャーを受けることで、熟練者の観測結果と遜色ない精

していく礎になると考える（図７）。

18）,19）度で観測できるようになった。単点の微動観測について

そのための方法（筆者の提案）として、防災科研がデは15分程度の観測で、その地点の揺れやすさの周期特性

ータセンターとしての機能を持ち、微動計の管理・デー（Tg）・増幅特性等がわかり観測データそのものもデー

タベース登録等システム・教育システムをパッケージとタベースに即時登録が可能なシステム「常時微動システ

して用意・提供し、機材・データの品質・情報管理を効

ム（先名・他（2012）） 」も開発されている（図６）。

率良く行っていけるものと考える。現在、上記微動シス

このシステムは、現場においてデータ登録・解析等が完テムを地方大学・研究機関等の希望の部署に観測レクチ

了するため、位置情報・写真情報・観測データ・時刻歴

ャー教育ともに提供し、大きな成果を得られつつある。

情報・解析結果等の登録時のヒューマンエラーが少なく先々収集される「ビックデータ」の管理を効率化するたなる仕組みとなっている。

め、今後において防災科研の微動システムを大規模に再20）構築する必要が出てくるであろう。

一方、常時微動観測は、研究自体は有名大学を中心に沢山実施されているものの、機材の高価さ・これまでに

図７

防災科研が目指す相互運用システムの仕組み

6

MSS技報・Vol.24

5．新しい微動探査法の構築と地盤情報サービスシステムの開発への取り組み
近年、筆者らが取り組んでいる新しい微動観測手法および地盤情報サービスシステムについて下記に簡単に紹
介する。
5.1

新しい微動探査法の活用（極小アレイによるＳ波速度断面の構築）
筆者らは、現在、新しい微動観測・解析手法を検討している。１つは、微動アレイ解析方法の１つである、
CCA 法（Cho et al.（2006））21）の手法をベースに、観測から解析までを自動的に結果を算出する（位相速度を
求 め、S波 速 度 構 造 に 変 換） よ う に 検 討 し た、Cho,
Senna and Fujiwara（2013）22），23）の手法である。アレ
イ半径に対する解析可能波長帯域、解析結果の信頼性、安定性の観点から、コヒーレンスCCA 法（nc-CCA法）
による極小アレイ（r＜１m）と、必要に応じて３点からなる不規則形の「小アレイ」（r＜5m）を別途実施することで従来よりも実用的な浅部（深さ数10m〜100 m 程度）構造探査が可能になるとしている（図８）。極小ア
レイは観測が簡単なため、測線に沿う複数点で観測して「チェーンアレイ」（原口・他、2010）24） と呼ばれるようなS波速度構造を2次元断面にて構造解析する手法も実施されている。この解析方法は、適当な変換式によって分散データをS波速度構造に変換する簡易逆解析である（図９）。もう１つは、H/Vスペクトル25）を利用した手法

図９である。過去の研究は微動の振幅情報が地下構造推定に

図８

微動観測結果から簡易逆解析を行う新手法（位相速度簡易変換・トラフ法・ピーク法）

極小アレイ（半径60cm）による各解析結果の比較（nc-CCAの位相速度が最新の研究によるもの）

7

MSS技報・Vol.24

有効なことを表しており、スペクトルのピーク周波数は（上図はインバージョン無し・中図はインバージョン有速度コントラストの強い軟弱層の厚さや速度に、ピークり）である。上記の結果は地盤図の地質構造断面（下高は層間の速度コントラストあるいはインピーダンス比図）と良く合っており、S波速度値も信頼できる値であに 関 連 す る と 考 え ら れ て い る（例 え ば、Konno andると考える。精度・有効性については、今後さらに他のOhmachi, 1998）26）。しかしH/V スペクトル比で地盤構地域についても検討し、研究論文として公表予定である。

造を推定しようとするすべての既存研究には、スペクト
5.2

ルのピークによって示される深さスケールの解釈には地新しい探査手法を利用した地盤情報サービス
システムの開発の取り組み

盤速度と層厚のトレードオフに起因する任意性があった
り、ピークが低く複数あるためにピークの読み取り自体

5.1節で示した手法の利用は、当初土木探査を想定しが難しい場合が多い等の問題が常につきまとってきた。

ていたが、高い周波数100Hz以上でも比較的S/Nが確保

本研究では３成分地震計で極小アレイと同時に実施するできるため、非常に浅い構造（深さ50cm〜1m程度）にことでトレードオフを解決している。

も対応出来る。現在、2000年の法改正により、住宅建設

浅部構造探査のために３成分単点観測とCCA法によ時に地盤調査が必須となり、SS（スウェーデン・サウ

る極小アレイ、小アレイを多点で展開して疑似S波速度

ンディング）方式27） による地耐力調査が行われている

断面とH/V断面（それぞれ疑似S波速度構造、H/V構造

が、その後、種々の調査方法（水平載荷試験・物理探の空間補間で得られる２次元プロット）を解析し、統合査）で大きく成長する地盤調査会社が出る市場となって的に解釈することで、これまでよりも多くの情報を簡易いる。本微動観測システムによる観測の簡便性や、解析に抽出できるようになった。CCA法による極小アレイ・の自動化による簡易性より、市場展開の可能性が十分に小アレイ観測は作業コストが低いので多点データの取得あるものと考える。更に、2011年東北地方太平洋沖地震に適する。またアレイサイズが小さいので、水平方向のによりハウスメーカーの調査方法に関する意識が変わり

空間分解能が高い。防災科学技術研究所と白山工業つつあり、地耐力だけでなく地震の卓越周期特性（Tg）（株）の共同開発による微動／強震観測キットJU215

や、観測地点で推定される応答スペクトルのニーズも出（先名・他、2006）

は、この目的に適している。図10てきている状況である。S波速度構造を比較的深い部分

は Ｉ 地 域 の 約 ３kmの 測 線 に 沿 う 疑 似S波 速 度 断 面 図まで得られる本手法は、基本的な要求を満たしつつ、さ

18）

図10

極小アレイによる解析結果の２次元断面図と地質断面図（下図）との比較（Ｉ地域）※真ん中のカラーの図が極小アレイによるＳ波速度構造断面図（上図は逆解析なし（位相速度簡易変換のみ）
・中図はピーク・トラフ法を使用した逆解析有り）

8

MSS技報・Vol.24

図11

地盤評価システムにおける結果のイメージ

らに他の手法と比較しても優位な結果が得られ、その結参考文献

果に基づいた新たなサービスの展開も可能となるものと考える。例えば、前述の浅部・深部統合地盤モデルと合

１）平成13年国土交通省告示第1347号（評価方法基準），

わせて利用することで、調査地点の地震基盤までの構造2001.

に基づいて長周期を含む広帯域の揺れの特性や、周辺と

２）平成13年国土交通省告示第1113号（基礎及び地盤の比べてどの程度揺れやすいのかという評価結果もサービ

耐震設計用許容応力等の算出方法），2001.

スすることが可能となる（図11）。また、防災科研の３）独立行政法人防災科学技術研究所：地震ハザードス

J-SHIS等に収められているデータとの相互運用を行うこ

テ ー シ ョ ン（J-SHIS）：http://www.j-shis.bosai.

とで、種々の地震動予測地図を利用したハザード情報もgo.jp.

出すことが可能となり、最先端で高付加価値のリスク情４）Koketsu, K., H. Miyake, Afnimar, and Y. Tanaka,

報サービスの提供も可能となる。

2009，A proposal for a standard procedure ofm o d e l i n g 3- D v e l o c i t y s t r u c t u r e s a n d i t s

6．まとめ

application to the Tokyo metropolitan area, Japan,

過去13年間、常駐者・契約研究員・客員研究員としてTectonophysics, 472，290〜300.

関わってきた防災科研研究業務における地盤構造に関す

５）纐纈一起：地下構造と長周期地震動，日本地震工学る取り組みについて、これまでの経緯に沿って記載し

会誌，No. 4，pp. 12-15，2006.

た。防災科研での地盤構造に関する取り組みについて

６）藤原広行・河合伸一・青井 真・森川信之・ 先名重は、通常は地震動予測のための調査・研究が一般的であ樹・工藤暢章・大井昌弘・はお憲生・早川 讓・遠るが、それ以外の用途も考慮した地下構造データのデー

山信彦・松山尚典・岩本鋼司・鈴木晴彦・劉瑛：

タベース化や、新しい探査手法の研究・開発等、その内強震動評価のための全国深部地盤構造モデル作成

容・ボリュームともに拡張しつつ、府省・官民連携によ手法の検討，防災科学技術研究所研究資料, 337，

る総合防災を意識した取り組みに大きく変化しつつあ2009.

る。その中で、現在の手法を駆使して、新しい地盤情報

７）Miyake, H., K. Koketsu, and T. Furumura, Source

サービスを提供していくことは、民間企業にとっては大modeling of subduction-zone earthquakes andきなチャンスであると考える。

long-period ground motion validation in the Tokyo

本稿をまとめるに当たり、原稿段階で査読を御願いし

metropolitan area, Eos Trans. AGU，Vol. 88，No.

た関係各機関の各位には御礼申し上げます。

52，Fall Meet. Suppl., Abstract S14C-04，2007.
８）文部科学省地震調査研究推進本部地震調査委員会強震動評価部会：警固断層（南東部）の地震を想定

9

MSS技報・Vol.24

した強震動評価について，http://www.jishin.go.jp/

18）先名重樹, 安達繁樹, 荒木恒彦, 飯澤清典, 藤原広行：

main/，2008．

微動探査観測システムの開発, 第115回物理探査学

９）防災科学技術研究所：統合化地下構造データベース

会予稿集, pp.120-122, 2006.

の構築＜地下構造データベース構築ワーキンググ

19）先名重樹, 藤原広行：微動探査システムの開発

そループ報告書＞，防災科学技術研究所研究資料,の１−常時微動解析ツールの開発−, 防災科学技術

361, 2011.

研究所研究資料, 第313号, 2008.

10）先名 重樹・藤原 広行・河合 伸一・青井 真・功刀

20）先名 重樹・安達 繁樹・浅香 雄太・藤原 広行：新卓・石井 透・早川 讓・森川 信之・本多 亮・小林

微動探査システムの開発，第126回物理探査学会予京子・大井 昌弘・八十島 裕・神野 達夫・奥村 直

稿集，pp.126-129，2012.

子：森本・富樫断層帯の地震を想定した地震動予21）Cho, I., T. Tada, and Y. Shinozaki（2006）, New

測地図作成手法の検討，防災科学技術研究所研究

method of microtremor exploration: the centerless

資料，第255号．

circular array method and the two-radius method,11）若松加寿江・松岡昌志・杉浦正美・久保純子・長谷3rd International Symposium on the Effects of

川浩一：日本の地形・地盤デジタルマップ ，東京

Surface Geology on Seismic Motion, Grenoble,

大学出版会，2005.11．

France, Paper Number 51.

12）松岡昌志・若松加寿江・藤本一雄・翠川三郎：日本

22）Cho, I., S. Senna, and H. Fujiwara（2013）,

全国地形・地盤分類メッシュマップを利用した地

Miniature array analysis, Geophysics, 78, KS13−

盤 の 平 均 Ｓ 波 速 度 分 布 の 推 定,土 木 学 会 論 文 集

KS23, 10.1190/GEO2012-0248.1.

No.794/I-72, 2005, pp.239-251.

23）Cho, I., S. Senna, and H. Fujiwara（2014）, Depth13）三宅弘恵・纐纈一起・古村孝志・稲垣賢亮・増田conversion of the H/V spectrum: A simple

徹・翠川三郎：首都圏の強震動予測のための浅層

method for acquiring information on the S wave

地盤構造モデルの構築，第12回日本地震工学シン

vel oci ty di sconti nui ty f rom mi crotre m o r s ,ポジウム，2006, pp.318-321.

Geological Survey of Japan（now printing）.

14） 先 名 重 樹・ 前 田 宜 浩・ 稲 垣 賢 亮・ 鈴 木 晴 彦・ 神

24）原口強、松岡達郎、南雄一郎、小野雅弘：チェーン

薫・宮本賢治・松山尚典・森川信之・河合伸一・

アレー微動探査法による上町断層地下構造のイメ

藤原広行：強震動評価のための千葉県・茨城県にージング、日本応用地質学会平成22年度研究発表

おける浅部・深部統合地盤モデルの検討，防災科会講演論文集、2010．

学技術研究所研究資料，第370号．

25）中村豊：常時微動計測に基づく表層地盤の地震動特

15）Senna, S., T. Maeda, Y. Inagaki, H. Suzuki, N.

性の推定, 鉄道総研報告, Vol.2, No.4, pp.18-27, 1988.

Matsuyama, and H. Fujiwara（2013b）, Modeling26）Konno, K., and T. Ohmachi（1998）, Ground-of the subsurface structure from the seismic

motion characteristics estimated from spectral

bedrock to the ground surface for a broadband

ratio between horizontal and vertical components

strong motion evaluation, submitted to Journal ofof microtremor, Bull. Seismol. Soc. Am., 88, 228-

Disaster, 2013, Vol8, No.5, pp.889-903.

241.

16）Aoi, S. and H. Fujiwara,3-D finite-difference27）日本工業規格：スウェーデン式サウンディング試験

method using discontinuous grids,Bull. Seismol.

方 法，http://www.jiban.or.jp/file/organi/bu/

Soc. Am., 89, 918-930, 1999.

kijyunbu/kouji/jis̲a̲1221̲kai̲201003 kaishikeisai.

17）Arai, H., and K. Tokimatsu（2004）, S-Wave

pdf

velocity profiling by inversion of microtremor H/
V Spectrum, Bull. Seismol. Soc. Am., 94, 53-63.

10

MSS技報・Vol.24