特許 6932614 亀裂連続性判定システム、判定方法及び判定プログラム、並びに、亀裂検出装置、検出方法及び検出プログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6932614

(24)【登録日】2021年8月20日

(45)【発行日】2021年9月8日

(54)【発明の名称】亀裂連続性判定システム、判定方法及び判定プログラム、並びに、亀裂検出装置、検出方法及び検出プログラム

(51)【国際特許分類】

   G01V 1/00 20060101AFI20210826BHJP

   G01V 8/02 20060101ALI20210826BHJP

   E02D 1/02 20060101ALI20210826BHJP

【ＦＩ】

   G01V1/00 C

   G01V8/02

   E02D1/02

【請求項の数】12

【全頁数】23

(21)【出願番号】特願2017-210823(P2017-210823)

(22)【出願日】2017年10月31日

(65)【公開番号】特開2019-82443(P2019-82443A)

(43)【公開日】2019年5月30日

【審査請求日】2020年4月28日

(73)【特許権者】

【識別番号】000001373

【氏名又は名称】鹿島建設株式会社

(74)【代理人】

【識別番号】110002468

【氏名又は名称】特許業務法人後藤特許事務所

(72)【発明者】

【氏名】升元 一彦

(72)【発明者】

【氏名】野中 隼人

(72)【発明者】

【氏名】中嶌 誠門

【審査官】 川瀬 正巳

(56)【参考文献】

【文献】 特開昭６４−０３９４９２（ＪＰ，Ａ）

【文献】 特開平０６−１４７８４６（ＪＰ，Ａ）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｖ １／００

Ｇ０１Ｖ ８／０２

Ｅ０２Ｄ １／０２

(57)【特許請求の範囲】

【請求項1】

地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して前記展開画像から亀裂に基づいた正弦波形を抽出し、抽出された前記正弦波形から、前記第１及び第２ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得する亀裂検出装置と、
前記亀裂検出装置により検出された前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定するための指標を算出する指標算出部と、
前記指標算出部により算出された前記指標に基づいて、前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定する連続性判定部と、を備える
亀裂連続性判定システム。

【請求項2】

前記指標は、前記第１ボアホールの前記亀裂の走向及び傾斜に基づいて前記第１ボアホールの前記亀裂を前記第２ボアホールに延長したときに前記第２ボアホールに形成される仮想亀裂と前記第２ボアホールの前記亀裂との深度差を含む
請求項１に記載の亀裂連続性判定システム。

【請求項3】

前記指標は、前記第１ボアホールの前記亀裂の走向及び傾斜から求められる法線ベクトルと、前記第２ボアホールの前記亀裂の走向及び傾斜から求められる法線ベクトルと、の交差角度を含む
請求項１又は２に記載の亀裂連続性判定システム。

【請求項4】

前記亀裂検出装置は、前記展開画像を解析して、前記亀裂の幅を取得する亀裂幅取得部を更に備え、
前記指標は、前記第１ボアホールにおける前記亀裂の幅と前記第２ボアホールにおける前記亀裂の幅との差を含む
請求項１から３のいずれか１項に記載の亀裂連続性判定システム。

【請求項5】

前記第１ボアホールにおける前記亀裂の幅と、前記第２ボアホールにおける前記亀裂の幅とは、前記傾斜を用いて補正された幅である
請求項４に記載の亀裂連続性判定システム。

【請求項6】

前記亀裂検出装置は、前記展開画像を解析して前記亀裂の色情報を取得する亀裂色情報取得部を更に備え、
前記指標は、前記第１ボアホールにおける前記亀裂の色情報と前記第２ボアホールにおける前記亀裂の色情報との差を含む
請求項１から５のいずれか１項に記載の亀裂連続性判定システム。

【請求項7】

前記指標は、前記第１ボアホールにおける前記亀裂の深度でのルジオン値と、前記第２ボアホールにおける前記亀裂の深度でのルジオン値と、の差を含む
請求項１から６のいずれか１項に記載の亀裂連続性判定システム。

【請求項8】

地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して、前記展開画像から亀裂に基づいた正弦波形を抽出し、
前記正弦波形から、前記第１及び第２ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得し、
検出された前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定するための指標を算出し、
算出された前記指標に基づいて、前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定する
亀裂連続性判定方法。

【請求項9】

地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して、前記展開画像から亀裂に基づいた正弦波形を抽出し、
前記正弦波形から、前記第１及び第２ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得し、
検出された前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定するための指標を算出し、
算出された前記指標に基づいて、前記第１及び第２ボアホールにおける前記亀裂どうしの連続性を判定する
処理をコンピュータに実行させるための亀裂連続性判定プログラム。

【請求項10】

地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、前記展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出部と、
前記正弦波形抽出部により抽出された前記亀裂に基づいた正弦波形から、前記ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得する走向傾斜取得部と、を備え、
前記正弦波形抽出部は、複数の仮の正弦波形を前記展開画像に重ね、前記複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を前記複数の仮の正弦波形毎に算出し、算出された前記代表値に基づいて、前記複数の仮の正弦波形のうちの１つを、前記亀裂に基づいた正弦波形として抽出する、
亀裂検出装置。

【請求項11】

地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、前記展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出ステップと、
前記亀裂に基づいた正弦波形から、前記ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得する走向傾斜取得ステップと、を備え、
前記正弦波形抽出ステップでは、複数の仮の正弦波形を前記展開画像に重ね、前記複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を前記複数の仮の正弦波形毎に算出し、算出された前記代表値に基づいて、前記複数の仮の正弦波形のうちの１つを、前記亀裂に基づいた正弦波形として抽出する、
亀裂検出方法。

【請求項12】

地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、前記展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出処理と、
前記亀裂に基づいた正弦波形から、前記ボアホールにおける前記亀裂を検出すると共に前記亀裂の走向及び傾斜を取得する走向傾斜取得処理と、
をコンピュータに実行させるための亀裂検出プログラムであって、
前記正弦波形抽出処理では、複数の仮の正弦波形を前記展開画像に重ね、前記複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を前記複数の仮の正弦波形毎に算出し、算出された前記代表値に基づいて、前記複数の仮の正弦波形のうちの１つを、前記亀裂に基づいた正弦波形として抽出する
処理を前記コンピュータに実行させるための亀裂検出プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、亀裂検出装置、亀裂連続性判定システム、亀裂検出方法及び亀裂検出プログラムに関する。

【背景技術】

【0002】

ダム工事では、基礎地盤の遮水性を高めるために、基礎地盤にボアホールを形成しグラウトを注入して基礎地盤中の亀裂を塞ぐグラウチングが行われる。特許文献１には、ボアホール中の亀裂の走向傾斜等のデータから亀裂の分布密度を算出し、算出された分布密度に基づいてグラウチングの仕様を決定することが提案されている。

【先行技術文献】

【特許文献】

【0003】

【特許文献1】特開２０１４−２１１３７４号公報

【発明の概要】

【発明が解決しようとする課題】

【0004】

特許文献１における亀裂の走向傾斜等のデータは、ボアホール内を撮影して得られた画像を測定者が見てその画像上の亀裂から手作業により計測しているのが現状である。そのため、亀裂のデータの取得に時間がかかると共に、測定者によって測定精度が異なるといった問題がある。

【0005】

本発明は、迅速かつ一定の精度で亀裂を検出することを目的とする。

【課題を解決するための手段】

【0006】

