地盤改良体の強度を推定するための強度推定プログラムであって、コンピュータを、前記地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得手段と、前記地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得手段と、前記地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得手段と、前記未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得手段と、前記供試体強度情報と、前記強度相関長情報と、前記未改良率情報と、前記未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる前記地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成手段と、前記地盤改良体モデル生成工程で生成された各々の前記三次元モデルの強度を算出するモデル強度算出手段と、として機能させることを特徴とする強度推定プログラム。地盤改良体の強度を推定するための強度推定装置であって、前記地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得手段と、前記地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得手段と、前記地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得手段と、前記未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得手段と、前記供試体強度情報と、前記強度相関長情報と、前記未改良率情報と、前記未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる前記地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成手段と、前記地盤改良体モデル生成工程で生成された各々の前記三次元モデルの強度を算出するモデル強度算出手段と、を備えたことを特徴とする強度推定装置。地盤改良体の強度を推定するための強度推定方法であって、前記地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得工程と、前記地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得工程と、前記地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得工程と、前記未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得工程と、コンピュータが、前記供試体強度情報と、前記強度相関長情報と、前記未改良率情報と、前記未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる前記地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成工程と、コンピュータが、前記地盤改良体モデル生成工程で生成された各々の前記三次元モデルの強度を算出するモデル強度算出工程と、を備えたことを特徴とする強度推定方法。前記地盤改良体モデル生成工程では、前記三次元モデルをなす各前記要素の各々の強度が、前記供試体強度情報と、前記強度相関長情報と、前記未改良率情報と、前記未改良部分サイズ情報と、に基づいて設定されることを特徴とする請求項3に記載の強度推定方法。前記強度相関長取得工程では、前記地盤改良体をボーリングして取得した当該地盤改良体のコアサンプルを採取し、前記コアサンプルの外周面に所定間隔で設定した試験位置に対して針貫入試験を行い、前記針貫入試験で得られた各試験位置の針貫入勾配値に基づいて前記相関長を算出することを特徴とする請求項3又は4に記載の強度推定方法。しかしながら、上記推定方法の考え方は、数m規模の地盤改良体の物性とコアサンプルの物性とを同一とみなすものであるので、得られた推定値が真の強度と大きく異なる場合もあると憶測される。このような地盤改良体の強度推定にあっては、更に合理的な手法が望まれていた。そこで、本発明は、地盤改良体の合理的な強度推定を可能にする強度推定プログラム、強度推定装置、及び強度推定方法を提供することを目的とする。本発明の強度推定プログラムは、地盤改良体の強度を推定するための強度推定プログラムであって、コンピュータを、地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得手段と、地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得手段と、地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得手段と、未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得手段と、供試体強度情報と、強度相関長情報と、未改良率情報と、未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成手段と、地盤改良体モデル生成工程で生成された各々の三次元モデルの強度を算出するモデル強度算出手段と、として機能させることを特徴とする。本発明の強度推定装置は、地盤改良体の強度を推定するための強度推定装置であって、地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得手段と、地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得手段と、地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得手段と、未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得手段と、供試体強度情報と、強度相関長情報と、未改良率情報と、未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成手段と、地盤改良体モデル生成工程で生成された各々の三次元モデルの強度を算出するモデル強度算出手段と、を備えたことを特徴とする。本発明の強度推定方法は、地盤改良体の強度を推定するための強度推定方法であって、地盤改良体からサンプリングされた複数の供試体の強度に関する情報を供試体強度情報として取得する供試体強度取得工程と、地盤改良体の強度に関する相関長を強度相関長情報として取得する強度相関長取得工程と、地盤改良体内で部分的に地盤改良がなされていない部分である未改良部分の体積の割合を未改良率情報として取得する未改良率取得工程と、未改良部分の1箇所当たりのサイズを未改良部分サイズ情報として取得する未改良部分サイズ取得工程と、コンピュータが、供試体強度情報と、強度相関長情報と、未改良率情報と、未改良部分サイズ情報と、に基づいて、三次元的に分割された複数の要素からなる地盤改良体の三次元モデルを複数生成する地盤改良体モデル生成工程と、コンピュータが、地盤改良体モデル生成工程で生成された各々の三次元モデルの強度を算出するモデル強度算出工程と、を備えたことを特徴とする。これらの強度推定プログラム、強度推定装置、及び強度推定方法によれば、地盤改良体の三次元モデルが複数生成され、各々の三次元モデルの強度がシミュレーションで算出される。その後、複数の三次元モデルの強度を統計処理することにより、地盤改良体全体としての強度を推定することができる。ここで、本発明者らは、地盤改良体内の各部位においてどの程度の距離範囲で互いに強度の関連性があるか地盤改良体の強度に関する相関長や、地盤改良体内において改良されていない部分未改良部分の空間的な分布状態が、地盤改良体全体としての強度に影響を及ぼすことに着目した。この知見に鑑み、上記三次元モデルの生成においては、地盤改良体からサンプリングされた複数の供試体の強度の情報に加えて、地盤改良体の強度に関する相関長、未改良部分の体積割合、及び未改良部分のサイズの情報を反映させた地盤改良体の三次元モデルが生成される。前述の知見によれば、地盤改良体の強度に関する相関長、未改良部分の体積割合、及び未改良部分のサイズは、何れも、地盤改良体全体としての強度に影響を及ぼすと考えられる。このように、地盤改良体の強度に影響を及ぼす要因を更に加えて導入することにより、現実の地盤改良体の状態がより詳細に反映された三次元モデルが生成される。そして、現実の地盤改良体の状態が詳細に反映された三次元モデルをシミュレーションのモデルとして採用するので、地盤改良体の強度の合理的な推定結果を得ることができる。また、具体的には、本発明の強度推定方法の地盤改良体モデル生成工程では、三次元モデルをなす各要素の各々の強度が、供試体強度情報と、強度相関長情報と、未改良率情報と、未改良部分サイズ情報と、に基づいて設定されることとしてもよい。この構成により、地盤改良体の各部位の強度の分布の統計的性質が再現された三次元モデルを生成することができる。また、本発明の強度推定方法の強度相関長取得工程では、地盤改良体をボーリングして取得した当該地盤改良体のコアサンプルを採取し、コアサンプルの外周面に所定間隔で設定した試験位置に対して針貫入試験を行い、針貫入試験で得られた各試験位置の針貫入勾配値に基づいて相関長を算出することとしてもよい。針貫入試験は、比較的短時間で実行可能であるので、多数の試験位置の針貫入勾配値のデータを採取することができる。また、コアサンプルの所定の試験位置の針貫入勾配値は、当該試験位置における部位の強度に関連するので、針貫入勾配値の分布の傾向を強度の分布の傾向とみなすことができる。従って、地盤改良体の強度に関する相関長の推定値を簡易に得ることができる。本発明の強度推定プログラム、強度推定装置、及び強度推定方法によれば、地盤改良体の合理的な強度推定が可能になる。以下、図面を参照しつつ本発明に係る強度推定プログラム、強度推定装置、及び強度推定方法の実施形態について詳細に説明する。図1は、本発明の一実施形態に係る強度推定装置1の構成を示す図である。強度推定装置1は、図2に示すように、例えば、ジェットグラウト高圧噴射攪拌工法や機械攪拌工法によって改良された地盤改良体100の強度例えば、一軸圧縮強度等を推定するための装置である。なお、以下では、強度推定装置1を用いて、地盤改良体100の一軸圧縮強度を推定する例について説明する。図1に示すように、強度推定装置1は、供試体強度入力部供試体強度取得手段5と、強度相関長取得部強度相関長取得手段7と、未改良率入力部未改良率取得手段9と、未改良部分サイズ入力部未改良部分サイズ取得手段11と、を備えている。更に、強度推定装置1は、モデル生成部改良地盤体モデル生成手段15と、モデル強度算出部モデル強度算出手段17と、統計処理部19と、を備えている。図3は、強度推定装置1のハードウエア構成を示す図である。図3に示されるように、強度推定装置1は、物理的には、CPUCentralProcessingUnit201、主記憶装置であるRAMRandomAccessMemory202、ROMReadOnlyMemory203、ハードディスクなどの補助記憶装置204、ネットワークカードなどのデータ送受信デバイスである通信モジュール205、外部記憶媒体からの情報を読み出す外部記憶媒体読取装置206、キーボードやマウスなどの入力デバイスである入力装置207、ディスプレイ装置などの出力デバイスである出力装置208などのハードウエアにより構成されるコンピュータである。図1に示す供試体強度入力部5、強度相関長取得部7、未改良率入力部9、未改良部分サイズ入力部11、モデル生成部15、モデル強度算出部17、及び統計処理部19等の各構成要素の機能は、強度推定装置1がコンピュータソフトウエアである強度推定プログラム210に従って動作することにより実現される。すなわち、上記各要素の機能は、図2に示すCPU201、RAM202などのハードウエア上に強度推定プログラム210を読み込ませることにより、CPU201の制御のもとで外部記憶媒体読取装置206、入力装置207、出力装置208などを動作させるとともに、RAM202や補助記憶装置204におけるデータの読み出しおよび書き込みを行うことで実現される。なお、例えばDVD等の外部記憶媒体211に電子情報として格納された強度推定プログラム210を、外部記憶媒体読取装置206を経由して補助記憶装置204やRAM202に読み込ませるインストールするようにしてもよい。続いて、図1に示す強度推定装置1の各構成要素の機能について説明する。供試体強度入力部5は、地盤改良体100図2参照からサンプリングされた供試体105の一軸圧縮強度供試体強度情報を入力させる部分である。供試体強度入力部5は、例えば、供試体強度情報の入力画面をディスプレイ装置等に表示し、ユーザによる入力操作を受け付ける。上記の供試体強度情報は、例えば、入力装置207を経由してユーザによるキー入力等で入力されてもよく、通信モジュール205を経由して外部ネットワークから入力されてもよく、外部記憶媒体読取装置206を経由して外部記憶媒体から入力されてもよい。供試体強度入力部5は、取得した供試体強度情報を後述のモデル生成部15に送出する。強度相関長取得部7は、針貫入勾配値入力部7aと強度相関長算出部7bとを有している。針貫入勾配値入力部7aは、地盤改良体100からサンプリングされたコアサンプル103の各試験位置における針貫入勾配値を入力させる部分である。強度相関長取得部7は、例えば、針貫入勾配値等の入力画面をディスプレイ装置等に表示し、ユーザによる入力操作を受け付ける。針貫入勾配値入力部7aにおいては、例えば、各試験位置の位置座標と針貫入勾配値とが関連付けられて入力される。この入力情報は、例えば、入力装置207を経由してユーザによるキー入力等で入力されてもよく、通信モジュール205を経由して外部ネットワークから入力されてもよく、外部記憶媒体読取装置206を経由して外部記憶媒体から入力されてもよい。強度相関長算出部7bにおいては、各試験位置の位置座標と、針貫入勾配値と、に基づいて、地盤改良体100における針貫入勾配値の相関長が算出される。強度相関長取得部7は、取得した針貫入勾配値の相関長を、強度相関長情報として後述のモデル生成部15に送出する。未改良率入力部9は、コアサンプル103の目視検査により得られる未改良率情報を入力させる部分である。ここで、一般に地盤改良体100内には、硬化材との混合不十分等の原因によって、部分的に地盤改良がなされていない部分未改良部分が存在している。地盤改良体100全体に含まれる未改良部分の体積割合を、未改良率を呼び、未改良率を示す情報を未改良率情報と呼ぶものとする。未改良率入力部9は、例えば、未改良率情報の入力画面をディスプレイ装置等に表示し、ユーザによる入力操作を受け付ける。未改良率情報は、例えば、入力装置207を経由してユーザによるキー入力等で入力されてもよく、通信モジュール205を経由して外部ネットワークから入力されてもよく、外部記憶媒体読取装置206を経由して外部記憶媒体から入力されてもよい。未改良率入力部9は、取得した未改良率情報を後述のモデル生成部15に送出する。未改良部分サイズ入力部11は、コアサンプル103の目視検査により得られる未改良部分サイズ情報を入力させる部分である。地盤改良体100全体に含まれる未改良部分の1箇所当たりのサイズ例えば、1箇所当たりの体積を、未改良部分サイズを呼び、未改良部分サイズを示す情報を未改良部分サイズ情報と呼ぶものとする。未改良部分サイズ入力部11は、例えば、未改良部分サイズ情報の入力画面をディスプレイ装置等に表示し、ユーザによる入力操作を受け付ける。未改良部分サイズ情報は、例えば、入力装置207を経由してユーザによるキー入力等で入力されてもよく、通信モジュール205を経由して外部ネットワークから入力されてもよく、外部記憶媒体読取装置206を経由して外部記憶媒体から入力されてもよい。未改良部分サイズ入力部11は、取得した未改良部分サイズ情報を後述のモデル生成部15に送出する。モデル生成部15、モデル強度算出部17及び統計処理部19は、モンテカルロシミュレーションにより、複数の供試体105の一軸圧縮強度等に基づいて地盤改良体100全体の一軸圧縮強度を求める。具体的には、モデル生成部15は、上述の供試体強度情報、強度相関長情報、未改良率情報、及び未改良部分サイズ情報に基づいて演算を行い、地盤改良体100の一軸圧縮強度に関する三次元モデルを生成する。図4に示されるように、生成される三次元モデル100Mは、均質の多数の要素101に分割されてなる。そして、各要素101の一軸圧縮強度は、供試体強度情報、強度相関長情報、未改良率情報、及び未改良部分サイズ情報に基づいて、地盤改良体100の強度分布の統計的性質が再現されるように設定される。モデル強度算出部17は、モデル生成部15で生成された三次元モデル100Mについてシミュレーションを行い、当該三次元モデル100M全体の一軸圧縮強度を算出する。三次元モデル100Mの一軸圧縮強度を求めるためのシミュレーション手法としては、例えば、有限要素法が用いられる。モデル生成部15による三次元モデル100Mの生成と、モデル強度算出部17による一軸圧縮強度の算出と、が複数回繰り返し行われ、多数の三次元モデル100Mの一軸圧縮強度が得られる。統計処理部19は、モデル強度算出部17で得られた多数の三次元モデル100Mの一軸圧縮強度を統計処理する。この統計処理では、例えば、多数の三次元モデル100Mの一軸圧縮強度の平均値と標準偏差とが算出される。統計処理部19は、例えば、ディスプレイ装置等に上記平均値と標準偏差とを表示する。このように得られた平均値は、地盤改良体100の一軸圧縮強度の推定値として用いることができる。続いて、上述の強度推定装置1を用いた地盤改良体100の強度推定方法の一例について、図5を参照しながら説明する。この強度推定方法は、以下に説明する強度相関長取得工程と、供試体強度取得工程と、未改良率取得工程と、未改良部分サイズ取得工程と、地盤改良体モデル生成工程と、モデル強度算出工程と、統計処理工程と、を備えている。強度相関長取得工程地盤改良体100の複数の箇所をボーリングし、複数の円柱状のコアサンプル103を採取する。コアサンプル103は例えば直径55mmであり、公知のコアリング装置を用いて採取される。続いて、得られたコアサンプル103の外周面に、円柱軸方向に所定間隔例えば、5cm間隔で試験位置を設定し、各試験位置に対して針貫入試験を行う図5のS501。針貫入試験は、例えば、地盤工学会基準針貫入試験方法に従って行い、コアサンプル103の所定間隔ごとの各試験位置における針貫入勾配値が得られる。ユーザは、入力装置207キーボード等により、各試験位置の位置座標と針貫入勾配値とを関連付けて強度推定装置1に入力する。なお、試験位置は、コアサンプル103の円柱軸方向に直線的に配列されてもよく、コアサンプル103の外周面上に二次元的に配置されてもよい。強度推定装置1の強度相関長算出部7bにおいては、各試験位置の位置座標と、針貫入勾配値と、に基づいて、地盤改良体100における針貫入勾配値の相関長が算出されるS503。ここで、各試験位置の針貫入勾配値に基づいて、相関長を算出する手法の一例について説明する。図6に示すように、地盤改良体100の各位置から採取された各コアサンプル103図6の例の場合、3つのコアサンプルNo1~No3について、前述の針貫入試験を行う。各コアサンプル103の外周面に設定された多数の試験位置同士について、2点の各試験位置同士の距離横軸と、2点の各試験位置の針貫入勾配値の相違度縦軸との関係をプロットしていき、プロットしたデータの曲線フィッティングを行う。そして、図6に示すように、曲線フィッティングで作成された曲線fが相違度=1に漸近する距離の値を、相関長とする。例えば、図6の例では、相関長=0.36mが求められる。一般に、地盤における針貫入勾配値と一軸圧縮強度との間には高い相関関係があるので、上記の針貫入勾配値の相関長は、地盤改良体100の一軸圧縮強度に関する相関長強度相関長情報に等しいものとみなしてよい。供試体強度入力部5は、取得した上記の強度相関長情報を後述のモデル生成部15に送出する。なお、この強度相関長取得工程では、針貫入試験の結果から一軸圧縮強度に関する相関長を別途途求めておき、この相関長のみを、強度相関長情報として強度推定装置1に入力してもよい。相関長は、地盤改良体100内の各部位においてどの程度の距離範囲で互いに強度の関連性があるかを示す指標である。供試体強度取得工程上記コアサンプル103からランダムに複数の供試体105を抜き取り、各供試体105について一軸圧縮試験を行い、供試体105の一軸圧縮強度を求めるS505。供試体105は、例えば、直径55mm、高さ100mmの円柱とする。一軸圧縮試験は、例えば、JISA1216に従って行う。続いて、強度推定装置1の供試体強度入力部5の入力処理に従って、ユーザは、入力装置207キーボード等により、各供試体の一軸圧縮強度の測定値供試体強度情報を強度推定装置1に入力する。供試体強度入力部5は、取得した供試体強度情報を後述のモデル生成部15に送出する。なお、この供試体強度取得工程では、供試体105の一軸圧縮強度の平均値と標準偏差とを別途求めておき、この平均値及び標準偏差のみを、供試体強度情報として強度推定装置1に入力してもよい。未改良率取得工程続いて、ユーザは、上記コアサンプル103の目視検査を行う。コアサンプル103の目視検査により、コアサンプル103に含まれる未改良部分を容易に確認することができる。すなわち未改良部分は、コアサンプル103のうち一部脆弱な部分、或いはコアサンプル103の一部分が採取時に既に崩壊した部分として目視確認可能である。ユーザは、コアサンプル103全体に対する未改良部分の体積割合を目測で判定し、地盤改良体100に含まれる未改良部分の体積割合の推定値未改良率情報とする。続いて、強度推定装置1の未改良率入力部9の入力処理に従って、ユーザは、入力装置207キーボード等により、未改良率情報を強度推定装置1に入力する。未改良率入力部9は、取得した未改良率情報を後述のモデル生成部15に送出するS507。未改良部分サイズ取得工程また、上記コアサンプル103の目視検査において、ユーザは、上記未改良部分の1箇所当たりのサイズを目測で判定し、地盤改良体100に含まれる未改良部分の1箇所当たりのサイズの推定値未改良部分サイズ情報とする。続いて、強度推定装置1の未改良部分サイズ入力部11の入力処理に従って、ユーザは、入力装置207キーボード等により、未改良部分サイズ情報を強度推定装置1に入力する。未改良部分サイズ入力部11は、取得した未改良部分サイズ情報を後述のモデル生成部15に送出するS507。地盤改良体モデル生成工程続いて、強度推定装置1のモデル生成部15は、上述の供試体強度情報、強度相関長情報、未改良率情報、及び未改良部分サイズ情報に基づいて演算を行い、地盤改良体100の一軸圧縮強度に関する三次元モデル100M図4参照を生成するS511。具体的には、三次元モデル100Mは、例えば、図4に示されるように、均質の多数の要素メッシュ101に分割されてなる。例えば、要素101は5cm立法の直方体とする。各要素101には、当該要素101の一軸圧縮強度と、この一軸圧縮強度に対応する剛性例えば、弾性係数等と、が配分される。すなわち、モデル生成部15は、供試体強度情報、強度相関長情報、未改良率情報、及び未改良部分サイズ情報に基づいて、公知の空間統計学的手法に基づく演算を行い、各要素101に対して、一軸圧縮強度と、剛性に関する物性値例えば、弾性係数と、を割り振る。このとき、地盤改良体100各部位の強度及び剛性の分布の統計的性質が、三次元モデル100Mに再現されるように、各要素101の一軸圧縮強度及び剛性が設定される。図4は、多数の直方体の要素101に分割された三次元モデル100Mを示し、各要素101に付された色の濃さが当該要素101の一軸圧縮強度及び剛性を表すものとする。すなわち、図4中では、色が濃い要素101ほど、大きい一軸圧縮強度と大きい剛性とを有し、色が薄い要素101ほど、小さい一軸圧縮強度と小さい剛性とを有している。例えば、未改良部分の存在割合やサイズが、強度相関長情報、未改良率情報及び未改良部分サイズ情報等で示される統計的性質に整合するように、三次元モデル100Mの要素として未改良部分要素101zが配置される。例えば、未改良部分要素101zの一軸圧縮強度及び剛性はゼロと設定される。例えば、未改良部分要素101zは、前述の未改良率情報で示される出現頻度で、三次元モデル100M中にランダムに配分される。例えば、未改良部分要素100zのサイズ未改良部分要素100zとして割り当てられた一塊の要素の個数は一様としてもよく、未改良部分のサイズに所定の分布例えば正規分布を持たせてもよい。なお、図4の例では、未改良部分要素100zは、要素101の1個分のサイズに設定している。モデル生成部15は、上述の処理を所定回数繰り返して、多数の異なる三次元モデル100Mを、所定のデータ数になるまで生成するS511~S515。ここで生成された多数の三次元モデル100Mは、すべて、実際の地盤改良体100各部位の強度及び剛性の分布の統計的性質を再現したものとなる。モデル強度算出工程強度推定装置1のモデル強度算出部17は、モデル生成部15で生成された三次元モデル100Mについて演算を行い、多数の均質の要素で構成された当該三次元モデル100M全体としての一軸圧縮強度を、有限要素法で算出するS513。その後再び、地盤改良体モデル生成工程S511が行われ、他の三次元モデル100Mが生成される。このように、モデル生成部15による三次元モデル生成S511と、モデル強度算出部17による一軸圧縮強度の算出S513と、が所定回数繰り返し行われることでS515、所定データ数の三次元モデル100Mの一軸圧縮強度が得られる。統計処理工程強度推定装置1の統計処理部19は、多数の三次元モデル100Mの一軸圧縮強度を統計処理するS517。この統計処理では、例えば、多数の三次元モデル100Mの一軸圧縮強度の平均値と標準偏差とが算出される。統計処理部19は、例えば、ディスプレイ装置等に上記平均値と標準偏差とを表示する。このように得られた平均値は、地盤改良体100全体としての一軸圧縮強度の推定値として用いることができる。続いて、上述した強度推定プログラム、強度推定装置1、及び強度推定方法による作用効果について説明する。上述した強度推定プログラム、強度推定装置1、及び強度推定方法によれば、地盤改良体100の三次元モデル100Mが複数生成され、各々の三次元モデル100Mの一軸圧縮強度が有限要素法を用いて算出される。その後、複数の三次元モデル100Mの一軸圧縮強度を統計処理することにより、地盤改良体100全体としての一軸圧縮強度を推定することができる。上記三次元モデル100Mの生成においては、地盤改良体100からサンプリングされた複数の供試体105の強度の情報に加えて、地盤改良体100の強度に関する相関長、未改良部分の体積割合、及び未改良部分のサイズの情報を反映させた三次元モデル100Mが生成される。地盤改良体100の強度に関する相関長、未改良部分の体積割合、及び未改良部分のサイズは、何れも、地盤改良体100全体としての強度に影響を及ぼすと考えられる。このように、地盤改良体100の強度に影響を及ぼす要因を更に加えて導入することにより、現実の地盤改良体100の状態がより詳細に反映された三次元モデル100Mが生成される。そして、現実の地盤改良体100の状態が詳細に反映された三次元モデル100Mを有限要素法のモデルとして採用するので、地盤改良体100の強度の合理的な推定結果を得ることができる。また、モンテカルロシミュレーションにより、現実の地盤改良体100の状態に近い多数の三次元モデルを用いて当該三次元モデル全体の一軸圧縮強度を多数算出するので、多数の供試体の一軸圧縮強度を単純に平均するといった従来の強度推定に比較して、標準偏差も小さくなり、精度が高い推定結果が得られると考えられる。また、強度相関長取得工程では、針貫入試験による各試験位置の針貫入勾配値に基づいて相関長を算出している。針貫入試験は、比較的短時間で実行可能であるので、多数の試験位置の針貫入勾配値のデータを採取することができる。また、コアサンプル103の所定の試験位置における針貫入勾配値は当該試験位置における部位の一軸圧縮強度に関連するので、針貫入勾配値の分布の傾向を一軸圧縮強度の分布の傾向とみなすことができる。従って、地盤改良体の一軸圧縮強度に関する相関長の推定値を簡易に得ることができる。また、針貫入試験は、コアサンプル103に与えるダメージも小さい。従って、針貫入試験後のコアサンプル103から供試体105を抜き取って、供試体強度取得工程における一軸圧縮試験の供試体としても再利用することができる。次に、上述した強度推定プログラム、強度推定装置1、及び強度推定方法による作用効果の確認のため、本発明者らが行った強度推定について説明する。前述の非特許文献1の304頁には、国土交通省の総合技術開発プロジェクト以下総プロによる試験結果が記載されている。総プロの試験では、直径1000mmの地盤改良体の一軸圧縮強度の平均と、当該地盤改良体から採取された直径55mmの供試体の一軸圧縮強度の平均及び標準偏差と、が実測されている。非特許文献1には、砂質土の地盤を改良してなる改良地盤体について、直径1000mmの地盤改良体の一軸圧縮強度の平均=3330kNm2直径55mmの供試体の一軸圧縮強度の平均=3858kNm2直径55mmの供試体の一軸圧縮強度の標準偏差=1165kNm2との実験結果が示されている。地盤改良体の一軸圧縮強度Quと供試体の一軸圧縮強度quとの比Qu/quは、約0.9である。一方、本発明者らの強度推定では、砂質土の地盤をジェットグラウト工法で改良した地盤改良体について、ボーリングを行いコアサンプルを採取した。このコアサンプルの外周面に、円柱軸方向に10cmごとの試験位置を設定し、針貫入試験を行って各試験位置の針貫入勾配値を得た。得られた針貫入勾配値のデータに基づいて針貫入勾配値の相関長を算出し、相関長0.086mを得た。この値を、強度相関長情報として強度推定装置1に入力した。また、一般的に、砂質土の地盤をジェットグラウト工法で改良した場合、未改良率は5程度であると推定されるので、この値を未改良率情報として強度推定装置1に入力した。また、未改良部分サイズは、三次元モデルのメッシュと同サイズ0.05m立方と仮定し、未改良部分サイズ情報として強度推定装置1に入力した。また、前述の総プロの試験結果に基づき、直径55mmの供試体の軸圧縮強度の平均=3858kNm2直径55mmの供試体の軸圧縮強度の標準偏差=1165kNm2の値を、供試体強度情報として強度推定装置1に入力した。以上の入力に基づいて強度推定装置1による地盤改良体の一軸圧縮強度の平均Quを求めたところ、3300kNm2との結果が得られた。従ってQu/quは約0.9との結果が得られ、総プロによる地盤改良体の一軸圧縮強度に近い値が得られた。よって、本実施形態の強度推定プログラム、強度推定装置1、及び強度推定方法による強度推定では、実際の地盤改良体の一軸圧縮強度に近い合理的な結果が得られることが確認された。以上、本発明の一実施形態について説明したが、本発明は、上記実施形態に限られるものではなく、各請求項に記載した要旨を変更しない範囲で変形したものであってもよい。