(57)【要約】【課題】誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することのできる地盤調査位置の決定方法、決定装置、地盤推定方法および地盤推定装置を提供する。【解決手段】事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップＳ１、Ｓ２と、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定するステップＳ３と、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定するステップＳ４と、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定するステップＳ５、Ｓ６とを備えるようにする。【選択図】図１
2019100011.tif 000002
【特許請求の範囲】
【請求項１】
地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する方法であって、
事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップと、
交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定するステップと、
クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定するステップと、
算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定するステップとを備えることを特徴とする地盤調査位置の決定方法。
【請求項２】
追加の地盤調査の位置を決定するステップは、
所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであることを特徴とする請求項１に記載の地盤調査位置の決定方法。
【請求項３】
地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する装置であって、
事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段と、
交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定する手段と、
クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定する手段と、
算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定する手段とを備えることを特徴とする地盤調査位置の決定装置。
【請求項４】
追加の地盤調査の位置を決定する手段は、
所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであることを特徴とする請求項３に記載の地盤調査位置の決定装置。
【請求項５】
請求項１または２に記載の地盤調査位置の決定方法を用いて、所定エリアの地盤特性の空間分布を推定する方法であって、
決定した追加の地盤調査の位置の地盤データを取得するステップと、
この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップとを備えることを特徴とする地盤推定方法。
【請求項６】
請求項３または４に記載の地盤調査位置の決定装置を用いて、所定エリアの地盤特性の空間分布を推定する装置であって、
決定した追加の地盤調査の位置の地盤データを取得する手段と、
この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段とを備えることを特徴とする地盤推定装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、地盤調査位置の決定方法、決定装置、これらを用いた地盤推定方法および地盤推定装置に関し、特に杭基礎構造物の支持層調査ボーリングの配置計画に好適な地盤調査位置の決定方法、決定装置、地盤推定方法および地盤推定装置に関するものである。
【背景技術】
【０００２】
従来、杭基礎建物の杭設計においては一般的に、ボーリングによる地質調査結果から杭の種類や配置を設計した後に、杭打設地点での支持基盤面の深さを推定して杭長の設計を行っている。全ての杭打設地点において調査をすることは現実的にできないため、調査をしていない地点の支持層深度は、限られた調査地点の支持層深度から地層の深度変化の傾向を読み取ることによって推定することになる。
【０００３】
離散的なボーリングデータから支持層深度を推定する手法として、地球統計学の一手法であるクリギング手法（例えば、特許文献１、２を参照）を用いる方法がある。この方法を用いれば、統計的に客観的な予測をすることができる。ただし、予測値に含まれる誤差を定量的に評価できないため、調査が杭設計をするに十分な精度であるかの判断が難しい。また、追加で調査を計画する場合でも、調査の位置や数を決める合理的な手法がないので、一般的には既存の調査位置から距離を離して、調査範囲（杭の配置範囲）に均等となるように配置している。調査数も杭本数に対する比率（一般的には杭本数の５〜１０％程度）で決めているのが現状である。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００４−３４６５７３号公報
【特許文献２】特開２０１５−１６１５９１号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
このため、離散的なボーリングデータから支持層深度をクリギング手法で推定する方法において、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することのできる技術が求められていた。
【０００６】
本発明は、上記に鑑みてなされたものであって、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することのできる地盤調査位置の決定方法、決定装置、地盤推定方法および地盤推定装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記した課題を解決し、目的を達成するために、本発明に係る地盤調査位置の決定方法は、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する方法であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップと、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定するステップと、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定するステップと、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定するステップとを備えることを特徴とする。
【０００８】
また、本発明に係る他の地盤調査位置の決定方法は、上述した発明において、追加の地盤調査の位置を決定するステップは、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであることを特徴とする。
【０００９】
また、本発明に係る地盤調査位置の決定装置は、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する装置であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段と、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定する手段と、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定する手段と、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定する手段とを備えることを特徴とする。
【００１０】
また、本発明に係る他の地盤調査位置の決定装置は、上述した発明において、追加の地盤調査の位置を決定する手段は、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであることを特徴とする。
【００１１】
また、本発明に係る地盤推定方法は、上述した地盤調査位置の決定方法を用いて、所定エリアの地盤特性の空間分布を推定する方法であって、決定した追加の地盤調査の位置の地盤データを取得するステップと、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップとを備えることを特徴とする。
【００１２】
また、本発明に係る地盤推定装置は、上述した地盤調査位置の決定装置を用いて、所定エリアの地盤特性の空間分布を推定する装置であって、決定した追加の地盤調査の位置の地盤データを取得する手段と、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段とを備えることを特徴とする。
【発明の効果】
【００１３】
本発明に係る地盤調査位置の決定方法によれば、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する方法であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップと、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定するステップと、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定するステップと、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定するステップとを備えるので、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができるという効果を奏する。
【００１４】
また、本発明に係る他の地盤調査位置の決定方法によれば、追加の地盤調査の位置を決定するステップは、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであるので、平均誤差指標が最も小さくなる候補位置を求めることで、追加調査の位置や数をより合理的に決定することができるという効果を奏する。
【００１５】
また、本発明に係る地盤調査位置の決定装置によれば、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する装置であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段と、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定する手段と、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定する手段と、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定する手段とを備えるので、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができるという効果を奏する。
【００１６】
また、本発明に係る他の地盤調査位置の決定装置によれば、追加の地盤調査の位置を決定する手段は、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであるので、平均誤差指標が最も小さくなる候補位置を求めることで、追加調査の位置や数をより合理的に決定することができるという効果を奏する。
【００１７】
また、本発明に係る地盤推定方法によれば、上述した地盤調査位置の決定方法を用いて、所定エリアの地盤特性の空間分布を推定する方法であって、決定した追加の地盤調査の位置の地盤データを取得するステップと、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップとを備えるので、杭の支持基盤面などの空間分布を十分な精度で推定することができるという効果を奏する。
【００１８】
また、本発明に係る地盤推定装置によれば、上述した地盤調査位置の決定装置を用いて、所定エリアの地盤特性の空間分布を推定する装置であって、決定した追加の地盤調査の位置の地盤データを取得する手段と、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段とを備えるので、杭の支持基盤面などの空間分布を十分な精度で推定することができるという効果を奏する。
【図面の簡単な説明】
【００１９】
【図１】図１は、本発明に係る地盤調査位置の決定方法、決定装置の実施の形態を示す追加調査の配置算定フロー図である。
【図２】図２は、想定事例を示す図であり、左図は支持基盤面調査対象エリアの平面図、右図は支持基盤面の真値を示すコンター図である。
【図３】図３は、支持基盤面レベル推定状況（ステップＳ１）を示す図である。
【図４】図４は、推定誤差分布の算定状況（ステップＳ２）を示す図である。
【図５】図５は、交差検証（ブラインドテスト）による誤差分布の算定状況（ステップＳ３）を示す図である。
【図６】図６は、追加調査配置のための誤差指標分布の算定状況（ステップＳ４）を示す図である。
【図７】図７は、誤差指標による追加調査位置の選定状況（ステップＳ５）（１カ所ごと逐次算出）を示す図である。
【図８】図８は、誤差指標による追加調査位置の選定状況（ステップＳ５）（１カ所ごと逐次算出／１カ所目）を示す図である。
【図９】図９は、誤差指標による追加調査位置の選定状況（ステップＳ５）（１カ所ごと逐次算出／２カ所目）を示す図である。
【図１０】図１０は、誤差指標による追加調査位置の選定状況（ステップＳ５）（１カ所ごと逐次算出／６０カ所目）を示す図である。
【図１１】図１１は、最適配置と予測される誤差指標分布（ステップＳ６）（６０カ所追加調査を実施する場合）を示す図である。
【図１２】図１２は、６０カ所の調査後の推定支持基盤面分布と誤差指標分布（検証）を示す図である。
【図１３】図１３は、真の分布（真の支持基盤面分布と真の誤差分布（検証））との比較を示す図である。
【図１４】図１４は、一般的な調査計画（クリギング推定誤差を評価指標とした６０カ所の配置計画）との比較を示す図である。
【発明を実施するための形態】
【００２０】
以下に、本発明に係る地盤調査位置の決定方法、決定装置、地盤推定方法および地盤推定装置の実施の形態を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。
【００２１】
本実施の形態は、地球統計学の一手法であるユニバーサル・クリギング手法を用いて、事前調査で取得されたボーリングデータ（地盤データ）から調査範囲（杭の配置エリア）における支持基盤面（地盤特性）の空間特性をモデル化して、そのモデルを基に合理的な追加調査位置の配置と調査後に予測される推定誤差分布を推定するものである。
【００２２】
（地盤調査位置の決定方法）
まず、本実施の形態に係る地盤調査位置の決定方法について説明する。
【００２３】
図１は、本実施の形態による追加調査の配置決定のフローを示したものである。具体的なフローを説明するため、本実施の形態では、図２の左図に示すような複数の杭基礎構造物を有する施設（例えば変電所など）の支持基盤面調査を事例として想定する。この事例では、２００ｍ×１６０ｍの敷地に５棟（Ａｒｅａ−０１〜０５）の杭基礎構造物（杭本数：計５３９本）を計画する場合を想定している。事前調査で取得済みのボーリングデータ（以下、事前ボーリングデータということがある。）はＡｒｅａ−０１〜０５内のＢ０１〜Ｂ１４の１４カ所である。なお、Ａｒｅａ−０１は左上の四角範囲、Ａｒｅａ−０２は左下の四角範囲、Ａｒｅａ−０３は中央下の四角範囲、Ａｒｅａ−０４は中央上右寄りの四角範囲、Ａｒｅａ−０５は右端の長方形範囲である。また、杭本数の内訳は、Ａｒｅａ−０１：１９２本、Ａｒｅａ−０２：４０本、Ａｒｅａ−０３：２１０本、Ａｒｅａ−０４：６４本、Ａｒｅａ−０５：３３本である。
【００２４】
本発明の有効性を示す上で、本事例では支持基盤面の真値として、図２の右図に示すように、敷地中央部分において右上から左下に向かって谷地形を有する分布を想定している。杭位置での支持基盤面深さは、平均値−２２．３（ＤＬｍ）、最浅値−１６．１（ＤＬｍ）、最深値−３０．２（ＤＬｍ）である。
【００２５】
さて、具体的なフローでは、まず図１のステップＳ１、Ｓ２において、事前ボーリングデータを基に、ユニバーサル・クリギング手法により支持基盤面の推定分布および推定誤差分布を算出する。ユニバーサル・クリギング手法では周知のように、空間補間により離散的なデータから任意の位置の推定値と同時に推定誤差を算出することができる。
【００２６】
図３および図４は、支持基盤面レベルをユニバーサル・クリギング手法により推定した結果の例を示したものである。この例では、ユニバーサル・クリギング手法によって、事前ボーリングデータから空間統計モデル（分散および空間相関性）を構築した後に、エリア内を一定間隔（この例では１ｍ×１ｍ）の格子に分割し、格子点座標における支持基盤面レベルを推定するとともに、推定誤差（以下、クリギング推定誤差ということがある。）を算定し、推定値と推定誤差それぞれを図化ソフトウェアによるコンター図として描画している。ユニバーサル・クリギング手法による推定値および推定誤差の計算は、データの空間分布を二次定常確率場（エリア内の平均値、分散および空間相関特性が一定である）と仮定して、事前ボーリングデータを条件とした条件付き確率場を推定するものである。ここで、ユニバーサル・クリギング手法による自己共分散モデルとして、Ｃ（ｈ）＝σ２・ｅｘｐ（−ｈ／ａ）の式で表される指数関数型モデルを用いた。ただし、μ（平均値）＝−２０．７（ｍ）、σ（標準偏差）＝３．３３（ｍ）、ａ（相関距離）＝７２．８ｍである。
【００２７】
その結果として図３、図４に示すような推定値と推定誤差分散の分布が求められる。図３の杭位置における支持基盤面は、平均値−２２．３（ＤＬｍ）、最浅−１６．７（ＤＬｍ）、最深０−２６．７（ＤＬｍ）である。図４の杭位置におけるクリギング推定誤差は、平均値１．８７（ｍ）、最大２．４２（ｍ）、最小０．０（ｍ）である。
【００２８】
なお、図４の推定誤差分布（クリギング推定誤差分布）は、条件付き確率場の分散分布を表しており、空間相関特性が一定であると仮定した下での、事前ボーリングデータの配置のみから決まる値であり、真の誤差とは異なる。しかしボーリングの配置だけから決まるものなので、仮に追加調査位置を配置した場合に、追加調査のボーリングデータが得られる前でも計算することが可能である。
【００２９】
次に、図１のステップＳ３では、事前ボーリングデータを１カ所ごと除いて、その位置をユニバーサル・クリギング手法により推定して真値との誤差を求める、いわゆる交差検証（ブラインドテスト）を実施する。これにより、その位置での真値との差を真の誤差（以下、ブラインド誤差ということがある。）として評価できる。この各位置の真の誤差を基に、エリア全体の誤差分布を描画したものが図５のコンター図（ブラインド誤差分布）である。誤差が大きい位置の事前ボーリングデータは重要度が高いと考えられ、その周辺の推定結果は変動が大きくなる可能性がある。なお、図５に示すように、Ｂ１１のボーリング位置が基盤面の谷部にあることで最も誤差が大きくなり、また、Ｂ０１やＢ１２のボーリング位置などは敷地の端に位置し、外挿領域となることで誤差が大きくなっている。
【００３０】
次に、図１のステップＳ４では、追加調査位置を選定するための指標として、クリギング推定誤差（ステップＳ２）とブラインド誤差（ステップＳ３）を掛けて平方根した値（合成誤差）を算定する。こうして算定された値を誤差指標と呼ぶことにする。この誤差指標は、本発明の根幹となる指標である。図６は、想定事例で求めた誤差指標の分布をコンター図で表したものである。これは、その調査段階での真の推定誤差に近い値を示していると考えられ、この指標を基に杭位置での誤差指標値の平均が最も小さくなるように追加の調査を計画する。なお、図６の杭位置における誤差指標は、平均値１．９６（ｍ）、最大２．５８（ｍ）、最小０．０（ｍ）である。
【００３１】
次に、図１のステップＳ５では、対象となる杭位置（候補位置）での誤差指標の平均（平均誤差指標）が最も小さくなる追加調査位置を１カ所ごと逐次決定（選定）する。具体的には、図７のように調査を追加する前の誤差指標分布から杭位置での平均誤差指標を算出する（この例では１．９６ｍ）。そして、図８に示すように、杭配置の中から１カ所を追加調査位置としたときのクリギング推定誤差分布（ステップＳ２）を再計算して、誤差指標の分布（ステップＳ４）を更新する。ただし、ブラインド誤差分布（ステップＳ３）は追加調査位置におけるボーリングデータが得られないと更新できないので現状の分布を利用して、ステップＳ４の誤差指標の分布を更新する。その際に杭位置での平均誤差指標を算出して、１カ所目の追加調査位置を全杭位置で探索して、最も平均誤差指標が小さくなる位置を１カ所目の最適位置として決定する。なお、この事例ではＡ０１が最適位置となり、平均誤差指標は１．９０ｍと０．０６ｍ減となる。
【００３２】
同様に、図９に示すように、２カ所目の追加調査位置Ａ０２の最適位置を決定する。これを逐次繰り返して、事前に設定した追加数もしくは平均誤差指標値が許容できる誤差に達するまで実施することで、最終的な追加調査計画を決定でき、また調査前の段階で想定される推定誤差分布も得ることができる（ステップＳ６）。図１０、図１１は６０カ所目まで追加調査配置を決定した例を示しており、最終的な平均誤差指標値は１．３２ｍとなり当初の１．９６ｍから０．６４ｍ誤差を低減できると想定される。ただし、この結果は、ブラインド誤差分布が更新されておらず、６０カ所の追加調査後にはブラインド誤差も小さくなっていると考えられるので、過大に評価されている。なお、６０カ所の追加調査位置は、Ａｒｅａ−０１：１９カ所、Ａｒｅａ−０２：６カ所、Ａｒｅａ−０３：２２カ所、Ａｒｅａ−０４：９カ所、Ａｒｅａ−０５：２カ所に分布している。杭本数に対する比率で考えると、Ａｒｅａ−０１、０５が少なめ、Ａｒｅａ−０４が多めに配置されている。
【００３３】
（本発明の効果とその検証）
上記の想定事例では、真の支持基盤面の分布がわかっているものとして実施している。そこで、この結果を基に本発明の効果を説明する。
【００３４】
上記のステップＳ１〜Ｓ６にしたがって、例えば６０カ所の追加調査配置を決定した後に、６０カ所の追加調査を実施した後を想定して、支持基盤面分布の推定と誤差指標分布を算出してコンター図として示したものが図１２である。推定結果のコンター図（左図）を見るとわかるように、真の分布（図２右図）にかなり近い形で、エリア中央部に右上から左下に向かって谷地形を有する分布が明瞭に推定されている。誤差指標分布（右図）では、追加調査前（図１１）と比べてかなり低減されている。なお、杭位置での予測誤差指標は、平均値０．５４ｍ、最大１．４４ｍ、誤差が１ｍ以上となった杭は２６カ所であった。杭位置での平均予測誤差は、１．３２ｍ→０．５４ｍと半分以下となっており、これは事前想定ではデータが得られていないために、ステップＳ３のブラインド誤差を更新できていないことによるもので、調査が進んだことによる低減と考えることができる。
【００３５】
真の支持基盤面分布がわかっているので、図１２の推定結果に含まれる本当の誤差を評価できる。図１３は真の支持基盤面分布と真値の誤差分布を示したコンター図である。杭位置での誤差の平均値は０．２２ｍと、予測値（０．５４ｍ）に比べてさらに半減している。誤差の最大値は１．４０ｍで、誤差が１ｍ以上となった杭は全部で９カ所であった。これは事前の平均誤差指標（０．５４ｍ）を正規分布に従う標準偏差σと考えれば、誤差１ｍ以内である範囲は約±２σとなるので、全体の９５％以上が１ｍ以下の誤差に抑えられることになる。５３９本中５３０本の杭位置で１ｍ以下の誤差に収まっていることから、それ以上の精度で推定できていると考えられる。また、真の誤差が大きくなった杭は、予測誤差分布（図１２）においても、比較的誤差の大きいエリアにあることがわかる。このように、本発明により得られる支持基盤面の推定結果と誤差指標は、実際に想定される結果と概ね一致していることがわかる。
【００３６】
次に、本発明による追加調査の配置による効果を検証するために、一般的な配置、ここではクリギング推定誤差（ステップＳ２）のみを用いて配置を決定する方法と比較した。クリギング推定誤差（ステップＳ２）を基に最適配置を行う場合、基本的には杭配置を考慮しながら距離が均等になるように配置されることになる。図１４は、杭配置を考慮しながら距離が均等になるように６０カ所を配置した結果であり、左図は支持基盤面の推定結果のコンター図、右図は真値との誤差分布のコンター図である。支持基盤面の推定結果は、この配置でも真の分布（図２右図）に近い形で推定されている。真値との誤差分布を見ると、平均値は０．２４ｍであり本発明の方法による誤差（０．２２ｍ）と比べてわずかに大きくなっている。また、誤差の最大値は１．９１ｍ、誤差１ｍ以上の杭は１２カ所と、いずれも本発明の方法に比べて大きいことがわかる。大きな差ではないものの配置の違いにより、推定結果に違いが生じており、本発明の方法の有効性が認められる。
【００３７】
したがって、本実施の形態の地盤調査位置の決定方法によれば、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができる。
【００３８】
（地盤調査位置の決定装置）
次に、本実施の形態に係る地盤調査位置の決定装置について説明する。
【００３９】
本実施の形態に係る地盤調査位置の決定装置は、上述した地盤調査位置の決定方法を装置として具現化したものである。この地盤調査位置の決定装置は、事前ボーリングデータを使用して、ユニバーサル・クリギング手法により杭の支持基盤面分布を推定するとともに、クリギング推定誤差分布を算定する手段と、交差検証によりブラインド誤差分布を算定する手段と、クリギング推定誤差とブラインド誤差から誤差指標分布を算定する手段と、算定した誤差指標分布に基づいて、追加調査の位置を決定する手段とを備える。この地盤調査位置の決定装置は、例えばＣＰＵを有するコンピュータ、メモリ、ディスプレイ、キーボード等のハードウェア、これらハードウェアを用いて実行されるコンピュータプログラム等のソフトウェアにより構成することができる。
【００４０】
本実施の形態の地盤調査位置の決定装置によれば、上述した地盤調査位置の決定方法の場合と同様に、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができる。
【００４１】
（地盤推定方法）
次に、本実施の形態に係る地盤推定方法について説明する。
【００４２】
本実施の形態に係る地盤推定方法は、上述した地盤調査位置の決定方法を用いて、杭の支持基盤面分布を推定する方法である。この方法は、上述した地盤調査位置の決定方法で決定した追加調査の位置のボーリングデータを取得するステップと、この追加調査で取得されたボーリングデータと、事前ボーリングデータを使用して、ユニバーサル・クリギング手法により杭の支持基盤面分布を推定するとともに、クリギング推定誤差分布を算定するステップとを備える。本実施の形態の地盤推定方法によれば、上述したように、支持基盤面分布を十分な精度で推定することができる。
【００４３】
（地盤推定装置）
次に、本実施の形態に係る地盤推定装置について説明する。
【００４４】
本実施の形態に係る地盤推定装置は、上述した地盤調査位置の決定装置を用いて、杭の支持基盤面分布を推定する装置である。この装置は、上述した地盤調査位置の決定装置で決定した追加調査の位置のボーリングデータを取得する手段と、この追加調査で取得されたボーリングデータと、事前ボーリングデータを使用して、ユニバーサル・クリギング手法により杭の支持基盤面分布を推定するとともに、クリギング推定誤差分布を算定する手段とを備える。本実施の形態の地盤推定装置によれば、上述したように、支持基盤面分布を十分な精度で推定することができる。
【００４５】
上記の実施の形態においては、地盤特性の空間分布として杭の支持基盤面を推定する場合を例にとり説明したが、本発明はこれに限るものではなく、例えば地質分布、地下水分布などの他の地盤特性の空間分布を推定する場合にも適用可能であり、この場合にも上記と同様の作用効果を奏することができる。
【００４６】
以上説明したように、本発明に係る地盤調査位置の決定方法によれば、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する方法であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップと、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定するステップと、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定するステップと、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定するステップとを備えるので、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができる。
【００４７】
また、本発明に係る他の地盤調査位置の決定方法によれば、追加の地盤調査の位置を決定するステップは、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであるので、平均誤差指標が最も小さくなる候補位置を求めることで、追加調査の位置や数をより合理的に決定することができる。
【００４８】
また、本発明に係る地盤調査位置の決定装置によれば、地盤調査により取得された離散的な地盤データに基づいて所定エリアの地盤特性の空間分布を推定する場合に、追加の地盤調査の位置を決定する装置であって、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段と、交差検証により、推定に使用した地盤データの取得位置の地盤特性を推定して、クリギング手法により推定された地盤特性との誤差の空間分布を算定する手段と、クリギング手法により算定された推定誤差と、交差検証により算定された誤差とに基づいて、誤差指標の空間分布を算定する手段と、算定した誤差指標の空間分布に基づいて、追加の地盤調査の位置を決定する手段とを備えるので、誤差を定量的に評価できるとともに、追加調査の位置や数を合理的に決定することができる。
【００４９】
また、本発明に係る他の地盤調査位置の決定装置によれば、追加の地盤調査の位置を決定する手段は、所定エリア内に複数の追加調査の候補位置を設定して、クリギング手法により算定された推定誤差の空間分布を候補位置の地盤データに応じて更新することにより誤差指標の空間分布を更新する一方で、全ての候補位置における誤差指標を平均した平均誤差指標を算定し、この平均誤差指標が最も小さくなる候補位置を全ての候補位置の中から探索して、探索された候補位置を追加の地盤調査の位置として決定するものであるので、平均誤差指標が最も小さくなる候補位置を求めることで、追加調査の位置や数をより合理的に決定することができる。
【００５０】
また、本発明に係る地盤推定方法によれば、上述した地盤調査位置の決定方法を用いて、所定エリアの地盤特性の空間分布を推定する方法であって、決定した追加の地盤調査の位置の地盤データを取得するステップと、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定するステップとを備えるので、杭の支持基盤面などの空間分布を十分な精度で推定することができる。
【００５１】
また、本発明に係る地盤推定装置によれば、上述した地盤調査位置の決定装置を用いて、所定エリアの地盤特性の空間分布を推定する装置であって、決定した追加の地盤調査の位置の地盤データを取得する手段と、この追加調査で取得された地盤データと、事前調査で取得された地盤データを使用して、クリギング手法により所定エリアの地盤特性の空間分布を推定するとともに、その推定誤差の空間分布を算定する手段とを備えるので、杭の支持基盤面などの空間分布を十分な精度で推定することができる。
【産業上の利用可能性】
【００５２】
以上のように、本発明に係る地盤調査位置の決定方法、決定装置、地盤推定方法および地盤推定装置は、杭基礎構造物の支持層調査ボーリングの配置計画に有用であり、特に、誤差を定量的に評価するとともに、追加調査の位置や数を合理的に決定するのに適している。
【図１】
2019100011.tif 000003
【図２】
2019100011.tif 000004
【図３】
2019100011.tif 000005
【図４】
2019100011.tif 000006
【図５】
2019100011.tif 000007
【図６】
2019100011.tif 000008
【図７】
2019100011.tif 000009
【図８】
2019100011.tif 000010
【図９】
2019100011.tif 000011
【図１０】
2019100011.tif 000012
【図１１】
2019100011.tif 000013
【図１２】
2019100011.tif 000014
【図１３】
2019100011.tif 000015
【図１４】
2019100011.tif 000016