本発明のある態様によれば、亀裂検出システムは、地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して展開画像から亀裂に基づいた正弦波形を抽出し、抽出された正弦波形から、第１及び第２ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得する亀裂検出装置と、亀裂検出装置により検出された第１及び第２ボアホールにおける亀裂どうしの連続性を判定するための指標を算出する指標算出部と、指標算出部により算出された指標に基づいて、第１及び第２ボアホールにおける亀裂どうしの連続性を判定する連続性判定部と、を備える。
また、本発明の別の態様によれば、亀裂連続性判定方法は、地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して、展開画像から亀裂に基づいた正弦波形を抽出し、正弦波形から、第１及び第２ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得し、検出された第１及び第２ボアホールにおける亀裂どうしの連続性を判定するための指標を算出し、算出された指標に基づいて、第１及び第２ボアホールにおける亀裂どうしの連続性を判定する。
また、本発明の別の態様によれば、亀裂連続性判定プログラムは、地盤に形成された第１及び第２ボアホールの内壁を撮影して得られた展開画像を解析して、展開画像から亀裂に基づいた正弦波形を抽出し、正弦波形から、第１及び第２ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得し、検出された第１及び第２ボアホールにおける亀裂どうしの連続性を判定するための指標を算出し、算出された指標に基づいて、第１及び第２ボアホールにおける亀裂どうしの連続性を判定する、処理をコンピュータに実行させる。
また、本発明の別の態様によれば、亀裂検出装置は、地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出部と、正弦波形抽出部により抽出された亀裂に基づいた正弦波形から、ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得する走向傾斜取得部と、を備え、正弦波形抽出部は、複数の仮の正弦波形を展開画像に重ね、複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を複数の仮の正弦波形毎に算出し、算出された代表値に基づいて、複数の仮の正弦波形のうちの１つを、亀裂に基づいた正弦波形として抽出する。

【0007】

また、本発明の別の態様によれば、亀裂検出方法は、地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出ステップと、亀裂に基づいた正弦波形から、ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得する走向傾斜取得ステップと、を備え、正弦波形抽出ステップでは、複数の仮の正弦波形を展開画像に重ね、複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を複数の仮の正弦波形毎に算出し、算出された代表値に基づいて、複数の仮の正弦波形のうちの１つを、亀裂に基づいた正弦波形として抽出する。

【0008】

また、本発明の別の態様によれば、亀裂検出プログラムは、地盤に形成されたボアホールの内壁を撮影して得られた展開画像を解析して、展開画像から亀裂に基づいた正弦波形を抽出する正弦波形抽出処理と、亀裂に基づいた正弦波形から、ボアホールにおける亀裂を検出すると共に亀裂の走向及び傾斜を取得する走向傾斜取得処理と、をコンピュータに実行させるための亀裂検出プログラムであって、正弦波形抽出処理では、複数の仮の正弦波形を展開画像に重ね、複数の仮の正弦波形の各々が通過する複数の画素における輝度の代表値を複数の仮の正弦波形毎に算出し、算出された代表値に基づいて、複数の仮の正弦波形のうちの１つを、前記亀裂に基づいた正弦波形として抽出する。

【発明の効果】

【0009】

本発明によれば、迅速かつ一定の精度で亀裂を検出することができる。

【図面の簡単な説明】

【0010】

【図1】本発明の実施形態に係る亀裂検出装置の構成を示す概略図である。

【図2】（ａ）は、２つの亀裂を貫通するボアホールを示す斜視図であり、（ｂ）は、（ａ）に示すボアホールの壁面を撮影して展開することにより得られる展開画像である。

【図3】亀裂の走向及び傾斜を説明するための図である。

【図4】本発明の実施形態に係る亀裂検出装置の構成を示すブロック図である。

【図5】展開画像から亀裂に近似する正弦波形を抽出する処理を示すフローチャートである。

【図6】（ａ）は、ボアホールの内壁の展開画像であり、（ｂ）は、代表Ｇｒａｙ輝度値を示すグラフである。

【図7】亀裂の走向及び傾斜を取得する処理を示すフローチャートである。

【図8】亀裂の幅を取得する処理を示すフローチャートである。

【図9】（ａ）は、ボアホールの内壁の展開画像の拡大図であり、（ｂ）は、Ｇｒａｙ輝度を示すグラフである。

【図10】亀裂の走向に平行にボアホールを透視した模擬拡大図である。

【図11】亀裂の色情報を取得する処理を示すフローチャートである。

【図12】本発明の実施形態に係る亀裂連続性判定システムの構成を示すブロック図である。

【図13】亀裂どうしの連続性を判定するための指標を説明する図である。

【図14】亀裂どうしの見かけの深度差を算出する処理を示すフローチャートである。

【図15】亀裂どうしの法線ベクトルの交差角度を算出する処理を示すフローチャートである。

【図16】亀裂どうしの幅差を算出する処理を示すフローチャートである。

【図17】亀裂どうしの色情報差を算出する処理を示すフローチャートである。

【図18】亀裂どうしのルジオン値差を算出する処理を示すフローチャートである。

【図19】亀裂どうしの連続性を判定する処理を示すフローチャートである。

【図20】亀裂どうしの連続性の判定結果の表示例を示す。

【図21】亀裂どうしの面間距離を説明する図である。

【発明を実施するための形態】

【0011】

以下、図面を参照して、本発明の実施形態に係る亀裂検出装置１００、亀裂連続性判定システム１０００、亀裂検出方法及び亀裂検出プログラムについて説明する。

【0012】

図１は、ダム等における基礎地盤１を示す断面図である。ダム工事では、基礎地盤１の遮水性を高めるために、基礎地盤１にボアホール２を複数形成しグラウトを注入して基礎地盤１中の亀裂３を塞ぐグラウチングが行われる。グラウチングでは、まず、１次ボアホール２を基礎地盤１に所定間隔で複数形成し、必要に応じて、１次ボアホール２どうしの間に２次ボアホール、３次ボアホールといった高次ボアホールを追加で形成する。高次ボアホールの位置や数は、基礎地盤１中における亀裂３の大きさや数に応じて決定されることが望ましい。

【0013】

本実施形態に係る亀裂連続性判定システム１０００（図１２参照）は、予め形成されたボアホール２の内壁における亀裂３の情報に基づいて亀裂３を推定するものである。亀裂連続性判定システム１０００により亀裂３を推定することにより、追加の高次ボアホール２の数や位置を容易に決定することが可能になる。また、亀裂検出装置１００は、基礎地盤１に形成されたボアホール２の内壁における亀裂３を検出すると共に亀裂３の走向及び傾斜等の亀裂３の情報を取得するものである。亀裂検出装置１００によりボアホール２における亀裂３の情報を取得することにより、亀裂連続性判定システム１０００における亀裂３の推定精度を向上させることができる。

【0014】

まず、亀裂検出装置１００及び亀裂検出方法について、図１〜図１１を参照して説明する。亀裂検出装置１００は、画像作成手段としてのボアホールカメラ１０により得られた展開画像を解析する解析手段としてのコンピュータ２０を備える。ボアホールカメラ１０は、ボアホール２の内壁を撮影して展開画像を作成する。コンピュータ２０は、演算処理を行うＣＰＵ（ｃｅｎｔｒａｌ ｐｒｏｃｅｓｓｉｎｇ ｕｎｉｔ）と、ＣＰＵにより実行されるプログラムが予め記憶されたＲＯＭ（ｒｅａｄ−ｏｎｌＹ ｍｅｍｏｒＹ）と、ボアホールカメラ１０から入力されるデータやＣＰＵの演算結果等を記憶するＲＡＭ（ｒａｎｄｏｍ ａｃｃｅｓｓ ｍｅｍｏｒＹ）と、を含む。

