1

日本情報地質学会シンポジウム2013 講演論文集, 1-19, 2013
地質情報等の三次元モデリングとＣＩＭについて

物理探査による河川堤防の 3 次元可視化とモデリング高橋 亨*

3D imaging and modeling of a river embankment with geophysical data
Toru Takahashi**公益財団法人深田地質研究所 Fukada Geological Institute,2-13-12 Honkomagome, Bunkyo-ku, Tokyo 113-0021, Japan. E-mail: takahashi@fgi.or.jp

キーワード：物理探査，河川堤防，３次元
Key words : Geophysical method, River embankment, Three dimensions

2

話の内容
河川堤防の被害
河川堤防調査と物理探査の役割
河川堤防の可視化（概要調査での利用）
河川堤防の可視化（詳細調査における３次元可視化）
物理探査データを用いた河川堤防のモデリングまとめ・近年，地震や局地的な集中豪雨等による河川堤防の被害の発生例が多く、その強靭化が求められている。
・特に、規模は小さくても将来の破堤につながるような前兆的な堤体の変状（漏水、亀裂等）の報告は極めて多い。

昨年7月の九州北部豪雨で決壊した矢部川－冠水した九州新幹線の筑後船小屋駅付近の状況－（NHK2012)3

国土交通省(2012)5国土交通省九州地方整備局筑後川河川事務所(2012)4国土交通省 (2012)

浸透が問題となる堤体と基礎地盤の土質構成

漏水被害の事例（国土技術研究センター、2002)（国土交通省近畿地方整備局、2006)

7

5スタート

河川堤防調査の流れ－平常時あるいは予防的な調査の流れ－

基礎調査（縦断方向）
自然・社会条件調査

堤防現況調査構造物現況調査

被災履歴調査堤防諸元等調査（堤防縦断）
土質調査（堤防縦断）
物理探査（堤防縦断）

一連区間の設定
一連区間の細分
浸透・浸食・地震に対する概略評価

概要調査

代表断面の選定

詳細調査

代表断面における詳細調査
土質調査（堤防横断）
物理探査（堤防横断）

代表断面の安全性照査

照査の基準

OK終了

NO
危険要因の分析

堤防強化対策設計・施工

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)終了

河川堤防調査の流れ

8

被害発生－被害発生時およびその後の調査の流れ－

堤防および構造物の被害調査

対策工設計のための緊急調査
土質調査（被害地点周辺重点調査）
物理探査（被害地点周辺詳細探査）

堤防の被害原因・メカニズムの解明

対策工の設計・施工
詳細調査
概要調査

被災堤防と類似堤防の選定のための広域調査
資料・文献調査
現地踏査
堤防・基礎地盤の調査
土質調査（堤防縦断）
物理探査（堤防縦断）対策の検討が必要な堤防の選定

平常時と同様の詳細調査と安全性照査（必要な場合は強化対策）

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

終了

96

河川堤防の調査対象と要件
調査対象
堤体

縦断方向

具体的な事例・堤防からの漏水・堤防表面の亀裂調査で把握すべき事項と深度・土質構成の連続性（弱点箇所の特定）

0-9m・堤体内異物の有無（亀裂、陥没、空洞、緩み、廃棄樋門など）

0-9m

横断方向

・堤防からの漏水

・土質構成
・築堤履歴

0-9m

堤防縦断方向・矢部川での堤防決壊

・要注意地形の分布範囲（旧河道、旧砂州など）・透水層の分布、層厚・液状化層の分布、層厚15m

堤防横断方向・矢部川での堤防決壊
・パイピング(漏水）・土質構成の連続性（特に透水層の連続性）
・液状化層の分布、層厚15m

河川構造物周辺・護岸破損
・樋門管内（底版下）の空洞・樋門等の構造物周辺におけるゆるみや空洞の有無・規模

0-9m

護岸・護岸背面空洞・護岸背面の空洞の有無・規模2m

基礎地盤

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

調査対象・段階に応じた探査法
調査段階

概要調査

詳細調査

要求事項

堤防を縦断方向に効率的(早く安く)に探査し、堤 ボーリングデータを補完して堤防横断面の詳体、基礎地盤の概略構造や土質構成および弱部、 細な構造や土質構成を推定する異常部（旧河道や空洞、緩み領域等）の概略位置を推定する

探査基本
仕様

・ベルト状の発受信機による牽引探査・リアルタイムに近いスピード解析
・成果：概略の構造や土質構成

・対象断面での横断あるいは３次元探査
・詳細解析
・成果：横断面の構造や土質構成

探査手法

・ランドストリーマ牽引型反射法
・ランドストリーマ牽引型表面波法
・牽引型電気探査・浅部電磁法探査
・地中レーダ探査
・複合探査（上記データの複合解析）

・３次元電気探査（２次元を含む）
・屈折法トモグラフィ
・複合探査（上記データの複合解析）

概念図

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)7

