Method and device for constructing SCR (Selective Catalytic Reduction) control simulation system

Abstract

The invention belongs to the technical field of SCR (Selective Catalytic Reduction) denitration, and discloses a method and a device for constructing an SCR control simulation system, and the method comprises the following steps: screening influence factors influencing the denitration efficiency of a reactor by constructing a related parametric variable thermodynamic diagram influencing the denitration efficiency; by calculating the correlation between the influence on the denitration efficiency of the reactor and the parameter variables, screening to obtain main influence factors of the simulation system; calculating the correlation between the valve action signal and the CEMS valve opening signal to obtain the delay characteristic of the CEMS instrument; and training historical data of main influence factors of the simulation system through an algorithm, establishing a corresponding mathematical model, and performing parameter adjustment on the mathematical model through CEMS instrument delay characteristics to obtain the SCR control simulation system. According to the method, the NOx concentration of the outlet can be stably controlled under the condition that the load of the unit is changed and other disturbances occur.

Description

Technical Field

[0001] The invention belongs to the technical field of SCR denitration, and in particular relates to a construction method and device of an SCR control simulation system. Background Art

[0002] For NOx emission 10mg / m 3 The control target has a higher standard for control accuracy. Under the conditions of load rate and coal quality changes, the NOx concentration and flue gas volume generated by the boiler vary greatly. The control object of the denitrification system exhibits complex nonlinear and servo characteristics, which increases the control difficulty. In order to ensure the stable control of the outlet NOx concentration of the unit under variable load and other disturbances, to keep the control system robust, and to ensure the input rate of the automatic system under the premise of safety and stability, it is necessary to design a control strategy based on simulation testing on the basis of process optimization and transformation and the development of more advanced instrumentation equipment.

[0003] The Chinese patent publication number is CN111044668A, and the name is a patent application for an online prediction method for the life of a catalyst in an SCR denitration system. The method includes: obtaining the NOx inlet concentration during the catalyst operation cycle and the NOx outlet concentration data of each detection point in different areas of the flue cross section; calculating the average denitration efficiency of different areas of the flue cross section according to the NOx inlet concentration and NOx outlet concentration data; establishing the relationship between the denitration efficiency and the catalyst activity according to the actual ammonia nitrogen molar ratio at the catalyst inlet; calculating the average activity of the catalyst in the cycle according to the average denitration efficiency and the relationship between the denitration efficiency and the catalyst activity; calculating the catalyst deactivation rate constant according to the average activity of the catalyst in the cycle, obtaining the relationship between the catalyst activity and the catalyst in-service time, establishing the catalyst life deactivation curve, and predicting the catalyst life. This patent application cannot achieve dynamic prediction and control of NOx concentration. Summary of the invention

[0004] In order to overcome the problems existing in the above-mentioned prior art, the purpose of the present invention is to provide a method and device for constructing an SCR control simulation system, and to design an SCR control simulation system by considering the parameters affecting the denitrification efficiency of the reactor and the delay characteristics of the instrument, so as to realize dynamic prediction and control of NOx concentration.

[0005] To achieve the above object, the technical solution adopted by the present invention is: In a first aspect, the present invention provides a method for constructing an SCR control simulation system, comprising the following steps: S1: Screening the factors affecting the denitrification efficiency of the reactor by constructing a heat map of parameter variables related to the denitrification efficiency; S2: By calculating the correlation between the denitrification efficiency of the reactor and the parameter variables, the main influencing factors of the simulation system are screened; S3: Obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the CEMS valve opening signal; S4: The algorithm is used to train the historical data of the main influencing factors of the simulation system, and the corresponding mathematical model is established. The parameters of the mathematical model are adjusted according to the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

[0006] Optionally, before step S1, the NOx fluctuation data at the reactor inlet, the NOx fluctuation data at the reactor outlet, the NOx fluctuation data at different load rates and the flow characteristic curve data of the ammonia injection regulating valves on both sides are obtained.

[0007] Optionally, in step S2, the Pearson correlation coefficient formula is used to calculate the correlation between the parameter variables.

[0008] Optionally, in step S3, based on the correlation between the valve action signal and the signal of the CEMS valve opening, the statistical variation concentration is smoothed based on the time series data to construct a correlation analysis data matrix.

[0009] Optionally, in step S4, a grey optimization algorithm is used to train the simulation system.

[0010] Optionally, after step S4, the dynamic process of the denitration reduction reaction in the reactor under different ammonia nitrogen molar ratios is tested based on a simulation system for verification.