【0015】

ボアホールカメラ１０は、ボアホール２の内壁を周方向に３６０度カラー撮影可能に形成される。ボアホールカメラ１０は、ウインチ１１に巻きつられたワイヤ１２の一端に取り付けられ、ウインチ１１の駆動によりボアホール２内を上下に移動する。ボアホールカメラ１０を上下に移動させながらボアホール２の内壁を撮影することにより、ボアホール２の全深度に渡るカラーの展開画像（図２（ｂ）参照）が得られる。

【0016】

図２（ａ）は、２つの亀裂３ａ，３ｂを貫通するボアホール２を示す斜視図であり、図２（ｂ）は、図２（ａ）に示すボアホール２の壁面を撮影して展開することにより得られる展開画像である。図２（ａ）に示すように、亀裂３ａ，３ｂは、ボアホール２の内壁では走向及び傾斜に応じた楕円に似た形状を描く。

【0017】

図３は、亀裂３の「走向」及び「傾斜」を説明するための図である。「走向」とは、亀裂３と仮想水平面ＨＳとの交差によりできる交線ＩＬが延びる方向である。ここでは、「走向」を、交線ＩＬと南北方向に延びる軸との間の角度θ［ｄｅｇ］として表し、北（Ｎ）を基準に西（Ｗ）回りを正方向とする。また、「傾斜」とは、交線ＩＬと直交する垂線ＰＬが亀裂３に沿って延びる方向である。ここでは、「傾斜」を、垂線ＰＬと仮想水平面ＨＳとの間の角度φ［ｄｅｇ］として表す。

【0018】

図２（ｂ）に示すように、展開画像では、亀裂３ａ，３ｂは、正弦波に似た形状を描く。具体的には、亀裂３ａの走向θａに応じて、展開画像での亀裂３ａの上側ピーク点ＰＵａと下側ピーク点ＰＬａの方位が定まり、亀裂３ａの傾斜φａとボアホール２の半径Ｒに応じて、展開画像での亀裂３ａの振幅Ａａが定まる。亀裂３ｂについても亀裂３ａと同様に、展開画像では亀裂３ｂの走向θｂ及び傾斜φｂに応じた正弦波に近い形状を描く。亀裂検出装置１００（図１参照）は、展開画像での亀裂３ａ，３ｂの形状に基づいて、亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂを取得する。

【0019】

図４に示すように、コンピュータ２０は、画像を取り込む画像取込部２１と、展開画像から亀裂３ａ，３ｂに近似する正弦波形を抽出する正弦波形抽出部２２と、亀裂３ａ，３ｂの走向θａ、θｂ及び傾斜φａ，φｂを取得する走向傾斜取得部２３と、を備える。コンピュータ２０は、信号線１３を介してボアホールカメラ１０に接続され、ボアホールカメラ１０により作成されたカラーの展開画像は、画像取込部２１によりコンピュータ２０に取り込まれる。以下、正弦波形抽出部２２及び走向傾斜取得部２３によって行われる処理について、具体的に説明する。

【0020】

まず、正弦波形抽出部２２によって行われる正弦波形抽出処理について、図２〜図６を参照して説明する。

【0021】

図５に示すように、まず、ステップＳ５０１にて、取り込まれたカラーの展開画像をグレースケール画像に変換する。具体的には、カラーの展開画像における各画素に対して赤（Ｒ）、緑（Ｇ）、青（Ｂ）の輝度の重み付け平均を算出する。この重み付け平均が、グレースケール画像における各画素のＧｒａｙ輝度となる。カラーの展開画像をグレースケール画像に変換することにより、各画素の輝度がＧｒａｙ輝度のみになり、後述する亀裂判定用閾値と展開画像の輝度との比較が容易になる。

【0022】

次に、ステップＳ５０２〜ステップＳ５０６にて、グレースケール画像から、亀裂３ａ，３ｂ（図２（ｂ）参照）に近似する正弦波形を求める。具体的には、次の式（１）の関数Ｆ（Ｙ）により表される正弦波形をグレースケール画像に重ね、この正弦波形が通過する複数の画素のＧｒａｙ輝度の平均値に基づいて、関数Ｆ（Ｙ）により表される正弦波形が亀裂３ａ，３ｂに近似しているか否かを判定する。

【0023】

Ｆ（Ｙ）＝Ｚ＋Ｒ×ｔａｎφ×ｓｉｎ（Ｙ−θ） ・・・・式（１）
ここで、Ｙ：北（Ｎ）を基準として任意の方位を表す角度［ｄｅｇ］
Ｚ：深度［ｍ］
Ｒ：ボアホール２の半径［ｍ］
θ：走向［ｄｅｇ］
φ：傾斜［ｄｅｇ］

【0024】

以下、ステップＳ５０２〜Ｓ５０６にて行われる亀裂３ａ，３ｂに近似する正弦波形を求める処理について、より詳細に説明する。

【0025】

図２には、関数Ｆ（Ｙ）における深度ＺをＺ１とし、正弦波形の振幅を表すＲ×ｔａｎφをＡ１とし、走向θをθ１としたときの正弦波形ＳＷ１を鎖線で示している。通常、亀裂３ａ，３ｂのＧｒａｙ輝度は、亀裂３ａ，３ｂ以外の部分のＧｒａｙ輝度よりも小さい。そのため、正弦波形ＳＷ１が通過する複数の画素のＧｒａｙ輝度の平均値が小さいほど、正弦波形ＳＷ１が亀裂３ａ，３ｂに近似していることになる。つまり、正弦波形ＳＷ１が通過する複数の画素のＧｒａｙ輝度の平均値に基づいて、正弦波形ＳＷ１が亀裂３ａ，３ｂに近似しているか否かを判定することが可能である。

【0026】

展開画像における亀裂３ａ，３ｂに近似する正弦波形ＳＷ１を得るために、まず、ステップＳ５０２にて、深度Ｚを一定値とし、走向θを−１８０〜＋１８０［ｄｅｇ］の範囲で変化させる（例えば１［ｄｅｇ］毎）と共に傾斜φを０〜９０［ｄｅｇ］の範囲で変化させる（例えば１［ｄｅｇ］毎）。このとき、走向θ及び傾斜φの各条件下での正弦波形をグレースケール画像に重ね、正弦波形が通過する複数の画素のＧｒａｙ輝度の平均値をそれぞれ算出する。

【0027】

次に、ステップＳ５０３にて、算出された複数のＧｒａｙ輝度平均値の最小値と、Ｇｒａｙ輝度平均値の最小値が算出される走向θ及び傾斜φと、を取得する。Ｇｒａｙ輝度平均値の最小値が算出される正弦波形は、その深度Ｚに亀裂がある場合にその亀裂に最も近似する正弦波形である。以下において、その深度ＺにおけるＧｒａｙ輝度平均値の最小値を、その深度Ｚにおける「代表Ｇｒａｙ輝度値」と称する。

【0028】

次に、ステップ５０４にて、予め定められた深度範囲に渡って代表Ｇｒａｙ輝度値を取得したか否かを判断する。予め定められた深度範囲は、例えば、ボアホール２の全深度である。予め定められた深度範囲は、操作者によって変更可能であってもよい。

【0029】

ステップＳ５０４にて、予め定められた深度範囲に渡って代表Ｇｒａｙ輝度値を取得していないと判断した場合には、ステップＳ５０５にて深度Ｚを例えば０．０１［ｍ］変化させ、ステップＳ５０２，Ｓ５０３を再び実行する。これにより、予め定められた深度範囲に渡って代表Ｇｒａｙ輝度値が取得される。ステップＳ５０４にて、予め定められた深度範囲に渡って代表Ｇｒａｙ輝度値を取得したと判断した場合には、ステップＳ５０６を実行する。