調査対象・段階に応じた探査法
調査段階

概要調査

詳細調査

要求事項

堤防を縦断方向に効率的(早く安く)に探査し、堤 ボーリングデータを補完して堤防横断面の詳体、基礎地盤の概略構造や土質構成および弱部、 細な構造や土質構成を推定する異常部（旧河道や空洞、緩み領域等）の概略位置を推定する

探査基本仕様・ベルト状の発受信機による牽引探査・リアルタイムに近いスピード解析
・成果：概略の構造や土質構成

・対象断面での横断あるいは３次元探査
・詳細解析
・成果：横断面の構造や土質構成

探査手法

・ランドストリーマ牽引型反射法
・ランドストリーマ牽引型表面波法
・牽引型電気探査・浅部電磁法探査
・地中レーダ探査
・複合探査（上記データの複合解析）

・３次元電気探査（２次元を含む）
・屈折法トモグラフィ
・複合探査（上記データの複合解析）

概念図

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

河川堤防の可視化（概要調査での利用）-宇治川での探査例-

測線を下流側（南側）から見た風景

宇治川左岸「旧巨椋池排水機場樋門」

100m

測線配置と土質調査位置（実線：概査測線500m，矩形範囲：精査範囲）

測線位置 国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)8

河川堤防の可視化（概要調査での利用）-宇治川での探査例-

表面波探査（S波速度分布）と牽引式電気探査（比抵抗分布）を実施し、複合的にデータを解析して、堤防のモデル化を行った。

探査範囲
国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

表面波探査の探査原理と解析法
震源

波長の短い波波長の長い波

表面波の分散特性
Phase-velocity(km/s)
21.9
1.8
1.7
1.6
1.5
1.4
1.3
1.2
1.1
10.9
0.8
0.7
0.6
0.5

Obserevd
Calculated

010

2030

Frequency(Hz)

4050

9

表面波探査の測定方法（牽引型）電気探査の探査原理
交流電源+ + + +- - - 絶縁体

- - - -+ + + +
地盤

電荷の移動

10

牽引型電気探査の測定方法

河川堤防の可視化（概要調査での利用）-宇治川での探査例0深度（ｍ）

ボーリングで推定された
土質断面図12深度（ｍ）

0

表面波探査で得られた
Ｓ波速度断面120

深度（ｍ）

牽引型電気探査により得られた比抵抗断面

12500m

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)11

調査対象・段階に応じた探査法
調査段階

概要調査

詳細調査

要求事項

堤防を縦断方向に効率的(早く安く)に探査し、堤 ボーリングデータを補完して堤防横断面の詳体、基礎地盤の概略構造や土質構成および弱部、 細な構造や土質構成を推定する異常部（旧河道や空洞、緩み領域等）の概略位置を推定する

探査基本仕様・ベルト状の発受信機による牽引探査・リアルタイムに近いスピード解析
・成果：概略の構造や土質構成

・対象断面での横断あるいは３次元探査
・詳細解析
・成果：横断面の構造や土質構成

探査手法

・ランドストリーマ牽引型反射法
・ランドストリーマ牽引型表面波法
・牽引型電気探査・浅部電磁法探査
・地中レーダ探査
・複合探査（上記データの複合解析）

・３次元電気探査（２次元を含む）
・屈折法トモグラフィ・複合探査（上記データの複合解析）

概念図

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

河川堤防の可視化（詳細調査での利用）-３次元電気探査の測定方法-

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

12河川堤防の可視化（詳細調査での利用）-３次元電気探査の探査結果-

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

河川堤防の可視化（詳細調査での利用）-屈折トモグラフィの測定方法-

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

