コークスの強度を表すDIを推定するための方法であって、前記コークスの気孔率Pに基づいて、下記式1を満たす強度基準値y1を算出する工程と、前記コークスの真比重Dtに基づいて、下記式2を満たす差分基準値y2を算出する工程と、強度基準値y1と差分基準値y2を加算することにより、DIの推定値を算出する工程とを含むことを特徴とするコークス強度の推定方法。式1y1=-a×P2b×P-c式2y2=d×Dt-ea、b、c、d、eは定数を示す。少量法により測定された前記コークスの見かけ比重Da1に基づいて、下記式3を満たす体積法による見かけ比重の推定値Da2を算出する工程と、推定値Da2及び真比重Dtに基づいて、下記式4を満たす気孔率Pを算出する工程とを含むことを特徴とする請求項1に記載のコークス強度の推定方法。式3Da2=f×Da1-g式4P=1-Da2/Dt×100f、gは定数を示す。コークスの強度を表すDIを推定するためのプログラムであって、前記コークスの気孔率Pに基づいて、下記式1を満たす強度基準値y1を算出する工程と、前記コークスの真比重Dtに基づいて、下記式2を満たす差分基準値y2を算出する工程と、強度基準値y1と差分基準値y2を加算することにより、DIの推定値を算出する工程とをコンピュータに実行させることを特徴とするコークス強度の推定プログラム。式1y1=-a×P2b×P-c式2y2=d×Dt-ea、b、c、d、eは定数を示す。少量法により測定された前記コークスの見かけ比重Da1に基づいて、下記式3を満たす体積法による見かけ比重の推定値Da2を算出する工程と、推定値Da2及び真比重Dtに基づいて、下記式4を満たす気孔率Pを算出する工程とをコンピュータに実行させることを特徴とする請求項3に記載のコークス強度の推定プログラム。式3Da2=f×Da1-g式4P=1-Da2/Dt×100f、gは定数を示す。しかしながら、石炭の性状を用いてコークス強度を推定する方法では、原料となる石炭の性状からの推定であるため、製造条件によって変化するコークス強度を精度よく推定することが困難である。また、原料に関する情報がない場合には、コークス強度を推定すること自体が不可能となる。コークスの性状を用いてコークス強度を推定する方法においても、上記のようにコークスを多孔質物質として評価する方法では、基質部分と気孔部分の性状をそれぞれ別々の試験で測定し破壊強度の測定等を行うため、手間や時間がかかるという問題がある。本発明は、上記実情に鑑みてなされたものであり、コークス強度を精度よく簡単に推定することができるコークス強度の推定方法及び推定プログラムを提供することを目的とする。また、本発明は、少量のコークスを用いて精度よくコークス強度を推定することができるコークス強度の推定方法及び推定プログラムを提供することを目的とする。本願発明者は、コークスの気孔率Pと真比重Dtに着目し、これらのパラメータとDIの実測値との関係を研究した結果、一定の関係が成り立つことを見出した。具体的には、コークスの気孔率Pを変数とする二次関数により、DIの実測値に近似した強度基準値y1を表すことができ、当該強度基準値y1とDIの実測値との差分を、コークスの真比重Dtを変数とする一次関数により表わされる差分基準値y2に近似することができることを見出した。下記式1、式2及び式3は、本願発明者の研究によって見出された式である。また、下記式4は、JIS規格JISK2151に規定された式に基づくものである。本発明に係るコークス強度の推定方法は、コークスの強度を表すDIを推定するための方法であって、前記コークスの気孔率Pに基づいて、下記式1を満たす強度基準値y1を算出する工程と、前記コークスの真比重Dtに基づいて、下記式2を満たす差分基準値y2を算出する工程と、強度基準値y1と差分基準値y2を加算することにより、DIの推定値を算出する工程とを含むことを特徴とする。式1y1=-a×P2b×P-c式2y2=d×Dt-ea、b、c、d、eは定数を示す。このような構成によれば、コークスの気孔率Pと真比重Dtを用いて、コークスの強度を表すDIの推定値を算出することができる。したがって、簡単に測定可能なコークスの気孔率Pと真比重Dtを用いて、コークス強度を精度よく簡単に推定することができる。また、コークスの気孔率Pと真比重Dtは、少量のコークスを用いて測定可能である。したがって、少量のコークスを用いて精度よくコークス強度を推定することができる。前記コークス強度の推定方法は、少量法により測定された前記コークスの見かけ比重Da1に基づいて、下記式3を満たす体積法による見かけ比重の推定値Da2を算出する工程と、推定値Da2及び真比重Dtに基づいて、下記式4を満たす気孔率Pを算出する工程とを含むことが好ましい。式3Da2=f×Da1-g式4P=1-Da2/Dt×100f、gは定数を示す。このような構成によれば、少量のコークスを用いた少量法によりコークスの見かけ比重Da1を測定するだけで、体積法による見かけ比重の推定値Da2を算出し、当該推定値Da2を用いてコークスの気孔率Pを算出することができる。このようにして算出したコークスの気孔率Pを用いてDIの推定値を算出することにより、少量のコークスを用いて精度よくコークス強度を推定することができる。本発明に係るコークス強度の推定プログラムは、コークスの強度を表すDIを推定するためのプログラムであって、前記コークスの気孔率Pに基づいて、下記式1を満たす強度基準値y1を算出する工程と、前記コークスの真比重Dtに基づいて、下記式2を満たす差分基準値y2を算出する工程と、強度基準値y1と差分基準値y2を加算することにより、DIの推定値を算出する工程とをコンピュータに実行させることを特徴とする。式1y1=-a×P2b×P-c式2y2=d×Dt-ea、b、c、d、eは定数を示す。前記コークス強度の推定プログラムは、少量法により測定された前記コークスの見かけ比重Da1に基づいて、下記式3を満たす体積法による見かけ比重の推定値Da2を算出する工程と、推定値Da2及び真比重Dtに基づいて、下記式4を満たす気孔率Pを算出する工程とをコンピュータに実行させることが好ましい。式3Da2=f×Da1-g式4P=1-Da2/Dt×100f、gは定数を示す。図1は、本発明の一実施形態に係るコークス強度の推定方法の一例を示すフローチャートである。本実施形態では、コークスの強度を表す指数としてJIS規格JISK2151に規定されたDIドラム強度指数を推定する方法について説明する。コークス強度を推定する際には、まず、コークスの真比重Dtを測定するステップS101。真比重Dtは、空隙を含まない比重を意味しており、少量のコークスをサンプルとして用いて、周知の方法により容易に測定することができる。真比重Dtの測定は、JIS規格JISK2151に準拠して行うことが可能である。ただし、コークスの真比重Dtを測定する工程ステップS101を行うことなく、予め定められた真比重Dtの値を用いてコークス強度を推定してもよい。その後、コークスの真比重Dtに基づいてコークスの気孔率Pを算出するステップS102。気孔率Pは、P=1-Da/Dt×100で表すことができる。ここで、Daはコークスの見かけ比重である。見かけ比重Daは、空隙を含む比重を意味しており、少量のコークスをサンプルとして用いて、周知の方法により容易に測定することができる。JIS規格JISK2151に規定された見かけ比重Daの測定方法としては、体積法と少量法とがある。体積法では約8kgのコークスが必要となるが、少量法では約250gのコークスで測定が可能である。本実施形態では、コークスの気孔率Pに基づいて、DIの推定値の基準となる強度基準値y1を算出するとともにステップS103、コークスの真比重Dtに基づいて、強度基準値y1とDIの実測値との差分の基準となる差分基準値y2を算出するステップS104。これらの工程ステップS103、S104は、どちらを先に行ってもよい。強度基準値y1としては、下記式1を満たす値を算出する。また、差分基準値y2としては、下記式2を満たす値を算出する。式1y1=-a×P2b×P-c式2y2=d×Dt-ea、b、c、d、eは定数を示す。上記式1及び式2について、表及び図面を用いて説明する。下記表1は、複数種類のコークスのサンプルを用いて、DIの実測値DI15015、体積法による見かけ比重、真比重及び気孔率などを測定した結果を示している。図2は、表1の結果から導き出されたDIの実測値と気孔率との関係を示すグラフである。表1における各サンプルのDIの実測値平均値と気孔率との関係は、図2のようにプロットすることができる。この結果から、図2に実線で示すように、気孔率を変数とする二次関数により、DIの実測値に対する近似値を表すことができる。この近似値が、上記式1に示した強度基準値y1となる。この例では、上記式1における各定数a、b、cの値が、a=3.9249、b=345.15、c=7500.4となっている。図3は、強度基準値y1とDIの実測値との差分に対する真比重の関係を示すグラフである。表1における各サンプルの真比重と上記差分との関係は、図3のようにプロットすることができる。この結果から、図3に実線で示すように、真比重を変数とする一次関数により、強度基準値y1とDIの実測値との差分に対する近似値を表すことができる。この近似値が、上記式2に示した差分基準値y2となる。この例では、上記式2における各定数d、eの値が、d=67.946、e=128.44となっている。各定数a、b、c、d、eの値は、複数種類のコークスのサンプルを用いて予め求めることができるが、上記値に限定されるものではない。算出された強度基準値y1と差分基準値y2は加算され、これによりDIの推定値が算出される図1のステップS105。図4は、算出されたDIの推定値とDIの実測値との関係を示すグラフである。表1における各サンプルについて、上記のようにして算出されたDIの推定値とDIの実測値平均値との関係は、図4のようにプロットすることができる。この結果から、図4に実線で示すように、算出されたDIの推定値とDIの実測値とがほぼ同等であることが分かる。このように、本実施形態では、コークスの気孔率Pと真比重Dtを用いて、コークスの強度を表すDIの推定値を算出することができる。したがって、簡単に測定可能なコークスの気孔率Pと真比重Dtを用いて、コークス強度を精度よく簡単に推定することができる。また、コークスの気孔率Pと真比重Dtは、少量のコークスを用いて測定可能である。したがって、少量のコークスを用いて精度よくコークス強度を推定することができる。図5は、気孔率Pの算出方法の一例を示すフローチャートである。この例では、少量法によりコークスの見かけ比重Da1を測定しステップS201、測定した見かけ比重Da1に基づいて、気孔率Pを算出する方法について説明する。まず、測定した見かけ比重Da1に基づいて、体積法による見かけ比重の推定値Da2を算出するステップS202。このとき、体積法による見かけ比重の推定値Da2としては、下記式3を満たす値を算出する。式3Da2=f×Da1-gf、gは定数を示す。その後、コークスの真比重Dtと、算出された体積法による見かけ比重の推定値Da2とに基づいて、下記式4を満たす気孔率Pを算出するステップS203。式4P=1-Da2/Dt×100上記式3について、表及び図面を用いて説明する。下記表2は、2つのサンプルを用いて、一方のサンプルの見かけ比重を少量法により測定するとともに、他方のサンプルの見かけ比重を体積法により測定し、それらの値から算出されたDIの推定値とDIの実測値との関係を比較した結果を示している。図6は、少量法による見かけ比重の実測値と体積法による見かけ比重の実測値との関係を示すグラフである。このグラフに示すように、少量法による見かけ比重の実測値を変数とする一次関数により、体積法による見かけ比重の実測値に対する近似値を表すことができる。この近似値が、上記式3に示した体積法による見かけ比重の推定値Da2となる。この例では、上記式3における各定数f、gの値が、f=1.1519、g=0.1022となっている。各定数f、gの値は、複数種類のコークスのサンプルを用いて予め求めることができるが、上記値に限定されるものではない。上記表2のサンプルLについては、少量法により測定した見かけ比重Da1の値1.036に基づいて、上記式3を満たす値1.09を体積法による見かけ比重の推定値Da2として算出した。このようにして算出した推定値Da2と、予め測定された真比重Dtとに基づいて、上記式4を満たす気孔率Pを算出し、当該気孔率Pを用いてDIの推定値を算出した。その結果、上記表2の通り、算出されたDIの推定値は、DIの実測値とほぼ同等であった。上記表2のサンプルMについては、体積法により見かけ比重Daを測定し、その見かけ比重と、予め測定された真比重Dtとに基づいて、P=1-Da/Dt×100の式により気孔率Pを算出した。当該気孔率Pを用いてDIの推定値を算出した結果、上記表2の通り、算出されたDIの推定値は、DIの実測値とほぼ同等であった。上述の式3及び式4を用いれば、少量のコークスを用いた少量法によりコークスの見かけ比重Da1を測定するだけで、体積法による見かけ比重の推定値Da2を算出し、当該推定値Da2を用いてコークスの気孔率Pを算出することができる。このようにして算出したコークスの気孔率Pを用いてDIの推定値を算出することにより、少量のコークスを用いて精度よくコークス強度を推定することができる。図1におけるステップS102~S105、及び、図5におけるステップS202、S203の各工程は、それぞれコンピュータがプログラムを実行することにより行うことができる。この場合、プログラムを実行するコンピュータは、コークス強度の推定装置として機能することとなる。ただし、上記工程の少なくとも一部を作業者が行うこともできる。また、図1及び図5に示す全ての工程を自動化することも可能である。