【0030】

図６（ｂ）に、ステップＳ５０２，Ｓ５０３を予め定められた深度範囲に渡って実行することによって得られる代表Ｇｒａｙ輝度値を示すグラフを示す。図６（ａ）及び（ｂ）に示すように、通常、代表Ｇｒａｙ輝度値は、亀裂３ａ，３ｂが存在する深度区間において、他の深度区間よりも小さくなる。そこで、亀裂検出装置１００では、代表Ｇｒａｙ輝度値に基づいて、亀裂３ａ，３ｂの有無を判定する。

【0031】

図５及び図６に示すように、ステップＳ５０６にて、各深度Ｚにおける代表Ｇｒａｙ輝度値と予め定められた亀裂判定用閾値とを比較し、代表Ｇｒａｙ輝度値が亀裂判定用閾値以下の場合に亀裂３ａ，３ｂが有ると判定する。また、亀裂判定用閾値以下の代表Ｇｒａｙ輝度値が算出される正弦波形ＳＷａ，ＳＷｂを、亀裂３ａ，３ｂに近似する正弦波形であると判定する。亀裂判定用閾値は、操作者によって変更可能であってもよい。

【0032】

以上により、展開画像から亀裂３ａ，３ｂに基づいて正弦波形が抽出され、亀裂抽出処理が終了する。

【0033】

次に、走向傾斜取得部２３により行われる走向傾斜取得処理について、図６及び図７を参照して説明する。

【0034】

図６（ａ），（ｂ）に示すように、亀裂３ａ，３ｂには、ボアホール２の深さ方向に幅がある。そこで、まず、ステップＳ７０１（図７参照）にて、亀裂３ａ，３ｂの代表の深度Ｚａ，Ｚｂを求める。具体的には、亀裂３ａ，３ｂがあると判定された各深度区間ＺＲａ，ＺＲｂにおいて、代表Ｇｒａｙ輝度値が最小となる深度Ｚを求め、その深度Ｚを、亀裂３ａ，３ｂの代表の深度Ｚａ，Ｚｂとする。以上により、亀裂３ａ，３ｂの代表の深度Ｚａ，Ｚｂが求められる。

【0035】

次に、ステップＳ７０２（図７参照）にて、代表の深度Ｚａ，Ｚｂにおいて亀裂３ａ，３ｂに近似すると判定された正弦波形から、走向θ及び傾斜φを求める。具体的には、ステップＳ５０３（図５）にて各深度Ｚに対して代表Ｇｒａｙ輝度値を算出する走向θ及び傾斜φが取得されているので、代表の深度Ｚａ，Ｚｂで取得された走向θ及び傾斜φを読み出す。

【0036】

前述のように、関数Ｆ（Ｙ）における走向θ及び傾斜φは、ボアホール２の内壁に描かれる楕円形状の走向及び傾斜に相当する。したがって、代表の深度Ｚａで取得された走向θ及び傾斜φを亀裂３ａの走向θａ及び傾斜φａ（図２参照）と見なすことにより、亀裂３ａの走向θａ及び傾斜φａを取得することができる。同様に、代表の深度Ｚｂで取得された走向θ及び傾斜φを亀裂３ｂの走向θｂ及び傾斜φｂ（図２参照）と見なすことにより、亀裂３ｂの走向θｂ及び傾斜φｂを取得することができる。

【0037】

以上により、正弦波形抽出部２２により抽出された正弦波形から、ボアホール２の内壁での亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂ（図２参照）が取得され、走向傾斜取得処理が終了する。

【0038】

なお、本実施形態では、ステップＳ５０２にて、深度Ｚを一定値とし、走向θを変化させると共に傾斜φを変化させてＧｒａｙ輝度の平均値をそれぞれ算出しているが、走向θ及び傾斜φを一定値とし、深度Ｚを変化させてＧｒａｙ輝度の平均値を算出してもよい。この場合、ステップＳ５０３では、Ｇｒａｙ輝度値の平均値が深度方向に所定の深度以上連続して閾値を下回る区間を亀裂候補区間として抽出し、亀裂候補区間におけるＧｒａｙ輝度平均値の最小値（代表Ｇｒａｙ輝度）と、代表Ｇｒａｙ輝度が算出される深度Ｚと、を取得する。そして、予め定められた走向範囲及び傾斜範囲に渡って亀裂候補区間及び代表Ｇｒａｙ輝度値を取得する。解析精度分異なる走向θ、傾斜φの亀裂候補区間を比較して、深度分布が所定の深度以上重なっておりかつ走向θ及び傾斜φの一方が同じで他方が解析精度分違う場合には、同一亀裂区間と判定する。同一亀裂区間の判定に用いられなかった亀裂候補区間は、全て個別の亀裂区間として判定する。このようにして、走向傾斜取得処理を行ってもよい。

【0039】

亀裂検出装置１００では、コンピュータ２０の処理によって、ボアホールカメラ１０により作成されたカラーの展開画像から、亀裂３ａ，３ｂに近似する正弦波形が抽出され、抽出された正弦波形から、亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂが取得される。したがって、迅速かつ一定の精度で亀裂３ａ，３ｂを検出することができると共に亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂを取得することができる。

【0040】

取得された亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂは、後述する亀裂連続性判定処理にて使用される。

【0041】

後述の亀裂連続性判定処理では、亀裂３ａ，３ｂの走向θａ，θｂ及び傾斜φａ，φｂに加え、亀裂３ａ，３ｂの幅及び色情報が使用される。亀裂３ａ，３ｂの幅及び色情報は、それぞれ、亀裂検出装置１００の亀裂幅取得部２４及び亀裂色情報取得部２５（図４参照）によって取得される。以下では、亀裂３ａ，３ｂの幅及び色情報を取得する処理についてそれぞれ説明する。

【0042】

まず、亀裂幅取得部２４により行われる亀裂３ａの幅を取得する処理について、図２、図８から図１０を参照して説明する。

【0043】

図２に示すように、展開画像における亀裂３ａの幅は、方位軸に関して一定ではない。そこで、ここでは、上側ピーク点ＰＵａにおける亀裂３ａの幅ＷＵａと、下側ピーク点ＰＬａにおける亀裂３ａの幅ＷＵｂと、に基づいて算出される幅を、亀裂３ａの代表幅として取得する。

【0044】

まず、ステップＳ８０１（図８参照）にて、亀裂３ａの代表の深度Ｚａにおいて亀裂３ａに近似する正弦波形を取得する。具体的には、ステップＳ５０３（図５参照）にて各深度Ｚに対して代表Ｇｒａｙ輝度値を算出する走向θ及び傾斜φが取得されているので、代表の深度Ｚａで取得された走向θ及び傾斜φを読み出し、関数Ｆ（Ｙ）に代入することにより正弦波形を取得する。

【0045】

次に、ステップＳ８０２（図８参照）にて、上側ピーク点ＰＵａにおける亀裂３ａの幅ＷＵａを、Ｇｒａｙ輝度に基づいて求める。亀裂３ａの幅ＷＵａを求める処理について、図９を参照してより具体的に説明する。

【0046】

図９（ａ）は、グレースケール画像の拡大図であり、亀裂３ａの上側ピーク点ＰＵａの近傍を示す。図９（ａ）では、代表の深度Ｚａにおいて亀裂３ａに近似する正弦波形ＳＷａが亀裂３ａに重ねられ、正弦波形ＳＷａの上側ピーク点Ｐを通りボアホール２の深さ方向に延びる直線Ｌが鎖線で示されている。また、図９（ｂ）は、直線Ｌに沿う深度毎のＧｒａｙ輝度を示すグラフである。

