Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JPH0582893B2 - - Google Patents
[go: Go Back, main page]

JPH0582893B2 - - Google Patents

Info

Publication number
JPH0582893B2
JPH0582893B2 JP25610085A JP25610085A JPH0582893B2 JP H0582893 B2 JPH0582893 B2 JP H0582893B2 JP 25610085 A JP25610085 A JP 25610085A JP 25610085 A JP25610085 A JP 25610085A JP H0582893 B2 JPH0582893 B2 JP H0582893B2
Authority
JP
Japan
Prior art keywords
laser beam
light
particles
particle
intensity
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.)
Expired - Lifetime
Application number
JP25610085A
Other languages
Japanese (ja)
Other versions
JPS62116224A (en
Inventor
Juichi Ide
Shohei Noda
Masami Okada
Kazutomo Ootake
Takeo Hirai
Naokazu Kimura
Yoshio Tanaka
Shigehiko Ueda
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electric Power Development Co Ltd
Furukawa Electric Co Ltd
Mitsubishi Heavy Industries Ltd
Original Assignee
Electric Power Development Co Ltd
Furukawa Electric Co Ltd
Mitsubishi Heavy Industries Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Electric Power Development Co Ltd, Furukawa Electric Co Ltd, Mitsubishi Heavy Industries Ltd filed Critical Electric Power Development Co Ltd
Priority to JP25610085A priority Critical patent/JPS62116224A/en
Publication of JPS62116224A publication Critical patent/JPS62116224A/en
Publication of JPH0582893B2 publication Critical patent/JPH0582893B2/ja
Granted legal-status Critical Current

Links

Landscapes

  • Radiation Pyrometers (AREA)

Description

【発明の詳細な説明】[Detailed description of the invention]

〔産業上の利用分野〕 本発明は、微粉炭だきボイラ火炉の燃焼監視等
に適用される粒子温度分布測定方法に関する。 〔従来の技術〕 微粉炭だきボイラ火炉の燃焼監視等の目的で従
来より行なわれている粒子温度分布測定方法を第
3図に示す。 円筒炉1の中心軸上に位置する微粉炭燃焼火炎
2内の微粉炭粒子3から放射された光はレンズ4
を透過し、ハーフミラー5により2系統に分けら
れる。ハーフミラー5で反射した光はレンズ4を
透過した後、波長λ1のモノクロメータ6で波長λ1
の単色光のみが取出され、検出器7及び電位計8
により微粉炭粒子3が放射する波長λ1の光強度
I〓1が検出される。一方、ハーフミラー5を透過
した光はレンズ4を通り、波長λ1に近い波長λ2
モノクロメータ9に入り、ここで波長λ2の単色光
のみが取出され、検出器7及び電位計8により微
粉炭粒子3が放射する波長λ2の光強度I〓2が検出
される。 微粉炭粒子3から放射される波長λにおける輻
射能I〓は、波長λの輻射率をε〓、粒子の温度をTp
とすると、Planekの式により I〓=ε〓・C1/λ5・1/eC2/TPp−1 ……(1) ここで、C1,C2は定数で、 C1=3.2179×103(Kcalμm4/m2h)、 C2=1.43388×104(μm6K)であるので、
eC2/Tpは103以上のオーダーとなる。従つて(1)式は
次のように変形される。 I〓=ε〓・C1/λ5・eC2/Tp ……(2) 従つて、波長λ1,λ2での輻射率をε〓1,ε〓2とす
ると、 I〓1/I〓2= 〔ε〓1・λ1 -5/(ε〓2・λ2 -5)〕・ e−C2/TP・λ2−λ1/λ1 λ2 ……(3) 従つて、ε〓1=ε〓2とみなせる範囲の2波長λ1
λ2では ln(I〓1/I〓2)= −5ln(λ1/λ2)−C2/Tp・λ2−λ1/λ1λ2……(
4) 故に、増幅器10によりln(I〓1/I〓2)を算出す
れば、微粉炭燃焼火炎2中の粒子温度Tpが求ま
る。 〔発明が解決しようとする問題点〕 微粉炭燃焼火炎中の濃度分布、粒子温度分布が
不均一な場合、従来の測定方法では検出器にて検
出される粒子からの光強度は光路上に存在する粒
子全部から検出器まで到達した輻射量の和である
し、またレンズにより求めようとしている粒子に
焦点を合せてもその粒子と検出器との間に存在す
る他の粒子に邪魔されて粒子からの真の光強度が
検出できなくなる等、粒子間の影響を受けるた
め、粒子濃度を考慮した測定方法を考えない限り
光路上の任意位置での粒子温度を知ることは不可
能である。 本発明はこのような問題を解消するためになさ
れたもので、その目的は、不均一な分布を有する
粒子の浮遊場の濃度分布や、粒子温度分布も計測
することができ、微粉炭燃焼場の燃焼解析を通じ
て火炉の高効率運転及び最適設計等に有効な粒子
温度分布測定方法を提供することにある。 〔問題点を解決するための手段〕 以上の目的を達成するために、本発明の粒子温
度分布測定方法では、粒子から放射される2波長
の光強度の情報の他に、レーザ光を粒子の浮遊場
に照射して浮遊場を透過した光強度の情報を取入
れる。さらにこれら3つの情報を多方向について
求め、それらのデータをCT(Computer
Tomography)演算処理する。 〔作用〕 したがつて、本発明の測定方法によれば、レー
ザ光を粒子の浮遊場に照射して浮遊場におけるレ
ーザ光の透過減衰特性のデータから光路上におけ
る全体の粒子濃度の情報、並びに上記浮遊場から
の放射光のうち2波長を選択して光路上における
全体の粒子温度の情報を2つの情報を多方向に取
入れ、それらをCT演算処理することにより任意
位置での粒子温度が求められる。 〔実施例〕 本発明の方法の一実施例を第1図及び第2図に
もとづき説明する。 空気によつて搬送された微粉炭はバーナ21か
ら円筒炉22内へ吹込まれ、そこで燃焼火炎を形
成する。円筒炉22の中心軸に対して垂直な断面
23(−断面)内における炉外の任意の1点
に位置するArレーザ光発振器24から同断面内
における炉中心すなわちθ0方向に向つて、円筒炉
22の炉壁に設けた孔25を通して入射された
Arレーザ光27は途中粒子28に邪魔されて光
強度を減衰し、上記孔25と対向する位置の炉壁
に設けた孔26より炉外へ出る。炉外へ出たAr
レーザ光27の一部はレーザ光選択ハーフミラー
29にて反射され、レンズ30を通り、波長λa
(=4880Å)のモノクロメータ31、検出器32、
電位計33によりその光強度I〓aが検出される。
この光強度I〓aと、Arレーザ光発振器24で発生
する光強度I′〓aとの間には次の関係がある。 I〓a=I′〓a・e−(∫L 0kπ/4dp2Ndl) ……(5) ここに、Kは減衰係数、dpは粒子の直径、N
は粒子濃度、Lは光路長である。ここで、I′〓a
既知であるため、次式にて強度減衰率Fjを求め、
計算機37のメモリにストアする。 Fj=−ln{(I〓a/I′〓aj} =∫L Okπ/4dp2Ndl =∫L Ofldl= Σi fi・Δli ……(6) 次に、Arレーザ光発振器24とArレーザ光選
択ハーフミラー29とをγ方向にΔγだけ平行移
動し、前記と同様にしてArレーザ光27の透過
強度I〓aを測定し、(6)式によりFjを計算し、これ
を計算機37のメモリにストアする。この操作を
繰返すことによりθ0方向のFj分布が求まる。θ0
向のFj分布が求まると、前記一対のArレーザ光
発振器24とArレーザ光選択ハーフミラー29
をθ方向にΔθだけ回転させ、前記と同様にして
(θ0+Δθ)方向のFj分布を求め、これらの操作を
(θ0+360°−Δθ)方向まで繰返し、あらゆる方向
のFj分布を求める。(6)式中のΔliはメツシユの大
きさで既知であるため、計算機37でCT演算処
理することにより、上記−断面内の光吸収係
数fi分布が求まる。 ここで、Fjの値を計算機37にストアすると
き、同一光路上に存在する全ての粒子からの輻射
光を、さらにレーザ光選択ハーフミラー29とは
別のハーフミラー34により反射光と透過光に分
け、反射光はレンズ30を通り、波長λ1(燃焼火
炎からの輻射光の極めて強い波長域である5890Å
前後を除く黄色から赤色レンジの波長)のモノク
ロメータ35、検出器32、電位計33にて波長
λ1の光強度I〓1を検出し、またハーフミラー34
からの透過光はレンズ30を通り、λ1に近い波長
λ2のモノクロメータ36、検出器32、電位計3
3にて波長λ2の光強度I〓2を検出する。 今、炉壁26よりlなる距離にある粒子からの
輻射エネルギI′〓が光路上に沿う他の粒子による減
衰のため検出器に到達するときのエネルギI〓は、 I〓=I′〓・e−(∫L lkπ/4dp2Ndl) =I′λ・e−(∫L 0fldl ……(7) 従つて、光路長全てに対して積分することによ
り(I〓1j,(I〓2jは次のようになる。
[Industrial Application Field] The present invention relates to a particle temperature distribution measuring method applied to combustion monitoring of a pulverized coal-fired boiler furnace. [Prior Art] Fig. 3 shows a particle temperature distribution measurement method that has been conventionally used for the purpose of monitoring combustion in pulverized coal-fired boiler furnaces. The light emitted from the pulverized coal particles 3 in the pulverized coal combustion flame 2 located on the central axis of the cylindrical furnace 1 is transmitted through a lens 4.
is transmitted and divided into two systems by a half mirror 5. After the light reflected by the half mirror 5 passes through the lens 4, it is converted to a wavelength λ 1 by a monochromator 6 with a wavelength λ 1.
Only the monochromatic light of is extracted, and the detector 7 and the electrometer 8
The light intensity of wavelength λ 1 emitted by pulverized coal particles 3 is
I〓 1 is detected. On the other hand, the light transmitted through the half mirror 5 passes through the lens 4 and enters the monochromator 9 with a wavelength λ 2 close to the wavelength λ 1 , where only the monochromatic light with the wavelength λ 2 is extracted, and the detector 7 and the electrometer 8 The light intensity I〓 2 of wavelength λ 2 emitted by the pulverized coal particles 3 is detected. The radiation intensity I〓 at the wavelength λ emitted from the pulverized coal particles 3 is given by the emissivity at the wavelength λ being ε〓, and the temperature of the particle being T p
Then, according to Planek's formula, I〓=ε〓・C 15・1/e C2/TP p−1 ……(1) Here, C 1 and C 2 are constants, and C 1 = 3.2179 ×10 3 (Kcalμm 4 /m 2 h), C 2 = 1.43388×10 4 (μm 6 K),
e C2/Tp is on the order of 10 3 or more. Therefore, equation (1) is transformed as follows. I〓=ε〓・C 15・e C2/Tp ...(2) Therefore, if the emissivity at wavelengths λ 1 and λ 2 are ε〓 1 and ε〓 2 , then I〓 1 / I〓 2 = [ε〓 1・λ 1 -5 / (ε〓 2・λ 2 -5 )]・e−C 2 /TP・λ 2 −λ 11 λ 2 ...(3) Therefore , 2 wavelengths λ 1 in the range that can be considered as ε〓 1 = ε〓 2 ,
At λ 2 , ln(I〓 1 /I〓 2 ) = −5ln(λ 12 )−C 2 /T p・λ 2 −λ 11 λ 2 ...(
4) Therefore, by calculating ln (I〓 1 /I〓 2 ) using the amplifier 10, the particle temperature T p in the pulverized coal combustion flame 2 can be found. [Problem to be solved by the invention] When the concentration distribution and particle temperature distribution in the pulverized coal combustion flame are uneven, in the conventional measurement method, the light intensity from the particles detected by the detector is on the optical path. It is the sum of the amount of radiation that has reached the detector from all the particles that are present, and even if you focus on the particle you are trying to find with a lens, other particles that exist between that particle and the detector will interfere with the radiation. It is impossible to know the particle temperature at any position on the optical path unless a measurement method that takes particle concentration into account is considered because the true light intensity from the particle cannot be detected. The present invention has been made to solve these problems, and its purpose is to be able to measure the concentration distribution of particles floating in a non-uniform distribution field and the particle temperature distribution, and to measure the particle temperature distribution. The purpose of this study is to provide a particle temperature distribution measurement method that is effective for high-efficiency operation and optimal design of furnaces through combustion analysis. [Means for Solving the Problems] In order to achieve the above object, in the particle temperature distribution measuring method of the present invention, in addition to information on the light intensity of two wavelengths emitted from particles, laser light is Information on the intensity of light that irradiates the floating field and passes through the floating field is taken in. Furthermore, these three pieces of information are obtained in multiple directions, and these data are
Tomography) calculation processing. [Function] Therefore, according to the measurement method of the present invention, a laser beam is irradiated onto a floating field of particles, and information on the overall particle concentration on the optical path can be obtained from data on the transmission attenuation characteristics of the laser beam in the floating field. Select two wavelengths of the emitted light from the floating field, take in the two pieces of information on the entire particle temperature on the optical path in multiple directions, and calculate the particle temperature at any position by performing CT calculation processing on them. It will be done. [Example] An example of the method of the present invention will be described based on FIGS. 1 and 2. The air-borne pulverized coal is blown from the burner 21 into the cylindrical furnace 22, where it forms a combustion flame. From the Ar laser beam oscillator 24 located at an arbitrary point outside the furnace in the cross section 23 (- cross section) perpendicular to the central axis of the cylindrical furnace 22, toward the furnace center in the same cross section, that is, in the θ 0 direction, the cylindrical Injected through the hole 25 provided in the furnace wall of the furnace 22
The Ar laser beam 27 is obstructed by particles 28 on the way and its light intensity is attenuated, and then exits from the furnace through a hole 26 provided in the furnace wall at a position opposite to the hole 25 . Ar that went out of the reactor
A part of the laser beam 27 is reflected by the laser beam selection half mirror 29, passes through the lens 30, and has a wavelength λ a
(=4880Å) monochromator 31, detector 32,
The light intensity I〓a is detected by the electrometer 33.
The following relationship exists between this light intensity I〓a and the light intensity I′〓a generated by the Ar laser beam oscillator 24. I〓 a =I′〓 a・e−(∫ L 0 kπ/4dp 2 Ndl) ……(5) Here, K is the attenuation coefficient, dp is the diameter of the particle, and N
is the particle concentration and L is the optical path length. Here, since I′〓 a is known, the intensity attenuation rate F j is calculated using the following formula,
It is stored in the memory of the computer 37. F j = −ln {(I〓 a /I′〓 a ) j } =∫ L O kπ/4dp 2 Ndl =∫ L O f l dl= Σ i fi・Δli ...(6) Next, Ar laser The optical oscillator 24 and the Ar laser beam selection half mirror 29 are translated in the γ direction by Δγ, the transmitted intensity I〓 a of the Ar laser beam 27 is measured in the same manner as above, and F j is calculated using equation (6). and stores it in the memory of the computer 37. By repeating this operation, the F j distribution in the θ 0 direction is determined. When the F j distribution in the θ 0 direction is determined, the pair of Ar laser beam oscillators 24 and the Ar laser beam selection half mirror 29
Rotate by Δθ in the θ direction, obtain the F j distribution in the (θ 0 + Δθ) direction in the same way as above, and repeat these operations up to the (θ 0 + 360° − Δθ) direction to obtain the F j distribution in all directions. demand. Since Δl i in equation (6) is known from the size of the mesh, the light absorption coefficient f i distribution within the above-mentioned cross section is determined by CT calculation processing by the computer 37. Here, when storing the value of F j in the computer 37, the radiation light from all the particles existing on the same optical path is further reflected and transmitted by a half mirror 34 that is different from the laser beam selection half mirror 29. The reflected light passes through the lens 30 and has a wavelength of λ 1 (5890 Å, which is the extremely strong wavelength range of the radiant light from the combustion flame).
A monochromator 35, a detector 32, and an electrometer 33 detect the light intensity I〓 1 with a wavelength λ 1 (wavelengths in the yellow to red range excluding the front and rear), and a half mirror 34
The transmitted light passes through a lens 30, and a monochromator 36 with a wavelength λ 2 close to λ 1 , a detector 32, and an electrometer 3
3, the light intensity I〓 2 of wavelength λ 2 is detected. Now, when the radiant energy I'〓 from a particle at a distance l from the furnace wall 26 reaches the detector due to attenuation by other particles along the optical path, the energy I〓 is I〓=I′〓・e−(∫ L l kπ/4dp 2 Ndl) = I′λ・e−(∫ L 0 f l dl ……(7) Therefore, by integrating over the entire optical path length, (I〓 1 ) j , (I〓 2 ) j is as follows.

【化】[ka]

【化】 そこで、Arレーザ光発振器24とArレーザ光
検出器を用いてFjを求めるとき同時に(I〓1j
(I〓2jも求め、あらゆる方向についての(I〓1j
(I〓2jの分布を求める。(8)式、(9)式中のΔliはメ

シユの大きさであり、また前にFjに関するCT演
算処理にてfiを求めておけば、
[C] Therefore, when finding F j using the Ar laser light oscillator 24 and the Ar laser light detector, (I〓 1 ) j ,
(I〓 2 ) j is also found, and (I〓 1 ) j ,
(I〓 2 ) Find the distribution of j . Δl i in equations (8) and (9) is the size of the mesh, and if f i is calculated by CT calculation processing regarding F j beforehand, then

〔発明の効果〕〔Effect of the invention〕

本発明により、従来の測定方法では求めること
ができなかつた不均一な分布を有する粒子の濃度
分布や、粒子温度分布も計測することができ、微
粉炭燃焼場の燃焼解析を通じて火炉の高効率運転
及び最適設計等に大きな効果をもたらすことがで
きる。
With the present invention, it is possible to measure particle concentration distribution and particle temperature distribution with non-uniform distribution, which could not be determined using conventional measurement methods, and to achieve high efficiency operation of furnaces through combustion analysis of pulverized coal combustion fields. This can have a great effect on optimal design and the like.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図は本発明の方法の一実施例を示す円筒炉
部分の縦断面図、第2図は第1図の−線に沿
う断面と全工程の概略を合せて示す図、第3図は
従来の方法を説明する概略説明図である。 22……円筒炉、24……Arレーザ発振器、
27……Arレーザ光、28……粒子、31,3
5,36……モノクロメータ、33……電位計、
37……計算機。
Fig. 1 is a vertical cross-sectional view of a cylindrical furnace section showing an embodiment of the method of the present invention, Fig. 2 is a cross-sectional view along the - line in Fig. 1 and an outline of the entire process, and Fig. 3 is FIG. 2 is a schematic explanatory diagram illustrating a conventional method. 22...Cylindrical furnace, 24...Ar laser oscillator,
27...Ar laser beam, 28...Particle, 31,3
5, 36... Monochromator, 33... Electrometer,
37...Calculator.

Claims (1)

【特許請求の範囲】[Claims] 1 粒子の浮遊している浮遊場を横断するように
レーザビームを多方向から照射し、このときに得
られたレーザビームの強度減衰率データを基に
CT演算処理を行なつて上記レーザビーム相互が
交わる位置における光吸収係数をそれぞれ算出す
るとともに、レーザビーム通路上の粒子から輻射
される輻射光のうちの2波長以上の輻射光強度を
ビーム軸延長線上においてそれぞれ検出し、検出
された輻射光強度データと前記光吸収係数データ
とを基にCT演算処理を行つて前記浮遊場の温度
分布データを得るようにしたことを特徴とする粒
子温度分布測定方法。
1 A laser beam is irradiated from multiple directions across a floating field where particles are suspended, and based on the laser beam intensity attenuation rate data obtained at this time,
CT calculation processing is performed to calculate the light absorption coefficients at the positions where the laser beams intersect, and the beam axis is extended by the intensity of the radiant light of two or more wavelengths of the radiant light emitted from the particles on the laser beam path. Particle temperature distribution measurement characterized in that temperature distribution data of the floating field is obtained by detecting each on the line and performing CT calculation processing based on the detected radiation light intensity data and the light absorption coefficient data. Method.
JP25610085A 1985-11-15 1985-11-15 Method for measuring particle temperature distribution Granted JPS62116224A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP25610085A JPS62116224A (en) 1985-11-15 1985-11-15 Method for measuring particle temperature distribution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP25610085A JPS62116224A (en) 1985-11-15 1985-11-15 Method for measuring particle temperature distribution

Publications (2)

Publication Number Publication Date
JPS62116224A JPS62116224A (en) 1987-05-27
JPH0582893B2 true JPH0582893B2 (en) 1993-11-22

Family

ID=17287888

Family Applications (1)

Application Number Title Priority Date Filing Date
JP25610085A Granted JPS62116224A (en) 1985-11-15 1985-11-15 Method for measuring particle temperature distribution

Country Status (1)

Country Link
JP (1) JPS62116224A (en)

Also Published As

Publication number Publication date
JPS62116224A (en) 1987-05-27

Similar Documents

Publication Publication Date Title
DE19611290C2 (en) Gas sensor
CN113820035B (en) Femtosecond laser filament remote non-contact temperature measurement device and measurement method
JPS62222144A (en) Apparatus for measuring particle size
EP0167272A2 (en) Particle size measuring apparatus
JPS63241336A (en) Particle size measuring apparatus
Shi et al. Reconstruction of 3D flame temperature and absorption coefficient field by the hybrid light-field imaging and laser extinction technique
JPH0263181B2 (en)
Lee et al. A new sensor for detection of coolant leakage in nuclear power plants using off-axis integrated cavity output spectroscopy
Tobiasson et al. An optical method for the measurement of combustion gas temperature in particle laden flows
JPH0582893B2 (en)
Linjewile et al. Optical probe measurements of the temperature of burning particles in fluidized beds
CN106769735B (en) A dust concentration measuring device
JPS62116225A (en) Method for measuring temperature distribution of gas
Shaw et al. An Instrument for In-Situ Temperature and Emissivity Measurements (INSITE) on Boiler Ash Deposits
CN104181128B (en) Trnaslucent materials based on time-correlated single photon counting t radiation property measurement method
Lewis et al. Insights from a new method providing single-shot, planar measurement of gas-phase temperature in particle-laden flows under high-flux radiation
CN105258808A (en) A method for determining the direction of a partially-coherent laguerre-gaussian light beam vortex
Woo et al. Measurement of gas temperature profile using spectral intensity from CO2 4.3 μm band
CN103499814B (en) A kind of high precision Doppler lidar Frequency Locking system
Bradley et al. Measurement of temperature PDFS in turbulent flames by the CARS technique
JPS6325522A (en) Temperature measurement method using anti-Stox Raman spectroscopy
JPS61240146A (en) Method for analyzing composition of article to be measured by x-rays
Kim et al. Measurement of gas temperature distributions in a test furnace using spectral remote sensing
JPH0310128A (en) Method for simultaneously measuring temperature and emissivity in high temperature furnace
Bates et al. Infrared monitoring of combustion