"時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するための非定常発塵源の探索方法であって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出するステップと、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出するステップと、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出するステップと、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定するステップと、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定するステップとを有することを特徴とする降下煤塵の非定常発塵源の探索方法。連続する2つ以上のTdit期間から構成されるTgk期間において、前記評価地点ixにおいて最大の降下煤塵量mを測定したTdit期間の風向平均値及び風速平均値を前記代表風向及び代表風速として用いることを特徴とする請求項1に記載の降下煤塵の非定常発塵源の探索方法。前記参照地点iyの総数をniyと定義し、前記評価関数として、下式を用いることを特徴とする請求項1又は2に記載の降下煤塵の非定常発塵源の探索方法。前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類するための方法として、捕集された降下煤塵粒子の放射線量を測定してその強度に基づいて降下煤塵粒子を分類することを特徴とする請求項1乃至3のいずれか1項に記載の降下煤塵の非定常発塵源の探索方法。前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類するための工程として、高炉法による製鉄プラントに由来する降下煤塵の煤塵種の判別を前提とし、捕集された降下煤塵試料に磁力を付与して降下煤塵を磁石に着磁する着磁性降下煤塵と磁石に付着しない非着磁性降下煤塵に分離する第1の工程と、前記着磁性降下煤塵と前記非着磁性降下煤塵それぞれに対して、降下煤塵をカメラで撮影して着磁性降下煤塵画像及び非着磁性降下煤塵画像を得る第2の工程と、前記着磁性降下煤塵画像及び前記非着磁性降下煤塵画像をそれぞれ画像処理して個々の降下煤塵粒子の代表明度を算出し、前記代表明度と所定の明度しきい値を比較して個々の粒子を暗色粒子と明色度粒子に粒子明度を区分し、さらに、前記粒子明度区分及び前記第1の工程で得られた前記個々の降下煤塵の着磁性の有無との組み合わせを用いて前記個々の降下煤塵粒子を着磁性暗色粒子、着磁性明色粒子、非着磁性暗色粒子、及び非着磁性明色粒子とに煤塵特性を区分する第3の工程と、前記製鉄プラントにおける代表的降下煤塵種毎の純粋な試料に対して前記第1~第3の工程を施して得られた煤塵特性の区分と個々の煤塵種の対応関係を求め、前記捕集された降下煤塵の前記第3の工程で得られた煤塵特性と前記純粋な煤塵種の煤塵特性を比較して最も煤塵種特性の類似する煤塵種を求め、前記捕集された降下煤塵の煤塵種を前記最も煤塵種特性の類似する煤塵種と判別する第4の工程とを有することを特徴とする請求項1乃至4のいずれか1項に記載の降下煤塵の非定常発塵源の探索方法。時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するための非定常発塵源の探索装置であって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出する手段と、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出する手段と、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出する手段と、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定する手段と、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定する手段とを備えたことを特徴とする降下煤塵の非定常発塵源の探索装置。時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するためのプログラムであって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出する手段と、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出する手段と、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出する手段と、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定する手段と、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定する手段としてコンピュータを機能させるためのプログラム。従来技術においては、以下の問題点があった。即ち、第1の問題点として、発生源を探索する対象の発生物が降下煤塵ではないことが挙げられる。例えば、特許文献1~4の技術においては、対象がガスであり、特許文献3の技術には、対象にSPMが含められているに過ぎない。SPMは、降下煤塵に比べて遥かに小さな粒子であり定義上、直径10μm以下の粒子、その大気中での拡散挙動は、微小な粒子沈降を生じることを除けば実質的にガスの挙動に等しい。一方、降下煤塵は、SPMに比べて遥かに大きな煤塵粒子であり直径約10μm超の粒子である、その落下速度が極めて大きいため、降下煤塵の大気中での拡散には、粒子降下速度の極めて大きな影響を受け、拡散挙動はガスとは大きく異なる。また、本発明が対象とする降下煤塵量とは、地表面への降下煤塵の沈着量のことである。特許文献1~4の技術では、降下煤塵評価地点におけるガス及びSPM濃度を観測及び管理対象としており、地表面への沈着速度を直接知ることはできない。確かに、数式2には、沈着速度Vdが記載されているので、Vdを正確に与えることができれば、降下煤塵評価地点上でのガス及びSPM濃度から地表面での沈着量に換算することが可能である。しかし、非特許文献1に記載されているように、SPMに対してVdは地表面の状態や大気乱流の影響を受けて大きく変動する。また、ガスに関しては、沈着速度を一般的に与える手法は開発されていない。従って、Vdの値を正確に与えることは実際には極めて困難であり、特許文献1~4の技術で降下煤塵を対象とすることは、少なくとも定量的には困難である。第2の問題点として、降下煤塵を対象とした発塵源探索手法は、従来、存在しなかった。これは、従来の発生源探索手法においては、特許文献3に代表されるように水平面内での発生源探索を前提としており、大きな粒子落下速度Vsを有し、かつ、地表面での沈着量を問題とする、降下煤塵の発生源を三次元的に取り扱うことが困難であったからである。特に、特許文献3に示されるような、降下煤塵評価地点から風上方向に発生源探索線を伸長する手法の場合、数式2における地表面反射項αexp[-Hez-Vsx/WS2/2σz2]の影響を定量的、かつ、一般的に取り扱うことが困難なことから、前記発生源探索線をプルーム式と関連付ける有効な手法は従来、提案されていない。第3の問題点として、従来技術においては、発生源の探索を行う際に発生源の位置及びそこでの概略の発生量を予め仮定する手順が必須であることが挙げられる。例えば、特許文献1及び2の技術においては、予め想定される全ての発生源及び全ての降下煤塵評価地点について、任意発生源での発生量と任意降下煤塵評価地点での濃度の関係を上記プルーム式等の気象条件の関数として予測し、次に、全ての降下煤塵評価地点における濃度実測値と、前記濃度予測値との差が最小となるように、前記関数のパラメータσyやQP等を最適化手法により調整している。従って、少なくとも全ての発生源位置を予め与える必要があり、また、最適化手法の計算過程の妥当性を確保するためには、各発生源での概略の発生量も初期条件として予め与えることが一般には望ましい。なぜならば、最適化問題においては、実情から極端に解離した初期条件を与えた場合、実情とは大きく異なる局所安定点に解の収束する場合があるからである。特許文献3の技術においては、図6に示すように、粉塵SPM発生箇所a、b、c、d、e等を予め仮定した上で、周辺の複数降下煤塵評価地点ix:2,ix:3,ix:4等でのSPM濃度を長期間、測定し、この期間内で各評価点で風向別SPM濃度平均値1の最も高い風向の風上方向に評価点ix:2,ix:3,ix:4からそれぞれ水平面内で発生源探索線3、4、5を伸長し、これらix:2,ix:3,ix:4発生源探索線の互いに交差した交点6、7、8の内、前記a、b、c、d、eのいずれかに合致した地点を、発生箇所を特に発生量の大きい発生箇所と判定している。特許文献4の技術においては、想定される発生源の近傍に計測機を設けることが前提であるので、発生源は、予め既知でなければならない。しかしながら、多数の発生源が存在する場合、これら全ての発生源位置と概略発生量を予め全て把握することは、実際には困難であり、もし、可能だとしても、多大な資源を必要とするため好適ではない。また、原子力発電所の事故地のように、そもそも発塵源に近寄ることのできない場合もある。従って、従来技術は、発生源の数が極めて少数であるか、あるいは、発生源の発生量を十分正確に把握し得る環境下でしか有効に適用することはできない問題がある。第4の問題点として、従来技術においては対象とする発生源は、発生量が時間的に変動しない定常発生源であるか、又は、発生量が時間平均値の近傍で僅かに時間変動するだけの準定常発塵源を基本的に対象としている。例えば、特許文献1及び2においては、この技術では最適化手法を適用するため、一般的には、降下煤塵評価地点数の数を、適用されるプルーム式等の関数の中で調整可能なパラメータの数よりも、多く設定しなければならない。もし、調整可能なパラメータの数が実質的に降下煤塵評価地点数よりも多ければ、得られる解は、一般に一意に定まらないので手法として破綻するからである。多数の発生源が存在する場合、経済性の観点から降下煤塵評価地点の数を発生源の数よりも少なく設定する場合が多い。このような場合でも、発生源を定常発生源に限定すれば即ち発生量QPを調整可能なパラメータとはしなければ、多数の異なる時刻での降下煤塵評価地点での測定値を用いることにより発生源数以上の測定データ数を確保することができ、最適化手法を適用することができる。一方、発生量QPが非定常的に大きく変動する、非定常発生源に対して特許文献1、2の技術を適用する際には、QPを調整可能なパラメータとせざるを得ないため、多数の発生源対象とする場合には、発生源数を超える極めて多数の降下煤塵評価地点を設ける必要があり、経済性の観点から現実的でない。特許文献3の技術においては、2カ月以上の期間内での離散的に採取された降下煤塵評価地点濃度データを平均化して発生源の探索を行う。従って、発生源は、定常発生源に限定される。特許文献4の技術においては、想定される発生源の近傍に降下煤塵評価地点を配置するので原理的には非定常発生源を探索することができる。しかし、この技術においては、複数の発生源からのガスが特定降下煤塵評価地点に同時に到達する場合にどの発生源が卓越するのかを判断する方法が開示されておらず、また、想定される全ての発生源の近傍に降下煤塵評価地点を設置するのかも記載されていない。従って、この技術で非定常発塵源を探索することが可能なのは、発生源間の距離が互いに影響を及ぼさない程度に遠い場合に限られる。即ち、この技術は、実質的に発生源と降下煤塵評価地点が一対一に特定できる対象にしか適用できない。しかしながら、現実の発生源では発生量が一般に大きく時間変動するので、定常発生源や、発生源と降下煤塵評価地点が一対一に特定できる発生源のみを対象とする従来技術では、現実の発生源探索に対して十分に適用できない問題があった。以上のように、捕集された降下煤塵試料を判別するための簡易な分析手法は、従来存在しなかった。本発明は上記のような点に鑑みてなされたものであり、発塵量が非定常的に変動する降下煤塵の発塵源を、発塵源周辺での降下煤塵測定値を基に探索する手法を提供することを目的とする。上記課題を解決するために、本発明者の研究の結果、以下の解決方法を発明するに至った。1時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するための非定常発塵源の探索方法であって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出するステップと、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出するステップと、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出するステップと、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定するステップと、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定するステップとを有することを特徴とする降下煤塵の非定常発塵源の探索方法。2連続する2つ以上のTdit期間から構成されるTgk期間において、前記評価地点ixにおいて最大の降下煤塵量mを測定したTdit期間の風向平均値及び風速平均値を前記代表風向及び代表風速として用いることを特徴とする1に記載の降下煤塵の非定常発塵源の探索方法。3前記参照地点iyの総数をniyと定義し、前記評価関数として、下式を用いることを特徴とする1又は2に記載の降下煤塵の非定常発塵源の探索方法。4前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類するための方法として、捕集された降下煤塵粒子の放射線量を測定してその強度に基づいて降下煤塵粒子を分類することを特徴とする1乃至3のいずれか一に記載の降下煤塵の非定常発塵源の探索方法。5前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類するための工程として、高炉法による製鉄プラントに由来する降下煤塵の煤塵種の判別を前提とし、捕集された降下煤塵試料に磁力を付与して降下煤塵を磁石に着磁する着磁性降下煤塵と磁石に付着しない非着磁性降下煤塵に分離する第1の工程と、前記着磁性降下煤塵と前記非着磁性降下煤塵それぞれに対して、降下煤塵をカメラで撮影して着磁性降下煤塵画像及び非着磁性降下煤塵画像を得る第2の工程と、前記着磁性降下煤塵画像及び前記非着磁性降下煤塵画像をそれぞれ画像処理して個々の降下煤塵粒子の代表明度を算出し、前記代表明度と所定の明度しきい値を比較して個々の粒子を暗色粒子と明色度粒子に粒子明度を区分し、さらに、前記粒子明度区分及び前記第1の工程で得られた前記個々の降下煤塵の着磁性の有無との組み合わせを用いて前記個々の降下煤塵粒子を着磁性暗色粒子、着磁性明色粒子、非着磁性暗色粒子、及び非着磁性明色粒子とに煤塵特性を区分する第3の工程と、前記製鉄プラントにおける代表的降下煤塵種毎の純粋な試料に対して前記第1~第3の工程を施して得られた煤塵特性の区分と個々の煤塵種の対応関係を求め、前記捕集された降下煤塵の前記第3の工程で得られた煤塵特性と前記純粋な煤塵種の煤塵特性を比較して最も煤塵種特性の類似する煤塵種を求め、前記捕集された降下煤塵の煤塵種を前記最も煤塵種特性の類似する煤塵種と判別する第4の工程とを有することを特徴とする1乃至4のいずれか一に記載の降下煤塵の非定常発塵源の探索方法。6時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するための非定常発塵源の探索装置であって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出する手段と、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出する手段と、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出する手段と、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定する手段と、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定する手段とを備えたことを特徴とする降下煤塵の非定常発塵源の探索装置。7時間周期Δtd毎のit番目の時刻であるTditにおいて、Tdit-1からTditまでの期間であるTdit期間に任意の降下煤塵観測地点iで降下煤塵を捕集して算出される降下煤塵量mi,itと、前記各降下煤塵観測地点iの近傍において、前記時間周期Δtdよりも長くない時間周期Δtwdにて連続的に計測された風向及び風速測定値を基に算出される代表風向WD及び代表風速WSと、1つ以上の連続する前記時刻Tditを含む時間周期Δtg毎の、k番目の時刻tgkを設けてtgk-1からtgkまでの期間であるTgk期間を設定し、任意の時刻tgkにおいて、Tgk期間に捕集された降下煤塵を他から区別して降下煤塵試料としてそれぞれ保存し、特定のTgk期間において、前記特定のTgk期間に捕集された降下煤塵試料を分析して算出される個々の降下煤塵粒子の粒子落下速度Vsを用いて個々の降下煤塵粒子を2つ以上の粒子落下速度区分ivのいずれかに分類すると共に、前記特定のTgk期間に捕集された降下煤塵試料を分析して個々の降下煤塵粒子を2つ以上の煤塵種区分isのいずれかに分類し、前記粒子落下速度区分iv及び前記煤塵種区分isを用いて算出される、前記特定のTgk期間に特定の降下煤塵観測点iで捕集された降下煤塵試料の区分毎の算出された質量又は体積構成率Comi,iv,isと、を用いて、発塵量が非定常に変動する降下煤塵発生源の位置を特定するためのプログラムであって、1つ以上の降下煤塵観測地点である評価地点ixを設けて、前記各評価地点ixにおいて降下煤塵量mix,itを測定してTgk期間における代表降下煤塵量Mixを算出する手段と、3つ以上の降下煤塵観測地点である参照地点iyを設けて、前記各参照地点iyにおいてそれぞれ降下煤塵miy,itを測定してTgk期間における代表降下煤塵量Miyを算出する手段と、2つ以上の、発塵源として仮定した仮定発塵源jを設け、前記仮定発塵源jに関して、単位発塵量を前提としてプルーム式を用いて前記評価地点ixにおける降下煤塵量計算値mcj,ix及び前記参照地点iyにおける降下煤塵量計算値mcj,iyを算出する手段と、前記特定のTgk期間で前記評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義し、Comiy,ivmaxix,ismaxixMiy/CommaxixMixで表現される規格化実績降下煤塵量Raiy及びmcj,iy/mcj,ixで表現される規格化計算降下煤塵量Rciyの比率Rciy/Raiyを用いた評価関数を設定する手段と、全ての前記仮定発塵源jについて、前記規格化計算降下煤塵量Rciyの前記規格化実績降下煤塵量Raiyに対する一致性を前記評価関数を用いて評価し、最も一致性の高い仮定発塵源jを前記特定のTgk期間における非定常発塵源と判定する手段としてコンピュータを機能させるためのプログラム。本発明によれば、少数の降下煤塵評価地点での降下煤塵計測によって、非定常発塵源を含む発塵源の探索を効率的、かつ、正確に実施することが可能になる。以下、添付図面を参照して、本発明の好適な実施形態について説明する。尚、本明細書及び図面において、実質的に同一の機能構成を有する構成要素については、同一の符号を付することにより重複説明を省略する。第1の実施形態降下煤塵量計測手段によって周期Δtd毎に降下煤塵量計測値の出力される時刻Tditにおいて、Tdit-1からTditまでの時間をTdit期間と定義する。itは、降下煤塵計測開始時刻を0とし、1ずつ増加する整数である。このTdit及びΔtdは、複数の降下煤塵の連続測定手段においてそれぞれ存在するが、全ての降下煤塵の連続測定手段において、Tdit及びΔtdを揃えることが、解析の便宜上、好ましい。このように揃えられたTdit期間において、各降下煤塵の連続測定手段で測定された降下煤塵量測定値をそれぞれmi,itiは、測定手段の通し番号に対応すると定義する。また、1以上の自然数であるnt個の、連続するTdit期間から構成される時間をTgk期間と定義する。ここで、Tgk期間の始点の時刻をTgk-1とし、このときのitを0とする。Tgk期間の終点の時刻をTgkとし、このときのitをntとする。kは、降下煤塵計測開始時刻を0とし、1ずつ増加する整数である。本実施形態は、個々のTgk期間における降下煤塵の発生源を特定するものであり、Δtg=ntΔtd以上の時間スケール即ち、発塵継続時間を有する発塵源を対象とする。降下煤塵量の測定及び降下煤塵試料の捕集のためには、例えば、特許文献6に記載される連続式降下煤塵計を用いることができる。この装置を用いる場合、例えば、時間周期Δtdを10分に、Δtgを60分に設定することができる。降下煤塵は、後に降下煤塵の分類を分析するために、捕集を行い、Tgk期間毎に、他の時刻の試料と分離してそれぞれ保管することができる。このように、降下煤塵の捕集時間単位をΔtdではなく、Δtgとする理由は、一般に、連続式降下煤塵測定手段において、降下煤塵量測定値の測定期間のΔtd毎の切り替えにはさほど時間を必要としない数秒程度でよいのに対し、Δtg毎の降下煤塵の捕集単位を切り替える際には、捕集フィルタを交換する等の物理的操作が必要なために通常数分単位の時間を必要とし、この交換時間中は発塵源探索を行うための測定値が得られないからである。本実施形態で対象とする非定常発塵源の時間スケールは、高々10分のオーダであるので、頻繁な捕集単位の切り替えは、解析可能な時間割合を著しく減少させるので好ましくない。このため、比較的頻繁に切り替え可能な降下煤塵量測定の周期Δtdを、捕集時間周期Δtgよりも短く設定することにより、時間スケールΔtg以上の非定常発塵源を効率的に特定することができる。また、降下煤塵の捕集を行ってTgk期間毎に、他の時刻の試料と分離してそれぞれ保管する方法としては、例えば、市販のテープ式煤塵捕集フィルタを用いて、前記捕集フィルタの特定個所にTgk期間の間に煤塵の捕集を行い、時刻Tgkにおいて、テープ状のフィルタを移動させて、次のTgk1期間用に、捕集フィルタ上の新たな箇所で煤塵捕集を行うと共に、Tgk期間に捕集された煤塵の付着した捕集フィルタ上を保護テープ等で覆うことにより、この箇所に捕集された煤塵の飛散や脱落を防止することによって実現することができる。一例として、図2の平面図例えば、地上高度15mの水平面の配置を対象とした実施形態に対して、特定の時刻TgkにおけるTgk期間における特定の降下煤塵評価地点ixにおける主要な発塵源jの特定手順として、本実施形態のフローの具体例を図1のフローチャートに沿って示す。図1のフローチャートは、非定常発塵源の探索装置以下、発塵源探索装置と称するによる発塵源探索処理の一例である。発塵源探索装置は、例えば図3に示すように、CPU等の演算装置301、メモリ302、HDD303、及び各種のインターフェース304を備えた情報処理装置例えば、市販のパーソナルコンピュータPCを用いることにより実現される。例えば、図1のフローチャートは、C言語等のプログラミング言語を用いて実行可能なコンピュータプログラムに翻訳され、予め、HDD303等に保存される。情報処理装置における発塵源探索処理の実行時には、CPU等の演算装置301によって、HDD303等に記憶された前記実行可能なコンピュータプログラムが読み出されて起動し、前記実行可能なコンピュータプログラムの指令に基づいた演算をCPU等の演算装置301が順次実行することにより実現される。発塵源探索処理の起動タイミングは、手入力で前記実行可能なコンピュータプログラムを起動してよく、また、定期的に自動的に起動するようにしてもよい。前述したように、本実施形態の発塵源探索装置は、ある時刻において、Tdit期間における降下煤塵の発塵源を探索する。発塵源探索装置において、評価地点座標点等の位置情報、降下煤塵量風向風速等の測定値や煤塵種に関する分析値等の必要入力情報は、情報処理装置に接続されたキーボードやコンソール画面等を用いて、予め人力で入力することができる。入力された前記入力情報は、HDD303等に保存され、発塵源探索処理実行の進行に応じて、適宜、読み出される。発塵源探索装置において、算出された特定座標点に対する非定常発塵源判定結果及び発塵量等の算出結果は、HDD303等に保存されると共に、コンソール画面等に表示することができる。尚、上記の発塵源探索装置の処理の一部又は全部を手計算等の他の手段に置き換えても何ら問題ない。図1のステップS101において、降下煤塵評価地点ix、降下煤塵参照地点iy、及び仮定発塵源jからなる検討系を設定する。図2の例において、これら検討系では、降下煤塵評価地点ix、降下煤塵参照地点iyは、それぞれ、地上2mの高さで降下煤塵の測定を連続的に測定する地点である。ここでは、特定の降下煤塵評価地点ix:1での発塵源を特定するために、iy:1~iy:8なる8点の降下煤塵参照地点iyを設けている。降下煤塵評価地点ix及び降下煤塵参照地点iyには、時系列的に連続に降下煤塵量を測定できる測定手段を設ける。仮定発塵源jは、任意の数を設定可能であるが、図2の検討系では一例として、2つの仮定発塵源jを配置する。これらの発塵源の地上高さhjとしてhj1及びhj2をそれぞれ設定する。以上の対象系に対して、仮定発塵源j:1、j:1の内いずれが降下煤塵評価地点ixに関する主要発塵源であるかを、後続のステップを通じて判定する。ステップS102において、発塵源を探索する降下煤塵分類、即ち、煤塵種is及び粒子落下速度粒径区分ivの組み合わせを設定する。煤塵の種類の分類には後述の分析方法のいずれかを用いて行う。発塵源特定の解析は、ここで分類された特定の煤塵の種類例えば、降下煤塵評価地点において最も構成率の高い煤塵の種類を採用することができる毎に実施する。全ての降下煤塵測定地点において、当該Tgk期間における降下煤塵分類毎の質量又は体積構成率Comix,iv,isを求める。当該降下煤塵評価地点ixにおいて最大の構成率Comix,iv,isをCommaxixと定義すると共に、このときの粒子落下速度区分iv及び煤塵種区分isをivmaxix及びismaxixと定義する。ステップS103において、当該Tgk期間におけるixに関する最大降下煤塵量mmaxix及び最大降下煤塵量を生じるTdit時刻Tditmaxを求める。尚、nit=1の場合、mmaxix=mix,it、かつ、itmax=1である。ステップS104において、各降下煤塵測定手段における代表降下煤塵量、並びに、当該Tgk期間における代表風向WD及び代表風速WSを求める。降下煤塵評価地点ixにおける代表降下煤塵量Mixとして、ステップS103で求めたmmaxixを用いることができる。各降下煤塵参照地点iyにおける代表降下煤塵量Miyとして、ステップS103で求めたitmaxを用いて、それぞれmmaxiy,itmaxを用いることができる。代表風向WD及び代表風速WSは、風向風速計測値を用いて次のように設定する。当該Tgk期間における代表風向及び代表風速値として、例えば、前記Tditmax期間に含まれる風向及び風速計測値の平均値を用いることができる。このように、代表風向値及び代表風速値を設定することによって、当該Tgk期間において当該降下煤塵評価地点に関する主要発塵源の方向についての情報を的確に得ることができる。あるいは、より単純に、Tgk期間の風向平均値及び風速平均値をそれぞれ代表風向値及び代表風速値として用いてもよい。代表風向及び代表風速が決まれば、代表風のベクトル1を算出できるので、このイメージを図2上に示す。風向及び風速は、降下煤塵評価地点、又は、降下煤塵評価地点の近傍において、Δtdよりも長くない時間周期Δtwdにて連続的に計測を行う。風向及び風速の計測手段としては、プロペラ式風向風速計等の市販の計測機を用いることができる。この装置の場合、例えば、Δtwdを1秒に、風向測定の空間分解能として、1°に設定できる。ここで、降下煤塵評価地点の近傍とは、風向風速が降下煤塵評価地点上空における風向風速と高い相関を示す範囲であればよく、例えば、降下煤塵評価地点から1km以内の水平距離とすることができる。地形が単調で風向風速分布の少ない地域ではこれ以上の水平距離であってもよい。また、風向風速測定点の高さは、気象庁の推奨する測定高さである、地表面から10mを採用することができる。想定される発塵源の高さが10mよりも十分に高い場合には、地表面とこの発塵源高さの中間の高さを測定点高さとしてもよい。ステップS105において、全ての降下煤塵参照地点iyに対して、数式6で規格化された降下煤塵量測定値Raiyを求める。Raiy=Comiy,ivmaxix,ismaxixMiy/CommaxixMix...6ステップS106において、特定の仮定発塵源jを発塵源とした場合の当該降下煤塵量評価地点ix及び全ての降下煤塵参照地点iyでの降下煤塵量の計算値mcj,ix及びmcj,iyを求める。これらmcj,ix及びmcj,iyの算出には、例えば、数式3~6で表現されるプルーム式を用いることができる。後述のように、発塵量の絶対値の影響は降下煤塵を規格化することによって除外されるので、これらの式で用いる発塵量Qpは任意の値を設定すればよい。例えば、発塵量Qp=1としてよい。Cx,y,z=QP/2πσyσzWSexp[-y2/2σy2]×exp[-He-z-Vsx/WS2/2σz2]...4mcx,y=VsCx,y,z...5ステップS107において、全ての降下煤塵参照地点iyに対して、数式7で規格化された降下煤塵量計算値Rciyを求める。Rciy=mcj,iy/mcj,ix...7ステップS108において、前記規格化計算降下煤塵量Rciyと前記規格化実績降下煤塵量Raiyの比率Rciy/Raiyを用いて、当該仮定発塵源jに関する後述の評価関数値Evjを算出する。全ての仮定発塵源jj=1~njについて評価関数値Evjを求めた後、ステップS109において、最も評価関数値Evjの評価の高い即ち、RciyとRaiyの一致性の元雄も高い仮定発塵源を当該Tgk期間における降下煤塵評価地点ixに関する主要な非定常降下煤塵発塵源と判定する。降下煤塵の分類方法:落下速度区分前記捕集された個々の降下煤塵粒子を分類する方法は、例えば、次のように実施することができる。まず、粒子の落下速度区分ivを判別する方法として、捕集された個々の煤塵粒子を顕微鏡等で観察測定して、寸法や形状を記録する。次に、個々の煤塵粒子の寸法形状密度を用いて当該粒子の等価粒径を算出し、この等価粒径に対して、所定の粒径しきい値を用いて分類することができる。一般に、同一の煤塵種に関して、価粒径が同一の粒子は、粒子落下速度Vsが同一になると考えられているので、等価粒径を用いた粒子の分類は、実質的に同一煤塵種におけるVsによる分類と等価である。等価粒径の算出には各種の手法が提案されており、例えば、個々の煤塵粒子の体積と等しい体積の球形の直径に、煤塵粒子の形状アスペクト比等に基づいて予め経験的に定められた補正係数を乗じて等価粒径を算出することができる。粒子の密度は、煤塵種毎に実測してもよいし、文献値が存在するならこれを用いてもよい。粒子落下速度Vsは、実測してもよいし、あるいは、価粒径をストークスの終末速度の式に適用して次の数式8で求めてもよい。Vs={4gDpρp-ρf/3ρfCR}12...8g:重力加速度Dp:等価粒子径ρP,ρf:粒子、流体の密度CR:抵抗係数降下煤塵の分類方法:煤塵種次に、粒子の煤塵種区分isを判別する方法として、放射性煤塵粒子の飛散を問題とする場合には、捕集した降下煤塵粒子の放射線量を測定し、その強さに応じて降下煤塵粒子を放射性粒子、非放射性粒子等に分類することができる。捕集した効果煤塵粒子の放射線量の測定方法には、特許文献9~11に記載の方法等の公知の手法を適用することができる。煤塵粒子の放射性を問題とせずに、他の性質によって粒子の煤塵種区分isを判別する方法として、当該Tgk期間に捕集された個々の降下煤塵粒子の物性を分析すればよい。粒子の物性分析方法としては、例えば、全ての粒子を樹脂等に埋め込み、研磨によって粒子断面を露出させた上でX線を走査しながら粒子断面に照射してその反射の特性から粒子中の元素を特定するEPMAを適用することができる。あるいは、EPMAのような手間と費用のかかる手法を避けるために、高炉法による製鉄プラントに由来する降下煤塵の煤塵種であることが事前に明らかな場合、簡易に煤塵種を判別する次の方法を用いてもよい。即ち、第1の工程において、捕集された前記降下煤塵試料に磁力を付与して前記降下煤塵を磁石に着磁する着磁性降下煤塵と磁石に付着しない非着磁性降下煤塵に分離する。ここで用いる磁石は、0.1T~0.4T程度の磁力を表面で保持できる電磁石又はネオジウム磁石等の永久磁石である。そして、捕集された前記降下煤塵試料に磁力を付与して前記降下煤塵を磁石に着磁する着磁性降下煤塵と磁石に付着しない非着磁性降下煤塵に分離する。このような磁力の範囲であれば、大部分の鉄鉱石や製鋼スラグ一般に鉄分を多く含んでいるを着磁させることができ、かつ、極微量な鉄分を含有する石炭等を着磁させないので、個々の煤塵を着磁性煤塵と非着磁性種煤塵に分離することができる。また、高炉法による製鉄プラントに由来する降下煤塵の煤塵種としては、鉄鉱石や鉄粉等の鉄系煤塵、石炭やコークス等のカーボン系煤塵、高炉スラグ煤塵、並びに、製鋼スラグ煤塵等が代表的なものである。着磁性の煤塵には鉄系煤塵及び製鋼スラグ煤塵が対応し、非着磁性の煤塵にはカーボン系煤塵及び高炉スラグ煤塵が対応する。第2の工程において、前記着磁性煤塵試料と前記非着磁性種煤塵試料のそれぞれに対して、煤塵粒子が互いに接触しないように二次元的に分散配置させた上で、これらを市販のデジタルカメラ等で撮影して粒子の画像を得る。前記着磁性煤塵試料の画像を着磁性煤塵試料の画像と、前記非着磁性煤塵試料の画像を非着磁性煤塵試料の画像と、以下に呼ぶことにする。煤塵試料の個々の粒子を互いに接触させないためには、煤塵粒子を高所から散布する等の方法を用いればよい。第3の工程において、前記着磁性降下煤塵画像及び前記非着磁性降下煤塵画像は、画像処理装置に入力されて、それぞれ画像処理がなされる。個々の画像になされる画像処理の内容は、次のとおりである。まず、画像中の各画素の位置及び明度情報を基に、独立した粒子と判別される連続する画素群を算出する粒子判別を行う。次に、個々の降下煤塵粒子の代表位置及び代表明度を算出する。代表位置としては、個々の粒子において対応する画素群の各画素位置の中心を用いることができる。代表明度としては、個々の粒子において対応する画素群の各画素明度の平均値を用いることができる。次に、前記代表明度と所定の明度しきい値を比較して個々の粒子を暗色粒子と明色度粒子に粒子明度を区分する。所定しきい値に関しては、予め、明色粒子群と暗色粒子群の代表試料を準備して上記と同様の画像処理を行い、得られたそれぞれの粒子群の明度平均値の中間の明度をしきい値として用いることができる。次に、前記粒子明度区分及び前記第1の工程で得られた前記個々の降下煤塵の着磁性の有無との組み合わせを用いて、前記着磁性粒子画像及び非着磁性粒子画像中の全ての粒子を着磁性暗色粒子、着磁性明色粒子、非着磁性暗色粒子、並びに、非着磁性明色粒子のいずれかに分類して煤塵特性を区分する。画像処理装置として、市販の画像処理ソフト例えば、ImageProPlusversion5等を組み込んだ市販のパーソナルコンピュータ等を用いることができ、上記の画像処理操作は、このようなソフトの標準機能画像の、二値化境界判別粒子計測機能等を用いて実現することができる。第4の工程において、前記着磁性粒子画像及び非着磁性粒子画像中の全ての粒子を、前記煤塵特性に基づいて、所定の煤塵種のいずれかに判定する。所定の煤塵種とは、代表的な同一種類鉄鉱石、石炭、コークス、鉄粉、高炉スラグ、転炉スラグ等の粒子試料に対して上記の磁力選別及び画像処理による粒子明度判別を行って、平均的な煤塵特性着磁性暗色粒子、着磁性明色粒子、非着磁性暗色粒子、並びに、非着磁性明色粒子のいずれかを求める。本発明者の調査結果から、鉄鉱石及び鉄粉は着磁性暗色粒子に分類されるので、着磁性暗色粒子を鉄系煤塵として所定の煤塵種の1つとして採用できる。また、石炭及びコークスは非着磁性暗色粒子に分類されるので、非着磁性暗色粒子をカーボン系煤塵として所定の煤塵種の1つとして採用できる。さらに、製鋼スラグと高炉スラグは、着磁性明色粒子と非着磁性暗色粒子にそれぞれ分類されるので、着磁性明色粒子を製鋼スラグ系煤塵として、非着磁性明色粒子を高炉スラグ系煤塵として、それぞれ所定の煤塵種の1つとして採用できる。煤塵質量の算出方法煤塵種毎に等価粒径を用いて捕集された全降下煤塵粒子を分類した後、各粒径区分に含まれる個々の粒子の密度及び体積を用いて算出された個々の粒子の質量を積算することにより、当該煤塵種に当該粒径区分における降下煤塵質量を求めることができる。さらに、降下煤塵の全分類に対して同様に求めた降下煤塵質量を積算して全降下煤塵質量を算出し、当該分類の降下煤塵質量と全降下煤塵質量の比を降下煤塵の当該分類の質量構成率Comiy,iv,isとして算出することができる。評価関数上記の規格化計算降下煤塵量Raiyと規格化計算降下煤塵量Rciyとの一致性を評価する評価関数として、例えば、次の関数を用いることができる。ここで、対数の底は、10やネイピア数eを用いることができる。また、降下煤塵量測定値に負値を生じない場合測定誤差の小さい場合には、数式9中で|Rc/Ra|の代わりに単にRc/Raを用いてもよい。数式9の関数において、全ての参照地点において、RcがRaと一致する場合即ち、降下煤塵の実測値分布と計算値分布の一致度が最も高い場合、評価関数値Evは0となる。RcとRaが異なれば、Ev>0となり、RcとRaの違いが大きいほど、評価関数値Evの値は大きくなる。従って、複数の仮定発塵源おいてそれぞれ算出された評価関数値Evの内、最も小さい値を示す仮定発塵源において、降下煤塵量分布の計算値と実測値の一致度が最も高いと判断することができる。また、前述のように、この評価関数を用いれば、特定の降下煤塵評価地点で仮にRcとRaとの間でオーダの異なる場合でも、評価関数値が過度にこの特定の降下煤塵評価地点の影響を受けることがない。このため、測定誤差や計算モデル誤差が比較的大きく、いずれの仮定発塵源での評価においてもRcとRaとの間でオーダの異なる降下煤塵評価地点を含むことが避けられないような検討系に対しても、仮定発塵源間での妥当性を安定的に比較評価することができ、好適である。この評価関数の場合、評価関数値Evの小さいものほど、降下煤塵の発塵源としての妥当性が高いことを意味する。この評価関数のバリエーションとして、例えば、Tgk期間毎にniyが変動する場合に、この評価関数をniyで除したものを評価関数として用いてもよい。あるいは、測定誤差や計算モデル誤差が十分に小さいと期待できる場合には、評価関数として、次の数式10に示す関数を用いてもよい。あるいは、次の数式11に示す評価関数として用いてもよい。ここで、max[a,b]は、数値a,bの内、大きい方の値を表す関数である。このように、評価関数にRc/Raを用いることにより、仮定発塵源の妥当性評価に際して、個々の降下煤塵量評価地点における降下煤塵量の多寡がその降下煤塵参照地点の評価関数に対する重みに影響しなくなるので、降下煤塵量の多寡によって特定の降下煤塵参照地点の影響が評価関数値算出時に支配的となるという構造的な問題を回避することができる。以上述べたような本発明を適用した降下煤塵の非定常発塵源の探索手法の第1の特徴は、降下煤塵評価地点における降下煤塵を直接、測定することによって降下煤塵の発塵源を探索することができる点である。第2の特徴は、本発明者が見出した降下煤塵用に修正したプルーム式を用いて同一時刻Tg期間における複数の降下煤塵測定地点における降下煤塵測定量間の整合性を評価することにより、最も整合性の高い仮定発塵源を主要な発塵源と判定することである。本発明で修正されたプルーム式について具体的に説明する。従来のプルーム式を用いた降下煤塵量の整合性あるいは妥当性の評価の際には、数式2における地表面反射項αexp[-Hez-Vsx/WS2/2σz2]の取り扱いが困難なため、高い精度での評価が困難であった。しかし、本発明者の調査の結果、この地表面反射項が問題となるのは、従来技術が主にガスやSPMを対象としていたためであることを突き止めた。降下煤塵の場合には、粒子の落下速度が大きいために、沈着速度Vd≒落下速度Vsとなるので、地表面での反射の影響は小さく、α=0と見做すことができる。そこで、修正した降下煤塵濃度に対する大気拡散式として、数式2にα=0を代入した次の数式4を用いる。Cx,y,z=QP/2πσyσzWSexp[-y2/2σy2]×exp[-He-z-Vsx/WS2/2σz2]...4ここで、プルーム拡散幅水平σy及びプルーム拡散幅鉛直σzは、Δtdを固定してこれを基準期間とする求められた数表化又は図表化された値として、非特許文献1に記載される、Pasquill-GiffordによるものやBriggsによるもの等を用い、Δtdの影響を経験式で補正して求められる。Δtdの影響を経験式で補正する方法は、非特許文献2に示されるように、σyに、[実際に使用するΔtd]/[基準純時間のΔtd]Pを乗じるものである。煤塵種と煤塵粒径が与えられれば、粒子落下速度Vsが終末速度として決まるので、降下煤塵量計算値mcは、濃度CにVsを乗じた次の数式5で表現できる。mcx,y=VsCx,y,z...5このように、気象地表状態煤塵種等の影響によって大きく変化してプルーム式の予測精度を悪化させ得る要因である経験常数αを用いずに降下煤塵量を表現できるので、プルーム式での予測誤差が小さい利点がある。以上のような降下煤塵量のプルーム式を変形した発塵源探索範囲の単純、かつ、定量的な表現は、従来のガスやSPMを前提としたプルーム式では実現し得なかったものであり、本発明者が降下煤塵ではVsが比較的大きいことに着目した上で行った一連の洞察によって、初めて可能になったものである。第3の特徴は、降下煤塵の評価地点ixにおける主要発塵源を評価する際に、仮定発塵源を設定する地域とixとの間に3つ以上の降下煤塵観測地点である参照地点iyを設けてiyでの降下煤塵量測定値を用いることである。従来技術では仮に複数地点での降下煤塵測定値を用いて発塵源の妥当性を評価する場合であっても、各測定地点の関係は本質的に同等であり、主たる評価地点ixに対する主要発塵源を評価するために、補助的に他の測定点iyを設けるという発想は、従来存在しなかった。このような測定地点レイアウトを設定することにより、仮定発塵源を設定する地域から評価地点ixに到達する煤塵プルームの中心位置を、ixの風上に存在する3つ以上の参照地点iyでの測定値を互いに比較することによって容易に判定することができる。もし、参照地点iyが2つしかない場合には、降下煤塵量測定値のより大きい測定地点とより小さい地点の内側に煤塵プルームの中止が存在するのか、あるいは、降下煤塵量測定値のより大きい測定地点の外側より小さい地点の反対方向に煤塵プルームの中心が存在するのかを判断することができない。3つ以上の参照地点iyを設ける場合には、煤塵プルームの中心が両端のiyの内側、外側いずれに存在するのかは自明に判断できる。即ち、両端の参照地点iyのいずれかで最大の降下煤塵量を測定値を観測すれば、煤塵プルームの中心は外側に存在し、両端以外の参照地点iyで最大の降下煤塵量を測定値を観測すれば、煤塵プルームの中心は内側に存在し、最大の降下煤塵量の測定値を観測すれば、煤塵プルームの中心は外側最大の測定値を観測した参照地点iyの近傍に煤塵プルームの中心が存在すると判定できる。このように、煤塵プルームの中心位置の内外を判定しないと、仮に、プルーム式を用いて測定値とよく一致する発塵源を見出したとしても、その発塵源が煤塵プルームの中心一意であるのか判断できない。尚、ここで、降下煤塵量がより多い少ないと表現しているのは、仮定発塵源と参照点間の距離の影響を補正して降下煤塵量を比較したものである。第4の特徴は、煤塵を捕集して粒径区分と煤塵種を判別し、粒径区分と煤塵種毎にプルーム式を設定して、複数の地点で測定された降下煤塵量間の整合性を評価することである。発塵源が単一である場合には、上述した第3の特徴における方法のみで発塵源の位置を特定することができる。しかし、そもそも本発明が目的とする発塵源の特定が必要となる対象系とは、発塵源が複数存在するために主要な発塵が自明でないものに限られる。このような系では、複数地点で降下煤塵を測定した場合に、個々の測定地点では、一般に複数の発塵源からの降下煤塵が捕集測定されている。従って、図4に示すように、主要発塵源探索の過程で仮定発塵源点源を設定してプルーム式での降下煤塵計算結果と実測値を比較しても、計算値でのように、プルーム端部で降下煤塵量が減少する傾向を実測値が示すとは限らないため、整合性の評価が困難となる前記第3の特徴の方法のみでは発塵源を特定できない。例えば、図4での仮定発塵源j:1、j:2のどちらがより整合性が高いかを評価できない。ここで、図4は、水平面上の配置図であり、降下煤塵評価地点ix及び降下煤塵参照地点iyにおいて、測定された降下煤塵量に比例する面積の円網かけ部分を各地点の周囲にそれぞれ記載したものである。尚、符号2は仮想発塵源の風下直下を示す。ところで、多くの発塵地帯において、発塵源で発生する煤塵種は、必ずしも同一ではなく、煤塵種毎に発塵源を分類できる場合が多い。そこで、本発明では、捕集された煤塵種を細かく煤塵種に分類することにより、各降下煤塵測定地点での煤塵種毎の降下煤塵量を算出し、当該煤塵種のみに関して、プルーム式による降下煤塵量分布の整合性を評価する。このようにすることで、個々の整合性評価において、近接する発塵源から対象を絞り込むことができるので、複数発塵源からの同一地点への降下煤塵飛来による、前記の整合性評価困難性を大幅に低減することができる。尚、前記の近接する発塵源から対象を絞り込む際に、個々の発塵源候補での主要煤塵種を予め把握しておく必要は必ずしもなく、近接する発塵源間で煤塵種が互いに異なりさえすれば、自動的に、前述のように本発明での発塵量分布整合性評価が高精度化するのである。さらに、同一発塵源であっても、発塵する煤塵の粒径には大きな分布が一般に存在する。発塵源毎に煤塵の粒径分布が予め把握できるとは限らないので、従来技術においては、全ての煤塵粒径を含んだ煤塵プルームを用いて妥当性の検討がなされてきた。しかし、煤塵が飛散する際には粒径の大きなものほど早く沈降するので、煤塵プルームが進行するに従って、プルーム中の平均粒径は減少する。正確な粒径分布の把握なしにこの影響を正確にプルーム式で評価することは極めて困難であるため、従来技術でのプルーム式を用いた降下煤塵量の妥当性判断では大きな誤差が避けられなかった。一方、本発明では、煤塵の粒径区分毎にプルーム式を設定するので、このような粒径分布に起因するプルーム式の誤差を回避することができ、降下煤塵量分布整合性評価の精度が高いという利点がある。図5は、図4と同様の対象系における本発明での降下煤塵量分布の計算値‐実績値間の対応を示している。図5では、実績値分布がプルーム式と相似に近い状態中心で高い値、端部で低い値を示しており、計算値との整合性評価が容易である。即ち、仮定発塵源j:1では整合性が低く、仮定発塵源j:2では整合性が高いと評価できる。第5の特徴は、計算降下煤塵量分布と実績降下煤塵量分布を比較する際に、参照地点での降下煤塵量を、評価地点での降下煤塵量によって規格化したものを比較することである。もし、規格化を行わずに、降下煤塵量の計算値-実測値間の比較を直接行う場合には、発塵源探索を行う際に設けられる仮定発塵源での発塵量を正確に与える必要がある。しかし、発塵源での発塵量を事前に把握できない場合が多いので問題である。一方、本発明では、降下煤塵量の規格化を行っているので、仮定発塵源に正確な発塵量を与える必要がない。第6の特徴は、非定常発塵源の特定を行うことができる点である。本発明の発塵源探索手法では、降下煤塵量測定値取得周期毎にその時間帯における主要な発塵源を特定することができる。従って、降下煤塵量測定値取得周期の数周期分以上の時間スケールで変動する非定常発塵源であれば、これを把握することができる。また、非定常発塵源を特定する際に必要な降下煤塵評価地点数は、潜在的な発塵源数よりも十分少なくてよい。第7の特徴は、任意の仮定発塵源に対して、複数の参照地点での規格化された降下煤塵量計算値と実測値を用いて、計算値‐実測値間の一致性を評価する方法として、各参照地点での降下煤塵量計算値と実測値の比の絶対値の対数を、全参照地点分、積算することによって与えられる評価関数を用いることである。仮定発塵源と各参照地点間の距離が大きく分布するときには、参照地点間で測定値間の値のオーダの異なる場合がある。このとき、各参照地点での誤差同一参照地点での降下煤塵量計算値と実測の差等で表現できるを単に積算するだけでは、発塵源に近く、降下煤塵量又は、前記規格化された降下煤塵量が元々、多い参照地点での誤差に対して一致性評価の際に極端に大きな重みが与えられるので、全評価地点での平均的な一致性を評価するのに適当でない場合がある。また、個々の参照地点での降下煤塵量計算値と当該参照地点での実測値の比を単純に積算することも考えられる。しかし、この場合、実績降下煤塵量が極端に小さい参照点ではこの比が極端に大きな値を示し、一致性評価の際にこの参照点に対して極端に大きな重みが与えられるので、全評価地点での平均的な一致性を評価するのに適当でない場合がある。一方、本発明では、規格化された降下煤塵量の計算値と実測値の比の絶対値の対数を用いるため、仮に参照地点間で降下煤塵量や降下煤塵量の比のオーダが異なる場合でも、一致性評価に際して各参照点に生じる重みを比較的一様に保つことができるので、一致性評価手法として好適である。第8の特徴は、対象とする降下煤塵種が放射性である場合、強い放射性発塵源から発生したと考えられる放射性降下煤塵の発生源位置を、放射性発塵源に近寄ることなく、遠方での降下煤塵計測データを用いて推定することができる。第9の特徴は、捕集された煤塵試料に対して強力磁石による磁力選別と画像処理による明度測定を用いて煤塵種の判別を行うことにより、安価かつ簡易に煤塵種分析を実施可能なことである。従来、このような技術が困難であったのは、以下の理由によるものである。代表的な鉄鉱石である赤鉄鉱は、代表的な永久磁石であるフェライト磁石の磁力では磁石に吸着しない。このため、一般の鉄鋼業技術者は、赤鉄鉱と石炭の分離に磁力を用いるという発想をそもそも持ち得なかった。また、本発明者の調査の結果、単に赤鉄鉱が着磁し得る強力な磁石を用いれば煤塵微粒子を前記の代表的な煤塵種毎に分離できる訳ではないことを本発明者は見出した。これは、強力過ぎる磁石を煤塵種に接触させた場合、僅かな鉄分鉄鉄鋼銑鉄酸化鉄を含む煤塵でも着磁してしまうので、単に、全く鉄分を含まない微粒子と鉄分を含む微粒子に選別できるのみである。高炉法による製鉄プラントで発生する煤塵は、多かれ少なかれ鉄分を含有することが普通なので石炭には灰分中に1程度の酸化鉄が含まれことが多い。また、高炉スラグにも通常、少なくとも03程度以上の鉄分が含有されている、実質的には全く鉄分を含まない微粒子とは製鉄所内土壌由来の石英系土砂等のような主要ではない煤塵種に限定される。従って、このような磁力選別は、少なくとも製鉄プラントでの煤塵種の判別にはあまり役立たない。そこで、本発明では、適度に強い磁石を用いることで、鉄分を含むが含有量の少ない煤塵種石炭系と高炉スラグ系を着磁させないことで、着磁性の煤塵種鉄系と製鋼スラグ系から分離することができる。さらに、煤塵微粒子の明度を用いて暗色煤塵種石炭系鉄系と明色煤塵種高炉スラグ系製鋼スラグ系に分離できる。その結果、高炉法による製鉄プラントで発生する前記主要な4つの煤塵種に限定すれば、着磁性と明度の組み合わせを用いることで、個々の煤塵を実用的、かつ、簡易に、全ての主要な煤塵種に分類できる。このことを本発明者は初めて見出した。これに対して、一般の粒子分析技術者特に電子顕微鏡を用いた分析技術者は、より汎用性の高い煤塵種分析手法を指向していた。このため、本発明法でのように、単に4つの主要な煤塵種に個々の煤塵を分類することだけでも高い実用性を有する分析を実現できることに、一般の粒子分析技術者は思い至らず、一般の粒子分析技術者は、徒に分析可能な元素数分子数を増やすことや、成分構成率の高精度化に労力を傾注していたため、本発明法に想至することはなかった。具体的値を用いて、本発明を適用した実施例を説明する。図1の処理フローに従って、非定常発塵源の探索を行った。まず、ステップS101において、図2に示す降下煤塵評価地点ix、降下煤塵参照地点iyの水平面内配置において、j:1及びj:2なる2つの仮定発塵源を図中に配置した。ステップS102において、最大構成率の煤塵特性として、粒径区分が粒径30~60μm、かつ、煤塵種として石炭煤塵カーボン系煤塵を採用した。以下、試料内のこのグループを当該煤塵特性と表記する。ここで、降下煤塵評価地点ix及び全降下煤塵参照地点iyにおける当該煤塵特性の当該地点での試料中の構成率を100とした。nit=6、Δtd=10分とし、特定のTgk期間において、ステップS103でixにおいて最大降下煤塵量を観測するTdit時刻として、itmax=3を採用した。ステップS104において、降下煤塵評価地点ix及び各降下煤塵参照地点iyでのTditmaxにおける降下煤塵量測定値をそれぞれの地点での代表降下煤塵量として採用し、その値は、表1に示されるものとした。風向及び風速は、超音波式風向風速計を用いて1秒周期で測定を行ったものとして、代表風向風速として、Tditmax期間の風ベクトル平均である、図2に示す風ベクトル1に対応する風向風速をそれぞれ採用した。また、同様に、ステップS106において、j:1及びj:2を発塵源と仮定した場合の降下煤塵評価地点ix及び各降下煤塵参照地点iyにおける降下煤塵量計算値を、数式3のプルーム式及び数式4を用いてそれぞれ表1に示す値とした。ステップS105において、各降下煤塵参照地点iyにおけるRaは、数式6を用いて表2のように算出された。ここで、降下煤塵評価地点ix及び各降下煤塵参照地点iyは、上記のように全て当該煤塵特性に含まれるものとしたので即ち、Comiy,ivmaxix,ismaxix=Commaxix=1、数式6は、各iyでの降下煤塵量測定値をixでの降下煤塵量測定値で除した数式に簡略化される。同様に、ステップS107において、各降下煤塵参照地点iyにおけるRcは、数式7を用いて表2のように算出された。ステップS108において、まず、上記の結果を用いて、log10[|Rc/Ra|]2を計算すると、表3のようになった。次に、上記log10[|Rc/Ra|]2の値を数式9に代入して、表3の最下段に示す評価関数値Evを求めた。この評価関数値Evは、j:1について1.487556559、j:2について19.68380755となる。ステップS109において、上記の結果から、j:1の評価関数値Evがj:2でのものよりも小さいので、仮定発塵源j:1をix:1に関する当該時刻における主要発塵源と判定した。"