【0047】

図９（ａ）及び（ｂ）に示すように、通常、直線Ｌに沿うＧｒａｙ輝度は、亀裂３ａが存在する深度区間において、他の深度区間よりも小さくなる。そこで、ステップＳ８０２（図８参照）では、直線Ｌに沿うＧｒａｙ輝度と予め定められた亀裂幅取得用閾値とを比較し、Ｇｒａｙ輝度が亀裂幅取得用閾値以下となる深度区間の長さを亀裂３ａの幅ＷＵａとして取得する。亀裂幅取得用閾値は、操作者によって変更可能であってもよい。

【0048】

以上の処理により、亀裂３ａの幅ＷＵａが求められる。

【0049】

次に、ステップＳ８０３（図８参照）にて、図２に示す下側ピーク点ＰＬａにおける亀裂３ａの幅ＷＬａを、Ｇｒａｙ輝度に基づいて求める。ステップＳ８０３における処理は、ステップＳ８０２における処理と略同じであるため、ここではその説明を省略する。

【0050】

次に、ステップＳ８０４にて、亀裂３ａの幅ＷＵａと幅ＷＬａとの平均値を算出する。

【0051】

ここで、展開画像から求められる亀裂３ａの幅ＷＵａ及び幅ＷＬａと、亀裂３ａの傾斜φａと、の関係について、図１０を参照して説明する。図１０は、ボアホール２における亀裂３ａを、その走向に平行に透視した模擬拡大図である。図１０には、亀裂３ａの傾斜φａを小さくした仮想の亀裂３ａ’を示している。

【0052】

ステップＳ８０３，Ｓ８０４（図８参照）にて求められる幅ＷＵａ，ＷＬａは、ボアホール２の内壁での亀裂３ａの幅、つまりボアホール２の深度方向における亀裂３ａの幅に相当する。そのため、亀裂３ａの幅ＷＵａと幅ＷＬａとの平均値は、ボアホール２の中心での深度方向における亀裂３ａの幅ＷＡａに相当する。亀裂３ａの幅ＷＡａは、亀裂３ａと仮想の亀裂３ａ’とを比較すると明らかなように、傾斜φが大きいほど大きくなる。そのため、後述する亀裂連続性判定処理において、亀裂３ａの幅ＷＡａを用いて亀裂の連続性を判定すると、亀裂３ａの傾斜φａの影響により判定精度が低下するおそれがある。

【0053】

そこで、亀裂検出装置１００では、Ｓ８０５（図８参照）にて、亀裂３ａの傾斜φａに基づいて、亀裂３ａの幅ＷＵａと幅ＷＬａとの平均値（幅ＷＡａ）を補正する。具体的には、幅ＷＡａに傾斜φａの余弦値（ｃｏｓφａ）を乗じて、幅ＷＡａを補正する。これにより、ボアホール２の中心での亀裂３ａの法線方向における幅ＷＮａが得られる。以下において、幅ＷＮａを「代表幅ＷＮａ」とも称する。

【0054】

以上により、展開画像から亀裂３ａの代表幅ＷＮａが取得され、亀裂幅取得処理が終了する。亀裂３ｂの代表幅を取得する処理は、亀裂３ａの代表幅ＷＮａを取得する処理と略同じであるため、ここではその説明を省略する。

【0055】

なお、本実施形態では、亀裂３ａの幅ＷＵａと幅ＷＬａとの平均値（幅ＷＡａ）を算出し更に亀裂３ａの傾斜φａの余弦値（ｃｏｓφａ）を乗じて幅ＷＡａを補正して代表幅ＷＮａを取得しているが、亀裂区間の区間長（見かけ亀裂幅）を算出し傾斜φａで補正して代表幅ＷＮａを取得してもよい。この場合、代表幅ＷＮａは、実際の割れ目幅ではなく、平均的な亀裂幅として算出される。

【0056】

亀裂検出装置１００では、コンピュータ２０の処理によって、展開画像から亀裂３ａの代表幅ＷＮａを取得する。したがって、亀裂３ａの代表幅ＷＮａを迅速かつ一定の精度で取得することができる。

【0057】

また、亀裂検出装置１００では、ボアホール２の深さ方向における亀裂３ａの幅ＷＡａを、傾斜φの余弦値を乗じて補正することにより、亀裂３ａの法線方向における代表幅ＷＮａを取得する。したがって、亀裂３ａの傾斜φａの影響を受けない幅を取得することができる。

【0058】

次に、亀裂色情報取得部２５により行われる亀裂色情報取得処理について、図１１を参照して説明する。図示を省略するが、亀裂３ａの色は、図２に示す亀裂３ａの幅と同様に、方位軸に関して一定ではない。そこで、ここでは、上側ピーク点ＰＵａにおける亀裂３ａの色情報と、下側ピーク点ＰＬａにおける亀裂３ａの色情報と、の平均値を、亀裂３ａの代表色情報として取得する。

【0059】

まず、図１１に示すように、ステップＳ１１０１にて、亀裂３ａの代表の深度Ｚａにおいて亀裂３ａに近似する正弦波形を取得する。具体的には、ステップＳ５０３（図５参照）にて各深度Ｚに対して代表Ｇｒａｙ輝度値を算出する走向θ及び傾斜φが取得されているので、代表の深度Ｚａで取得された走向θ及び傾斜φを読み出し、関数Ｆ（Ｙ）に代入することにより正弦波形を取得する。

【0060】

次に、ステップＳ１１０２にて、図９（Ａ）に示す上側ピーク点ＰＵａにおける亀裂３ａのＲＧＢ輝度を、カラーの展開画像から抽出する。具体的には、直線Ｌに沿うＲＧＢ輝度を、ステップＳ８０２（図８参照）にて求められる幅ＷＵａの深度区間で抽出する。

【0061】

次に、ステップＳ１１０３にて、ＲＧＢ表色系からＬ＊ａ＊ｂ＊表色系へ変換する。これにより、上側ピーク点ＰＵａにおけるＬ＊値、ａ＊値及びｂ＊値が算出される。

【0062】

次に、ステップＳ１１０４にて、図２に示す下側ピーク点ＰＬａにおける亀裂３ａのＲＧＢ輝度を、カラーの展開画像から抽出する。次に、ステップＳ１１０５にて、下側ピーク点ＰＬａにおけるＲＧＢ輝度に基づいて、Ｌ＊値、ａ＊値及びｂ＊値を算出する。ステップＳ１１０４，Ｓ１１０５における処理は、ステップＳ１１０２，Ｓ１１０３における処理と略同じであるため、ここではその説明を省略する。

【0063】

次に、ステップＳ１１０６にて、上側及び下側ピーク点ＰＵａ，ＰＬａにおけるＬ＊値、ａ＊値、ｂ＊値の平均値をそれぞれ算出し、これらの平均値を、亀裂３ａ，３ｂの代表色情報として取得する。

【0064】

以上により、展開画像から亀裂３ａの代表色情報が取得され、亀裂色情報取得処理が終了する。亀裂３ｂの代表色情報を取得する処理は、亀裂３ａの代表色情報を取得する処理と略同じであるため、ここではその説明を省略する。

【0065】

