(57)【特許請求の範囲】
【請求項１】
路面下に埋設した光ファイバ温度レーダにより計測した地中温度計測値に基づき地中側から路面に流入する熱量（Ｗｅ）を求め、
道路近傍に設置した気象計（気温計、湿度計、日射計、雪雨量計、風速計、放射収支計を含む）により計測した大気の温度、日射強度などの気象量計測値と、前記地中温度計測値から算出し、あるいは前記光ファイバ温度レーダなどの温度計測器で計測した路面温度（Ｔｅ）とに基づき大気側から路面に流入する熱量（Ｗａ）を求め、
前記路面温度（Ｔｅ）と、前記気象量計測値の気温計測値、相対湿度計測値、風速計測値とに基づき路面上の水分の蒸発による蒸発熱（Ｗｅｖ）を求め、
地中側から路面に流入する熱量（Ｗｅ）と大気側から路面に流入する熱量（Ｗａ）の差と、前記蒸発熱（Ｗｅｖ）とから湿潤判定値（Ｙ）を求めて、路面の湿潤状態を判定し、前記地中温度計測値及び前記気象量計測値と、前記湿潤判定値（Ｙ）とに基づき路面の状態を推定することを特徴とする路面状態推定方法。
【請求項２】
路面への薬剤散布がない場合は、前記路面温度（Ｔｅ）とその時間当たりの変化量（ΔＴｅ）とから路面の積雪状態を判定し、路面への薬剤散布がある場合は、前記路面温度（Ｔｅ）と前記変化量（ΔＴｅ）と上記蒸発熱（Ｗｅｖ）と積雪の融解または凝固による潜熱量（δＷ）とから路面の積雪状態を判定することを特徴とする請求項１に記載の路面状態推定方法。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、湿潤、積雪、凍結等の路面状態を把握するための方法に係り、特に、低コストで、路面状態推定を自動化し、かつ局地的な判断が可能な路面状態推定方法に関するものである。
【０００２】
【従来の技術】
道路交通において、湿潤、積雪、凍結等の路面状態が生じると、スリップ事故が発生しやすい。このため湿潤、積雪、凍結等の危険な路面状態を速やかに検知して自動車の運転手に注意を促す必要がある。しかし、長距離に及ぶ道路を常に巡回監視するのは、人手、コストの面で困難であるため、このような路面状態を把握する手段として路面状態検出センサが用いられる。
【０００３】
従来の路面状態検出センサは、主に２つに分けられる。
【０００４】
第１のグループは、道路表面（路面）に赤外線または電磁波を照射し、路面からの反射率を計測するものである。このセンサでは、道路近傍に赤外線照射装置、受光装置および信号処理装置を設置し、道路への入射光と反射光との強度比が路面状態によって異なることを用いて路面状態を推定する。
【０００５】
第２のグループとしては、道路近傍の気象量及び路面温度を計測するセンサが挙げられる。路面温度計測には、熱電対、赤外放射温度計等の単点計測用のセンサに加えて、多点計測が可能で長距離の温度を計測するのに適した光ファイバ温度レーダを用いる手法も提案されている。このセンサは、道路管理担当者が過去の知見に基づいて路面状態を推定するときの判断材料として用いられる。
【０００６】
【発明が解決しようとする課題】
赤外線または電磁波を照射する方式のセンサを用いて路面状態を道路の長手方向に沿って連続的に把握するためには、多数のセンサを設置する必要があり、コストがかかる。
【０００７】
一方、気象量と路面温度とを計測する場合は、気象量観測地がある程度の広域に対して適用可能であること、および光ファイバ温度レーダが長距離に敷設可能であることから低コストである。しかし、路面状態を判断する道路管理担当者に経験が求められること、及び危険と判断する地域を限定することが難しいという問題がある。
【０００８】
そこで、本発明の目的は、上記課題を解決し、低コストで、路面状態推定を自動化し、かつ局地的な判断が可能な路面状態推定方法を提供することにある。
【０００９】
【課題を解決するための手段】
上記目的を達成するために本発明は、路面下に埋設した光ファイバ温度レーダにより計測した地中温度計測値に基づき地中側から路面に流入する熱量（Ｗｅ）を求め、道路近傍に設置した気象計（気温計、湿度計、日射計、雪雨量計、風速計、放射収支計を含む）により計測した大気の温度、日射強度などの気象量計測値と、前記地中温度計測値から算出し、あるいは前記光ファイバ温度レーダなどの温度計測器で計測した路面温度（Ｔｅ）とに基づき大気側から路面に流入する熱量（Ｗａ）を求め、前記路面温度（Ｔｅ）と、前記気象量計測値の気温計測値、相対湿度計測値、風速計測値とに基づき路面上の水分の蒸発による蒸発熱（Ｗｅｖ）を求め、地中側から路面に流入する熱量（Ｗｅ）と大気側から路面に流入する熱量（Ｗａ）の差と、前記蒸発熱（Ｗｅｖ）とから湿潤判定値（Ｙ）を求めて、路面の湿潤状態を判定し、前記地中温度計測値及び前記気象量計測値と、前記湿潤判定値（Ｙ）とに基づき路面の状態を推定するものである。
【００１０】
路面への薬剤散布がない場合は、前記路面温度（Ｔｅ）とその時間当たりの変化量（ΔＴｅ）とから路面の積雪状態を判定し、路面への薬剤散布がある場合は、前記路面温度（Ｔｅ）と前記変化量（ΔＴｅ）と上記蒸発熱（Ｗｅｖ）と積雪の融解または凝固による潜熱量（δＷ）とから路面の積雪状態を判定してもよい。
【００１４】
【発明の実施の形態】
以下、本発明の一実施形態を添付図面に基づいて詳述する。
【００１５】
本発明の原理は、２つの異なる方法で求めた路面への流入熱量から路面状態を求めるものである。まず、２つの熱量計算方法について説明する。
【００１６】
第１の方法は、道路近傍に気象計（気温計、湿度計、日射計、雪雨量計、風速計、放射収支計を含む）を設置して、これらの気象量計測値から路面に流入する熱量を求めるものである。
【００１７】
ここで、路面の熱収支関係は、概念的には図１の内容で表される。路面の熱収支は、大気からの熱伝達、大気赤外放射と路面赤外放射、日射エネルギの吸収、及び道路内部熱伝達に分類することができる。本発明では、第１の方法を用いて、道路内部熱伝達を除いた道路の上側（大気側）から流入する熱量の総和Ｗａを計算する。実際には、濡れ路面の乾燥、積雪の融解等の路面上水分の相変化に伴う熱移動（潜熱伝達）も存在するが、この熱量を検出することが本発明の目的であるため、ここでは潜熱伝達は考慮しない。
【００１８】
Ｗａ＝［大気からの熱伝達］＋［大気赤外放射と路面赤外放射］＋［日射エネルギの吸収］
上記の各熱量を気象量計測値から求める計算式を以下に示す。
（１）大気からの熱伝達
大気からの熱伝達Ｑ１は、大気と路面との接触による熱移動であり、次式で求める。
【００１９】
Ｑ１＝ｈT ×（Ｔａ−Ｔｅ）×Ａ      （式１）
ただし、  ｈT ：熱伝達係数
Ｔａ：大気の温度
Ｔｅ：路面温度
Ａ：熱が伝わる面積
である。
（２）大気赤外放射と路面赤外放射
大気から路面に対する放射Ｑ２と、路面から大気に対する放射Ｑ３との計算は、それぞれ次の式で表される。
【００２０】
Ｑ２＝εａ×σ×Ｔａ4           （式２）
Ｑ３＝εｒ×σ×Ｔｅ4           （式３）
ただし、  εａ：大気の放射率
Ｔａ：大気の温度
εｒ：路面の放射率
Ｔｅ：路面温度
σ：ステファン・ボルツマン定数
である。
（３）日射エネルギの吸収
日射量と路面の日射吸収率とから、路面に吸収される日射エネルギＱ４を算出する。
【００２１】
Ｑ４＝（１−ｒ）×Ｉ            （式４）
ただし、  ｒ：路面のアルベト数
Ｉ：日射強度
である。
【００２２】
よって、
Ｗａ＝Ｑ１＋Ｑ２＋Ｑ３＋Ｑ４
により、大気側から路面に流入する熱量の総和Ｗａを求めることができる。
【００２３】
第２の方法は、地中に埋設した光ファイバ温度レーダによる地中温度計測値に基づき、地中側から路面に流入する熱量（第１の方法で除外した道路内部熱伝達による流入熱量）を求めるものである。
【００２４】
大気側から路面への流入熱量ΔＷがΔＷ* の時、道路に埋設した光ファイバで計測した地中温度Ｔf がΔｔ時間当たりΔＴf だけ変化したとする。このとき、ΔＴf は、ΔＷ＝０の時の光ファイバ埋設位置ｆでの温度変化量δＴf0と、大気側から路面へ熱量ΔＷ* が流入してきたときの温度変化量δＴf との和で表すことができる。
【００２５】
ΔＴf ＝δＴf0＋δＴf           （式５）
式５において、ΔＷ＝０の時の光ファイバ埋設位置ｆでの温度変化量δＴf0は、次の方法により求める。地中から路面への流入熱Ｗ2 は、地中内部の温度分布から求まる。ここでは、地中から路面への流入熱Ｗ2 を基底温度Ｔ0 から求める方式で説明する。
【００２６】
基底温度Ｔ0 を利用する方式は、地中深いところ（数１０ｍ程度の深さ）の温度が年間を通してほぼ一定になるということを利用するものであり、この温度を基底温度Ｔ0 とすると、路面から地中深いところまでの熱等価回路は図２で表すことができる。即ち、地表の測定点から所定の深さごとに地中温度の測定点があるとき、各測定点間には熱抵抗が存在し、各測定点には熱容量が存在し、地表の測定点には大気側から路面への流入熱Ｗ1 が入力される。
【００２７】
熱量計算は、計算刻み時間Δｔごとに行う。地中を図２に示すような階層構造と見なす場合、地表を１番目としたときｋ番目の階層における地中温度Ｔk は式６であらわされる。
【００２８】
0003809519.tif 000002
ここで、    Ｒ＝Δｈ／λT
Ｃ＝ρ・Ｃp ・ΔｈD
ただし、
Ｒ：当該階層の熱抵抗
Ｃ：当該階層の熱容量
ΔｈD ：道路下深さ方向の刻み幅
λT ：当該階層における熱伝導率
ρ：当該階層における地中物質の密度
Ｃp ：当該階層における面積当たり熱容量
k ：階層を表す添字変数（1 ，2 ，3 ，…）
である。
【００２９】
ここで、境界条件は式７、式８で表される。
【００３０】
（Ｔ1 −Ｔ2 ）／Ｒ1 ＝Ｗ1           （式７）
Ｔn ＝Ｔ0                           （式８）
式５から式８までの式を解くことで、光ファイバ埋設位置ｆにおける温度変化量δＴf0を求めることができる。
【００３１】
次に、図３に示すように、路面から光ファイバ埋設位置ｆへ流入する熱量と、それに伴う光ファイバ埋設位置ｆの温度変化量が直線で近似できる範囲内で路面に対する流入熱の変化量ΔＷを任意の値ΔＷz として与える。このとき、光ファイバ埋設位置ｆでの仮の温度変化量ΔＴfzは、式７のＷ1 にΔＷz を代入して式６から式８までの式を解くことで求めることができる。温度変化量ΔＴfzは、ΔＷ＝０の時の光ファイバ埋設位置ｆでの温度変化量δＴf0と、任意に定めた路面への流入熱量ΔＷz によって与えられる温度変化量δＴfzとの和で表すことができる。
【００３２】
ΔＴfz＝δＴf0＋δＴfz            （式９）
よって、Ｗ＝０の時の光ファイバ埋設位置ｆでの温度変化量δＴf0と仮の温度変化量ΔＴfzとが上記までに求められるので、温度変化量δＴfzが求まる。
【００３３】
また、図３に示す直線関係を仮定しているので、式１０の関係が成り立つ。
【００３４】
δＴf ／ΔＷ* ＝δＴfz／ΔＷ      （式１０）
式９、式１０から、路面に流入する熱量ΔＷ* を求める式１１が得られる。
【００３５】
ΔＷ* ＝（ΔＷz ／δＴfz）・δＴf         （式１１）
式１１において、δＴf は実際の光ファイバ埋設位置ｆでの温度変化量δＴfzは、式９から得られ、また、ΔＷz は任意に定めることのできる値であることから、実際の熱量の変化量ΔＷ* を算出することができる。このようにして求めた熱量ΔＷ* が地中側から路面への流入熱Ｗｅである。
【００３６】
上記の計算に必要となる地中の熱定数は、地中の物質の組成によって一定なので、事前に熱流計を用いて計測しておけば、地中側から路面への流入熱Ｗｅを正確に算出することができる。
【００３７】
次に、第１の方法で求めた熱量Ｗａと第２の方法で求めた熱量Ｗｅとを用いて路面状態を推定する手法を以下に述べる。
【００３８】
まず、第１の方法で求めた熱量Ｗａ及び熱量Ｗｅを用いて、路面の湿潤状態を判定する方法を述べる。路面上の水分の厚みをｈW （ｔ）とすると、水分の厚みｈW （ｔ）と蒸発熱Ｗｅｖ（ｔ）との関係は式１２で表すことができる。
【００３９】
ｈW （ｔ）＝ｈW ０−Ｗｅｖ（ｔ）／λW     （式１２）
ここで、Ｗｅｖ（ｔ）は式１３により気象量より求められる値である。
【００４０】
0003809519.tif 000003
ただし、
ｈW ０：降雨直後の路面に溜った水分の厚み
λW ：水の蒸発潜熱
α，β：定数
ｖ：風速
Ｈ：相対湿度
Ｔｅ：路面温度
Ｔａ：気温
Ｐ（Ｔｅ）：温度Ｔｅにおける水分の飽和水蒸気圧
Ｐ（Ｔａ）：温度Ｔａにおける水分の飽和水蒸気圧
である。
【００４１】
このように、本発明では、路面上の水分の蒸発による蒸発熱Ｗｅｖを路面水分厚ｈW ０、風速ｖ、路面温度Ｔｅ、相対湿度Ｈ、気温Ｔａの測定値（一部推定値も含む）と既知の諸量とから求める。
【００４２】
湿潤状態では、大気側から路面への流入熱Ｗａ、大地（地中側）から路面への流入熱Ｗｅ、蒸発熱Ｗｅｖの間には、以下の関係式が成立する。
【００４３】
Ｗａ−Ｗｅ＝Ｗｅｖ              （式１４）
ここで、湿潤判定値Ｙとして式１５を導入する。
【００４４】
Ｙ＝（Ｗａ−Ｗｅ）／Ｗｅｖ      （式１５）
湿潤判定値Ｙは、理論的に、湿潤状態では１、乾燥状態及び凍結状態では０となるので、この値をモニタしていると湿潤の有無が判定できる。即ち、本発明では、大気側から路面への流入熱Ｗａ、地中側から路面への流入熱Ｗｅ、路面上の水分の蒸発による蒸発熱Ｗｅｖを用いて湿潤判定値Ｙを求め、この湿潤判定値Ｙにより路面の湿潤状態を判定する。
【００４５】
次に積雪状態の判定方法を述べる。積雪状態の判定条件は、路面への薬剤散布がある場合とない場合とで異なる。路面への薬剤散布がない場合は、次の判定条件によって容易に積雪状態判定が可能である。
【００４６】
Ｔｅ＝０℃  かつ  ΔＴｅ＝０℃
ただし、
ΔＴｅ：計測時間当たりの路面温度の変化量
である。
【００４７】
このように、路面への薬剤散布がない場合は、路面温度Ｔｅを用い、路面温度Ｔｅが時間経過によらず０℃を維持していることで積雪状態が判定できる。
【００４８】
薬剤散布が有り得る場合は、まず、積雪判定値Ｚを式１６で定義する。
【００４９】
Ｚ＝ΔＴｅ／ΔＴｃ              （式１６）
ただし、
ΔＴｅ：計測時間当たりの路面温度の上昇値（測定値）
ΔＴｃ：大気側からの流入熱Ｗａが流入したときの路面温度の上昇値（計算値）
である。路面温度の変化量であるΔＴｅ及びΔＴｃは、路面への流入熱Ｗｅ及びＷａにそれぞれほぼ比例する。これらの比例関係（ΔＴｅがＷｅに比例、ΔＴｃがＷａに比例）によって、ほぼΔＴｅ／ΔＴｃ＝Ｗｅ／Ｗａとなる。従って、
Ｚ＝Ｗｅ／Ｗａ                  （式１６´）
と表すことができる。
【００５０】
また、湿潤状態の判定に用いた湿潤判定値Ｙについて考えると、積雪の融解または凝固に使われる潜熱量をδＷとすれば、
0003809519.tif 000004
となり、湿潤判定値Ｙは、積雪状態では
Ｙ＝１以上    融解過程（δＷが負）のとき
Ｙ＝１以下    凝固過程（δＷが正）のとき
となる。
【００５１】
次に、大気側から路面への流入熱Ｗａと地中側から路面への流入熱Ｗｅとの差をΔＷとすると、
0003809519.tif 000005
となる。ここで、δＷは、融解過程において正、凝固過程において負となる。
【００５２】
式１６´、式１８から、積雪判定値Ｚは、次のように表すことができる。
【００５３】
0003809519.tif 000006
式１９から以下の条件によって、積雪状態判定を行うことができる。
【００５４】
0003809519.tif 000007
の条件が考えられる。結局、積雪が継続する条件は、
条件１＝Ｗｅ＜０かつＴｆ＜ＴｅかつδＷ＞０かつＷｅｖ＞０かつＺ≧１、
条件２＝Ｗｅ＜０かつＴｆ＞ＴｅかつδＷ＞０かつＷｅｖ＜０かつ｜δＷ｜＞｜Ｗｅｖ｜かつ（Ｗｅｖ＋δＷ）／Ｗｅ＞０かつＺ＜１、
条件３＝Ｗｅ＞０かつＴｆ＜ＴｅかつδＷ＞０かつＷｅｖ＞０かつＺ＜１、
条件４＝Ｗｅ＞０かつＴｆ＞ＴｅかつδＷ＜０かつＷｅｖ＞０かつ｜δＷ｜＞｜Ｗｅｖ｜かつ（Ｗｅｖ＋δＷ）／Ｗｅ＜０かつＺ＞１、
の４つの条件であり、条件１〜４のいずれかが満たされている場合は、積雪状態と判定する。
【００５５】
このように、本発明では、路面への薬剤散布がない場合は、路面温度Ｔｅを用い、路面温度Ｔｅが時間経過によらず０℃を維持しているかどうか積雪状態を判定し、路面への薬剤散布がある場合は、路面温度の変化方向（熱量Ｗｅの正負）と、水の凝固点Ｔｆに対する路面温度Ｔｅの高低と、蒸発熱Ｗｅｖの正負と、積雪の融解または凝固による潜熱量δＷの正負とから路面の積雪状態を判定する。
【００５６】
次に、凍結状態の判定方法を述べる。熱的には凍結状態と乾燥状態とでは、ほぼ同じ挙動が生じる。即ち、凍結状態でも乾燥状態でも、積雪判定値Ｚ＝ほぼ１、かつ湿潤判定値Ｙ＝ほぼ０となる。凍結状態と乾燥状態とを区別するためには、現時点に至るまでの過去の路面状態の時系列を用いる。この時系列で前回の状態を基準にする。即ち、積雪判定値Ｚ＝ほぼ１、かつ湿潤判定値Ｙ＝ほぼ０であるとき、前回の状態が湿潤状態であれば、今回は乾燥状態に移行すると考える。また、積雪判定値Ｚ＝ほぼ１、かつ湿潤判定値Ｙ＝ほぼ０であるとき、路面温度Ｔｅが前回、前々回の路面温度Ｔｅとほぼ同じであれば遷移状態と考え、さらに路面への流入熱ΔＷが負の値であれば冷却過程であると考え、冷却過程の遷移状態であれば、次回は凍結状態に移行すると考える。
【００５７】
様々の路面状態において計算した判定値Ｙ，Ｚをそれぞれ横軸、縦軸にとりプロットすると、各路面状態に対応した領域を図４のように示すことができる。この図に即して各路面状態の特長を次に述べる。
１）乾燥状態、凍結状態（定常状態）
Ｙ＝０，Ｚ＝１のポイントで示される。乾燥状態か凍結状態かの区別はトレンド（過去の路面状態の時系列）から判定することになる。
２）湿潤状態
Ｙ＝１，Ｚ＝０〜１（温度上昇過程のとき）の範囲か、Ｙ＝１，Ｚ＝１以上  （温度下降過程のとき）の範囲で示される。
３）積雪状態（薬剤なし）、凍結状態（遷移状態であって薬剤なし）
Ｚ＝０，Ｙ＝１以下（凝固過程のとき）の範囲か、Ｚ＝０，Ｙ＝１以上（融解過程のとき）の範囲で示される。
４）積雪状態（薬剤あり）、凍結状態（遷移状態であって薬剤あり）
薬剤を含んだ水の凝固点Ｔｆと路面温度Ｔｅとの組み合わせによって、Ａ〜Ｄの４つの領域に区分される。
【００５８】
この４つの領域は、路面温度Ｔｅの時間的変化に対して図５のように区分される。即ち、路面温度Ｔｅが正の温度から負の温度に降下した後、凝固点Ｔｆに至るまでは領域Ａ、路面温度Ｔｅが凝固点Ｔｆより下で降下すると領域Ｂ、路面温度Ｔｅが凝固点Ｔｆより下で上昇すると領域Ｃ、路面温度Ｔｅが凝固点Ｔｆより上で上昇し０度に至るまでは領域Ｄとなる。
【００５９】
これまでに述べた路面状態の判定方法を実際に行った結果を図６に示す。対象期間は３日間で、その間に路面状態を目視観測した結果を図の上部に示してある。縦線縞模様を描いた期間は湿潤状態であり、その後の空白を描いた期間は乾燥状態であった。これに対して、潜熱相当の熱量比である湿潤判定値Ｙは、実線グラフで示したとおりに変化した。なお、図中ＷｗはＷｅｖと同義である。この湿潤判定値Ｙに適切な閾値（ここでは０．１程度）を適用することで湿潤状態を判定することが可能である。
【００６０】
本発明の路面状態推定方法を実施するためのシステム構成について図７を用いて説明する。このシステムは、対象となる道路４の近傍に設置された気象センサ１、路面下に埋設された光ファイバ２（光信号を温度分布値に変換する光ファイバ温度レーダを含む）、光ファイバ２により計測されたデータを信号線６を介して伝送する信号伝送装置３、気象センサ１により計測されたデータを信号線６を介して伝送する信号伝送装置５、計測データを受信する情報収集装置７、受信したデータに基づいて路面状態推定計算を行う路面状態推定装置８、その推定結果を出力する表示部９により構成される。
【００６１】
気象センサ１は、電気的な出力が得られるものであれば、気象諸量の計測方式は特に限定されない。気象センサ１は、気温計、湿度計、日射計、雪雨量計、風速計、放射収支計などを総称したものであり、大気の温度Ｔａ、相対湿度Ｈ、日射強度Ｉ、風速ｖなどを計測することができる。
【００６２】
図示しないが路面上の塩分濃度を測定する塩分量計が設けられている。塩分量計で測定される塩分濃度によって、路面上の水分の凝固点Ｔｆが計算される。また、塩分濃度によって、薬剤散布の有無を判定することができる。
【００６３】
また、水膜厚ｈW （ｔ）については、予めレーザレーダ等により路面粗さ（路面凹凸深さの平均値）を計測しておき、降雨時の雪雨量計の計測値が路面粗さ以下であれば雪雨量計の計測値を初期水膜厚ｈW ０とし、降雨時の雪雨量計の計測値が路面粗さを超えるときには、路面粗さを初期水膜厚ｈW ０とし、ある時間が経過した後の水膜厚ｈW （ｔ）は式１２で求める。
【００６４】
路面温度Ｔｅは、前述した第２の方法で求めた路面への流入熱Ｗｅを用いて路面温度の変化量ΔＴｅが算出可能であるので、前回の路面温度ＴｅにΔＴｅを加算することで今回の路面温度Ｔｅの値が求まる。ただし、システム起動時に限り、路面温度Ｔｅの初期値は路面温度が直接計測可能な温度計測器を使用して実測する。
【００６５】
光ファイバ２の埋設位置は、図２の温度Ｔ１，Ｔ２，Ｔ３…の測定点のいずれかであり、実際の深さ位置は施工上の都合で決定することができる。
【００６６】
信号伝送装置３、５は、無線伝送装置であってもよい。信号伝送装置３、５が無線伝送装置である場合の構成を図８に示す。図７と符号の同じ部材は同一部材である。図７との相違として、信号線６の代わりに受信アンテナ６が設置される。
【００６７】
図７、図８のいずれの場合も、情報収集装置７は、気象センサ１及び光ファイバ２からの情報を受信し、データの種類や遅延時間等を考慮して必要なデータを路面状態推定装置８に与える。
【００６８】
路面状態推定装置８は、これら入力データを用い、これまでに説明した計算式、手順に従い路面状態を判定する。
【００６９】
【発明の効果】
本発明は次の如き優れた効果を発揮する。
【００７０】
（１）長距離、広範囲での路面状態推定が従来より低コストで可能となる。
【図面の簡単な説明】
【図１】本発明の基本となる路面の熱収支の概念図である。
【図２】本発明において地中側から路面に流入する熱量を求めるための地中熱伝導の等価回路図である。
【図３】本発明で使用する地中側から路面に流入する熱量と地中温度変化量との関係図である。
【図４】本発明による判定値Ｙ，Ｚと路面状態との関係図である。
【図５】本発明による路面温度変化と凝固点温度とによる路面状態判定区分を示した路面温度の時間変化図である。
【図６】本発明で実際に路面状態を判定したときの、熱量比及び路面状態の時間変化図である。
【図７】本発明の路面状態推定方法を実施するためのシステム構成の一例を示すブロック図である。
【図８】本発明の路面状態推定方法を実施するためのシステム構成の一例を示すブロック図である。
【符号の説明】
１  気象センサ（気象計）
２  光ファイバ（光ファイバ温度レーダ）
４  道路
８  路面状態推定装置
【図１】
0003809519.tif 000008
【図２】
0003809519.tif 000009
【図３】
0003809519.tif 000010
【図４】
0003809519.tif 000011
【図５】
0003809519.tif 000012
【図６】
0003809519.tif 000013
【図７】
0003809519.tif 000014
【図８】
0003809519.tif 000015