日本情報地質学会シンポジウム 2016 講演論文集, 75-90, 2016
i-Const./CIM 時代の地形・地質三次元モデリングWeb 版ジオモデラーを使用したサーフェスモデルの推定
中田 文雄*・土屋 彰義*：安藤潤*：坂森 計則*：和田 里絵*

Presumption of 3D Geologic Model using WEB Modeling System
Fumio NAKADA,* Akiyoshi TSUCHIYA*, Jun ANDOH*,
Kazunori SAKAMORI* and Rie WADA** 一般社団法人 全国地質調査業協会連合会 Japan Geotechnical Consultants Association,
1-5-13, Uchikanda, Chiyoda-ku, Tokyo, 101-0047, Japan E-mail: nakadaf@gupi.jp

キーワード：三次元地質モデル，システム開発，ウェブ技術
Key words : Three-dimensional geological model, System development,75

3.CIM 対応三次元地盤モデルの構築を支援するウェブサイト

2★URL：https://geonews.zenchiren.or.jp/cim3d/index.html★開発理念： 誰でも公開用ウェブサイトにアクセスするだけで利用できる。

メインメニュー

3.1 電子納品・電子成果品の二次利用を支援するウェブサイト
3.3 三次元地盤モデル作成支援用ウェブサイト
3.4 三次元地盤モデルのデモサイト★研究報告書・プレゼンテーション等資料閲覧ウェブサイト

3.1 電子納品・電子成果品二次利用支援ウェブサイト

3

サブメニュー

3.1.1 位置座標の読取り・確認処理(パソコン)
3.1.2 位置座標の読取り・確認処理(スマートフォン・タブレット)
3.1.3 地質・土質成果，ボーリング交換用データ 位置座標確認処理
3.1.4 ボーリング交換用データ 位置座標の測地系変換処理
3.1.5 ボーリング交換用データ → メタデータ抽出処理
3.1.6 原位置試験データ・土質試験データ抽出処理
3.1.7 土質試験データシート表示処理76

3.1.1 位置座標の読取り・確認処理(パソコン)

4

主な機能★緯度・経度(度単位，度・分・秒単位，地図メッシュコード・六次)★中心点の標高(国土地理院)，住所(GoogleMaps)★緯度・経度(度，度・分・秒)と住所での地図移動とマーキング★産総研シームレス地質図など，8種類の地図タイルのオーバーレイ機能

3.1.2 位置座標の読取り・確認処理(スマートフォンなど)

単発測定

連続測定

標高と住所検索

主な機能★内蔵のGPSによる位置測定とデータの表示(単発と連続)★地図の中心を測定位置に移動させる機能★その他の機能は，パソコン版とほぼ同じ

地図タイル

77

5

2.5 サーフェスモデル (2)形状データの作成手順 (再掲)6

3.3 三次元地盤モデル作成支援用ウェブサイト

7

三次元地盤モデル＝サーフェスモデル

サブメニュー
3.3.2 (1)地層境界面の形状推定用入力データの作成(入力データ)
3.3.3 (2)地層の論理モデルの作成
3.3.4 (3)標高データの取得
3.3.5 (4)Terramod_BS による地層境界面形状の推定
3.3.6 (5)サーフェスモデルの推定
3.3.7 (6)パネルダイアグラムの推定
3.3.8 (7)三次元可視化処理(Web-GL)
3.3.9 (8)三次元表示用ボーリングモデル(柱状図)の作成
3.3.10 (9)パラメータデータの作成

78

三次元地盤モデルの作成(推定)を支援するウェブサイト

8

ウェブ版OCU GeoModeller
地質技術者が三次元地盤モデルを初めて作成する入門用処理システム★誰でも自由(無料)に利用できる。
★地層の統一凡例と，それに準拠したボーリング交換用データ(XML)を用意すれば，専用のウェブサイトにアクセスするだけで「サーフェスモデル」と「パネルダイアグラム」を推定できる。
注 Terramod_BSのみ，プログラムのダウンロードが必要★推定処理に必要な，全てのパラメータファイルも作成できる。
★サーフェスモデルの形状データのファイル形式は「3D-DXF」。
⇒ 目視確認と再利用が容易。
☆初心者用のため「地層数は10層まで」，地層の論理モデルは「堆積層と侵食・堆積層の2種類のみ」という制限を設けた。
注 論理モデルを自作すればこの制限は無い。