なお、本実施形態では、上側ピーク点ＰＵａにおける亀裂３ａのＲＧＢ輝度と、下側ピーク点ＰＬａにおける亀裂３ａのＲＧＢ輝度と、の平均値を求めているが、正弦波上の全ての画素でのＲＧＢ輝度の平均値を求めてもよい。

【0066】

このように、亀裂検出装置１００では、コンピュータ２０の処理によって、カラーの展開画像から亀裂３ａ，３ｂの色情報を取得する。したがって、亀裂３ａ，３ｂの色情報を迅速かつ一定の精度で取得することができる。

【0067】

また、亀裂検出装置１００では、ＲＧＢ輝度を、Ｓ８０２（図８参照）にて求められる深度区間で抽出する。そのため、亀裂３ａ，３ｂ以外の部分での色情報を抽出するのを防止することができ亀裂３ａ，３ｂにおける色情報をより正確に取得することができる。

【0068】

亀裂検出装置１００の以上の処理により、亀裂３ａ，３ｂの深度、走向、傾斜、亀裂幅及び色情報が取得される。

【0069】

次に、亀裂連続性判定システム１０００について、図１２から図２０を参照して説明する。亀裂連続性判定システム１０００は、亀裂検出装置１００により取得される亀裂情報（亀裂の深度、走向、傾斜、幅及び色情報）を用いて算出される指標や、いわゆるルジオン試験を実施することによって得られるルジオン値（Ｌｕ値）を用いて算出される指標に基づいて、亀裂どうしの連続性を判定し、その結果を表示する。

【0070】

なお、各指標について連続性判定用閾値を予め設定し、連続性判定用閾値を越える指標又は連続性判定用閾値よりも小さい指標をエラー値として排除してもよい。

【0071】

図１２に示すように、亀裂連続性判定システム１０００は、亀裂検出装置１００から亀裂情報を取り込む亀裂情報取込部３１と、亀裂どうしの連続性を判定するための指標を算出する指標算出部３３と、亀裂どうしの連続性を判定する連続性判定部３４と、を備える。指標算出部３３により算出される指標は、亀裂どうしの見かけの深度差、法線ベクトルの交差角度、幅差、色情報差、及びルジオン値差と、を含む。連続性判定部３４は、これらの指標を用いて算出される連続性判定スコアに基づいて、亀裂どうしの連続性を判定する。

【0072】

以下、亀裂どうしの連続性を判定するための指標を算出する処理、及びこれらの指標を用いて算出される連続性判定スコアに基づいて亀裂どうしの連続性を判定する処理について、図１３から図１８を参照して具体的に説明する。図１３は、第１及び第２ボアホール２ｃ，２ｄを示す斜視図である。ここでは、第１ボアホール２ｃにおいて検出される複数の亀裂のうちの１つの亀裂３ｃと、第２ボアホールにおいて検出される複数の亀裂のうちの１つの亀裂３ｄと、の連続性を判定する指標の算出について説明する。

【0073】

まず、指標として亀裂３ｃ，３ｄどうしの見かけ深度差を算出する場合について説明する。

【0074】

ここで、「見かけの深度差」について、図１３を参照して説明する。「見かけの深度差」とは、例えば、第１ボアホール２ｃにおける亀裂３ｃを第２ボアホール２ｄまで直線的に延長した仮想亀裂３ｃ’が第２ボアホール２ｄと交差する交点Ｐｃでの深度と、亀裂３ｄの深度と、の差Ｄｃである。亀裂は、通常、基礎地盤１（図１参照）中において略直線状に延びるので、見かけの深度差Ｄｃが小さいほど亀裂３ｃ，３ｄどうしは連続している可能性が高くなる。換言すれば、亀裂３ｃ，３ｄどうしの見かけの深度差を指標として用いることにより亀裂３ｃ，３ｄどうしの連続性を判定することが可能である。以下、亀裂３ｃ，３ｄどうしの見かけの深度差を算出する場合について、図１３及び図１４を参照して具体的に説明する。

【0075】

まず、ステップＳ１４０１にて、亀裂３ｃの代表深度、走向及び傾斜に基づいて、亀裂３ｃを第２ボアホール２ｄまで直線的に延長する。次に、ステップＳ１４０２にて、延長された仮想亀裂３ｃ’と第２ボアホール２ｄの中心軸との交点Ｐｃの深度を算出する。次に、ステップＳ１４０３にて、交点Ｐｃと亀裂３ｄとの深度差（見かけの深度差）Ｄｃを算出する。

【0076】

ステップＳ１４０１，Ｓ１４０２，Ｓ１４０３は、第２ボアホール２ｄにおける亀裂３ｄに対しても行なわれる。ステップＳ１４０１，Ｓ１４０２，Ｓ１４０３が両方の亀裂３ｃ，３ｄに対して行われたかは、ステップＳ１４０４にて判断される。ステップＳ１４０１，Ｓ１４０２，Ｓ１４０３を亀裂３ｄに対して実行することにより、仮想亀裂３ｄ’と第１ボアホール２ｃの中心軸との交点Ｐｄの深度が算出され、交点Ｐｄと亀裂３ｃとの深度差（見かけの深度差）Ｄｄが算出される。

【0077】

次に、ステップＳ１４０５にて、２つの見かけの深度差Ｄｃ，Ｄｄの平均値を算出する。次に、ステップＳ１４０６にて、見かけの深度差Ｄｃ，Ｄｄの平均値を、亀裂３ｃ，３ｄ間の距離で除す。以上により、正規化された見かけの深度差が得られる。

【0078】

本実施形態では、深度差Ｄｃ、Ｄｄの平均値を見かけの深度差として算出しているが、見かけの深度差をいわゆる面間距離で算出してもよい。面間距離の算出方法を、図２１を参照して具体的説明する。まず、亀裂３ｃの代表深度、走向及び傾斜に基づいて、亀裂３ｃを平面的に延長する。次に、延長された仮想亀裂平面３ｃ’’に、亀裂３ｄの中心から垂線ＮＬｄを延ばし、垂線ＮＬｄの長さＬｄを算出する。同様に、亀裂３ｄを平面的に延長した仮想亀裂平面３ｄ’’に、亀裂３ｃの中心から垂線ＮＬｃを延ばし、垂線ＮＬｃの長さＬｃを算出する。垂線ＮＬｃ，ＮＬｄの長さＬｃ，Ｌｄの平均値を算出し、亀裂３ｃ，３ｄ間の距離で除す。以上により、正規化された面間距離が得られる。

【0079】

次に、指標として、亀裂３ｃの単位法線ベクトルと、亀裂３ｄの単位法線ベクトルと、の交差角度を算出する場合について説明する。亀裂は、通常、基礎地盤１（図１参照）中において略直線状に延びる。そのため、亀裂３ｃ，３ｄの法線ベクトルの方向が近いほど、つまり亀裂３ｃ，３ｄどうしの単位法線ベクトルの交差角度が０に近いほど亀裂３ｃ，３ｄどうしは連続している可能性が高くなる。換言すれば、亀裂３ｃ，３ｄどうしの単位法線ベクトルの交差角度を指標として用いることにより亀裂３ｃ，３ｄどうしの連続性を判定することが可能である。以下、亀裂３ｃ，３ｄどうしの単位法線ベクトルの交差角度を算出する場合について、図１３及び図１５を参照して具体的に説明する。

【0080】

まず、ステップＳ１５０１にて、亀裂３ｃの走向及び傾斜に基づいて、亀裂３ｃの平面の方程式を算出する。次に、ステップＳ１５０２にて、平面の方程式から亀裂３ｃにおける単位法線ベクトルＶｃを算出する。

【0081】