13河川堤防の可視化（詳細調査での利用）-屈折トモグラフィの探査結果-

Ｓ波速度断面

Ｐ波速度断面

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

センサー
震源

３次元堤防可視化装置の開発
電探器

弾性波探査器

PC名称

3次元可視化装置仕様
ベルトスリング型ランドストリーマ
70mm×50m
12測点/本受振器ベース＆電極一体型
測点間隔4m（最大）上下動受振器・水平動受振器 交互配置
受振器：MarkProduct L40-A
電極：ステンレス製

備

考

4セット
(8セット予定)

高橋ほか (2012)14

３次元堤防可視化装置の評価試験

旧堤防
現在の堤防

評価試験作業状況
３次元電気探査解析結果（比抵抗分布図）
高橋ほか (2012)

物理探査データを用いた河川堤防のモデリング－物理探査で得られる物性と河川堤防調査で必要な特性の関係－

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

15

物理探査データを用いた河川堤防のモデリング－特性間の相関関係、経験式を利用する方法－

Ｓ波速度とＮ値の相関関係をもとに、表面波探査によるＳ波速度からＮ値断面を推定した結果

比誘電率（電磁波速度）と体積含水率の経験式をもとに、地中レーダ探査による電磁波速度から体積含水率を推定した結果
国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

物理探査データを用いた河川堤防のモデリング－物理モデルを利用する方法－

0深度（ｍ）

ボーリングで推定された
土質断面図

120深度（ｍ）

表面波探査で得られた
Ｓ波速度断面

120深度（ｍ）

牽引型電気探査により得られた比抵抗断面

12500m

国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)

16

物理探査データを用いた河川堤防のモデリング－物理モデルを利用する方法－

350

K2/G dryZ

G2/ 0( 4 / 3) G 2

1
Ks

1/ 0Gs Z0

Z/ 0( 4 / 3) G 2

14G231

Z

Hashin-Shtrikman ModelG 2 9 K 2 8G 2
6 K 2 2G 22K2n 2 (1 0 ) 2 GsP
18 2 (1 s ) 2G25 4 s5(2 s )S波速度 (m/s)

Unconsolidated sand model
K dry

Model (C=0%)
Molde (C=30%)
Model (C=60%)
Model (C=100%)
Clayey gravel
Sand
Clayey sand

400

Sand
Clay

300
250
200
150

1/ 3

23n 2 (1 0 ) 2 GsP
2 2 (1 s )2

Hertz-Mindlin Model

1/ 310020Ks{(1 C ) K sand

CK clay 1 /((1 C ) / K sandC / K clay )} / 2{(1 C )Gsand

CGclay 1 /((1 C ) / Gsand

C / Gclay )} / 2ｎφ0Kclay
Gclay
Ksand
GsandHill average

K dry

G sat

Gdry

(1 K dry / K1 ) 2/Kf(1) / K1

4050607080

間隙率 (%)

GsK sat30

K dry / K 1

250.8
21 GPa7 GPa36.6 GPa
45GPa

Coordination number
Critical porosity
Bulk modulus of clay
Shear modulus of clay
Bulk modulus of Quartz
Shear modulus of Quartz

Gassmann’s Relation

Takahashi & Yamamoto (2010)

物理探査データを用いた河川堤防のモデリング－物理モデルを利用する方法－

300

Model (C=0%)
Molde
Molde(C=30%)(C=30%)
Model(C=60%)
Model(C=60%)
Model(C=100%)
Model(C=100%)
clayeygravel(RI-1)
clayly sand(RI-1)
sand(RI-1)
sand(RI-1)
clayey sand(RI-1)
clayey gravel(RI-1)

350
350
S波速度
S波速度(m/s)
(m/s)

350

250
200

300
300
250
250

Reclaimed
soil

200
200

150

150
150

100

100100

Sand/Clayey gravel
Sand/Clay Mixture
Clayey sand2030

405060
間隙率 (%)

7000

80

50
50

100 150150 200
200 250250 300
300 350350 400
400100比抵抗 (ohm-m)(ohm-m)
比抵抗RI-1
RI-2
RI-3
RI-4300
比抵抗 (ohm-m)

S波速度 (m/s)

400
400

Model (C=0%)
Molde (C=30%)
Model (C=60%)
Model (C=100%)
Clayey gravelSand
Clayey sand