[0011] In a second aspect, the present invention provides a system for constructing an SCR control simulation system, comprising: A screening module is used to screen the factors affecting the denitration efficiency of the reactor by constructing a thermal map of parameter variables related to the denitration efficiency; The first calculation module is used to screen out the main influencing factors of the simulation system by calculating the correlation between the denitration efficiency of the reactor and the parameter variables; The second calculation module is used to obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the signal of the CEMS valve opening; The system building module is used to train the historical data of the main influencing factors of the simulation system through algorithms, establish the corresponding mathematical model, and adjust the parameters of the mathematical model through the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

[0012] In a third aspect, the present invention provides an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the SCR control simulation system when executing the computer program.

[0013] In a fourth aspect, the present invention provides a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the SCR control simulation system method is implemented.

[0014] In a fifth aspect, the present invention provides a computer program product comprising a computer-readable medium, on which a computer-readable program code is contained, and the program code executes the SCR control simulation system method.

[0015] Compared with the prior art, the present invention has the following beneficial effects: The basic structure of the simulation system constructed by the construction method of an SCR control simulation system of the present invention is based on data analysis software as the bottom layer, modularized bottom layer functions, and a high-fidelity simulation system is constructed by combining mechanism, experience, and data model to verify the reliability of the control strategy in advance. The system of the present invention guarantees the investment rate under the premise of safety and stability. The present invention designs advanced control strategies based on intelligent data analysis and simulation tests. The present invention can ensure the stable control of the outlet NOx concentration of the unit under variable load and other disturbances, so that the control system maintains robustness.

[0016] The method of the present invention can provide customized simulation system construction and combine the unit's own operating parameters (load, coal volume change, coal mill combination, catalyst service life) training model, which can more realistically reflect the operating status of the unit's denitrification system. At the same time, it will also benefit the overall life management of the denitrification system. BRIEF DESCRIPTION OF THE DRAWINGS

[0017] The drawings described herein are for illustration purposes only and are not intended to limit the scope of the present disclosure in any way.

[0018] In the attached picture: Figure 1 The NOx fluctuation at the inlet of the reactor of the present invention; Figure 2 The NOx fluctuation at the reactor outlet of the embodiment of the present invention; Figure 3 The NOx fluctuations under different load rates in the embodiment of the present invention; Figure 4 This is a heat map of the parameters and variables related to the impact on denitration efficiency in the embodiment of the present invention; Figure 5 This is a curve diagram of the delay characteristic analysis of the CEMS according to an embodiment of the present invention; Figure 6 This is a comparison diagram of the effect of the SCR reaction dynamic simulation model parameter optimization algorithm according to an embodiment of the present invention; Figure 7 A flowchart of the steps of constructing a method according to an embodiment of the present invention. DETAILED DESCRIPTION

[0019] In order to enable those skilled in the art to better understand the technical solutions in the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of the present invention.

[0020] Unless otherwise defined, all technical and scientific terms used herein have the same meaning as those commonly understood by those skilled in the art to which the present invention belongs. The terms used herein in the specification of the present invention are only for the purpose of describing specific embodiments and are not intended to limit the present invention. The term "and / or" used herein includes any and all combinations of one or more of the related listed items.

[0021] It should be noted that in the claims, any reference signs placed between brackets shall not be construed as limiting the claims. The word "comprising" does not exclude the presence of components or steps not listed in the claims. The word "a" or "an" preceding a component does not exclude the presence of a plurality of such components. The present application may be implemented by means of hardware comprising several different components and by means of a suitably programmed computer. In a unit claim enumerating several means, several of these means may be embodied by the same item of hardware. The use of the words first, second, and third etc. does not indicate any order. These words may be interpreted as names. In addition, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of the indicated technical features. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of this application, the meaning of "plurality" is two or more, unless otherwise clearly and specifically defined. The present invention is described in detail below with reference to the accompanying drawings.

[0022] A method for constructing an SCR control simulation system of the present invention comprises the following steps: S1: Screening the factors affecting the denitrification efficiency of the reactor by constructing a heat map of parameter variables related to the denitrification efficiency; S2: By calculating the correlation between the denitrification efficiency of the reactor and the parameter variables, the main influencing factors of the simulation system are screened; S3: Obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the signal measured by the CEMS; S4: The algorithm is used to train the historical data of the main influencing factors of the simulation system, and the corresponding mathematical model is established. The parameters of the mathematical model are adjusted according to the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