ステップＳ１５０１，Ｓ１５０２は、亀裂３ｄに対しても行われる。ステップＳ１５０１，Ｓ１５０２が両方の亀裂３ｃ，３ｄに対して行われたかは、ステップＳ１５０３にて判断される。ステップＳ１５０１，Ｓ１５０２を亀裂３ｄに対して実行することにより、亀裂３ｄにおける単位法線ベクトルＶｄが算出される。

【0082】

次に、ステップＳ１５０４にて、亀裂３ｃ，３ｄどうしの単位法線ベクトルＶｃ，Ｖｄの交差角度を算出する。以上により、亀裂３ｃ，３ｄの単位法線ベクトルＶｃ，Ｖｄの交差角度が得られる。

【0083】

次に、指標として、亀裂３ｃ，３ｄどうしの幅差を算出する場合について説明する。連続する１つの亀裂においては、通常、第１ボアホール２ｃでの幅と、第２ボアホール２ｄでの幅と、は略一致する。したがって、亀裂３ｃ，３ｄどうしの幅差が小さいほど亀裂３ｃ，３ｄどうしは連続している可能性が高くなる。換言すれば、亀裂３ｃ，３ｄどうしの幅差を指標として用いることにより亀裂３ｃ，３ｄどうしの連続性を判定することが可能である。以下、亀裂３ｃ，３ｄどうしの幅差を算出する場合について、図１３及び図１６を参照して具体的に説明する。

【0084】

まず、ステップＳ１６０１にて、亀裂３ｃの代表幅ＷＮｃと、亀裂３ｄの代表幅ＷＮｄと、を取得する。次に、Ｓ１６０２にて、代表幅ＷＮｃと代表幅ＷＮｄとのうち大きい方の幅を小さい方の幅で除す。以上により、亀裂３ｃ，３ｄどうしの相対的な幅差が得られる。

【0085】

ステップＳ１６０１で取得される代表幅ＷＮｃ，ＷＮｄは、それぞれ、亀裂３ｃ，３ｄの傾斜φにより補正されている。したがって、亀裂３ｃ，３ｄの傾斜φの影響を受けない幅差を取得することができ、後述する亀裂連続性の判定精度が向上する。

【0086】

次に、指標として、亀裂３ｃ，３ｄどうしの色情報差を算出する場合について説明する。連続する１つの亀裂においては、通常、第１ボアホール２ｃでの亀裂の色情報と、第２ボアホール２ｄでの亀裂の色情報と、は略一致する。したがって、亀裂３ｃ，３ｄどうしの色情報差が小さいほど亀裂３ｃ，３ｄどうしは連続している可能性が高くなる。換言すれば、亀裂３ｃ，３ｄどうしの色情報差を指標として用いることにより亀裂３ｃ，３ｄどうしの連続性を判定することが可能である。以下、亀裂３ｃ，３ｄどうしの色情報差を算出する場合について、図１３及び図１７を参照して具体的に説明する。

【0087】

まず、ステップＳ１７０１にて、亀裂３ｃ，３ｄの色情報（Ｌ＊値、ａ＊値及びｂ＊値）を取得する。次に、ステップＳ１７０２にて、亀裂３ｃのＬ＊値と、亀裂３ｄのＬ＊値とのうち大きい方のＬ＊値を小さい方のＬ＊値と比較して、Ｌ＊値の差を算出する。ａ＊値及びｂ＊値についても同様に、ステップＳ１５０３，Ｓ１５０４にて、ａ＊値及びｂ＊値の差をそれぞれ算出する。次に、Ｓ１５０５にて、Ｌ＊値、ａ＊値及びｂ＊値の差の平均値を算出する。以上により、亀裂３ｃ，３ｄの色情報差が得られる。

【0088】

なお、色情報差を相対差とした場合は、正の値のみならず、負の値をとることもあるので、便宜、相対差、絶対差を選択してもよい。

【0089】

次に、指標として、亀裂３ｃ，３ｄどうしのルジオン値差を算出する場合について説明する。ルジオン値とは、地盤が高い水圧の作用下にあるときの水の通しやすさを示す指標であり、ルジオン値が高いほど水の通りやすい（透水性の高い）亀裂があることを意味する。ルジオン値は、いわゆるルジオン試験により求められる。ルジオン試験とは、地盤に形成されたボアホールに１［ＭＰａ］の水圧をかける試験であり、１［ＭＰａ］の水圧をかけた状態で１分間に地盤に浸透する水の量がルジオン値である。ルジオン試験は、ボアホールの所定の深度毎（例えば１ｍ毎）に行われ、ルジオン値は、ボアホールの所定の深度毎に求められる。ルジオン試験により求められたルジオン値は、ルジオン値取込部３２（図１２参照）により亀裂連続性判定システム１０００に取り込まれる。

【0090】

連続する１つの亀裂においては、通常、第１ボアホール２ｃでのルジオン値と、第２ボアホール２ｄでのルジオン値と、は略一致する。したがって、亀裂３ｃ，３ｄどうしのルジオン値差が小さいほど亀裂３ｃ，３ｄどうしは連続している可能性が高くなる。換言すれば、亀裂３ｃ，３ｄどうしのルジオン値差を指標として用いることにより亀裂３ｃ，３ｄどうしの連続性を判定することが可能である。以下、亀裂３ｃ，３ｄどうしのルジオン値差を算出する場合について、図１８を参照して具体的に説明する。

【0091】

まず、ステップＳ１８０１にて、亀裂３ｃ，３ｄの代表深度を取得する。次に、ステップＳ１８０２にて、亀裂３ｃ，３ｄの代表深度におけるルジオン値を取得する。次に、ステップＳ１８０３にて、亀裂３ｃのルジオン値と、亀裂３ｄのルジオン値とを比較し、大きい方のルジオン値を小さい方のルジオン値で除す。これにより、ルジオン値の相対差が算出される。

【0092】

以上により、亀裂３ｃ，３ｄどうしの見かけの深度差、法線ベクトルの交差角度、幅差、色情報差及びルジオン値差が指標として算出される。

【0093】

次に、亀裂３ｃ，３ｄどうしの見かけの深度差、法線ベクトルの交差角度、幅差、色情報差及びルジオン値差を用いて亀裂どうしの連続性を判定する処理について、図１９を参照して説明する。連続性を判定する処理は、連続性判定部３４（図１２参照）にて行われる。

【0094】

まず、ステップＳ１９０１にて、亀裂３ｃ，３ｄどうしの見かけの深度差、法線ベクトルの交差角度、幅差、色情報差及びルジオン値差の累積相対度数分布を作成する。次に、ステップＳ１９０２にて、各亀裂３ｃ，３ｄの各指標について累積度を算出する。次に、ステップＳ１９０３にて、累積度の重み付け平均値を算出し、連続性判定スコアとして取得する。

【0095】

次に、ステップＳ１９０４にて、連続性判定スコアに基づいて、亀裂３ｃ，３ｄどうしの連続性を判定する。具体的には、連続性判定スコアが予め定められた連続性判定用閾値以上の場合には、亀裂３ｃ，３ｄどうしが連続していると判定し、連続性判定用閾値未満の場合には、亀裂３ｃ，３ｄどうしは連続していないと判定する。

【0096】

以上により、亀裂３ｃ，３ｄどうしの連続性を判定する処理が終了する。

【0097】