OCU GeoModeller：升本 眞二・野々垣 進・サラウット ニンサワット・岩村 里美・
櫻井健一・生賀 大之・ベンカテッシュ ラガワン・塩野 清治：Web-GISを用いた3次元地質モデル構築システム，情報地質，vol.20，no.2，pp.94-95，2009
注 OCU：Osaka City University（公立大学法人 大阪市立大学）

ウェブ版 OCU GeoModeller を構成するプログラム群(3)メッシュ座標の地形データは，国土地理院のウェブサイトから直接取得できるよう，専用のウエブサイトを開発した。
(4)サーフェスモデルの基となる 地層境界面のメッシュデータは，ボーリングデータからTerramod_BSで推定する。
★地層の論理モデルが違うと，推定結果に不具合が生じるため，手順(2)の判断が最も重要。

799

ウェブ版 OCU GeoModeller を構成するプログラム群
名称

主な機

能

(1)入力データの作成

・ボーリングデータ(XML)から，地層境界面の三次元形状を推定する・ツールに必要な入力データ群を作成する。
・ボーリングの経緯度から直交座標値を計算する機能がある。

(2)地層の論理モデルの作成

・地層境界のサーフェスモデルを作成する際に必要となる論理モデル(堆積，浸食・堆積)データを作成する。

(3)標高データの取得

・国土地理院のウェブサイトから，(4)で設定するメッシュ座標の標高を自動的に入手する。

(4)地層境界面の推定(Terramod_BS)

・(1)のデータに対して，Bスプライン関数を利用した推定法により，地層境界面の三次元形状を推定する。
・出力は，地層境界面のメッシュ(四角)標高データである。

(5)サーフェスモデルの推定・(4)で得られた地層境界面のメッシュ標高データと，(2)で得られた各境界面の論理モデルからサーフェスモデルを推定し，3D-DXF形式で出力する。

(6)パネルダイアグラムの推定・(4)で得られた地層境界面のメッシュ標高データと，(2)で得られた各境界面の論理モデルからサーフェスモデルを推定し，任意断面のパネルダイアグラム(イメージ)を出力する。

(7)三次元の可視化(WebGL)

・Web-GL技術を利用して(3)と(4)の出力データを三次元で可視化する。
・ボーリングモデルのイメージも合わせて可視化できる。
(現時点で公開準備中)(8)ボーリングモデル(イメージ)の作成

・柱状図を仮想的な三次元円柱として描画し，その外周に貼り付ける柱状図イメージ(png)を作成する。

(9)パラメータデーの作成

・(1)(5)(6)(8)で使用するパラメータファイルを作成する。
・事前に，地層名の統一凡例とカラーコードを決めておく必要がある。

【参考】ボーリングデータからサーフェスモデルへ
統 一 凡 例

ボーリング交換用データ

三次元地質凡例
年代

地層完
新世

完
新統

更新世

更新統

新第三紀 中新世

地質名記号

境界面

順序

色

埋土・盛土層

B

b5

S5

侵食面

D9D9D9

上部層

AU

b4

S4

堆積面

FFFF99

火山灰層

AV

b3

S3

堆積面

FF99FF

下部層

AD

b2

S2堆積面

3399FF

上部砂礫層

DG

b1

S1

初期地表面

339933

上部層

DU--0066FF

下部層

DD--009900

軟岩層

R--

996633-は出現せず地層パラメータ<地層岩体区分>タグ★緯度・経度 ⇒平面直角座標系への変換が必要
地層データ群

801011

【参考】ボーリングデータからサーフェスモデルへ

12

地層境界面の形状推定(ワイヤーフレームモデル)
地層の論理モデル

b5層

b4層

サーフェスモデル

サンプルデータを使用したサーフェスモデルの作成

13

サンプルデータ※n：1～5～8
・bn.txt：「 X，Y，Z，l」形式で保存した地層の上面標高・BEDSample000n.png：ボーリングイメージ
・BEDSample000n.XML：ボーリングデータ
・ChisoHanrei.txt：当該地域の統一凡例とカラーコード
・color_tbl.txt：地層の色コードの指定ファイル(パネルダイアグラム)・fs.txt：地層境界面の形状データファイル群の記載ファイル。 S1(基盤層) ⇒ S5(地表面)・logic.txt：地層の論理モデルファイル
・PanelDiagram-X5400.png：X位置5400mのパネルダイアグラム
・Sn.grd 地層境界面の形状を 記載したファイル群
・SurfaceModel.dxf：サーフェスモデルの例
・SurfaceModel.jpg：同上の確認結果※赤字：デフォルトファイル名(変更可)81サンプルデータを使用したサーフェスモデルの作成
ボーリング位置(架空)とモデル作成範囲