[0023] The basic structure of the simulation system constructed by the construction method of an SCR control simulation system of the present invention is based on data analysis software as the bottom layer, modularized bottom layer functions, and a high-fidelity simulation system is constructed by combining mechanism, experience, and data models to verify the reliability of the control strategy in advance. It can ensure the stable control of the outlet NOx concentration of the unit under variable load and other disturbances, so that the control system maintains robustness.

[0024] Example 1 A method for constructing an SCR control simulation system in this embodiment includes the following steps: First, data is obtained, including: NOx fluctuation data at the reactor inlet, NOx fluctuation data at the reactor outlet, NOx fluctuation data at different load rates, and flow characteristic curve data of the ammonia injection regulating valves on both sides. Figures 1 to 3 shown.

[0025] Construct a heat map of parameter variables related to denitrification efficiency to find out the main factors affecting the denitrification efficiency of the reactor.

[0026] The Pearson correlation coefficient formula is used to calculate the correlation between variables:

[0027] Among them, x i and i are the observed values ​​of the variables, and x‾ and y‾ are their means.

[0028] In this embodiment, an absolute value of the correlation coefficient greater than 0.8 is considered to be a strong correlation, 0.5 to 0.8 is a moderate correlation, 0.2 to 0.5 is a weak correlation, and less than 0.2 is considered to be no significant correlation.

[0029] The statistical change concentration is smoothed based on the time series data, and the correlation analysis data matrix is ​​constructed to obtain the approximate delay characteristics of the original CEMS instrument. The change concentration refers to the NOx concentration measured by the CEMS.

[0030] Specifically, for data such as reactor inlet NOx concentration and outlet NOx concentration, the moving average method is used to smooth the time series data to eliminate accidental changes, periodic fluctuations and irregular changes in the time series, thereby more clearly showing its long-term trend or seasonal characteristics.

[0031] Taking the smoothed outlet NOx concentration as the reference variable, the time window of the inlet NOx concentration and other related variables is shifted, and the Pearson correlation coefficients between the related variables and the reference variables under different offsets are calculated to establish a correlation coefficient matrix. For the established correlation coefficient matrix, the time point difference when the correlation coefficient between the related variables and the reference variables reaches the maximum value is identified. This time point difference can be regarded as the approximate delay time of the CEMS instrument.

[0032] Matlab was used for data analysis. The construction of the reaction simulation model of the denitrification reactor is based on the denitrification reaction mechanism, and the dynamic process of the catalytic reduction reaction mechanism is based on the Eley-Rideal mechanism as the model structure.

[0033] The following formula is the Eley-Rideal mechanism formula for NO and NH3 on the catalyst surface. When the concentrations of NO and NH3 are constant, the reaction rate can be calculated:

[0034] Where, θNH3: coverage of NH3 on the catalyst surface (ratio of adsorbed molecules to total active sites).

[0035] CNO: NO concentration in the gas phase. ka: NH3 adsorption rate constant. kd: NH3 desorption rate constant. kr: Chemical reaction rate constant under ER mechanism: The ammonia adsorption process is ; The ammonia desorption process is .

[0036] Perform optimization calculations based on historical operation data.

[0037] Specifically, the optimization process includes the following steps: 1) Initialization: Set the basic parameters of the algorithm (such as the number of iterations, population size, etc.) and generate the initial solution set; 2) Evaluation: Use the fitness function to evaluate the quality of the current solution set; 3) Search and generate new solutions: Based on the current solution, generate a new set of candidate solutions through strategies such as mutation, crossover, and neighborhood search; 4) Selection: Selectively replace or update the solution set according to the fitness function and retain the better solution; 5) Iteration: Repeat steps 2 to 4 until a termination condition is met; the termination condition is that the maximum number of iterations or the improvement of the solution is less than a threshold; 6) Output results: Output the optimal solution set found when the algorithm terminates.

[0038] The reactor inlet nitrogen oxide and reactor outlet nitrogen oxide data in the historical data include the instrument measurement lag time, which is used as a new dimension for comprehensive optimization and iterative calculation. Through empirical analysis, the accuracy of the SCR control simulation system is verified.

[0039] After the simulation system was constructed, the dynamic generation process of ammonia escape in the reactor under different ammonia-nitrogen molar ratios was tested based on the simulation system. As the ammonia-nitrogen molar ratio increased, ammonia escape would increase.