400

250
200
150
100
50304050
間隙率 (%)

6070

Takahashi & Yamamoto (2010)

17

物理探査データを用いた河川堤防のモデリング－物理モデルを利用する方法により推定された土質断面－400

S波速度 (m/s)

Molde (C=30%)

350

Model (C=60%)
Model (C=100%)

300

clayly sand(RI-1)
sand(RI-1)
clayey gravel(RI-1)

250

Reclaimed
soil

200

Sand/Clayey gravel
Sand/Clay Mixture

150

Clayey sand

100
050

100

150

200

250

比抵抗 (ohm-m)

300

350

400

Takahashi & Yamamoto (2010)

物理探査データを用いた河川堤防のモデリング－物理モデルを利用する方法により推定された土質断面－

S波速度 (m/s)

400

Molde (C=30%)

350

Model (C=60%)
Model (C=100%)

300

clayly sand(RI-1)
sand(RI-1)
clayey gravel(RI-1)

250

Reclaimed
soil

200

Sand/Clayey gravel
Sand/Clay Mixture

150

Clayey sand

100
050

100

150

200250比抵抗 (ohm-m)

300

350

400

33
Takahashi & Yamamoto (2010)

18

物理探査データを用いた河川堤防のモデリング－Ｓ波速度と比抵抗を利用した透水係数の推定法－
探査データから透水係数の推定（飽和地盤での適用）（Ｋｏｓｅｎｙ-Ｃａｒｍａｎの式＋ Rock Physicsに基づく方法）

砂礫：1mm
砂 ：0.1mm
粘土：0.001mm

3B(1)2d2

高橋・稲崎 (2013)

物理探査データを用いた河川堤防のモデリング－複数の物理探査データを利用した河川堤防の３次元モデリング－

土質構成、間隙比、透水係数、強度等の３次元画像

３次元物理探査結果

堤防特性の３次元モデル

Takahashi et al. (2013)を修正

19

まとめ
河川堤防内部を非破壊的に可視化する技術として物理探査が利用されている。
特に、横断を対象にした詳細調査において３次元探査が適用され、探査物性による３次元可視化が行われている。
３次元探査を効率的に実施するための複合探査装置の開発や探査で得られた弾性波速度や比抵抗値を複合的に解析し、堤防の土質構成、さらには透水係数等のモデル化も試みられている。
近い将来、複数の３次元物理探査による河川堤防の３次元モデル化が可能となると考えられる。

謝辞
本講演で使用させていただいたデータの多くは，以下の研究で取得，解析されたデータです．関係者の皆様にここに記して感謝いたします。
京都大学河川堤防の内部構造調査と探査機器の開発研究（2006年～2008年）
国土交通省近畿地方整備局近畿技術事務所・京都大学河川堤防の内部構造調査と探査機器の開発研究委員会 (2009)；河川堤防調
ける物理探査適用の手引き-新しい技術の適用を中心に-.

査におTakahashi T. and Yamamoto T., 2010, An attempt at soil profiling on a river embankment using geophysical data, Exploration Geophysics, 41,
102-108.

国土交通省建設技術開発助成事業－被災堤防緊急対応のための３次元堤防可視化ツール及び対策工設計システムの開発－（実施中）高橋亨・相澤隆生・西尾英貴・松岡俊文(2012)；被災堤防緊急対応のための３次元可視化ツール及び対策設計支援システムの開発－研究開発の全体計画と2011年度成果の概要－、物理探査学会第127回学術講演会講演論文集、63-66.
村田和則・相澤隆生・高橋亨・松岡俊文(2012)；被災堤防緊急対応のための３次元可視化ツール及び対策設計支援システムの開発（その２）－３次元可視化ツールの概要－、物理探査学会第127回学術講演会講演論文集、67-69.
Takahashi T., Aizawa T., Nishio H. and Matsuoka, T., 2013, 3D geophysical soil profiling system for supporting countermeasure design fordamaged river embankments, Proc. of Near-Surface Geophysics Asia Pacific Conference.

河川整備基金助成事業－統合物理探査による河川堤防の安全性評価技術の開発－（実施中）高橋亨・稲崎富士(2013)；物理探査データを用いた河川堤防の浸透性の推定、物理探査学会第128回学術講演会講演論文集、166-169.

20