14

三次元地盤モデル

S5

S1

実際は地図を使用
統一凡例
年代

完
新世

更新世

地層

完
新統

更新統

地質名

記号

地層

境界面

堆積/ 侵食カラーコード

埋土・盛土層B

b5

S5堆積/ 侵食

217,217,217

上部層

AU

b4

S4

堆積

255,255,153

火山灰層

AV

b3

S3

堆積

255,153,255

下部層

AD

b2

S2

堆積

51,153,255

上部砂礫層

DG

b1

S1

堆積

51,153,51

3.3.2 (1)地層境界面の形状推定用入力データの作成(入力データ)15

前提条件(パラメータファイル)☆統一凡例の例：ChisoHanrei.txt
DG,AD,AV,AU,B,4
rgba(51,153,51,1):rgba(51,153,255,1):rgba(255,153,255,1):rgba(255,255,153,1):
rgba(217,217,217,1):
改行なし☆地層境界面の形状データファイル群の例：fs.txt☆地層の論理モデルファイルの例：logic.txt☆地層の色コード指定ファイル(パネルダイアグラム) の例：color_tbl.txt
S1.grd
S2.grd
S3.grd
S4.grd
S5.grd

5,6,S1,S2,S3,S4,S5
b1,-1,0,0,0,-1
b2,1,-1,0,0,-1
b3,1,1,-1,0,-1
b4,1,1,1,-1,-1
b5,1,1,1,1,-1alpha,0,0,0,0,1

0 0 0 01 51 153 51
2 51 153 255
3 255 153 255
4 255 255 153
5 217 217 217
6 255 255 255

fs.txtlogic.txt

color_tbl.txt82

3.3.2 (1)地層境界面の形状推定用入力データの作成(入力データ)

16☆ボーリングデータの留意点(必須登録データ)<経度_度>133</経度_度><経度_分>33</経度_分><経度_秒>30.4077</経度_秒><緯度_度>33</緯度_度><緯度_分>34</緯度_分><緯度_秒>2.8401</緯度_秒>

必須登録データ★位置情報(緯度・経度・標高)★地層岩体区分<孔口標高>0.4</孔口標高><地層岩体区分><地層岩体区分_上端深度>0.00</地層岩体区分_上端深度><地層岩体区分_下端深度>1.10</地層岩体区分_下端深度><地層岩体区分_地層岩体名>B</地層岩体区分_地層岩体名></地層岩体区分><地層岩体区分><地層岩体区分_上端深度>1.10</地層岩体区分_上端深度><地層岩体区分_下端深度>17.80</地層岩体区分_下端深度><地層岩体区分_地層岩体名>AU</地層岩体区分_地層岩体名></地層岩体区分

3.3.2 (1)地層境界面の形状推定用入力データの作成(入力データ)

ChisoHanrei.txt

BEDSample000n.XML

83

17

3.3.2 (1)地層境界面の形状推定用入力データの作成(入力データ)

18

OCU Geomodeller 用入力データ(Horizon 形式(X，Y，Z，l 形式))の出力

3.3.3 ( 2 )地層の論理モデルの作成

19

チェックで「侵食面」へ論理モデルは，以下のみ可能。
「侵食面」と「侵食面」☆地表面は，侵食面☆際下境界面は，初期地表面★断層などは非対応。

84

20

3.3.4 ( 3 )標高データの取得★国土地理院の公開API
により，直交座標値 ⇒緯度・経度値の自動計算★標高値の自動取得★出力データ形式は，(4)地層境界面の推定と同一である☆(4)地層境界面の推定と同じメッシュデータが必要☆南西端座標値(m 単位)☆X軸(南北)とY 軸(東西)の間隔(Dx，Dy，いずれもm 単位)☆区分数(Nx，Ny)

21

3.3.4 ( 3 )標高データの取得

留意点
表示データには，確認しやすいように・側線の開始座標・区切りは「,」が付加されている。
起点

データの取得方向

85