[0040] The simulation system of the present invention can optimize control strategy parameters, verify the robustness of the control method, and construct a target ammonia nitrogen molar ratio function model under different denitration efficiencies and a target ammonia nitrogen molar ratio function model under different flue gas temperatures. The change in the NOx concentration at the reactor inlet can be predicted based on the change in the NOx concentration at the reactor inlet and the change trend. The change in the NOx concentration at the reactor outlet can be predicted based on the change in the NOx concentration at the reactor outlet and the change trend.

[0041] Among them, the theoretical ammonia injection amount is calculated as (NOx concentration at reactor inlet - NOx concentration at reactor outlet + NOx concentration prediction variable at reactor inlet + NOx concentration prediction variable at reactor outlet) * flue gas flow rate * ammonia nitrogen molar ratio corrected by denitrification efficiency * ammonia nitrogen molar ratio corrected by flue gas temperature.

[0042] Based on the NOx concentration feedback value at the reactor outlet, the theoretical ammonia injection amount is corrected by the adaptive feedback controller. The adaptive feedback controller parameters constitute a fuzzy control function set based on the denitration efficiency, the deviation between the NOx concentration at the reactor outlet and the set value, and the rate of change of the deviation to achieve online dynamic correction.

[0043] Example 2 This embodiment is based on the construction method of the SCR control simulation system in Embodiment 1, wherein the delay inertia time of a simple CEMS instrument is greater than 4 minutes, and the system exhibits a large delay characteristic.

[0044] like Figure 1 The figure shows the NOx fluctuation at the reactor inlet.

[0045] In this embodiment, based on the collection and analysis of the actual historical big data of a power plant's operation, data cleaning and feature engineering are carried out, a big data analysis algorithm is constructed, the data regularity characteristics are mined, the reaction characteristics of the actual reactor of the original unit and the delay characteristics of the CEMS instrument are analyzed, and the initialization conditions are used as guidance for the scheme design.

[0046] Based on the Pearson correlation coefficient, a variable correlation analysis model is constructed, such as Figure 4This is the result of the variable correlation thermodynamic diagram based on the analysis of the historical operation data of the power plant. It can be seen that the load rate, ammonia nitrogen molar ratio, and reactor inlet flue gas temperature are the main factors affecting the reactor denitrification efficiency.

[0047] like Figure 5 As shown, based on time series data smoothing, a correlation analysis data matrix is ​​constructed, and the approximate delay characteristics of the original CEMS instrument are calculated. It is found that the delay inertia time of the pure CEMS instrument is greater than 4 minutes, and the system exhibits a large delay characteristic.

[0048] The basic structure of the simulation system is based on data analysis software, with modular underlying functions. It combines mechanisms, experience, and data models to build a high-fidelity simulation system to verify the reliability of the control strategy in advance, while providing certain verification and guidance for the transformation design.

[0049] The construction of the reaction simulation model of the denitrification reactor is the core. The dynamic process of the catalytic reduction reaction mechanism is based on the Eley-Rideal mechanism as the model structure. However, the dynamic reaction process of each reactor will change with the accumulation of usage time, the transformation of process equipment, the aging and replacement of catalysts, and other model parameters, such as Figure 6 As shown, the reaction mechanism model parameters are optimized and calculated. The nitrogen oxide data at the reactor inlet and the nitrogen oxide data at the reactor outlet in the historical data include the measurement lag time of the instrument, which is used as a new dimension for comprehensive optimization and iterative calculation.

[0050] Finally, the Gray Wolf Optimization Algorithm was selected as the parameter optimization method for the dynamic process model of the denitrification reaction. The prediction model in the control logic needs to consider the adaptive characteristics of the parameters. After data accumulation, the simulation model parameters can be iterated in the simulation system to predict the control effect.

[0051] After the simulation system was built, the dynamic process of denitrification reduction reaction in the reactor under different ammonia nitrogen molar ratios was tested based on the simulation system. The NOx concentration at the SCR outlet decreased with the increase of the ammonia nitrogen molar ratio, and the influence on the denitrification efficiency weakened with the increase of the ammonia nitrogen molar ratio.

[0052] Example 3 Based on the method for constructing the SCR control simulation system of embodiment 1, a system for constructing the SCR control simulation system is disclosed, including: A screening module is used to screen the factors affecting the denitration efficiency of the reactor by constructing a thermal map of parameter variables related to the denitration efficiency; The first calculation module is used to screen out the main influencing factors of the simulation system by calculating the correlation between the denitration efficiency of the reactor and the parameter variables; The second calculation module is used to obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the signal of the CEMS valve opening; The system building module is used to train the historical data of the main influencing factors of the simulation system through algorithms, establish the corresponding mathematical model, and adjust the parameters of the mathematical model through the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

