JP2011117806A - 粒度測定方法および粒度測定システム - Google Patents
粒度測定方法および粒度測定システム Download PDFInfo
- Publication number
- JP2011117806A JP2011117806A JP2009274760A JP2009274760A JP2011117806A JP 2011117806 A JP2011117806 A JP 2011117806A JP 2009274760 A JP2009274760 A JP 2009274760A JP 2009274760 A JP2009274760 A JP 2009274760A JP 2011117806 A JP2011117806 A JP 2011117806A
- Authority
- JP
- Japan
- Prior art keywords
- particle size
- displacement data
- measured
- measuring
- granular material
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 238000005259 measurement Methods 0.000 title claims description 11
- 238000006073 displacement reaction Methods 0.000 claims abstract description 77
- 239000008187 granular material Substances 0.000 claims abstract description 33
- 238000005070 sampling Methods 0.000 claims abstract description 20
- 239000002245 particle Substances 0.000 claims description 69
- 238000012545 processing Methods 0.000 claims description 8
- 238000000691 measurement method Methods 0.000 claims description 7
- 239000000843 powder Substances 0.000 abstract description 13
- 238000012546 transfer Methods 0.000 abstract description 3
- 239000003245 coal Substances 0.000 description 39
- 238000007873 sieving Methods 0.000 description 7
- 239000000571 coke Substances 0.000 description 6
- 238000000197 pyrolysis Methods 0.000 description 2
- 230000008021 deposition Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000011802 pulverized particle Substances 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000010025 steaming Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
【課題】搬送装置によって搬送される粉粒体の粒度をオンラインであっても容易に精度よく求めることができる粒度測定方法および粒度測定システムを提供する。
【解決手段】搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法であって、搬送される粉粒体表面の変位を所定のサンプリング周期で計測するステップと、計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理するステップと、数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出するステップと、を備える。
【選択図】図1
【解決手段】搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法であって、搬送される粉粒体表面の変位を所定のサンプリング周期で計測するステップと、計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理するステップと、数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出するステップと、を備える。
【選択図】図1
Description
本発明は、搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法および粒度測定システムに関するものである。
コークスは、石炭をコークス炉で乾留(蒸し焼き)することにより製造される。原料となる石炭は、粉砕機で細かく粉砕された後、ベルトコンベア等によりコークス炉へ搬送される。搬送された石炭は、コールビンに一旦貯蔵された後、装入車でコークス炉に装入され、乾留される。
コークス炉に装入される石炭の粒度は、装入する際の装入量や、製造されるコークスの品質などに影響する。そのため、装入される石炭の粒度を測定し、管理する必要がある。
下記特許文献1には、ベルトコンベアで搬送されている粉粒体の粒度および粒度分布をリアルタイムに自動測定する方法として、ベルトコンベア上に載せられて移動している粉粒体の表面の変位を所定のサンプリング時間で測定し、こうして得られた変位データを時系列にプロットした図形から粒子間の接触もしくは重なりを示す部分を不連続点として検出し、これより粒径を測定して粒度分布を求める方法が開示されている。
しかしながら、引用文献1に開示された方法は、粒子がほぼ完全な球形をしている場合には、粒子間の接触もしくは重なりを示す不連続点を検出可能であるが、粉砕された粒子がほぼ完全な球形をしていることは稀であり、粒子間の接触を示す不連続点のみを検出することは通常は困難である。そのため、不連続点を検出し、これより粒径を測定し、粒度分布を容易に精度よく求めることは困難である。
また、粉粒体の表面の変位を計測する場合、粉粒体の粒径の差に基づく変位に、重力方向の上下変動が重なることは、避けることができない。この上下変動が生じる理由として、回転機械特有の機械的な振動や、複数のベルトコンベアに粉粒体が乗り継ぐ時に発生する落下衝撃による揺れ、ベルトのエンドレス加工部に見られる凹凸面の存在などが考えられる。これら避けられない上下変動によって、粉粒体の表面の変位を計測しただけでは、粉粒体の粒度を測定することは不可能である。
本発明は上記実情に鑑みてなされたものであり、その課題は、搬送装置によって搬送される粉粒体の粒度をオンラインであっても容易に精度よく求めることができる粒度測定方法および粒度測定システムを提供することである。
上記の課題を解決するために、本発明の粒度測定方法は、
搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法であって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するステップと、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理するステップと、
数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出するステップと、を備えることを特徴とするものである。
搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法であって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するステップと、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理するステップと、
数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出するステップと、を備えることを特徴とするものである。
かかる構成による粒度測定方法の作用効果を説明する。搬送装置によって搬送される粉粒体表面の変位を所定のサンプリング周期で計測し、計測された粉粒体表面の計測変位データを時系列にプロットすると、大きな上下変動と小さな上下変動が組み合わされた形状となる。本発明者らが鋭意研究した結果、大きな上下変動は、粉粒体が載せられた搬送装置表面の変動と、粉粒体の表面プロフィールの変動とが組み合わさったものであり、小さな変動は、粉粒体を構成する個々の粒子の大小による変動であることが分かった。本発明によれば、計測変位データを数値化処理することにより、計測変位データから小さな上下変動が取り除かれ、大きな上下変動のみからなる数値化処理変位データを得ることができる。この数値化処理された数値化処理変位データと計測変位データとの差は、小さな上下変動に相当し、上述のように小さな上下変動は粉粒体を構成する個々の粒子の大小による変動である。さらに、本発明者らは、数値化処理された数値化処理変位データと計測変位データとの差から求められる誤差分散が、粉粒体を構成する粒子の大きさの割合である粒度と高い相関関係があることを見出した。よって、数値化処理された数値化処理変位データと計測変位データとの差から誤差分散を算出することで、粉粒体の粒度をオンラインであっても容易に精度よく求めることができる。
本発明に係る粒度測定方法において、予め作成しておいた誤差分散と粒度の関係から、前記算出された誤差分散に対応する粒度を求めるステップを備えることが好ましい。
上記のように誤差分散と粒度との間には高い相関関係があるので、誤差分散と粒度の関係を予め作成しておくことができる。この誤差分散と粒度の関係に、数値化処理変位データと計測変位データとの差から算出された誤差分散を当てはめ、対応する粉粒体の粒度を求めることができる。
本発明に係る粒度測定方法において、前記数値化処理は、時系列にプロットした計測変位データを非線形のカーブフィッティングにより数値化することが好ましい。
時系列にプロットした計測変位データを非線形のカーブフィッティングにより数値化することで、計測変位データから小さな上下変動を容易に取り除くことができる。
上記の課題を解決するために、本発明の粒度測定システムは、
搬送装置によって搬送される粉粒体の粒度を測定する粒度測定システムであって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するレーザー変位計と、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理し、数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出する演算処理装置と、を備えることを特徴とするものである。
搬送装置によって搬送される粉粒体の粒度を測定する粒度測定システムであって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するレーザー変位計と、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理し、数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出する演算処理装置と、を備えることを特徴とするものである。
かかる構成による粒度測定システムの作用効果は、すでに述べた通りであり、本発明に係る粒度測定システムによれば、搬送装置によって搬送される粉粒体の粒度をオンラインであっても容易に精度よく求めることができる。
以下、本発明の実施の形態について、詳細に説明する。図1は、本発明の粒度測定システムの概要構成を示している。粒度測定システムは、搬送装置によって搬送される粉粒体の粒度を測定することができる。以下、粉粒体の一例として石炭を挙げて説明する。
石炭Cは、ホッパー1からベルトコンベア2(搬送装置に相当)のベルト20上に落とされる。ベルト20上に落とされた石炭Cは、ベルトコンベア2により図1の左方向へと搬送されていく。ベルト20上の石炭Cは、均し手段3によって表面が平滑化された石炭層C1となる。
表面が平滑化された石炭層C1は、さらにベルトコンベア2で搬送されていく。レーザー変位計4は、移動する石炭層C1の表面の変位を所定のサンプリング周期で計測し、出力信号を出力する。出力された出力信号は、演算処理装置5に取り込まれる。演算処理装置5は、取り込んだ出力信号から変位データを算出するとともに、後述する数値化処理などを行なうことができる。
<粒度測定方法>
石炭層C1の表面の変位を計測した変位データ(以下、計測変位データ)の一例を図2に示す。図2では、計測変位データを時系列にプロットしている。ここでは、3mm篩下65%の石炭を用いた結果を示す。なお、本実施形態において、3mm篩下は、3mm篩を通過する粉体の質量から計算され、3mm未満の大きさの粉体の重量百分率で表される。すなわち、3mm篩下65%とは、3mm篩を通過した粉体の重量が全体の重量の65%であることを示している。
石炭層C1の表面の変位を計測した変位データ(以下、計測変位データ)の一例を図2に示す。図2では、計測変位データを時系列にプロットしている。ここでは、3mm篩下65%の石炭を用いた結果を示す。なお、本実施形態において、3mm篩下は、3mm篩を通過する粉体の質量から計算され、3mm未満の大きさの粉体の重量百分率で表される。すなわち、3mm篩下65%とは、3mm篩を通過した粉体の重量が全体の重量の65%であることを示している。
図2では、計測開始時におけるレーザー変位計4から石炭層C1の表面までの計測距離を0とし、この計測距離より短くなったときをマイナス、この計測距離より長くなったときをプラスとしている。また、図2には、石炭Cを載せない状態でのベルト20表面の変位も同時に記載している。
本実施形態で使用したレーザー変位計(キーエンス社製、LK−G150)は、20,50,100,200,500,1000μsecのサンプリング周期にて計測可能である。例えば、ベルトコンベア2の搬送速度が160m/minの場合、サンプリング周期を100μsecとすると約0.26mmの間隔で変位を計測可能である。このとき、粒径が3mmであれば、一つの粒子に対して約11点のデータをサンプリング可能となり、石炭粒子の表面プロフィールを十分に判別できる。仮にベルトコンベア2の搬送速度がもっと速い場合、サンプリング周期を短くすれば、一つの粒子に対するサンプリング数が少なくなることはない。反対にベルトコンベア2の搬送速度がもっと遅い場合には、サンプリング周期をもっと長くしてもよい。すなわち、サンプリング周期は、ベルトコンベアの搬送速度、必要とするサンプリング数、計測する粒径などによって適宜選択可能である。
図のように、計測された石炭層C1表面の計測変位データには、大きな上下変動と、小さな上下変動の2つの変動があることが分かる。一方、ベルト20表面のみの計測変位データの場合、大きな上下変動しか見られない。すなわち、石炭層C1表面の計測変位データは、ベルト20自体のうねりと、石炭層C1のベルト上堆積プロフィールとが合わさった大きな上下変動に、石炭粒子の大小による小さな上下変動が加わったものと考えられる。よって、石炭層C1表面全体の計測変位データから大きな上下変動を差し引き、小さな上下変動を抽出することで、石炭粒子の大小による変位が分かる。
そのため、初めに、石炭粒子の大小による小さな上下変動を全体の計測変位データから取り除くことで、大きな上下変動のみを抽出した。小さな上下変動を取り除くため、全体の計測変位データを数値化処理した。図3に、石炭層C1表面全体の計測変位データ、およびこれを数値化処理した変位データ(数値化処理変位データに相当)を合わせて示す。この実施形態では、数値化処理は、非線形のカーブフィッティングにより行なった。非線形のカーブフィッティングは、数学ソフト(ヒューリンクス社製、Flex Pro7)を用いて行なった。
非線形のカーブフィッティングは、計測変位データを例えば、
計算値i = p0+p1×exp(−X/p2)
の式に最小二乗法でフィッティングし、パラメータp0,p1,p2を導出することにより行なわれる。ここで、Xはデータの先頭位置からの番号である。
計算値i = p0+p1×exp(−X/p2)
の式に最小二乗法でフィッティングし、パラメータp0,p1,p2を導出することにより行なわれる。ここで、Xはデータの先頭位置からの番号である。
誤差分散は、計測変位データである計測値と上記の数値化処理によって求められた計算値との誤差を計算し、この誤差の二乗値の全総和を自由度で割ることで求められる。すなわち、誤差分散は、下記数式1で求められる。なお、自由度はX−1である。
誤差分散を求めるために用いるデータ数Xは、25点以上であって、好ましくは500点、より好ましくは2000点である。
例えば、本実施形態ではベルトコンベア2の搬送速度を24m/minとしており、最小のサンプリング周期である20μsecで計測すると、0.008mmのサンプリング間隔で変位を計測可能であり、最大のサンプリング周期である1000μsecで計測すると、0.4mmのサンプリング間隔で変位を計測可能である。最大のサンプリング間隔である0.4mmで計測する場合、25点未満のデータ数では、石炭が粉砕されて大きさ10mm未満の一粒のプロファイルは、計測できないことになるため、データ数は25点以上必要である。また、乾留工業の分野では、10〜50mmの大きさの石炭を取り扱うことが多いが、50mmの石炭の一粒をプロファイルするためには、0.1mmのサンプリング間隔で、500点のデータ数とすることが好ましい態様である。
一方、500点よりデータ数が多くなると、工場の操業では解析プログラムや制御系の機械への負荷が多くなってしまう。そのため、520mm=(0.26mm間隔×2000点)のデータであれば、32.5〜39kgの石炭を評価することになり、重量基準での粒度分布を評価するに匹敵する量を処理できる。
上記の3mm篩下65%の石炭のほか、3mm篩下75%、3mm篩下85%の石炭についても同様に石炭層C1表面の変位を計測し、上記の方法により誤差分散を算出した。
3mm篩下割合と誤差分散の関係について図4に示す。図のように、誤差分散が大きければ3mm篩下割合が小さく、誤差分散が小さければ3mm篩下割合が大きくなっており、誤差分散と3mm篩下割合には高い相関関係があることが分かる。これにより、計測した石炭層C1表面の計測変位データから誤差分散を求め、予め作成した誤差分散と3mm篩下割合との関係式から、この誤差分散に基づいて3mm篩下割合、すなわち粒度を推定することができる。
<別実施形態>
上記の実施形態では、粉粒体の一例として石炭を挙げて説明したが、本発明の粒度測定方法および粒度測定システムは、石炭以外の様々な粉粒体の粒度測定にも使用可能である。
上記の実施形態では、粉粒体の一例として石炭を挙げて説明したが、本発明の粒度測定方法および粒度測定システムは、石炭以外の様々な粉粒体の粒度測定にも使用可能である。
時系列にプロットした計測変位データを数値化処理する方法としては、上記の式を用いた非線形のカーブフィッティングに限定されない。
粒度と誤差分散の関係は、図4のような関係式を予め作成しておく方法のほか、誤差分散と対応する粒度の関係テーブルを予め作成しておく方法でもよい。
2 ベルトコンベア
3 均し手段
4 レーザー変位計
5 演算処理装置
20 ベルト
C 石炭
C1 石炭層
3 均し手段
4 レーザー変位計
5 演算処理装置
20 ベルト
C 石炭
C1 石炭層
Claims (4)
- 搬送装置によって搬送される粉粒体の粒度を測定する粒度測定方法であって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するステップと、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理するステップと、
数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出するステップと、を備える粒度測定方法。 - 予め作成しておいた誤差分散と粒度の関係から、前記算出された誤差分散に対応する粒度を求めるステップを備えることを特徴とする請求項1に記載の粒度測定方法。
- 前記数値化処理は、時系列にプロットした計測変位データを非線形のカーブフィッティングにより数値化することを特徴とする請求項1又は2に記載の粒度測定方法。
- 搬送装置によって搬送される粉粒体の粒度を測定する粒度測定システムであって、
搬送される粉粒体表面の変位を所定のサンプリング周期で計測するレーザー変位計と、
計測された粉粒体表面の計測変位データを時系列にプロットした後、数値化処理し、数値化処理された数値化処理変位データと前記計測変位データとの差から誤差分散を算出する演算処理装置と、を備える粒度測定システム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009274760A JP2011117806A (ja) | 2009-12-02 | 2009-12-02 | 粒度測定方法および粒度測定システム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009274760A JP2011117806A (ja) | 2009-12-02 | 2009-12-02 | 粒度測定方法および粒度測定システム |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2011117806A true JP2011117806A (ja) | 2011-06-16 |
Family
ID=44283314
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009274760A Pending JP2011117806A (ja) | 2009-12-02 | 2009-12-02 | 粒度測定方法および粒度測定システム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2011117806A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2020203255A1 (ja) * | 2019-04-02 | 2021-04-30 | Jfeスチール株式会社 | 粒度分布監視装置、粒度分布監視方法、コンピュータプログラム、炉、高炉、炉の制御方法、及び高炉操業方法 |
-
2009
- 2009-12-02 JP JP2009274760A patent/JP2011117806A/ja active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2020203255A1 (ja) * | 2019-04-02 | 2021-04-30 | Jfeスチール株式会社 | 粒度分布監視装置、粒度分布監視方法、コンピュータプログラム、炉、高炉、炉の制御方法、及び高炉操業方法 |
JP7020561B2 (ja) | 2019-04-02 | 2022-02-16 | Jfeスチール株式会社 | 粒度分布監視装置、粒度分布監視方法、コンピュータプログラム、炉、高炉、炉の制御方法、及び高炉操業方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3605064B1 (en) | Raw material particle size distribution measuring device, particle size distribution measuring method, and void ratio measuring device | |
JP6519034B2 (ja) | 粉率測定装置および粉率測定システム | |
Le et al. | Experimental analysis of powder layer quality as a function of feedstock and recoating strategies | |
Ren et al. | Determination of particle size distribution by multi-scale analysis of acoustic emission signals in gas-solid fluidized bed | |
Teffo et al. | Determination of the coefficients of restitution, static and rolling friction of Eskom-grade coal for discrete element modelling | |
Król et al. | Laboratory tests of idlers rotational resistance–selected issues | |
JP2011117806A (ja) | 粒度測定方法および粒度測定システム | |
Qin et al. | Study on physical and contact parameters of limestone by DEM | |
Grima et al. | Discrete element modelling: Trouble-shooting and optimisation tool for chute design | |
Kukreja et al. | Development of transfer chute design through Discrete Element Modelling for using Refused Derived fuel in Indian cement plants | |
Rahman et al. | Sensor-based control in eddy current separation of incinerator bottom ash | |
CN108627107A (zh) | 管道内沉积粉尘厚度监测装置及方法 | |
Wangchai et al. | The simulation of particle flow mechanisms in dustiness testers | |
CN108693072B (zh) | 料种辨识方法及装置 | |
CN113646446B (zh) | 粉率测定方法及装置 | |
CN206074113U (zh) | 斗轮机配煤称量修正系统 | |
Väyrynen et al. | Mass-flow estimation in mineral-processing applications | |
Van Laarhoven et al. | Toward a desktop attrition tester; validation with dilute phase pneumatic conveying | |
Žídek et al. | Effective use of DEM to design chain conveyor geometry | |
KR20180045943A (ko) | 로드셀 채널 분리를 이용한 동적 계량 시스템의 정확도 향상 장치 | |
JP6866240B2 (ja) | 改質材投入量算出方法および改質材投入量算出装置 | |
JP7427607B2 (ja) | 配合物の複数の画分の定量的割合を決定する方法及び装置 | |
JPH06167438A (ja) | 粒状性物体の粒度分布測定装置 | |
JP2018205271A (ja) | 改質材投入量算出方法および改質材投入量算出装置 | |
Heyduk | Machine vision monitoring and particle size feed analysis |