223.3.5 (4)地層境界面の推定(Terramod_BS)★本プログラムは，当該処理の管理用ウェブページからダウンロードし，Windows PC にコピーして使用する(インストール必要なし)。
★CIM 対応として，・位置データ形式：地層区分(層序)の位置情報(X，Y)と標高情報(Z)・不等式を含む位置データ形式(位置データ＋不等式(l))★原プログラムとは違い，X=縦(南北)，Y=横(東西)に改良した。
⇒平面直角座標系に合わせてある。

23

3.3.5 (4)地層境界面の推定(Terramod_BS)

平面直角座標系に留意

Y=横(東西)X=縦(南北)No.

YXZ8l6

1,5426.2854,62931.5093,-27.5,02,5521.8685,63006.8672,-27.7,03,5647.3262,62955.3699,-27,04,5521.7946,63135.7371,-36.5,05,5699.0645,63068.4111,-36.8,06,5774.7599,63171.5966,-27.5,07,5348.4870,63094.0504,-31.1,08,5396.2115,63250.6729,-27.3,09,9e9,9e9,9e9,97

186

24

3.3.5 (4)地層境界面の推定(Terramod_BS)
境界面の推定に使用するパラメータ
地表面の標高を取得したメッシュと同じ☆始点と終点の各座標値は，パネルダイアグラムの推定時に必要なので，別途記録しておくと良い。

コンター描画用のパラメータ設定の説明は省略

X=縦(南北)Y=横(東西)25

3.3.5 (4)地層境界面の推定(Terramod_BS)
推定結果の保存最新版は，小数点3桁に固定★ダンプリストを左90゜回転すると，等値線の向きに一致する。

87

3.3.6 (5)サーフェスモデルの推定

26

3.3.7 (6)パネルダイアグラムの推定

27★サーフェスモデルの作成処理に同じ・論理モデルファイル
・地層境界面ファイル名一覧ファイル
・地層境界面データファイル★主なパラメータ☆断面サイズ幅・高：ピクセルで指定するイメージの大きさ☆断面Z座標最低値・最高値：標高(m)で指定する断面高。
☆始点・終点座標：以下を参照
Terramod_BSのパラメータ★境界線の有無：有/無から選択★VRMLの出力：有/無から選択

88

28

3.3.8 (7)ボーリングモデル(イメージ)の作成★前提条件
・ボーリング交換用データの「<地層岩体区分>タグ」に，所定の地層情報を入力
・「地層データ(CSV)」が必要

ChisoHanrei.txt

CANVAS上のイメージ

BEDSample000n.XML＊.png

イメージの編集作業

29

3.3.9 (8 )パラメータデータの作成★作成可能データ・ChisoHanrei.txt：当該地域の統一
凡例とカラーコード
・color_tbl.txt：地層の色コードの指定
ファイル(パネルダイアグラム用)・fs.txt：地層境界面の形状データファイル群の記載ファイル。
S1(基盤層) ⇒ S5(地表面)・logic.txt：地層の論理モデルファイル

地表面

R，G，B
0～255で指定
平面直角座標系
・系番号( 1 ～19 )基盤層

テキストエディタに
コピ・ペ，保存でもよい

カラーコード確認用

89

地表面データは，地図(イメージ)にする場合でも，必ず作成のこと。

サンプルデータを使用したサーフェスモデルの作成(実践)30(9)パラメータデータの作成(1)地層境界面の形状推定用入力データの作成(入力データ)(2)地層の論理モデルの作成(3)標高データの取得(4)Terramod_BS による地層境界面形状の推定(5)サーフェスモデルの推定(6)パネルダイアグラムの推定
デモ［(7)三次元可視化処理(Web-GL)］(8)三次元表示用ボーリングモデル(柱状図)の作成

31

おわりに★簡単な三次元モデルの形状データを作成できる。
⇒ より複雑な地層の論理モデルが必要な場合は，(有料)の専門サーフェスモデラーを利用されたい。
★下流工程を勘案して，最も利用されているファイル形式(の一つ)での形状データのファイル形式を採用。
★三次元地質モデルの特徴
・地層境界面の形状を忠実に復元するためには，極めて多くのボーリングデータを必要とする。
⇒ 設計に必要な精度を持つモデルの推定は困難。
・三次元地質モデル作成のために，ボーリング箇所などを増やしてくれる保証は無い。
⇒ VRの中に三次元モデルを表示する程度の利用しか望めないだろう。
◆地元説明(建設同意)，議会説明(予算獲得)

90