• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    A Novel Sensing Imaging Equipment Under Extremely Dim Light for Blast Furnace Burden Surface:Starlight High-Temperature Industrial Endoscope

    2024-04-15 09:36:50ZhipengChenXinyiWangWeihuaGuiJilinZhuChunhuaYangandZhaohuiJiang
    IEEE/CAA Journal of Automatica Sinica 2024年4期

    Zhipeng Chen , Xinyi Wang , Weihua Gui ,,, Jilin Zhu ,Chunhua Yang ,,, and Zhaohui Jiang ,,

    Abstract—Blast furnace (BF) burden surface contains the most abundant, intuitive and credible smelting information and acquiring high-definition and high-brightness optical images of which is essential to realize precise material charging control,optimize gas flow distribution and improve ironmaking efficiency.It has been challengeable to obtain high-quality optical burden surface images under high-temperature, high-dust, and extremelydim (less than 0.001 Lux) environment.Based on a novel endoscopic sensing detection idea, a reverse telephoto structure starlight imaging system with large field of view and large aperture is designed.Combined with a water-air dual cooling intelligent self-maintenance protection device and the imaging system, a starlight high-temperature industrial endoscope is developed to obtain clear optical burden surface images stably under the harsh environment.Based on an endoscope imaging area model, a material flow trajectory model and a gas-dust coupling distribution model, an optimal installation position and posture configuration method for the endoscope is proposed, which maximizes the effective imaging area and ensures large-area, safe and stable imaging of the device in a confined space.Industrial experiments and applications indicate that the proposed method obtains clear and reliable large-area optical burden surface images and reveals new BF conditions, providing key data support for green iron smelting.

    I.INTRODUCTION

    ENERGY saving and emission reduction are the most important themes in the steel industry development [1].As the most critical equipment in steel smelting, the energy consumption and cost of the blast furnace (BF), accounts for about 60%-70% of the steel manufacturing process, and its CO2emission accounts for more than 90% of that, which has significant high-carbon characteristics [2].Therefore, controlling energy consumption of BF is the key to realize green production in the steel industry [3].

    The BF top is the core window for the operation and fine control of the blast furnace.As the most critical parameter for furnace top control, the burden surface topography information can effectively reflect the operation conditions and distribution of the flow field in BF top, which are important references for guiding the refinement of burden distribution,improving gas flow and dust distribution, ensuring smooth operation of BF and saving energy [4].However, the complex environment inside the BF severely limits the detection range and performance of traditional equipment, making it difficult to obtain clear and bright burden surface images, specifically in the following aspects: First, to prevent the leakage of highpurity toxic gases, BF is enclosed without light, and its internal burden layer has a strong light-absorbing effect, making the illumination inside the BF top reach a starlight level of 0.001 Lux, so it is difficult to acquire bright images of burden surface under this extremely dim light environment [5]; Second, the BF top temperature typically ranges from 200 °C to 400 °C.Under abnormal operating conditions, it even reaches as high as 800 °C to 1000 °C.In such a high-temperature environment, conventional detection equipment faces challenges in maintaining long-term stable operation, making it extremely difficult to perform close-range burden surface detection [6]; Finally, there is high dust concentration within the BF top ranging from approximately 40 to 100 g/m3.Under the influence of pressure difference, the suspended dust in the BF top space forms a layer of dust.This dust layer not only erodes and damages the equipment shell but also blocks the detection hole and interferes the detection accuracy, severely limiting the value of detection equipment in BF top [7], [8].

    Therefore, traditional detection equipment in BF top has given up the direct detection of the most abundant, intuitive,and reliable optical image information of the BF burden surface.Instead, utilize the mechanical probe, industrial radar,laser scanner and infrared camera to indirectly estimate the burden surface morphology, gas flow field, and dust field status based on limited point-level information or flow field information [9], [10].These methods commonly suffer from low detection accuracy and lack of furnace condition information, making it difficult to achieve precise charging control and leading to frequent abnormal operating conditions.

    To obtain high-precision and real-time morphological imaging information of the BF top burden surface, a starlight hightemperature industrial endoscope (SHIE) detection method is proposed.This method consists of three parts: the design of a starlight optical imaging system with large depth of field and wide field of view, development of an intelligent self-maintenance protection device with water-air dual cooling system and the optimal installation posture configuration of SHIE.Firstly, to improve the imaging quality of endoscope in starlight environment, a starlight optical imaging system with large depth of field and wide field of view is designed.Specifically, a reverse telephoto structure and the high-transmission Hopkins relay lens are employed to increase the field of view and minimize light loss of the optical system, enabling the endoscope to have an extremely wide imaging range with large field of view characteristics, achieving starlight-level imaging performance.Secondly, to ensure long-term stable operation of SHIE in the harsh environment of BF top, the intelligent self-maintenance protection device with water-air dual cooling system is developed.Thirdly, to capture the complete burden surface morphology information, the maximum burden-surface imaging optimization model for the endoscope is established.Based on the material flow trajectory model and the dust distribution model in BF top, the constraints for the model are determined.By utilizing the particle swarm optimization (PSO) algorithm, the optimal installation scheme for the endoscope is obtained, which enables the endoscope to capture the complete burden surface images while avoiding the material flow and high-concentration dust inside the BF top, providing sufficient smelting status information for the onsite operation of BF.

    To summarize, the main contributions of this paper are as follows:

    1) A starlight optical imaging system with large depth of field and wide field of view is designed to obtain high-definition and high-brightness images of whole burden surface.

    2) An intelligent self-maintenance protection device with water-air dual cooling system is developed to ensure the safe operation of equipment.

    3) The optimal installation posture configuration of the equipment is obtained, which enables the endoscope to capture the complete burden surface morphology image and to avoid the material flow and high-concentration dust inside the furnace.

    4) Both the qualitative and quantitative results on the captured burden surface images demonstrate the superiority of our approach over existing methods.

    The rest of this paper is organized as follows: Related works are discussed in Section II.Section III describes the principle and structure of the developed industrial endoscope; Section IV introduces the installation scheme of the endoscope; Section V discusses the simulation and field test results; Finally,the work is concluded in Section VI.

    II.RELATED WORK

    This section gives a brief introduction of existing methods on the detection of BF burden surface information.There are two types of methods: the mathematical methods and the physical methods.

    A. The Mathematical Methods

    The mathematical methods utilize detection equipment to obtain the height information of each point on the burden surface and then adopt surface fitting methods to indirectly reconstruct the burden surface morphology.These methods include mechanical stock rod measurement, indirect judgment,laser scanning, and radar imaging.

    The mechanical stock rod measurement method uses a mechanical probe between the inner walls of the BF throat to get 2- or 3-points’ height information of the BF stockline [11],which suffers from long-time measurement and limited information, and cannot measure continuously.The indirect judgment method combines the worker experience with the temperature data of the burden surface measured by a cross-temperature measurement device to estimate the approximate distribution of the burden surface.This method is subjective and lacks a quantitative description of the burden surface [12].The laser scanning method utilizes cameras to measure the twodimensional coordinates of the laser spot projected on the burden surface, thereby reconstructing the appearance of the burden surface [13].However, under the influence of dust and water vapor, the accuracy of this method is reduced and the laser head is easily scabbed and blocked.For the radar imaging method, using the phase difference between the transmitted signal and the received signal of the mechanical swing radar, the height information on the burden surface is measured point by point to reconstruct the shape of the burden surface.This method is greatly affected by the dust in the furnace, and multiple detection holes must be opened on the BF top, which easily affects the airtightness of BF and brings safety hazards to the steel production [14], [15].

    To sum up, the mathematical methods can only obtain limited point-level height information and indirectly estimate the burden surface morphology, which are unable to provide accurate and reliable burden surface morphology information.

    B. The Physical Methods

    The physical methods include the video capturing method,the infrared thermal imaging method and the high temperature industrial endoscope method, which directly obtain image information of the BF burden surface from the detection equipment.

    The video capturing method collects the burden surface information through the camera, but due to the dark and dusty environment, the expected image cannot be captured [16].The infrared thermal imaging method estimates the burden surface shape by obtaining the temperature distribution of the burden surface through infrared images.This method cannot accurately obtain the burden surface shape and is easily influenced by the high-temperature gas [17].The high temperature industrial endoscope method is a pioneering idea of us to use endoscopic optical imaging to detect the burden surface in BF top.The high-temperature industrial endoscope got far into the hinterland of the BF top, placing high-temperature-resistant optical lens inside the BF through the guiding fibers with the imaging chip outside.This method can obtain burden surface images, breaking the “black box” state on the BF top.However, the obtained burden surface images suffer from small imaging field of view, insufficient depth of field, low signalto-noise ratio and limited information in the starlight environment, therefore, the equipment requires to be optimized[18]-[20].

    Fig.1.Structure and prototype of starlight high-temperature industrial endoscope with large depth of field and wide field of view.

    To intuitively reflect the difference between different burden surface detection methods, data detected by different detection methods are compared in detail under different evaluation indicators.The comparison results are listed in Appendix A.As shown in Appendix A, the mechanical probe,the mechanical swing radar and the laser scanner can only obtain the height of material level at limited points on the burden surface, and indirectly estimate the trend of burden surface topography; The resolution of the images obtained by infrared camera is low, so it is impossible to see the burden surface morphology clearly; The images captured by the hightemperature industrial endoscope can observe the outline of the burden surface, but cannot reflect the specific details of the coke particles; While SHIE can directly obtain high-definition images of the topography of the burden surface, which contains higher quality data and larger amount of information.Thereby the proposed method can provide intuitive information on the distribution of the burden surface for the onsite workers.

    III.THEORY AND STRUCTURE OF THE STARLIGHT HIGH-TEMPERATURE INDUSTRIAL ENDOSCOPE

    The SHIE is the basis of the proposed burden surface detection method.With the purpose of obtaining clear and bright burden surface images in the harsh environment of BF, the developed equipment is mainly composed of two functional modules: a starlight optical imaging system with large depth of field and wide field of view and an intelligent self-maintenance protection device with water-air dual cooling system.The endoscope structure is shown in Fig.1(a), and its prototype is shown in Fig.1(b).The following section will introduce the specific principle and structure of each module.

    A. Starlight Imaging System With Large Depth of Field and Wide Field of View

    For the starlight imaging system, its performance indices are designed according to the relevant data of the BF, and based on the optical principle, the optical path and structure are developed.Then the imaging quality of the designed optical system is evaluated, to ensure that it satisfies the performance indicators.The specific steps are as follows.

    1) Analysis of Design Indices of Starlight Imaging System

    To ensure the BF airtightness, the installation aperture of the endoscope is limited, so the outer diameter of the starlight imaging system is set to 24 mm through the field test.The field of view 2 ω of the starlight imaging system is set to 103°according to technological conditions, so that the optical system can obtain the image of the whole burden surface and avoid the interference of dust and material flow on the BF top.Moreover, the starlight imaging system requires a large depth of field.The depth of field is calculated by

    where ?L1and ?L2are the front and rear depths of field of the starlight imaging system respectively,Fis the aperture value,δis the diameter of the permissible circle of confusion,fis the focal length of the system,Lis the distance between the focal object surface and the system.According to the definition of large depth of field, the focus is set at the hyperfocal distance point, that is, the front depth of field isL/2, and the rear depth of field is ∞, thenf2=FδL.Based on the latest data of high temperature industrial endoscope, the shortest object distance between the burden surface and the industrial endoscope is about 0.77 m, so ?L1=0.77 m, thenLis 1.54 m.To make the imaging results of the starlight optical imaging system clearer, the diameter of the permissible circle of confusion is equal to the pixel sizeIdand the diameter of Airy diskAd, thus

    where λ ( = 0.5876 μm) denotes the primary wavelength,Dis the diameter of the entrance pupil of the optical imaging system,Idis the pixel size, and the P2160 CMOS chip of 9.7 μm is used here.Then, substituting (2) into (1) obtains

    Through the above analysis, the design indices of the starlight imaging system are listed in Table I.

    TABLE I DESIGN INDICES OF OPTICAL IMAGING SYSTEM

    2) Structure of Starlight Imaging System

    With the above design indices, the designed starlight imaging system is composed of an objective lens group, three relay lens groups and a zoom lens group.

    a) Structure and principle of the objective lens group

    The objective lens group is the imaging element of the starlight imaging system.To satisfy the large field of view and short focal length goals in the design indices, the reverse telephoto objective structure is adopted as shown in Fig.2, which includes six lenses.The first and second lenses form the front lens groupL1, which has the negative power φ1and has capability of large field of view.The third to sixth lenses form the rear lens groupL2, which has the positive power φ2and larger relative aperture, and can realize strong light transmission ability.As illustrated in the optical path structure in Fig.2, the off-axis large-angle α beam enters the front lens group and diverges, and the angle decreases to α′after passing the rear lens group, then the image is formed on the focal plane, which is the Principle of large field imaging.

    b)Structure and principle of the relay lens group

    To realize long-distance transmitting the acquired image to the back-end image processing unit without damage, the developed system adopts three identical relay lens groups.As displayed in Fig.2, the structure of each lens group adopts two symmetrical Hopkins rod lenses, which is cemented by a thick biconvex lens and two symmetrical negative lenses.The Hopkins rod lens reduces the Spherical aberration, vertical axis aberration and Petzval sum.And its refractive index is greater than that of the air, thus the light energy transfer efficiency is better.

    The optical path structure of the relay lens group is shown in Fig.2, the focal length of the front rod lens is the same as the focal lengthf′of the objective lens group, so that the light emitted by the image plane of the objective lens group becomes parallel light after passing through the front rod lens.After the parallel light passes through the rear rod lens, it is focused on the focal plane again, so as to proportionally transmit the image.

    c) Structure and principle of the zoom lens group

    Since the BF stockline keeps changing in the whole distribution cycle, the distance between the burden surface and the starlight imaging system also changes, so the system utilizes a zoom lens group to realize the automatic focusing during the shooting process and guarantee the real-time clarity of the video.The zoom lens group consists of a front fixed lens, a movable lens, and a rear fixed lens as illustrated in Fig.2.According to Gaussian optics theory, the optical power of an optical system composed of two components is

    where φ1, φ2are the optical powers of the two components respectively and cannot be changed.Therefore, the zoom lens group moves the lens through the stepping motor and changes the component intervaldto realize the change of the system’s focal length.

    3)Design and Evaluation of Optical System

    With the above design indices and structural parameters, the starlight imaging system with large depth of field and wide field of view is designed by using the Zemax software, and iterative optimization is performed.The structure of the system is shown in Fig.3.To evaluate the imaging performance of the optical imaging system, simulation tests were conducted on the system, as depicted in Fig.4.Fig.4(a) illustrates the modulation transfer function (MTF) performance chart of the designed optical imaging system.At the optical system resolutionN=1/(2Nd)=50 lp/mm, the MTF values in both the tangential and sagittal directions across the entire field of view are greater than 0.25, approaching the diffraction limit, thus meeting the imaging quality requirements.Fig.4(b) exhibits the spot diagrams of the optical system at different field of views, with the maximum root mean square radius for each field of view measuring 9.910 μm, which is smaller than the pixel size of 10 μm.This indicates that the optical system effectively controls the aberrations across different field of views.Fig.4(c) displays the field curvature and distortion curves of the optical system.The field curvature in both the tangential and sagittal directions within the entire field of view is less than 0.25 mm, and the maximum distortion is -40% of barrel distortion.In a wide-field optical system, distortion is inevitable, but it does not affect the imaging quality and can be corrected during the image processing stage.Fig.4(d) shows the wave aberration of the optical system, with the maximum wave aberration being less than 0.0418 wavelengths, which is far below 1/4 wavelength.Thus,according to the Rayleigh criterion, the imaging quality of the optical system is close to that of an ideal system.

    Compared with other wide-field objective optical systems,as shown in Table II, it can be observed that the designed optical system has a larger field of view and can obtain a wider range of images.And the system pixel size is smaller, the system resolution is higher, and the captured image can reflect more details of objects.The aperture valueFis large, so the system has a large amount of light input, and can image in a dark environment.Besides, the ratio of wave aberration to wavelength is small, thus the imaging quality is relatively high, which can satisfy the clear imaging in starlight illumination environment.

    Fig.2.Optical path and lens structure of starlight-level optical imaging system with large depth of field and wide field of view.

    Fig.3.Structure of starlight optical imaging system.

    Fig.4.The performance indicators of the starlight optical imaging system.

    B. Intelligent Self-Maintenance Protection Device With Waterair Dual Cooling System and the Overall Operation Flow

    To ensure long-term stable operation of industrial endoscopes in the harsh environment, the intelligent self-maintenance protection device with water-air dual cooling system is designed, and its structure is shown in Fig.5.The front and rear sleeve shells, the fixed sleeve ring, the water-cooling circuit and the air-cooling circuit constitute the water-air dual cooling protection device.The intelligent self-maintenance device consists of the cleaning nozzles, actuator, shading cover.

    The overall workflow of the SHIE is as follows: Firstly, the water-air dual cooling protection device utilizes cooling water and industrial nitrogen to cool the internal components and increase the internal pressure of the equipment, preventing dust in BF top from entering the device.Secondly, the starlight imaging system derives the obtained image of burden surface to the low temperature area at the back end of the equipment, and the charge-coupled device (CCD) chip image digitally with the imaging drive circuit.Finally, the video signal line interface transmits the acquired digital video information.In addition, when the endoscope lens becomes blurry, the intelligent self-maintenance device will control the actuator to rotate the shading cover, opening the light-passing hole, and automatically sprays high-pressure water to clean the dust scab at the front end of the lens, restoring clear imaging and ensuring the long-term stable operation of the equipment.

    IV.OPTIMAL CONFIGURATION OF INSTALLATION POSITION AND POSTURE OF THE INDUSTRIAL ENDOSCOPE

    The installation posture of the industrial endoscope is essential to be designed.The installation requires to overcome three major challenges: First, the equipment must be installed in the no-material-flow area on the BF top to prevent the material flow particles from causing damage to the endoscope; Second,the equipment should be settled in the low dust area to avoid the interaction between dust and gas flow.Third, the vertical distance between the starlight imaging system and the standard stockline of the BF must ensure that the imaging area of the equipment is large enough and the information of burden surface is complete.

    Based on these, experiments were carried out on 2650 m3BF: Firstly, based on the obtained image area model of burden surface, a maximum burden-surface imaging optimiza-tion model of the industrial endoscope is established.Secondly, in order to obtain a safe installation region with few material flow and low dust on the top of the BF, the trajectory model of material flow and the coupling distributionk-Sε-upmodel of dust and gas are established.Finally,based on the field data, the optimal installation position and posture of the SHIE are calculated.

    TABLE II COMPARISON OF DESIGN INDICES BETWEEN DIFFERENT OPTICAL SYSTEMS

    Fig.5.Intelligent self-maintenance protection device.

    A. Maximum Burden-Surface Imaging Optimization Model

    To obtain the complete detail information of the burden surface, the maximum imaging area model of the SHIE should be constructed.Since the imaging area is determined by the image ellipse of SHIE, it is necessary to calculate the cone equation of the field of view for the endoscope first, and then calculate the device’s imaging ellipse equation to acquire the imaging area model.The model is shown in Fig.6(a).The field-of-view cone equation can be regarded as the generatrixABof the field-of-view coneA-BO′Crotates once around the axis of rotationDA.The height of the burden surfaceqis the distance from the burden level to the standard stockline,and its coordinates are generally between -0.5 m and 0.7 m relative to the standard stockline, that isq∈[-0.5,0.7].Using the parallel latitude circle method, the cone equation of the field of view can be calculated as

    Fig.6.The field of view model of designed industrial endoscope.

    where

    The imaging area is the middle shadow part of the crosssectional circle formed by the image ellipse and the BF wall on the two-dimensional plane in Fig.6(b).To calculate the imaging area of the endoscope, (6) is rewritten as

    To maximize the imaging area of the burden surface in the whole distribution period, according to the height of the burden surfaceq, the integral (7) can be calculated to obtain the imaging area of the endoscope in the entire material distribution cycle, thus the maximum burden-surface imaging optimization model of the industrial endoscope is established as

    wherexR1,xR2,xR3andxR4are thex-coordinate of the intersections of the image ellipse and the section circle, respectively;Ris the radius of cross-sectional circle;C1-C10are model parameters; α and 2 ω are the shooting angle ∠WAYand the field of view of the industrial endoscope;dandzAare the horizontal distance and vertical distance respectively between the front probe of the endoscope and the center of burden surface; λ0, λ1, λ2and λ3are constraint parameters which need to be established.

    B. Determination of Constraints for the Optimization Model

    1)No-Material-Flow Safety Region on the Top of the BF

    For the purpose of finding the no-material-flow region, the trajectory model of material flow in BF is established according to Newton’s equation of motion as follows [18]:

    In Fig.7(a),ndenotes the distance from the landing point of the material flow to the central line,Lxrepresents the horizontal distance on thexaxis from the landing point to the end of the rotating chute,K1is the velocity of the material, β is the angle between the chute and the horizontal direction,l0,lβ,ω0,eand μ describe the length, the effective working length,rotating speed, chute tilted holding, and coefficient of the sliding friction of the chute, respectively,gis gravity acceleration,andhis the distance from the bottom end of the chute to the top end of the stock heap.

    Fig.7.Material flow trajectory and low dust area on the top of the BF.

    Based on this model and the relevant parameters of BF distribution, the theoretical calculation value of the distance between the drop point of the material flow on the 1.5 m standard stockline and BF wall is 0.54 m.Thus, the distance between the front-end probe of the industrial endoscope and BF wall should not exceed 0.54 m to keep a certain safety margin.

    2)Low Dust Installation Region on the Top of the BF

    To determine the low dust installation region, the turbulent flow mixed with gas flowk-Sε model and the dust particle motionupmodel on the BF top are established based on the engineering model theory and Newton’s theorem of motion[23].And based on the interaction between phases, the coupled distributionk-Sε-upmodel is established as (9).

    There are four independent equations in the model, whereKˉ , εˉ are used to describe the motion of the gas flow field,τpdescribes the interaction between phases,upidenotes the motion of dust particles.Andui,ujare the speeds of air in different directions.μ, ρ,Kˉ and εˉ are dynamic viscosity parameters, density, turbulent kinetic energy and turbulent dissipation rate of gas flow respectively.Prepresents pressure on top of BF,fand τpare the resistance correction coefficient and relaxation time of the dust particles respectively.spi,av,upianduiare the displacement, acceleration, velocity of dust particle, velocity of the gas flow, respectively.Cε1,Cε2,Cε3,Cεε,Cεk,CT,CDare the key model constants parameters to be determined, which can be obtained by relevant engineering manuals.

    According to the model, the installation position of the industrial endoscope can be determined as the center of the interface between the cylindrical surface and the conical surface of BF roof as shown in Fig.7(b).Compared with the installation position of the traditional equipment, the structure of this installation position is simple, with only one layer of furnace shell, which is easy to open holes, and does not affect the safety of BF operation.And the installation position is far from the top gas flow riser pipe.Therefore, the proposed installation is more reasonable.

    3)Determination of Optimal Configuration Constraints

    Based on the field data, the radius of the section circleRis 4.15 m, the coordinate of installation pointDis(0, 4.64, 2.46)and that ofAis (xA,yA,zA).According to the material flow trajectory model, the distance betweenAand BF wall is less than 0.54 m.To prevent particle flow from impacting the endoscope, the distancedbetween the front end of the endoscope and the BF should be greater than 3.68 m and less thanR,therefore, 3.68 ≤d≤4.15.As shown in Fig.6(a), due to the angle ∠DNMbetween the line formed by pointDand the edge pointNof BF wall and the horizontal line is 60°, the projectionof∠DNMontheyOzp√lane cannot be bigger than 60°,then

    Based on above analysis, the model parameters and the constraints parameters are listed in Appendix B.

    C. Solution of the Maximum Imaging Optimization Model

    The calculation results of the PSO algorithm, the simulated annealing algorithm and the genetic algorithm are shown in Table III.It can be seen that the global optimal value calculated by the PSO algorithm is larger, the number of iterations is less, and the calculation time is shortest, indicating that the efficiency of PSO algorithm is higher, and it has a strong global search ability, which is easier to escape from the local optimum.Therefore, the PSO algorithm is more suitable for solving the proposed optimization model.

    TABLE III CALCULATION RESULTS OF PSO AND OTHER OPTIMIzATION ALGORITHMS

    From the table, the coordinate of pointAthat maximizes the imaging area of the burden surface is (1.9,3.6,0.7), that of pointBis (2.45,2.72,0.84).Thus, the projection angle β on thexOzplaneofαis 42.8?,and theprojectionangleγontheyOzplaneis59.4?,maxS dq=36.77.Finally, theequipment installation method for safe and close-range imaging of BF burden surface is shown in Fig.8(a).

    V.RESULTS AND DISCUSSION

    A. Application of the Designed Industrial Endoscope

    Fig.8.Installation method and results of starlight industrial endoscope.

    This section introduces the actual installation structure and operation of SHIE on 2650 m3BF.Fig.8(b) shows the opening hole of the BF.The diameter of the hole is 110 mm, and the shape of the burden surface inside can be clearly seen through it.Fig.8(c) shows the rear sleeve shell of the device outside, which is installed on the fixed flange and weld at the opening hole, to guarantee the airtightness of the BF.And the back end of rear sleeve shell is connected to the water-cooling circuit, the industrial nitrogen pipe, the power and video signal lines.Fig.8(d) indicates that the industrial endoscope is installed under the cross-temperature measurement device with proper installation posture.Fig.8(e) displays the distribution box of the equipment, which contains the power supply,video signal module and temperature measurement module,where the temperature measurement module shows the realtime temperature in the front and rear sleeves of the device.

    The images captured by the developed endoscope during the entire BF operation are shown in Fig.9.The device can not only capture the burden surface images during normal operation, but also can reveal new abnormal BF conditions, which could not be obtained by the existing equipment.For example,Figs.9(a)-9(c) indicate the slow sinking and sudden drop of the burden surface.From Figs.9(a) and 9(b), the burden surface only decreases a little in two minutes, while from Figs.9(b) and 9(c), the burden surface drops rapidly within five seconds, and the abnormal furnace condition of collapsing burden appears.Figs.9(d)-9(f) show the abnormal furnace condition of splashing burden, the splashing burden is the bright spot falling on the edge area.Due to the splashing burden, more dust is generated, reducing the resolution of the images.Figs.9(g)-9(i) show the images when splashing water on the BF top, where the water mist and dust are shrouded over the BF top, but the shape of the water column and the burden surface can still be observed.After splashing water,the highlight of the image becomes dark, that is, the flame of the splashed burden gradually goes out.

    B. Qualitative Analysis of the Imaging Performance of the Starlight High-Temperature Industrial Endoscope

    Fig.9.The obtained image by the starlight industrial endoscope.

    Fig.10.Burden surface images obtained by different devices.

    To further verify the imaging performance of the designed device, the captured images by the device with the high-temperature industrial endoscope and the infrared camera are compared.Figs.10(a)-10(c) are pictures when starting BF by these devices, Figs.10(d) and 10(e) show the images of obtained burden surface during the idle period.To be specific,Figs.10(a) and 10(d) are images taken by the infrared camera,in which the burden surface can hardly be seen, and only the shape of central highlight area and cross-temperature measurement device can be seen blurred.Figs.10(b) and 10(e) are images from the high-temperature industrial endoscope, where the approximate shape of burden surface and furnace wall can be observed, but the image brightness is low.While the images obtained by the developed device are distinct and bright.The shape of the coke particles and the outline of the burden surface from concave to convex are clear.Based on the data measured by the south and north detection points, the inclination angle and inclination direction of the burden surface can be determined, the device can thus provide intuitive and reliable information of the burden surface.

    C. Quantitative Analysis of the Imaging Performance of the Starlight High-Temperature Industrial Endoscope

    The accuracy, reliability and image quality of the acquired images are evaluated to quantitatively analyze the performance of the developed equipment.The correct evaluation must be based on a large number of images.Thus, the dataset used in this study consists of 97 730 burden surface images collected from the developed BF burden surface detection system as shown in Fig.11.It includes 19 546 original burden surface images captured by different devices installed on the BF top, namely SHIE, the high-temperature industrial endoscope, the mechanical swinging radar, the laser scanner and the infrared camera.Based on the dataset, the image accuracy,the recognition rate for abnormal furnace conditions and the image quality of different equipment are quantitatively compared from the following aspects.

    Fig.11.3D monitoring system of BF burden surface.

    1)Comparison Results of the Accuracy and Reliability of the Images

    Generally, in industrial production sites, the ideal empirical curve of the burden surface is usually used as the gold standard for judging the material distribution of the burden surface.Considering the current smelting situation, workers always plot an ideal experience curve based on their own experiences and the measured data of the mechanical probe.Therefore, to compare the similarity between the burden surface images collected by the above detection equipment and the gold standard, and judge the accuracy of the images, we use Canny edge extraction operator to extract the burden surface curves in the images collected by the SHIE, infrared camera, laser scanner and mechanical swing radar, and set them in the same coordinate system, as shown in Fig.12.

    As shown in Fig.12(a), the burden surface heights measured by the north and south mechanical probes are 1.22 m and 1.55 m, respectively, and the ideal experience curve drawn by workers is a typical M-shaped burden surface.It can be seen that the detection points of the two mechanical probes are located on the burden surface contour curve of SHIE, and the contour curve is in good agreement with the ideal empirical curve, especially in the high light area in the BF center, the concave trend of burden surface shape is obvious.The contour curve of burden surface detected by the laser scanner is consistent with the ideal curve shape, but it does not include the detection point of southern mechanical probe, and there is a certain error.The burden surface contour curve detected by the mechanical swing radar includes the detection point of the north mechanical probe but does not include the detection point of the south mechanical probe, and the curve exhibits a significant deviation from the ideal empirical curve as it protrudes at the high light position in the center of the blast furnace.The burden surface contour curve of the infrared camera is a closed curve, which does not include the detection points of the north-south mechanical probe, which is completely different from the ideal empirical curve.Therefore,when detecting the normal burden level of BF, the burden surface image captured by SHIE is more accurate, the burden surface profile detected by the laser scanner has a certain error, and the burden surface profile obtained by the mechanical swing radar has a large error, while the infrared camera cannot obtain the correct burden surface profile.

    Fig.12.The burden surface contour curves under normal and inclined burden condition.

    When the inclined burden condition occurs as shown in Fig.12(b), the measured heights of the burden surface by the north and south mechanical probes are 0.45 m and 1.75 m,respectively.Both the north and south mechanical probe points locate on the contour curve of the burden surface captured by the SHIE, which closely aligns with the trend of the ideal empirical curve and exhibits a significant depression at the center.The contour curve of the burden surface detected by the laser scanner exhibits a resemblance to the ideal empirical curve in terms of overall shape but does not include the measurement points from the north and south mechanical probes, resulting in noticeable errors.The contour curve obtained from the mechanical swing radar does not include the measurement points from the north and south mechanical probes and exhibits an opposite trend compared to the ideal empirical curve.The curve obtained from the infrared camera approximates a circle-shape, deviating greatly from the ideal empirical curve.Therefore, in the presence of abnormal furnace conditions such as inclined burden, the detection accuracy of the SHIE remains high, while the laser scanner experiences significant declines in accuracy.The mechanical swing radar and the infrared camera are unable to accurately identify abnormal furnace conditions.

    The above analysis compares the accuracy of the burden surface image at a specific time.The correctness evaluation must also be based on a large set of images.Therefore,according to the golden standard, the correctness of 97 730 images captured by the system was evaluated.For the 19 546 burden surface images collected by the SHIE, 17 277 are correct and 2269 are wrong.For the 19 546 burden surface images acquired by the high temperature industrial endoscope,16 648 are correct and 2898 are wrong.For the 19 546 burden surface images acquired by the mechanical swing radar,15 000 are correct and 4546 are wrong.For the 19 546 burden surface images acquired by the laser scanner, 13 996 are correct and 5550 are wrong.For the 19 546 burden surface images acquired by the infrared camera, 317 are correct and 19 299 are wrong.All these images are checked by workers and manually sorted into five classes.Class 1 to Class 5 represent absolutely correct images, probably correct images, indeterminate images, probably wrong images, and absolutely wrong images, respectively.The classification results are shown in Table IV.

    TABLE IV EVALUATION RESULTS OF 19 546 IMAGES

    According to Table IV, the ROC (receiver operating characteristic) curves are drawn in Fig.13.From Fig.13, the best working point of SHIE is the first point, while that of the high-temperature industrial endoscope and mechanical swing radar is the second point, and the best working point of the laser scanner and infrared camera is the third point.That is, in practical industrial application, according to the gold standard judgment, the infrared camera, laser scanner, mechanical swing radar and high-temperature industrial endoscope capture the definite correct burden surface images corresponding to the worker’s subjective recognition of uncertain images,uncertain images, possible correct images and possible correct images, respectively, which makes the collected images have uncertainty and affects the reliability of the device.As for the SHIE, the determined image is very consistent with the absolutely correct image, which means that the reliability of the equipment is greatly improved.On the other hand, the ROC curve of SHIE is closer to the upper left corner of the coordinate system than the ROC curves of the infrared camera, laser scanner, mechanical swing radar, and high-temperature industrial endoscope, which reflects that the image of BF burden surface obtained by SHIE has higher correctness.

    Fig.13.The ROC curves of different devices.

    The AUC (area under curve), S.E.(standard error of area)and C.I.(confidence interval) at a 95% confidence level are calculated by ROC curve in Table V, to further verify the performance of SHIE.For AUC, the closer the value of AUC to 1, the higher the precision of the device.For the S.E.and C.I.indicator at the 95% level confidence, the smaller the S.E.value or the narrower the confidence interval, the higher the reliability of the device.In Table V, all indicators of SHIE are superior to the other devices, denoting that the burden surface image captured by the developed equipment has higher accuracy and reliability.Additionally, according to industry standards, images with AUC in the range of 0.95-0.99, 0.90-0.95,0.85-0.90, 0.80-0.95 and below 0.80 are defined as “completely credible images”, “credible images”, “basically credible images”, “uncertain images” and “untrustworthy images”,respectively.Therefore, the evaluations of images obtained by SHIE, the high-temperature industrial endoscope, the mechanical swing radar, the laser scanner and the infrared camera are completely credible, credible, basically credible, uncertain and untrustworthy, respectively.

    TABLE V INDICATORS OF ROC CURVES

    2)Recognition Rate of Abnormal Furnace Conditions by Different Detection Equipment

    The abnormal furnace conditions of BF refer to the abnormal physical and chemical phenomena that occur during the production process, such as hanging burden, collapsing burden, partial burden surface, low stockline, etc.These abnormal furnace conditions not only affect the stability, production, quality, and energy consumption of the BF but also bring a high risk of safety accidents.Therefore, it is important to correctly judge the abnormal furnace conditions, so that the onsite workers can control the BF reasonably in time and ensure the smooth operation.Based on the image dataset, the recognition rates of the current burden surface detection equipment under abnormal furnace conditions were calculated and are illustrated in Fig.14.

    Fig.14.The box plot of recognition rates of different equipment for abnormal furnace conditions.

    Fig.14 shows the box plot of the recognition rates of SHIE,mechanical swinging radar, high-temperature industrial endoscope, infrared camera, and laser scanner under abnormal furnace conditions at the BF top.The line within the box represents the median of the data, and the box contains 50% of the data, reflecting the fluctuation of the data.The upper and lower edges of the box represent the maximum and minimum values of the data set, respectively.As shown in Fig.14, the recognition rates of SHIE under four abnormal furnace conditions are generally between 0.8 and 0.95, much higher than the other devices.The box height of SHIE is relatively low,concentrated around 0.9, and the lengths of the upper and lower edges are similar, indicating its relatively stable overall recognition rate.Therefore, compared to the other existing burden surface detection equipment, SHIE exhibits a high recognition rate and good stability for abnormal furnace conditions.It can accurately identify abnormal furnace conditions most of the time and even detect abnormal phenomena that other devices may ignore, greatly contributing to the safe and stable operation and the improvement of smelting quality of BF.

    3)Comparison Results of the Image Quality

    To evaluate the quality of images captured by the three different devices, the image entropy, sum of image power spectral components, gray mean gradient, edge intensity, point sharpness and NIQE (the natural image quality evaluator) are calculated in Table VI [24].The standard image database is from the image quality evaluation database.Compared with the other detection equipment, the value of the image entropy and the sum of image power spectral components of SHIE are closer to the standard value, that means the images captured by developed industrial endoscope have more information than those captured by others.Also, the gray mean gradient,edge intensity and point sharpness of developed equipment aremuch larger than those of the infrared camera and high-temperature industrial endoscope, but slightly smaller than the standard values, which indicates that the image clarity of the developed device has been significantly improved compared with existing equipment.The NIQE indicator of SHIE is smaller than that of the other devices, slightly larger than the standard value, indicating that its image has better perceptual quality objectively.In summary, the burden surface images captured by SHIE have high clarity, rich information content,and good perceptual quality, which is far superior to the images obtained by the other inspection devices, but there is still room for improvement.

    TABLE VI COMPARISONS OF THE IMAGE QUALITY ASSESSMENT INDICATORS

    VI.CONCLUSIONS

    In this paper, a BF burden surface detection method based on the SHIE is proposed.The designed starlight imaging system of the method can acquire high-brightness, high-resolution burden surface images under extremely dim light and run steadily in the high-temperature and high-pressure environment of BF.On this basis, this paper discovers the optimum configuration of the installation position of the industrial endoscope, which ensures the equipment to obtain the maximum burden surface image and prevents the damage of highspeed material flow and high-concentration dust.The developed industrial endoscope has been installed on the top of a 2650m3BF and runs stably.Compared with the burden surface images obtained by the other detection equipment, it is verified that the images captured by the designed equipment are more reliable, where the contours and details of the burden surface can be clearly observed.The developed optical sensing imaging device provides a novel and beneficial method for detecting the burden surface of BF.

    APPENDIX A

    The comparison of data indicators of different burden surface detection methods is listed in Table VII.

    TABLE VII COMPARISON OF DATA INDICATORS OF THE BURDEN SURFACE DETECTION METHODS

    TABLE VIII OPTIMIzATION MODEL PARAMETERS

    APPENDIX B

    The important parameters and constants in the optimal model are arranged in Table VIII.

    国产精品乱码一区二三区的特点| 男女下面进入的视频免费午夜| 91九色精品人成在线观看| 国产成人一区二区三区免费视频网站| 性色av乱码一区二区三区2| 宅男免费午夜| 亚洲精品在线观看二区| 人人妻人人看人人澡| 特大巨黑吊av在线直播| 男人的好看免费观看在线视频| 亚洲精品美女久久久久99蜜臀| 国产av在哪里看| 黄色丝袜av网址大全| 亚洲无线在线观看| 日韩人妻高清精品专区| 国产精品一区二区精品视频观看| 久久精品亚洲精品国产色婷小说| 精品久久久久久久人妻蜜臀av| 一级毛片精品| 国产精品乱码一区二三区的特点| 亚洲精品在线美女| 一进一出抽搐gif免费好疼| 日本熟妇午夜| 国产97色在线日韩免费| 国产精品日韩av在线免费观看| 999久久久国产精品视频| 最近最新中文字幕大全电影3| 男女之事视频高清在线观看| 99国产极品粉嫩在线观看| 日韩欧美在线乱码| 免费看日本二区| 麻豆av在线久日| 久久中文字幕人妻熟女| 丁香欧美五月| 国产毛片a区久久久久| 美女扒开内裤让男人捅视频| 啪啪无遮挡十八禁网站| 男人舔女人下体高潮全视频| 亚洲午夜精品一区,二区,三区| 法律面前人人平等表现在哪些方面| 亚洲精品久久国产高清桃花| 十八禁网站免费在线| 国产精品1区2区在线观看.| 国产成人影院久久av| 黄片小视频在线播放| 午夜福利高清视频| 看黄色毛片网站| 精品午夜福利视频在线观看一区| 丝袜人妻中文字幕| 亚洲av成人av| 国产一区二区在线观看日韩 | 无人区码免费观看不卡| 久久亚洲真实| 午夜福利在线观看吧| 成人午夜高清在线视频| 国产精品自产拍在线观看55亚洲| 亚洲av第一区精品v没综合| 观看美女的网站| 久久久久久久精品吃奶| 亚洲七黄色美女视频| 男女之事视频高清在线观看| 免费看美女性在线毛片视频| 色综合站精品国产| 一级a爱片免费观看的视频| 国产精品影院久久| 亚洲 国产 在线| 精品国产亚洲在线| 精品国产亚洲在线| 欧美日韩综合久久久久久 | 热99re8久久精品国产| 欧美日韩福利视频一区二区| 午夜精品在线福利| 1024手机看黄色片| 久99久视频精品免费| 少妇的丰满在线观看| 好男人在线观看高清免费视频| 最近在线观看免费完整版| 免费观看的影片在线观看| 亚洲精华国产精华精| 亚洲国产精品成人综合色| 一个人免费在线观看电影 | 亚洲在线观看片| 99热这里只有是精品50| 成人亚洲精品av一区二区| 欧美一级毛片孕妇| 国产单亲对白刺激| 桃色一区二区三区在线观看| 久久久久性生活片| 在线免费观看不下载黄p国产 | 99在线视频只有这里精品首页| 亚洲人与动物交配视频| 香蕉国产在线看| 久久久久性生活片| 久久久国产成人精品二区| 中出人妻视频一区二区| 一本精品99久久精品77| 特大巨黑吊av在线直播| 日本在线视频免费播放| 精品一区二区三区av网在线观看| 国产久久久一区二区三区| 国产精品一区二区免费欧美| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区三区四区五区乱码| 成年人黄色毛片网站| 欧美xxxx黑人xx丫x性爽| 免费人成视频x8x8入口观看| 这个男人来自地球电影免费观看| 中亚洲国语对白在线视频| 淫秽高清视频在线观看| 国产精品精品国产色婷婷| 国产精品日韩av在线免费观看| 99在线人妻在线中文字幕| 国产av不卡久久| 麻豆国产97在线/欧美| 美女黄网站色视频| 国产免费av片在线观看野外av| 久久久成人免费电影| 色噜噜av男人的天堂激情| 五月伊人婷婷丁香| 久久亚洲精品不卡| 久久精品夜夜夜夜夜久久蜜豆| 久久国产乱子伦精品免费另类| 悠悠久久av| 成人av一区二区三区在线看| 三级国产精品欧美在线观看 | 国产免费av片在线观看野外av| 啦啦啦免费观看视频1| 午夜a级毛片| 久久久成人免费电影| 久久这里只有精品中国| 精品久久久久久,| 成人国产一区最新在线观看| 久久久国产欧美日韩av| 婷婷精品国产亚洲av在线| 国产三级中文精品| 亚洲自拍偷在线| 美女免费视频网站| 成人三级黄色视频| 精华霜和精华液先用哪个| 九九久久精品国产亚洲av麻豆 | 琪琪午夜伦伦电影理论片6080| 极品教师在线免费播放| 国模一区二区三区四区视频 | 国产亚洲精品久久久久久毛片| 大型黄色视频在线免费观看| 国产一区在线观看成人免费| 国产精品久久久久久精品电影| 亚洲av免费在线观看| 听说在线观看完整版免费高清| 嫩草影院精品99| 综合色av麻豆| 亚洲aⅴ乱码一区二区在线播放| 黄色视频,在线免费观看| 岛国在线免费视频观看| 一个人免费在线观看电影 | 亚洲av中文字字幕乱码综合| 欧美日韩瑟瑟在线播放| 一二三四在线观看免费中文在| 国产精品久久久久久久电影 | 熟女电影av网| 日韩欧美国产一区二区入口| 男人舔奶头视频| 日本 av在线| 97超视频在线观看视频| 色综合亚洲欧美另类图片| 少妇丰满av| 丰满的人妻完整版| 亚洲精华国产精华精| 桃色一区二区三区在线观看| 精品一区二区三区av网在线观看| 欧美精品啪啪一区二区三区| 国产高潮美女av| 欧美黄色片欧美黄色片| 亚洲av熟女| 午夜免费激情av| 一卡2卡三卡四卡精品乱码亚洲| 九色国产91popny在线| 成年版毛片免费区| 欧美+亚洲+日韩+国产| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩卡通动漫| 中亚洲国语对白在线视频| 午夜免费激情av| 在线观看免费午夜福利视频| 在线观看免费午夜福利视频| av女优亚洲男人天堂 | 亚洲成av人片在线播放无| 亚洲 欧美一区二区三区| av黄色大香蕉| 天堂影院成人在线观看| 男女那种视频在线观看| 最新美女视频免费是黄的| 中文字幕最新亚洲高清| 精品久久久久久久久久免费视频| 在线国产一区二区在线| 欧美xxxx黑人xx丫x性爽| 成人av在线播放网站| 99久久综合精品五月天人人| 欧美国产日韩亚洲一区| 一级毛片高清免费大全| 小说图片视频综合网站| 色在线成人网| 一区福利在线观看| 老汉色av国产亚洲站长工具| 九九热线精品视视频播放| 99国产综合亚洲精品| 18禁美女被吸乳视频| 国产精华一区二区三区| 欧美日韩一级在线毛片| www国产在线视频色| 亚洲第一电影网av| 免费在线观看日本一区| 久久精品国产综合久久久| 国内精品一区二区在线观看| 女同久久另类99精品国产91| 色在线成人网| 琪琪午夜伦伦电影理论片6080| 99riav亚洲国产免费| 成人鲁丝片一二三区免费| 国内久久婷婷六月综合欲色啪| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成人免费电影在线观看| 精品人妻1区二区| 一区福利在线观看| 国产一区二区激情短视频| 一边摸一边抽搐一进一小说| 欧美日韩一级在线毛片| 99热这里只有精品一区 | 欧美乱色亚洲激情| 国产精品99久久99久久久不卡| 91av网一区二区| 97超视频在线观看视频| cao死你这个sao货| 国产精品美女特级片免费视频播放器 | 久久久久久久久久黄片| 国产欧美日韩精品亚洲av| 制服丝袜大香蕉在线| 最近最新免费中文字幕在线| 1024手机看黄色片| 女人被狂操c到高潮| 搡老熟女国产l中国老女人| 国产99白浆流出| 美女午夜性视频免费| 色综合站精品国产| 午夜福利在线观看吧| 久久久久国产一级毛片高清牌| 国产精品电影一区二区三区| 在线观看免费视频日本深夜| 婷婷精品国产亚洲av在线| 色吧在线观看| 国产一区二区三区在线臀色熟女| 欧美在线一区亚洲| 在线观看免费午夜福利视频| 制服人妻中文乱码| 一本综合久久免费| 香蕉久久夜色| 国内少妇人妻偷人精品xxx网站 | a在线观看视频网站| 亚洲人与动物交配视频| 国产精品 欧美亚洲| 欧美乱码精品一区二区三区| 老汉色∧v一级毛片| 久久婷婷人人爽人人干人人爱| 日本精品一区二区三区蜜桃| 91av网站免费观看| 亚洲五月婷婷丁香| 美女高潮的动态| 国产高清三级在线| netflix在线观看网站| www.自偷自拍.com| 老汉色av国产亚洲站长工具| 欧美xxxx黑人xx丫x性爽| 精品午夜福利视频在线观看一区| 一级毛片高清免费大全| 精品国产超薄肉色丝袜足j| 热99在线观看视频| 色播亚洲综合网| 色综合欧美亚洲国产小说| 在线观看舔阴道视频| 日本五十路高清| 偷拍熟女少妇极品色| 美女 人体艺术 gogo| 99热这里只有是精品50| 亚洲欧美日韩东京热| av天堂在线播放| 哪里可以看免费的av片| 免费搜索国产男女视频| 免费在线观看成人毛片| 国产在线精品亚洲第一网站| 天天一区二区日本电影三级| 此物有八面人人有两片| 91麻豆精品激情在线观看国产| 我要搜黄色片| 精品无人区乱码1区二区| 91老司机精品| 身体一侧抽搐| 亚洲熟妇中文字幕五十中出| 国产精品久久久久久人妻精品电影| aaaaa片日本免费| 中文字幕av在线有码专区| 国产99白浆流出| 国产精品久久久久久精品电影| 黄频高清免费视频| 三级毛片av免费| 精品一区二区三区av网在线观看| 一区二区三区激情视频| 一进一出好大好爽视频| 欧美xxxx黑人xx丫x性爽| 高清在线国产一区| 国产精品 国内视频| 99久国产av精品| 日韩欧美在线乱码| 12—13女人毛片做爰片一| 亚洲国产欧美人成| 国产美女午夜福利| 国产精品电影一区二区三区| 精品福利观看| 黄色女人牲交| 亚洲九九香蕉| 中文字幕精品亚洲无线码一区| 脱女人内裤的视频| 久久欧美精品欧美久久欧美| 亚洲熟妇中文字幕五十中出| 99国产精品一区二区蜜桃av| 在线免费观看不下载黄p国产 | 又大又爽又粗| 特大巨黑吊av在线直播| 色综合婷婷激情| 人妻夜夜爽99麻豆av| 天堂动漫精品| 99久久国产精品久久久| 亚洲无线观看免费| 亚洲自拍偷在线| 国产高清有码在线观看视频| 欧美国产日韩亚洲一区| 免费av不卡在线播放| 香蕉av资源在线| 一个人观看的视频www高清免费观看| 亚洲人成网站在线播| 麻豆成人av视频| 亚洲精品,欧美精品| 久久99热这里只频精品6学生 | 国产淫语在线视频| 插逼视频在线观看| 久久久成人免费电影| 欧美xxxx黑人xx丫x性爽| 免费观看精品视频网站| 日本黄大片高清| eeuss影院久久| 日本免费a在线| 中文字幕免费在线视频6| 伊人久久精品亚洲午夜| 日本与韩国留学比较| 嫩草影院精品99| 中文字幕精品亚洲无线码一区| 大又大粗又爽又黄少妇毛片口| 亚洲中文字幕日韩| 久久这里只有精品中国| 91精品伊人久久大香线蕉| 建设人人有责人人尽责人人享有的 | 免费观看a级毛片全部| 欧美激情在线99| 老女人水多毛片| av.在线天堂| 天堂√8在线中文| 久99久视频精品免费| 国产午夜精品一二区理论片| 美女大奶头视频| 免费搜索国产男女视频| 亚洲精品乱码久久久v下载方式| 在线观看66精品国产| av又黄又爽大尺度在线免费看 | 国产乱人偷精品视频| 小说图片视频综合网站| 国产一区亚洲一区在线观看| 最近2019中文字幕mv第一页| 国产成人freesex在线| 亚洲av.av天堂| 久久久国产成人免费| 国产成人午夜福利电影在线观看| 亚洲精品久久久久久婷婷小说 | 精品无人区乱码1区二区| 国产伦一二天堂av在线观看| 深夜a级毛片| 一区二区三区乱码不卡18| 看片在线看免费视频| 2022亚洲国产成人精品| 亚洲av熟女| 九九爱精品视频在线观看| 日本黄色视频三级网站网址| 国产av不卡久久| 国产精品蜜桃在线观看| 亚洲激情五月婷婷啪啪| 午夜老司机福利剧场| 国产伦在线观看视频一区| 少妇的逼好多水| 亚洲精品一区蜜桃| 亚洲精品乱码久久久久久按摩| 淫秽高清视频在线观看| 亚洲一区高清亚洲精品| 国产精品乱码一区二三区的特点| 午夜a级毛片| 亚洲国产最新在线播放| 欧美变态另类bdsm刘玥| 亚洲美女视频黄频| 成人性生交大片免费视频hd| 精品人妻一区二区三区麻豆| 国产免费福利视频在线观看| 乱人视频在线观看| 中文字幕免费在线视频6| 国产精品久久久久久久电影| 国产精品日韩av在线免费观看| 桃色一区二区三区在线观看| 精品酒店卫生间| 男女边吃奶边做爰视频| 欧美xxxx性猛交bbbb| 国产精品国产高清国产av| 99热网站在线观看| 好男人在线观看高清免费视频| 国产精品不卡视频一区二区| 日本五十路高清| 偷拍熟女少妇极品色| 一卡2卡三卡四卡精品乱码亚洲| 老司机影院毛片| 欧美精品国产亚洲| 国产亚洲一区二区精品| 黄色日韩在线| 国产国拍精品亚洲av在线观看| 日韩一区二区三区影片| 看十八女毛片水多多多| 久久人人爽人人爽人人片va| 99热6这里只有精品| 搞女人的毛片| 国产精品一及| 免费在线观看成人毛片| 小说图片视频综合网站| 亚洲精品国产成人久久av| 国产精品久久电影中文字幕| 国产色婷婷99| 日产精品乱码卡一卡2卡三| 美女内射精品一级片tv| 级片在线观看| 亚洲精品aⅴ在线观看| 亚洲av电影在线观看一区二区三区 | 精品无人区乱码1区二区| 看免费成人av毛片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久精品久久久久久久性| 亚洲久久久久久中文字幕| 最近中文字幕2019免费版| 有码 亚洲区| 嘟嘟电影网在线观看| 人体艺术视频欧美日本| 99在线人妻在线中文字幕| 成人av在线播放网站| 欧美日韩国产亚洲二区| 国产伦精品一区二区三区视频9| 大香蕉97超碰在线| 精品99又大又爽又粗少妇毛片| 少妇熟女aⅴ在线视频| 久久久久久久久大av| 久久久色成人| 黄色一级大片看看| 欧美一区二区国产精品久久精品| 18禁动态无遮挡网站| 亚洲中文字幕一区二区三区有码在线看| 国产高清国产精品国产三级 | 欧美成人精品欧美一级黄| 国产精品1区2区在线观看.| 偷拍熟女少妇极品色| 成人美女网站在线观看视频| eeuss影院久久| a级毛色黄片| 成人毛片60女人毛片免费| 成人无遮挡网站| 91狼人影院| 国产精品日韩av在线免费观看| 国产在线一区二区三区精 | 国产私拍福利视频在线观看| 啦啦啦啦在线视频资源| 亚洲在久久综合| 国产精品一及| 久久精品夜色国产| 国产精品.久久久| 99在线人妻在线中文字幕| 欧美3d第一页| 国产极品精品免费视频能看的| 青青草视频在线视频观看| 久久久a久久爽久久v久久| 国产极品天堂在线| 2021少妇久久久久久久久久久| 高清视频免费观看一区二区 | 久久国产乱子免费精品| 热99re8久久精品国产| 亚洲中文字幕日韩| 69av精品久久久久久| 黄色欧美视频在线观看| 精品国产三级普通话版| 亚洲人与动物交配视频| 成年版毛片免费区| 青春草亚洲视频在线观看| 老师上课跳d突然被开到最大视频| 两个人的视频大全免费| 亚洲av电影在线观看一区二区三区 | 亚洲av免费在线观看| 中文字幕免费在线视频6| 色噜噜av男人的天堂激情| 成人二区视频| 日韩欧美精品v在线| 精品一区二区三区视频在线| 日韩,欧美,国产一区二区三区 | 亚洲美女搞黄在线观看| 又粗又爽又猛毛片免费看| 能在线免费看毛片的网站| 91狼人影院| 欧美性猛交黑人性爽| av黄色大香蕉| 精品人妻视频免费看| 亚洲人与动物交配视频| 亚洲国产高清在线一区二区三| 嫩草影院入口| 国产精品99久久久久久久久| 国国产精品蜜臀av免费| 最近最新中文字幕大全电影3| 色噜噜av男人的天堂激情| 超碰97精品在线观看| 熟女电影av网| 18禁在线播放成人免费| 又黄又爽又刺激的免费视频.| 男人舔女人下体高潮全视频| 中文字幕精品亚洲无线码一区| 国产一级毛片七仙女欲春2| 好男人视频免费观看在线| 特级一级黄色大片| 久久久久久国产a免费观看| videos熟女内射| 男女下面进入的视频免费午夜| 麻豆乱淫一区二区| 中文天堂在线官网| 久久国内精品自在自线图片| av天堂中文字幕网| 一二三四中文在线观看免费高清| 日日摸夜夜添夜夜添av毛片| 精品久久久久久久久av| 高清毛片免费看| 亚洲自偷自拍三级| 成人性生交大片免费视频hd| 午夜福利在线观看免费完整高清在| 六月丁香七月| 国产成人午夜福利电影在线观看| 亚洲国产欧美在线一区| av黄色大香蕉| 少妇丰满av| 中文字幕制服av| 日本色播在线视频| 国产69精品久久久久777片| 夫妻性生交免费视频一级片| 男女视频在线观看网站免费| 欧美日韩国产亚洲二区| 91久久精品国产一区二区三区| 亚洲三级黄色毛片| 国产一区亚洲一区在线观看| 亚洲av电影不卡..在线观看| 亚洲欧美精品专区久久| 三级国产精品片| 午夜福利高清视频| 亚洲国产色片| 一区二区三区高清视频在线| 中文乱码字字幕精品一区二区三区 | 国产成人a∨麻豆精品| 高清视频免费观看一区二区 | 亚洲一级一片aⅴ在线观看| 别揉我奶头 嗯啊视频| 少妇被粗大猛烈的视频| 嘟嘟电影网在线观看| 久久99蜜桃精品久久| 国产精品久久久久久精品电影小说 | 欧美极品一区二区三区四区| 男女下面进入的视频免费午夜| 在线播放国产精品三级| 国产精品,欧美在线| 国产在线一区二区三区精 | 超碰av人人做人人爽久久| 中国美白少妇内射xxxbb| 听说在线观看完整版免费高清| 亚洲精华国产精华液的使用体验| 舔av片在线| 看黄色毛片网站| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲av涩爱| 国产黄片视频在线免费观看| 亚洲最大成人中文| 岛国在线免费视频观看| 少妇高潮的动态图| 国产精品久久久久久av不卡| h日本视频在线播放| 国模一区二区三区四区视频| 夫妻性生交免费视频一级片| 18禁裸乳无遮挡免费网站照片| 欧美区成人在线视频| 一本久久精品| 午夜福利在线观看免费完整高清在| 国产亚洲av嫩草精品影院| 寂寞人妻少妇视频99o| 久久久久网色| 日韩人妻高清精品专区| 欧美xxxx性猛交bbbb| 久久久午夜欧美精品| 国产视频内射| 亚洲精品日韩av片在线观看| 老司机福利观看| 亚洲av成人精品一二三区| 高清在线视频一区二区三区 | 我要搜黄色片| 日本熟妇午夜|