"内部にガスが貯留されるライナーと、該ライナーの外周を覆うシェルとを備えた圧力容器の内部温度の推定値を、演算装置を用いて算出する圧力容器の内部温度推定方法であって、前記演算装置が、ライナー温度センサにより検出された前記ライナーの温度を、予め設定された前記ライナーの温度と前記圧力容器の内部温度との相関データに適用して、該検出された前記ライナーの温度に対応する前記圧力容器の内部温度である、対応内部温度を算出する対応内部温度算出工程と、前記演算装置が、前記対応内部温度算出工程で算出された前記対応内部温度を、前記圧力容器の内部圧力を検出する圧力センサの検出圧力と、前記圧力容器の周囲温度を検出する周囲温度センサの検出圧力とのうちの少なくともいずれか一方を用いて補正して、前記圧力容器の内部温度の推定値を算出する内部温度推定値算出工程とを含むことを特徴とする圧力容器の内部温度推定方法。請求項1記載の圧力容器の内部温度推定方法において、前記内部温度推定値算出工程における前記内部圧力を用いた補正は、前記演算装置が、所定期間の開始時点での前記圧力センサの検出圧力と、該所定期間の終了時点での前記圧力センサの検出圧力と、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度とに基づいて、該所定期間での前記圧力容器の内部圧力の変化に応じた圧力対応補正値を算出し、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度に、該圧力対応補正値を加算するものであることを特徴とする圧力容器の内部温度測定方法。請求項1又は請求項2記載の圧力容器の内部温度推定方法において、前記内部温度推定値算出工程における前記周囲温度を用いた補正は、前記演算装置が、所定期間における前記周囲温度センサの検出温度と、該所定期間の開始時点での前記ライナー温度センサの検出温度との温度差、又は、該所定期間における前記周囲温度センサの検出温度と、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度との温度差に基づいて、前記圧力容器の周囲温度と内部温度との温度差に応じた周囲温度対応補正値を算出し、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度に、該周囲温度対応補正値を加算するものであることを特徴とする圧力容器の内部温度測定方法。請求項1から請求項3のうちのいずれか1項記載の圧力容器の内部温度測定方法において、前記内部温度推定値算出工程における前記内部圧力と前記周囲温度とを用いた補正は、前記演算装置が、前記圧力センサの検出圧力を用いて算出した圧力対応補正値と、前記周囲温度センサの検出圧力を用いて算出した周囲温度対応補正値とを、所定の重み付けを行って加算した総補正値を、前記対応内部温度算出工程で算出された前記対応内部温度に加算するものであることを特徴とする圧力容器の内部温度推定方法。本発明は、圧力容器のライナーの検出温度に基づく内部温度の推定精度を、さらに高めた圧力容器の内部温度推定方法を提供することを目的とする。本発明は上記目的を達成するためになされたものであり、内部にガスが貯留されるライナーと、該ライナーの外周を覆うシェルとを備えた圧力容器の内部温度の推定値を、演算装置を用いて算出する圧力容器の内部温度推定方法に関する。そして、前記演算装置が、ライナー温度センサにより検出された前記ライナーの温度を、予め設定された前記ライナーの温度と前記圧力容器の内部温度との相関データに適用して、該検出された前記ライナーの温度に対応する前記圧力容器の内部温度である、対応内部温度を算出する対応内部温度算出工程と、前記演算装置が、前記対応内部温度算出工程で算出された前記対応内部温度を、前記圧力容器の内部圧力を検出する圧力センサの検出圧力と、前記圧力容器の周囲温度を検出する周囲温度センサの検出圧力とのうちの少なくともいずれか一方を用いて補正して、前記圧力容器の内部温度の推定値を算出する内部温度推定値算出工程とを含むことを特徴とする第1発明。第1発明によれば、前記対応内部温度算出工程において、前記演算装置が、前記ライナーの検出温度に対応する前記圧力容器の内部温度を算出する。そして、前記圧力容器の内部圧力が変化すると、この変化に応じて前記圧力容器の内部温度も変化する。また、前記圧力容器の内部温度と周囲温度との差があると、前記圧力容器の内部と外部との間の熱交換により、前記圧力容器の内部温度が変化する。そこで、前記内部温度推定値算出工程において、前記演算装置は、前記圧力容器の内部圧力を検出する圧力センサの検出圧力と、前記圧力容器の周囲温度を検出する周囲温度センサの検出圧力とのうちの少なくともいずれか一方を用いて補正して、前記圧力容器の内部温度の推定値を算出する。これにより、前記圧力容器の内部温度の推定値を前記ライナー温度センサの検出温度に基づいて算出するときの推定精度を高めることができる。また、第1発明において、前記内部温度推定値算出工程における前記内部圧力を用いた補正は、前記演算装置が、所定期間の開始時点での前記圧力センサの検出圧力と、該所定期間の終了時点での前記圧力センサの検出圧力と、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度とに基づいて、該所定期間での前記圧力容器の内部圧力の変化に応じた圧力対応補正値を算出し、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度に、該圧力対応補正値を加算するものであることを特徴とする第2発明。第2発明によれば、前記所定期間における前記圧力容器の内部圧力の変化に応じた前記圧力対応補正値を用いることにより、内部圧力の変化に伴う前記圧力容器の内部温度の上昇/下降分を補正して、前記圧力容器の内部温度の推定精度を高めることができる。また、第1発明又は第2発明において、前記内部温度推定値算出工程における前記周囲温度を用いた補正は、前記演算装置が、所定期間における前記周囲温度センサの検出温度と、該所定期間の開始時点での前記ライナー温度センサの検出温度との温度差、又は、該所定期間における前記周囲温度センサの検出温度と、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度との温度差に基づいて、前記圧力容器の周囲温度と内部温度との温度差に応じた周囲温度対応補正値を算出し、該所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度に、該周囲温度対応補正値を加算するものであることを特徴とする第3発明。第3発明によれば、前記所定期間における前記周囲温度センサの検出温度と、前記所定期間の開始時点での前記ライナー温度センサの検出温度との温度差、又は、前記所定期間における前記周囲温度センサの検出温度と前記所定期間の開始時点で前記対応内部温度算出工程により算出した前記対応内部温度との温度差に基づく、前記周期温度対応補正値を用いることにより、前記圧力容器の内部と外部間の熱交換に伴う前記圧力容器の内部温度の上昇/下降分を補正して、前記圧力容器の内部温度の推定精度を高めることができる。また、第1発明から第3発明のいずれかにおいて、前記内部温度推定値算出工程における前記内部圧力と前記周囲温度とを用いた補正は、前記演算装置が、前記圧力センサの検出圧力を用いて算出した圧力対応補正値と、前記周囲温度センサの検出圧力を用いて算出した周囲温度対応補正値とを、所定の重み付けを行って加算した総補正値を、前記対応内部温度算出工程で算出された前記対応内部温度に加算するものであることを特徴とする第4発明。第4発明によれば、前記圧力容器の構造や形状、前記圧力容器内のガスの熱容量、前記圧力容器周囲の空気の流れ等に応じて、前記重み付けを行うことにより、前記圧力容器の内部温度の推定精度をさらに高めることができる。本発明の実施形態の一例について、図1~図5を参照して説明する。図1を参照して、水素ガスタンク10本発明の圧力容器に相当するは、例えば燃料電池車両に搭載されて、燃料電池に燃料ガスとして水素ガスを供給するものである。水素ガスタンク10は、ガス遮蔽性を有するライナー11の外側を、耐圧性を有するシェル12で覆って構成されている。ライナー11は、高圧の水素ガスを密封状態で保持するために、例えば、厚さ10mm程度の高密度ポリエチレン等の樹脂材で形成されている。シェル12は、水素ガスタンク10の耐圧性及び機械的強度を確保している。水素ガスタンク10には、水素ガスタンク10への水素ガスの充填と水素ガスタンク10からの水素ガスの供給を行うためのインタンク電磁弁13と、ライナー11の外表面の温度を検出するライナー温度センサ14とが設けられている。インタンク電磁弁13はガス供給管20と接続され、水素ガスタンク10とガス供給管20間の連通と遮断とを切替える。ガス供給管20には、水素ガスタンク10からの水素ガスの供給圧を減少させる減圧弁21と、ガス供給管20内の圧力を検出する圧力センサ22とが設けられている。なお、インタンク電磁弁13の開弁時には、ガス供給管20内の圧力は、水素タンク10の内部圧力と等しくなる。ライナー温度センサ14は、サーミスタ等を用いて構成され、シェル12内に、ライナー11の外表面に接触するように埋設されている。また、水素ガスタンク10の近傍に、水素ガスタンク10の周囲温度を検出する周囲温度センサ30が配置されている。演算装置50は、CPU,メモリ,入出力インターフェース回路等、或いはこれらを集約したマイクロコンピュータ等により構成された電子ユニットであり、メモリに保持された水素ガスタンク10の内部温度推定用のプログラムを、CPUで実行することにより、水素ガスタンク10の内部温度の推定値を算出する。ライナー温度センサ14、圧力センサ22、及び周囲温度センサ30は、演算装置50に接続されている。演算装置50は、所定の制御周期本発明の所定期間に相当する毎に、圧力センサ22の圧力検出信号Ps、ライナー温度センサ14の温度検出信号TLs、及び、周囲温度センサ30の温度検出信号TCsを入力する。そして、演算装置50は、各制御周期毎に、これらの検出信号に基づいて水素ガスタンク10の内部温度の推定値を算出する。演算装置50は、対応内部温度算出工程と、内部温度推定値算出工程とを実行して、水素ガスタンク10の内部温度の推定値を算出する。以下、演算装置50による水素タンク10の内部温度の推定値の算出方法について説明する。[対応内部温度算出工程]図2は、制御周期間における、水素タンク10のライナー11の温度変化量ΔTaと、ライナー11の温度に対する水素タンク10の内部温度の補正量ALとの相関マップを示したものであり、横軸がライナー10の温度の変化量ΔTaに設定され、縦軸が水素タンク10の内部温度の補正量ALに設定されている。図2の相関マップ関係は、実験或いはコンピュータシミュレーション等により決定されたものであり、図2の相関マップのデータ本発明の相関データに相当するは、演算装置50のメモリに保持されている。なお、相関マップに変えて、水素タンク10のライナー11の温度変化量ΔTaと、ライナー11の温度に対する水素タンク10の内部温度の補正量ALとの相関式を用いてもよい。演算装置50は、各制御周期毎に、ライナー温度変化量ΔTaを算出して、対応する内部温度の補正量ALを求め、以下の式1により対応内部温度Tinrを算出する。Tinr=TL0AL1但し、Tinr:対応内部温度、TL0:制御周期の開始時点でのライナー11の検出温度、AL:ライナー温度変化に基づく内部温度の補正量。次に、演算装置50は、以下の[内部温度推定値算出工程]を実行する。[内部温度推定値算出工程]図3aは、水素ガスタンク10の内部と周囲との熱交換がない完全理想断熱場合の、水素ガスタンク10の内部温度の変化を示したものであり、この場合は、今回の制御周期における水素ガスタンク10の内部温度T2は、以下の式2で求めることができる。T2=T1[P2/P1]κ2但し、T2:今回の制御周期における水素ガスタンク10の内部温度、T1:前回の制御周期における水素ガスタンク10の内部温度、P2:今回の制御周期における水素ガスタンク10の内部圧力、P1:前回の制御周期における水素ガスタンク10の内部圧力、κ:断熱指数H2の断熱指数は1.409。また、図3bは、水素ガスタンク10の内部と周囲との熱交換が十分に行われる場合を示している。この場合は、熱交換により水素ガスタンク10の内部温度は、水素ガスタンク10の周囲温度に等しくなるので、T2=T1となる。また、図4は、水素ガスタンク10の内部と周囲との熱交換により、水素ガスタンク10の内部温度がポリトロープ変化する場合を示している。この場合は、圧力比P2/P1毎に異なるポリトロープ指数を用いて、今回の制御周期における水素ガスタンク10の内部温度を、以下の式3により算出することができる。T2=T1ΔTb=T1[P2/P1]n3但し、T2:今回の制御周期における水素ガスタンク10の内部温度、T1:前回の制御周期における水素ガスタンク10の内部温度、ΔTb:前回の制御周期における水素ガスタンク10の内部温度と今回の制御周期における水素ガスタンク10の内部温度との温度差、P2:今回の制御周期における水素ガスタンク10の内部圧力、P1:前回の制御周期における水素ガスタンク10の内部圧力、n:ポリトロープ指数1<n<κn-1}=Ap4そこで、演算装置50は、上記式4で算出したΔTbを、圧力対応補正値Apとする。図5は、ポリトロープ指数nと圧力比P2/P1との相関マップであり、縦軸がポリトロープ指数nに設定され、横軸が圧力比P2/P1に設定されている。演算装置50は、今回の制御周期における水素ガスタンク10の内部圧力P2と、前回の制御周期における水素ガスタンク10の内部圧力P1との比P2/P1を、図5の相関マップに適用して、対応するポリトロープ指数nを求める。例えば、制御周期が10秒であって、水素ガスタンク10からの水素ガスの放出により、水素ガスタンク10の内部温度と内部圧力が、20°C,35MPaから19.5°C,34MPaに変化したときには、圧力比P2/P1が約0.97、ポリトロープ指数nは約1.06となる。なお、図5の相関マップは、実験やコンピュータシミュレーション等により決定され、この相関マップのデータは演算装置50のメモリに保持されている。また、相関マップではなく、ポリトロープ指数nと圧力比P2/P1間の相関式を用いてもよい。次に、演算装置50は、今回の制御周期における周囲温度センサ30の検出温度TC2と、前回の制御周期における水素ガスタンク10の内部温度推定値Tine1との温度差ΔTc、又は今回の制御周期における周囲温度センサ30の検出温度TCs2と、今回の制御周期の開始時点におけるライナー11の検出温度TLs1との温度差ΔTdに基づいて、以下の式5又は式6により周囲温度対応補正値Atを算出する。At=TCs2-Tine1×α=ΔTc×α5At=TCs2-TLs1×β=ΔTd×β6但し、At:周囲温度対応補正値At、TCs2:今回の制御周期における周囲温度センサ30の検出温度、Tine1:前回の制御周期における水素ガスタンク10の内部温度対応内部温度算出工程により算出された対応内部温度、TLs1:今回の制御周期の開始時点におけるライナー温度センサ14の検出温度、α,β:係数。なお、α,βは、水素ガスタンク10の構造や形状により異なるため、予め実験やコンピュータシミュレーション等を行って決定する。また、水素ガスタンク10内の水素ガスの熱容量や、水素ガスタンク10の周辺の空気の流れによっても、ΔTb,ΔTcは変化する。そこで、水素ガスタンク10内の水素ガスの量や、水素ガスタンク10が搭載された車両の走行速度水素ガスタンク10の周辺の空気の流速に比例に応じて、α,βを決定するマップテーブルを用意して、α,βを適宜変更するようにしてもよい。演算装置50は、上記式4により算出した圧力対応補正値Apと、上記式5又は上記式6により算出した周囲温度対応補正値Atとを用いて、以下の式7により、総補正値Aalを算出する。そして、演算装置50は、以下の式8により、今回の制御周期における水素ガスタンク10の内部温度の推定値Tineを算出する。Aal=Ap×γAc×θ7但し、Aal:総補正値、Ap:圧力対応補正値、Ac:周囲温度対応補正値、γ,θ:係数。Tine=TinrAal8但し、Tine:水素ガスタンク10の内部温度の推定値、Tinr:上記式1により算出された対応内部温度、Aal:総補正値。γとθを、実験やコンピュータシミュレーションにより設定することによって、ライナー11の検出温度に基づく水素ガスタンク10の内部温度の推定精度を高めることができる。なお、上記式7におけるγ,θの設定は、本発明の圧力対応補正値と周囲温度設定補正値の重み付けに相当する。なお、本実施の形態では、内部温度推定値算出工程において、水素ガスタンク10の制御周期における内部圧力比P2/P1の基づく圧力対応補正値Apと、水素ガスタンク10の周囲温度TCsに基づく温度対応補正値Acとを用いて、上記式7により、総補正値Aalを算出した。しかし、圧力対応補正値Ap又は温度対応補正値Atのいずれか一方のみを用いて、対応内部温度Tineの補正値を算出し、この補正値を対応内部温度Tineに加算して、水素ガスタンク50の内部圧力の推定値を算出する場合にも、本発明の効果を得ることができる。"