[0053] Example 4 The purpose of this embodiment is to provide an electronic device, comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor implements the SCR control simulation system when executing the computer program.

[0054] Example 5 The purpose of this embodiment is to provide a computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the SCR control simulation system method is implemented.

[0055] Example 6 The purpose of this embodiment is to provide a computer program product including a computer-readable medium, on which a computer-readable program code is contained, and the program code executes the SCR control simulation system method.

[0056] The steps involved in the devices of the above embodiments 3, 4, 5 and 6 correspond to those of the method embodiment 1. For the specific implementation methods, please refer to the relevant description part of embodiment 1.

[0057] It should be understood by those skilled in the art that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes. The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks. These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes. These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0058] Unless otherwise specified, the working methods or control methods involved in the above embodiments are all conventional working methods or control methods in the art.

[0059] Although the preferred embodiments of the present application have been described, those skilled in the art may make other changes and modifications to these embodiments once they have learned the basic creative concept. Therefore, the appended claims are intended to be interpreted as including the preferred embodiments and all changes and modifications falling within the scope of the present application. Finally, it should be noted that the above embodiments are only used to illustrate the technical solution of the present invention rather than to limit it. Other modifications or equivalent substitutions made to the technical solution of the present invention by ordinary technicians in the field should be included in the scope of the claims of the present invention as long as they do not depart from the spirit and scope of the technical solution of the present invention.

Claims

1. A method for constructing an SCR control simulation system, characterized in that: The following steps are involved: S1: Screening the factors affecting the denitrification efficiency of the reactor by constructing a heat map of parameter variables related to the denitrification efficiency; S2: By calculating the correlation between the denitrification efficiency of the reactor and the parameter variables, the main influencing factors of the simulation system are screened; S3: Obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the CEMS valve opening signal; S4: The algorithm is used to train the historical data of the main influencing factors of the simulation system, and the corresponding mathematical model is established. The parameters of the mathematical model are adjusted according to the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

2. The method for constructing an SCR control simulation system according to claim 1, characterized in that: Before step S1, the NOx fluctuation data at the reactor inlet, the NOx fluctuation data at the reactor outlet, the NOx fluctuation data at different load rates and the flow characteristic curve data of the ammonia injection regulating valves on both sides are obtained.

3. The method for constructing an SCR control simulation system according to claim 1, characterized in that: In step S2, the correlation between the parameter variables is calculated using the Pearson correlation coefficient formula.

4. The method for constructing an SCR control simulation system according to claim 1, characterized in that: In step S3, based on the correlation between the valve action signal and the signal of the CEMS valve opening, the statistical change concentration is smoothed based on the time series data to construct a correlation analysis data matrix.

5. The method for constructing an SCR control simulation system according to claim 1, characterized in that: In step S4, the Grey Wolf Optimization Algorithm is used to train the simulation system.

6. The method for constructing an SCR control simulation system according to claim 1, characterized in that: After step S4, the dynamic process of the denitration reduction reaction in the reactor under different ammonia nitrogen molar ratios is tested based on the simulation system for verification.

7. A construction system for an SCR control simulation system, characterized in that: include: A screening module is used to screen the factors affecting the denitration efficiency of the reactor by constructing a thermal map of parameter variables related to the denitration efficiency; The first calculation module is used to screen out the main influencing factors of the simulation system by calculating the correlation between the denitration efficiency of the reactor and the parameter variables; The second calculation module is used to obtain the delay characteristic of the CEMS instrument by calculating the correlation between the valve action signal and the signal of the CEMS valve opening; The system building module is used to train the historical data of the main influencing factors of the simulation system through algorithms, establish the corresponding mathematical model, and adjust the parameters of the mathematical model through the delay characteristics of the CEMS instrument to obtain the SCR control simulation system.

8. An electronic device, characterized in that: The invention comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements the SCR control simulation system according to any one of claims 1 to 6 when executing the computer program.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the SCR control simulation system method according to any one of claims 1 to 6 is implemented.

10. A computer program product comprising a computer readable medium, characterized in that The computer-readable medium contains computer-readable program code, and the program code executes the SCR control simulation system method according to any one of claims 1 to 6.

Images