判定結果出力部３５（図１２参照）は、亀裂どうしの連続性を判定した結果を、モニターや紙といった媒体に表示する。図２０は、表示の一例を示す。図２０に示す例では、ボアホール２ｅにおいて検出された亀裂３ｅ１，３ｅ２，３ｅ３，３ｅ４，３ｅ５と、ボアホール２ｆにおいて検出された亀裂３ｆ１，３ｆ２，３ｆ３，３ｆ４の連続性の判定結果を表示している。具体的には、亀裂３ｅ２，３ｆ３どうしは連続していると判定され、亀裂３ｅ１，３ｅ３，３ｅ４，３ｅ５，３ｆ１，３ｆ２，３ｆ４は他のどの亀裂とも連続していないと判定された結果を表示している。

【0098】

まず、連続していないと判定された亀裂３ｅ１，３ｅ３，３ｅ４，３ｅ５，３ｆ１，３ｆ２，３ｆ４に関する表示について、亀裂３ｅ１を代表して説明する。判定結果出力部３５は、亀裂３ｅ１に関する連続性判定の結果として、亀裂３ｅ１の代表深度において円板状の亀裂が存在すると仮想して、円板状の仮想亀裂４ｅを表示する。仮想亀裂４ｅの中心、走向及び傾斜は、亀裂３ｅ１の中心、走向及び傾斜と一致するように設定される。仮想亀裂４ｅの半径は、仮想亀裂４ｅの外周がボアホール２ｅの内壁を超えボアホール２ｆに達しない大きさに設定される。

【0099】

亀裂３ｅ３，３ｅ４，３ｅ５，３ｆ１，３ｆ２，３ｆ４に関する連続性判定の結果についても、仮想亀裂４ｅと同様に表示される。

【0100】

次に、互いに連続すると判定された亀裂３ｅ２，３ｆ３に関する表示について説明する。判定結果出力部３５は、亀裂３ｅ２，３ｆ３に関する連続性判定の結果として、ボアホール２ｅからボアホール２ｆに渡る円板状の亀裂が存在すると仮想して、円板状の仮想亀裂４ｅｆを表示する。仮想亀裂４ｅｆの中心は、亀裂３ｅ２の中心と亀裂３ｆ３の中心を結ぶ線分の中点に設定される。仮想亀裂４ｅｆの走向及び傾斜は、それぞれ、亀裂３ｅ２と亀裂３ｆ３の走向及び傾斜のそれぞれの平均に設定される。仮想亀裂４ｅｆの半径は、仮想亀裂４ｅｆの外周がボアホール２ｅ，２ｆを超え他のボアホール（図示省略）には達しない大きさに設定される。

【0101】

このように、互いに連続すると判定された亀裂３ｅ２，３ｆ３に関しては、ボアホール２ｅからボアホール２ｆに渡る仮想亀裂４ｅｆが表示されるので、連続する亀裂の大きさ及び位置を容易に把握することができる。

【0102】

図示を省略するが、３つ以上のボアホールにおいて検出される３つ以上の亀裂が連続すると判定された場合には、これらの亀裂を含むように円板状の仮想亀裂を表示する。仮想亀裂は、円板状に限られず、矩形曲面として表示してもよい。

【0103】

以上のように、亀裂連続性判定システム１０００では、コンピュータ２０の処理によって、亀裂検出装置１００により検出された亀裂どうしの連続性が判定される。したがって、ボアホール２における亀裂どうしの連続性を迅速かつ一定の精度で検出することができ、基礎地盤１（図１参照）中の亀裂を推定することができる。これにより、より確実に基礎地盤１中の亀裂を塞ぐようにグラウチングの仕様を決定することが可能になる。

【0104】

例えば、図２０に示す例において、ボアホール２ｅ，２ｆの間に新たにボアホールを形成してグラウトを注入することによって、仮想亀裂４ｅｆに相当する亀裂を塞ぐといったグラウチングの仕様を提案することができるようになる。これにより、ボアホール２ｅ，２ｆの間に存在すると推定される透水性の高い仮想亀裂４ｅｆに相当する亀裂をグラウトにより塞ぐことができ、基礎地盤１の遮水性を高めることができる。

【0105】

以上の実施形態によれば、以下に示す作用効果を奏する。

【0106】

亀裂検出装置１００では、展開画像から亀裂に基づいて正弦波形が抽出され、ボアホール２の内壁での亀裂の走向及び傾斜が正弦波形から取得される。そのため、迅速かつ一定の精度で亀裂を検出することができると共に亀裂の走向及び傾斜を取得することができる。

【0107】

また、亀裂連続性判定システム１０００では、亀裂どうしの連続性を判定するための指標に基づいて、亀裂どうしの連続性が判定される。したがって、迅速かつ一定の精度で亀裂の連続性を判定することができる。これにより、グラウチングの仕様を容易に決定することができる。

【0108】

また、亀裂連続性判定システム１０００では、亀裂どうしの連続性を判定するための指標に、亀裂検出装置１００により取得された走向及び傾斜に基づいて算出される見かけの深度差及び法線ベクトルの交差角度が含まれる。亀裂の走向及び傾斜を迅速かつ一定の精度で検出できるため、見かけの深度差及び法線ベクトルの交差角度を迅速かつ一定の精度で算出することが可能になる。したがって、亀裂の連続性を迅速に判定できると共にその精度を高めることができる。

【0109】

また、亀裂連続性判定システム１０００では、亀裂幅取得部２４により取得される亀裂の幅に基づいて算出される幅差が指標に含まれる。亀裂の幅を亀裂検出装置１００により迅速かつ一定の精度で検出できるため、幅差を迅速かつ一定の精度で算出することが可能になる。したがって、亀裂の連続性を迅速に判定できると共にその精度を高めることができる。

【0110】

また、亀裂連続性判定システム１０００では、亀裂幅取得部２４により取得される亀裂の幅は、走向傾斜取得部２３により取得される傾斜を用いて補正された幅である。そのため、亀裂の傾斜の影響を受けない幅差を取得することができ、亀裂連続性の判定性の精度が向上する。

【0111】

また、亀裂連続性判定システム１０００では、亀裂色情報取得部２５により取得される亀裂の色情報の差が指標に含まれる。亀裂の色情報を亀裂検出装置１００により迅速かつ一定の精度で検出できるため、色情報の差を迅速かつ一定の精度で算出することが可能になる。したがって、亀裂の連続性を迅速に判定できると共にその精度を高めることができる。

【0112】

また、亀裂連続性判定システム１０００では、亀裂の深度におけるルジオン値差が指標に含まれる。したがって、連続しかつ透水性の高い亀裂を容易にかつ高い精度で判定することができる。

【0113】

以上、本発明の実施形態について説明したが、上記実施形態は本発明の適用例の一部を示したに過ぎず、本発明の技術的範囲を上記実施形態の具体的構成に限定する趣旨ではない。

【0114】

亀裂連続性判定システム１０００では、亀裂どうしの見かけの深度差、法線ベクトルの交差角度、幅差、色情報差、及びルジオン値差の全てを用いて亀裂の連続性を判定しているが、これらの指標のうち少なくとも１つに基づいて亀裂の連続性を判定してもよい。

【符号の説明】

【0115】

１・・・地盤
２，２ｃ，２ｄ，２ｅ，２ｆ・・・ボアホール
３,３ａ,３ｂ、３ｃ，３ｄ，３ｅ，３ｅｃ，３ｅｄ，３ｆ・・・亀裂
１０・・・ボアホールカメラ（画像作成部）
２２・・・亀裂抽出部
２３・・・走向傾斜取得部
２４・・・亀裂幅取得部
２５・・・亀裂色情報取得部
３３・・・指標算出部
３４・・・連続性判定部
１００・・・亀裂検出装置
１０００・・・亀裂連続性判定システム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】