Intelligent micro-area gas circulation control method and system for disinfection cabinet

Through the intelligent micro-region gas circulation control method, abnormal areas in the disinfection cabinet are identified and optimized, and uniform distribution and stable circulation of air flow are achieved, which solves the problem of low gas circulation control efficiency in the prior art, and improves the disinfection effect and system reliability.

CN120195994AActive Publication Date: 2025-06-24SHENZHEN FUYUBAO KITCHEN EQUIP CO LTD
View PDF 6 Cites 0 Cited by

Patent Information

Application Number
CN202510656043.3
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2025-05-21
Publication Date
2025-06-24
Estimated Expiration
2045-05-21

AI Technical Summary

Technical Problem

During the disinfection process, existing disinfection cabinets cannot flexibly adjust the disinfection conditions in each micro-area according to actual needs, resulting in excessive gas concentration or improper flow rate, affecting the disinfection effect and possibly causing damage to the items. At the same time, global control strategies are difficult to identify and deal with abnormal situations, such as temperature and humidity exceeding the standard or gas leakage, resulting in inefficient gas circulation control.

Method used

The intelligent micro-region gas circulation control method is adopted to obtain disinfectant concentration data and airflow velocity data, conduct regional state analysis, identify abnormal areas, and perform airflow optimization to generate optimization control parameters to achieve uniform distribution and stable circulation of airflow. When the airflow equalization is lower than the preset threshold, secondary optimization is triggered to further adjust the airflow distribution plan.

Benefits of technology

It improves the efficiency and accuracy of gas circulation control of the disinfection cabinet, ensures that each micro-area achieves the best disinfection effect, protects the items from damage, and enhances the reliability and adaptability of the system. It is suitable for occasions with high requirements for disinfection quality.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN120195994A_ABST
    Figure CN120195994A_ABST
Patent Text Reader

Abstract

The invention relates to the technical field of disinfection cabinet control, and discloses an intelligent micro-area gas circulation control method and system for a disinfection cabinet, and the method comprises the steps: obtaining disinfectant concentration data and airflow velocity data; performing regional state analysis according to the disinfectant concentration data and the airflow velocity data to obtain an abnormal region distribution diagram; performing airflow optimization according to the abnormal region distribution diagram to obtain an optimization control parameter; a control instruction is generated according to the optimization control parameters, and the gas circulation state is obtained after gas flow circulation is stable; performing balance degree calculation according to the gas circulation state to obtain a gas flow balance degree; when the airflow balance degree is lower than a preset balance degree threshold value, secondary optimization is triggered, and a distribution scheme after secondary optimization is generated. The method has the following effect that the gas circulation control efficiency of the local area in the disinfection cabinet can be improved.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the technical field of disinfection cabinet control, and particularly to an intelligent micro-region gas circulation control method and system for a disinfection cabinet. Background Art

[0002] As an important device used in the household, medical, and catering industries, the main function of a disinfection cabinet is to sterilize and disinfect items through methods such as high temperature, ultraviolet light, or chemical gases. With the improvement of people's requirements for the quality of life and the increasing strictness of hygiene standards, traditional disinfection methods have become difficult to meet modern needs. Especially for different types of items, such as glassware, plastic products, and medical devices, they have different requirements for the disinfection environment, which has promoted the development of intelligent micro-region gas circulation control technology. This technology aims to precisely regulate the gas flow and temperature and humidity conditions in each region inside the disinfection cabinet to achieve the best disinfection effect while protecting the items from damage.

[0003] Existing disinfection cabinets adopt an integrated gas circulation system, and its operation process is roughly as follows: First, place the items to be disinfected into the disinfection cabinet and close the cabinet door; then, start the corresponding program according to the set disinfection mode, including steps such as heating and injecting disinfection gas (such as ozone); during the disinfection process, the built-in fan will promote the overall circulation of the gas inside the cabinet to ensure the uniform distribution of the disinfection gas.

[0004] However, the entire disinfection process of the existing technical solutions relies on unified gas circulation and concentration settings, and it is impossible to flexibly adjust the disinfection conditions in each micro-region according to actual needs, resulting in the disinfection effect being affected in some regions due to excessive gas concentration or improper flow rate, and even causing damage to the items. Further, the global control strategy is also difficult to effectively identify and handle the emerging abnormal situations. For example, in a relatively large disinfection space, problems such as excessive temperature and humidity or gas leakage may occur in a small area. If these problems are not discovered and solved in time, not only will the disinfection effect of this local area be reduced, but it may also lead to the imbalance of the internal environment of the entire disinfection cabinet. However, due to the effect of overall data averaging, these abnormalities in a small range are often masked and cannot be fully emphasized and processed in time. In summary, the global control strategy is difficult to identify abnormal regions, resulting in low control efficiency of gas circulation. Summary of the Invention

[0005] The present invention provides an intelligent micro-region gas circulation control method and system for a disinfection cabinet, in particular, to improve the gas circulation control efficiency of the disinfection cabinet.

[0006] In a first aspect, to solve the above technical problems, the present invention provides an intelligent micro-region gas circulation control method for a disinfection cabinet, including: Obtaining disinfectant concentration data and air flow velocity data; Perform regional status analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal area distribution map; Perform air flow optimization based on the abnormal area distribution map to obtain optimized control parameters; Generate a control instruction according to the optimized control parameters. After the air flow circulation is stable, obtain the gas circulation status; Perform equilibrium degree calculation based on the gas circulation status to obtain the air flow equilibrium degree; When the air flow equilibrium degree is lower than a preset equilibrium degree threshold, trigger secondary optimization and generate a secondary optimized distribution plan.

[0007] In an alternative embodiment, the performing regional status analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal area distribution map includes: Perform data cleaning on the disinfectant concentration data to obtain the cleaned cleaning concentration data; Perform variance calculation on the cleaned cleaning concentration data to obtain a concentration deviation; Mark the area where the concentration deviation is higher than a preset concentration deviation threshold as a concentration abnormal area; Perform data cleaning on the air flow velocity data to obtain the cleaned cleaning flow velocity data; Perform standard deviation calculation on the cleaned cleaning flow velocity data to obtain the air flow uniformity; Mark the area where the air flow uniformity is lower than a preset uniformity standard value as a flow velocity abnormal area; Mark the overlapping area of the concentration abnormal area and the flow velocity abnormal area as a comprehensive abnormal area; Generate a heat map based on the concentration abnormal area, the flow velocity abnormal area and the comprehensive abnormal area to obtain an abnormal area distribution map.

[0008] In an alternative embodiment, the performing air flow optimization based on the abnormal area distribution map to obtain optimized control parameters includes: Divide the abnormal area distribution map into grids according to a preset grid size to obtain micro areas; Perform particle swarm optimization on each micro area, including: Initialize the particle swarm size, inertia weight, learning factor and maximum number of iterations, where each particle represents a set of disinfection cabinet control parameters; Randomly initialize the current position and current velocity of each particle; Perform fitness calculation according to the current position and the current velocity to obtain the current fitness of the current particle; Update the individual optimal solution and the global optimal solution according to the current fitness; Multiply the inertial weight by the current velocity to obtain the inertial term; Calculate the cognitive term according to the learning factor, the individual optimal solution, and the current position; Calculate the social term according to the learning factor, the global optimal solution, and the current position; Perform weighted summation on the inertial term, the cognitive term, and the social term to obtain the updated velocity; Calculate the updated position according to the updated velocity; Use the updated velocity and the updated position as the current velocity and the current position respectively, repeat the iterative process, and record the number of iterations; When the number of iterations reaches the maximum number of iterations or fitness convergence is detected, complete the iterative process and output the global optimal solution as the optimized control parameter.

[0009] In an alternative embodiment, generating control instructions according to the optimized control parameter, and after the air flow circulation is stable, obtaining the gas circulation state includes: Extract the air flow distribution ratio of each micro-region from the optimized control parameter; Generate a control instruction set for the gas circulation device according to the air flow distribution ratio; Drive the intake valve and the outlet valve of each micro-region according to the control instruction set, and record the gas deviation; When the gas deviation is less than a preset stability threshold, it is determined that the air flow circulation is stable, and the gas circulation state is obtained; Wherein the gas circulation state includes the disinfectant concentration value and the air flow velocity value.

[0010] In an alternative embodiment, calculating the air flow balance degree according to the gas circulation state to obtain the air flow balance degree includes: Calculate the air flow balance degree through the following formula: Wherein, represents the air flow balance degree, represents the total number of micro-regions, represents the micro-region number, represents the th micro-region's disinfectant concentration value, represents the average value of the disinfectant concentration values of all micro-regions.

[0011] In an alternative embodiment, when the air flow balance degree is lower than a preset balance degree threshold, secondary optimization is triggered to generate a secondary optimized distribution plan, including: Initialize the population according to the preset population size, where the gene sequence of each individual is a three-dimensional vector, including three dimensions: fan speed, air valve opening, and disinfectant cycle; Calculate the fitness according to the following formula: Among them, represents the fitness, represents the air flow balance, represents the normalization factor, represents the anti-zero constant; Arrange in descending order according to the magnitude of the fitness, and retain individuals with a preset ratio from large to small to obtain the remaining individual population; Perform arithmetic crossover and non-uniform mutation on the remaining individual population to obtain a new generation of population, repeat the iterative process, and record the number of iterations; When the number of iterations reaches the preset upper limit of the number of iterations, obtain the final population; Extract the gene sequence of the individual with the highest fitness in the final population, and after boundary verification, output the allocation plan.

[0012] In an optional implementation manner, the calculating the current fitness of the current particle according to the current position and the current speed includes: Calculate the current fitness according to the following formula: Among them, represents the current fitness, represents the maximum flow rate, represents the minimum flow rate, represents the average flow rate, represents the target concentration, represents the measured concentration, represents the normalized concentration, represents the comprehensive energy consumption, represents the normalized energy consumption, represents the index weight coefficient.

[0013] In a second aspect, the present invention provides an intelligent micro-region gas circulation control system for a disinfection cabinet, including: A data acquisition module for acquiring disinfectant concentration data and air flow velocity data; A regional status module for performing regional status analysis according to the disinfectant concentration data and the air flow velocity data to obtain an abnormal region distribution map; An air flow optimization module for performing air flow optimization according to the abnormal region distribution map to obtain optimized control parameters; The initial control module is used to generate a control instruction according to the optimized control parameter, and after the air flow circulation is stable, obtain the gas circulation state; The equilibrium calculation module is used to calculate the equilibrium degree according to the gas circulation state to obtain the air flow equilibrium degree; The secondary optimization module is used to trigger secondary optimization and generate a secondary optimized distribution plan when the air flow equilibrium degree is lower than a preset equilibrium degree threshold.

[0014] In a third aspect, the present invention also provides an electronic device, including a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor. When the processor executes the computer program, it implements the intelligent micro-region gas circulation control method of the disinfection cabinet described in any one of the above.

[0015] In a fourth aspect, the present invention also provides a computer-readable storage medium. The computer-readable storage medium includes a stored computer program. When the computer program runs, it controls the device where the computer-readable storage medium is located to execute the intelligent micro-region gas circulation control method of the disinfection cabinet described in any one of the above.

[0016] Compared with the prior art, the intelligent micro-region gas circulation control method of the disinfection cabinet of the present invention has the following beneficial effects: (1) The process of obtaining the disinfectant concentration data and the air flow velocity data ensures a comprehensive understanding of the internal environment state of the disinfection cabinet. By installing high-precision sensors inside the disinfection cabinet to monitor the disinfectant concentration and air flow velocity in real time, accurate data support is provided for subsequent analysis. This step not only improves the accuracy of data collection but also lays a foundation for achieving precise control, which helps to improve the disinfection effect.

[0017] (2) Perform regional state analysis according to the disinfectant concentration data and the air flow velocity data to obtain an abnormal area distribution map. The disinfectant concentration and air flow velocity data are respectively processed by using variance and standard deviation calculations to identify abnormal areas with insufficient disinfectant concentration or uneven air flow inside the disinfection cabinet. This method can quickly locate the problem areas, provide a clear direction for subsequent optimized control, and ensure that the ideal disinfection effect can be achieved in every corner. It improves the recognition efficiency and accuracy for local areas.

[0018] (3) Perform air flow optimization according to the abnormal area distribution map to obtain optimized control parameters. Based on the identified abnormal area distribution map, the system automatically adjusts the working parameters (such as air speed, air direction, etc.) of the air supply device to improve the air flow distribution and eliminate the abnormal areas. This dynamic optimization mechanism can significantly improve the uniformity of the air flow, enabling the disinfectant to be evenly distributed throughout the disinfection cabinet, thereby improving the control efficiency of the overall air flow circulation.

[0019] (4) Generate a control instruction according to the optimized control parameter. After the air flow circulation stabilizes, obtain the gas circulation state. Generate a specific control instruction by using the optimized control parameter and send it to the control system of the disinfection cabinet for execution. After a period of time, when the air flow circulation tends to be stable, collect the gas circulation state data again to evaluate the current air flow condition. This process ensures the effective implementation of the optimization measures and provides a basis for further adjustment.

[0020] (5) Calculate the balance degree according to the gas circulation state to obtain the air flow balance degree. By analyzing the gas circulation state data in detail, calculate the balance degree index of the air flow, which quantitatively represents the uniformity of the air flow in the disinfection cabinet. This step helps to objectively evaluate the current air flow distribution and provides a scientific basis for judging whether further optimization is needed.

[0021] (6) When the air flow balance degree is lower than the preset balance degree threshold, trigger secondary optimization and generate a secondary optimized distribution plan. If the calculated air flow balance degree does not meet the preset standard, the system will start the secondary optimization process to readjust the parameters of the air supply device until the requirements are met. This method ensures the uniformity and stability of the air flow through multiple rounds of iterative optimization, and finally realizes the consistency of the disinfection effect in each micro-region inside the disinfection cabinet, improving the disinfection quality and efficiency.

[0022] In summary, through a series of precise data collection, analysis and control steps, the present invention realizes the intelligent management of the air flow inside the disinfection cabinet, not only improving the disinfection effect, but also enhancing the reliability and adaptability of the system, and is particularly suitable for occasions with high requirements for disinfection quality. Description of the Drawings

[0023] Figure 1 is a schematic flow chart of an intelligent micro-region gas circulation control method for a disinfection cabinet provided by the first embodiment of the present invention; Figure 2 is a schematic structural diagram of an intelligent micro-region gas circulation control system for a disinfection cabinet provided by the second embodiment of the present invention. Detailed Embodiments

[0024] The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0025] Refer to Figure 1, the first embodiment of the present invention provides an intelligent micro-region gas circulation control method for a disinfection cabinet, including the following steps: S11, obtain disinfectant concentration data and air flow velocity data; S12, perform regional state analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal region distribution map; S13, perform air flow optimization based on the abnormal region distribution map to obtain optimized control parameters; S14, generate a control instruction according to the optimized control parameters. After the air flow circulation is stable, obtain the gas circulation state; S15, perform balance calculation according to the gas circulation state to obtain the air flow balance degree; S16, when the air flow balance degree is lower than a preset balance degree threshold, trigger secondary optimization and generate a secondary optimized distribution plan.

[0026] In step S11, obtain disinfectant concentration data and air flow velocity data.

[0027] In one implementation, the disinfectant concentration data is collected in real time by electrochemical sensors or infrared spectroscopy sensors distributed on the top, bottom, and side walls of the disinfection cabinet. The concentration unit is ppm. The air flow velocity data is obtained by hot-wire anemometers or ultrasonic velocimeters installed at the inlet, outlet, and key micro-regions of the circulation air duct. The velocity unit is uniformly m / s (meters per second). After the two types of data are converted by analog-to-digital conversion, they are transmitted to the main control module through the CAN bus or RS-485 communication protocol, stored as a structured data table indexed by time stamps, and recorded in CSV format for the spatial coordinates, concentration values (accurate to 0.1 ppm), wind speed values (accurate to 0.01 m / s), and acquisition time (millisecond-level accuracy) of each sensor node. At the same time, the CRC check code is used to ensure the integrity of data transmission, and all numerical values are attached with unit identifiers to avoid the risk of dimension confusion.

[0028] In step S12, perform regional state analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal region distribution map.

[0029] In one implementation, perform data cleaning on the disinfectant concentration data to obtain the cleaned concentration data; Perform variance calculation on the cleaned concentration data to obtain the concentration deviation; Mark the region where the concentration deviation is higher than the preset concentration deviation threshold as a concentration abnormal region; Perform data cleaning on the air flow velocity data to obtain the cleaned flow velocity data; Perform standard deviation calculation on the cleaned flow velocity data to obtain the air flow uniformity; Marking the area where the airflow uniformity is lower than a preset uniform standard value as an abnormal flow rate area; marking the overlapping area of ​​the abnormal concentration area and the abnormal velocity area as a comprehensive abnormal area; A thermal map is generated according to the abnormal concentration area, the abnormal velocity area and the comprehensive abnormal area to obtain an abnormal area distribution map.

[0030] In one embodiment, in the data cleaning implementation, the disinfectant concentration data is first subjected to sliding average filtering (the window width is set to 5 sampling points) to eliminate sensor noise, and outliers are identified by binning: the 5%-95% percentile interval of the concentration data is calculated, and data points outside the range are eliminated, and zero-value sequences with a duration of more than 3 seconds are marked as sensor failure states; for the airflow velocity data, a median filter (window width of 3 points) is applied to suppress pulse interference, combined with physical constraint filtering (such as when the flow rate exceeds the maximum design value of the circulation system of 3m / s, it is considered abnormal), and cubic spline interpolation is used to fill in the single point missing caused by sensor jitter; both types of data are normalized after preprocessing (the concentration value is mapped to the dynamic range of 0-50ppm, and the flow rate value is mapped to the interval of 0-2.5m / s), and finally the concentration and flow rate data sets are aligned at the millisecond level through the timestamp synchronization module, and the sampling time deviation is compensated by linear interpolation. The cleaned data is stored in a relational database with version control, and the original data and the cleaned data are bidirectionally traceable through the MD5 checksum.

[0031] In a specific embodiment, for the disinfectant concentration data (5 micro-area values ​​after cleaning: 12.3ppm, 14.8ppm, 11.5ppm, 13.2ppm, 25.6ppm), the mean μ=15.48ppm and the variance σ²=28.7 ppm² are calculated. According to the GB / T 18883-2022 standard, 1.5 times the standard deviation (8.04ppm) is set as the concentration deviation threshold, where the concentration deviation of area 5 reaches 10.12ppm (25.6ppm) to trigger an abnormal mark; for the air flow velocity data (measurement point values ​​after cleaning: 1.8m / s, 1.9m / s, 2.0m / s, 2.1m / s, 0.5m / s), its standard deviation σ=0.62m / s (μ=1.66m / s) results in a coefficient of variation CV=0.37, which is significantly higher than ISO The Grade II uniformity standard (CV≤0.25) specified in 5801:2017, especially the 0.5m / s low-speed airflow in area 5, is the main reason for the uniformity exceeding the standard.

[0032] In one implementation, the generation process of the abnormal area distribution map includes the following steps: First, based on the concentration detection data, the concentration abnormal areas deviating from the mean by 1.5 times the standard deviation (e.g., >8.04 ppm in the example) are screened out through variance calculation, and at the same time, the flow velocity abnormal areas are identified according to the standard deviation and coefficient of variation of the air flow velocity (e.g., CV>0.25); Subsequently, through spatial overlay analysis, the vector layers of the two types of abnormal areas are overlaid in a geographic information system (the GeoPandas library in Python can be used), and the spatial intersection is extracted as the comprehensive abnormal area; Finally, using heat map visualization technology, in Python's Matplotlib, different color codings are assigned to the three types of abnormal areas respectively (e.g., red gradient indicates concentration abnormality, blue gradient indicates flow velocity abnormality, and purple indicates comprehensive abnormality), a heat distribution map with transparency overlay effect is generated, and spatial coordinate grids and legend descriptions are embedded, and finally output as a high-resolution vector map or an interactive Web map (implemented by Leaflet.js).

[0033] In step S13, the air flow is optimized according to the abnormal area distribution map to obtain the optimized control parameters.

[0034] In one implementation, the abnormal area distribution map is divided into grids according to a preset grid size to obtain micro areas; Particle swarm optimization is performed on each of the micro areas, including: Initialize the particle swarm size, inertia weight, learning factor, and maximum number of iterations, where each particle represents a set of disinfection cabinet control parameters; Randomly initialize the current position and current velocity of each particle; Calculate the fitness according to the current position and the current velocity to obtain the current fitness of the current particle; Update the individual optimal solution and the global optimal solution according to the current fitness; Multiply the inertia weight by the current velocity to obtain the inertia term; Calculate the cognitive term according to the learning factor, the individual optimal solution, and the current position; Calculate the social term according to the learning factor, the global optimal solution, and the current position; Perform weighted summation on the inertia term, the cognitive term, and the social term to obtain the updated velocity; Calculate the updated position according to the updated velocity; Use the updated velocity and the updated position as the current velocity and the current position respectively, repeat the iteration process, and record the number of iterations; When the number of iterations reaches the maximum number of iterations, or fitness convergence is detected, the iteration process is completed, and the global optimal solution is output as the optimized control parameter.

[0035] It should be noted that in the grid division implementation for optimizing the air flow in the disinfection cabinet, taking the typical disinfection chamber size of 2.4m×1.8m as an example, the uniform grid division method can be adopted: the preset grid size is 0.3m×0.3m. The chamber is longitudinally divided into 8 columns (2.4m / 0.3m) and horizontally divided into 6 rows (1.8m / 0.3m), forming 48 rectangular micro-regions (numbered G01 - G48). The physical coordinate range corresponding to each grid is, for example, the G15 grid corresponding to the region of the X-axis from 1.2 to 1.5m and the Y-axis from 0.6 to 0.9m. When dividing, it is preferred to equally divide at equal intervals along the length and width directions of the chamber. If the size of the boundary region is less than 0.3m (such as the remaining 0.2m in the end grid), it is extended to the adjacent grid for combined processing to ensure full coverage without gaps. During anomaly detection, if a certain micro-region simultaneously satisfies the standard deviation of flow velocity ≥ 0.15m / s and the concentration deviation ≥ 8ppm, it is marked as an abnormal grid. For example, if the measured flow velocity fluctuation of the G22 grid is 0.18m / s and the concentration deviates from the target value by 12ppm, the particle swarm optimization process of this grid is triggered.

[0036] In one implementation, the current fitness is calculated according to the following formula: where, represents the current fitness, represents the maximum flow velocity, represents the minimum flow velocity, represents the average flow velocity, represents the target concentration, represents the measured concentration, represents the normalized concentration, represents the comprehensive energy consumption, represents the normalized energy consumption, represents the index weight coefficient.

[0037] It should be noted that in the fitness calculation formula, the units and physical meanings of each parameter are as follows: the units of the maximum flow velocity, minimum flow velocity, and average flow velocity are meters per second (m / s), the units of the target concentration and measured concentration are ppm, the normalized value is 1 ppm for dimensional normalization, the unit of the comprehensive energy consumption is watt-hour (Wh), the normalized energy consumption value is 100 watt-hour for dimensional normalization, and the weight coefficient is a unitless scaling factor (which can take values of 0.4, 0.3, 0.3 and satisfies the sum of 1). The fitness function is the core evaluation index of the optimization objective in the particle swarm algorithm. Its role is to guide the particle swarm to search for the optimal solution in the parameter space by quantifying the comprehensive performance of flow velocity uniformity, concentration control accuracy and energy consumption efficiency. For example, when a certain set of control parameters results in small flow velocity fluctuations (the first term approaches 0), low concentration deviation (the second term approaches 0) and low energy consumption (the third term approaches 0), the fitness value f will be minimized. At this time, the parameter combination corresponding to this particle is the optimization direction, and the algorithm converges the population to the global optimal solution through iteration.

[0038] It should be noted that in the implementation of particle swarm optimization for the control parameters of the disinfection cabinet, the algorithm parameters are first set as follows: the particle swarm size is set to 20 - 50 (e.g., 30), the inertia weight w ∈ [0.4, 1.2] (typical value 0.8), the individual and social learning factors are set to 1.5 and 2.0 respectively, and the maximum number of iterations is set to 100; each particle is encoded by a three-dimensional vector, representing wind speed (0.5 - 3.0 m / s), wind direction (0 - 360°) and injection volume (50 - 200 ml / min) respectively. For example, the initial parameters of randomly generated particle B are [2.1 m / s, 45°, 180 ml / min], and its velocity vector is initialized to [-0.1 m / s, +15°, -10 ml / min].

[0039] It should be noted that in each iteration, CFD (Computational Fluid Dynamics) software is used for dynamic simulation: ANSYS Fluent is used to establish a three-dimensional disinfection chamber model (size 2.4m×1.8m×1.2m), the multiphase flow coupling (air-disinfectant mixture) and the k-epsilon turbulence model are set, and the transient solver (time step 0.1 second) is used to simulate the air flow field distribution. The parameter combination of each particle (such as [2.3m / s, 75°, 170ml / min]) is input as the boundary condition, and the software automatically outputs the extreme values of the flow velocity, concentration and total system power consumption in each micro-region; among them, the wind speed is controlled by the inlet velocity boundary condition, the concentration is calculated through the component transport model, and the energy consumption is based on the integral of the fan power curve, all of which are automatically calculated within the simulation software. After the simulation data is extracted by a Python script, it is substituted into the fitness formula for calculation (such as f = 0.445), and the particle update logic is triggered. For example, when the simulation parameters of COMSOL Multiphysics are [2.2m / s, 90°, 180ml / min], the standard deviation of the flow velocity of the G22 grid is detected to be 0.12m / s (lower than the threshold of 0.15m / s) and the concentration deviation is 6ppm (meeting the standard) through the post-processing module, then this parameter combination is marked as a candidate optimal solution. The entire simulation process runs in parallel on the HPC cluster, and the time-consuming for a single iteration is about 3 minutes (based on 480,000 grids and 32 CPU cores), and the final output parameters are transmitted to the virtual PLC for verification through the OPC UA protocol. The obtained flow velocity, concentration and power consumption data are substituted into the fitness formula (the smaller the value, the better) for calculation; in the update rule, if the fitness of particle C is 0.38 at the 20th iteration, which is better than its historical best of 0.42, then the parameter [2.3m / s, 75°, 170ml / min] is updated to the individual best, and at the same time, if this value is lower than the global best of 0.40, the global best is replaced.

[0040] It should be noted that the cognitive term calculation formula is: Among them, represents the cognitive term, represents the individual learning factor, represents the random number, represents the individual optimal solution, represents the current position; The social term calculation formula is: Among them, represents the social term, represents the social learning factor, represents the random number, represents the global optimal solution, represents the current position.

[0041] It should be noted that the random number takes values from 0 to 1 to prevent premature convergence.

[0042] It should be noted that in the particle swarm optimization algorithm, the velocity update is achieved by integrating three factors: First, the original motion trend of the particle is retained (inertia term, weight 0.6), second, the attraction of the particle's own historical best position is introduced (cognitive term, weight 1), and finally, the guiding effect of the global optimal position of the group is superimposed (social term, weight 1).

[0043] It should be noted that in the particle swarm optimization algorithm, the calculation of the new position is achieved by vectorially superimposing the updated velocity vector and the current position. Specifically, the position update of each particle can be decomposed into the following steps: First, algebraically add the current parameter value (such as the wind speed value in the disinfection cabinet case) and the velocity increment of the corresponding dimension (i.e., the component of the velocity vector in that dimension). For example, when the current wind speed is 2.0 m / s and the velocity increment is +0.3 m / s, the new position wind speed is updated to 2.3 m / s. Second, in the multi-dimensional parameter space, such addition operations are performed independently for each dimension (such as wind speed, temperature, disinfection time, etc.) to form a new parameter combination. Finally, the system will detect whether the new position exceeds the preset parameter range (such as the wind speed is limited between 0.5 - 3.0 m / s). If it exceeds, the boundary handling mechanism will be triggered - for example, when 3.2 m / s is calculated, it is directly truncated to 3.0 m / s (hard boundary) or reflected back to 2.8 m / s according to the exceeding ratio (soft boundary). This process ensures that the movement of particles in the solution space follows both the dynamic laws and meets the engineering feasibility through the step-by-step superposition of parameter dimensions and boundary constraints.

[0044] In step S14, a control instruction is generated according to the optimized control parameter. After the air flow circulation is stable, the gas circulation state is obtained.

[0045] In one implementation, the air flow distribution ratio of each micro-region is extracted from the optimized control parameter; A control instruction set for the gas circulation device is generated according to the air flow distribution ratio; According to the control instruction set, the intake valves and outlet valves of each micro-region are driven, and the gas deviation amount is recorded; When the gas deviation amount is less than the preset stability threshold, it is determined that the air flow circulation is stable, and the gas circulation state is obtained; Wherein the gas circulation state includes the disinfectant concentration value and the air flow velocity value.

[0046] It should be noted that the optimized control parameters are encoded by a three-dimensional vector, which respectively represent the wind speed (0.5 - 3.0 m / s), wind direction (0 - 360°), and injection volume (50 - 200 ml / min). A control instruction set including parameters such as valve opening and timing is generated, and then the intake and exhaust valves of each micro-region are driven to perform dynamic adjustment. After the valve action, the system real-time monitors the gas distribution deviation in each region (such as pressure difference, flow difference, etc.). When the deviation is lower than the preset stability threshold for 3 consecutive sampling periods (for example, ±5% relative error), it is determined that the air flow circulation reaches a stable state. At this time, the measured values of the disinfectant concentration sensor and the wind speed sensor are synchronously collected to form gas circulation state data including the spatial concentration gradient distribution and the real-time flow velocity matrix, providing input for subsequent disinfection efficiency evaluation. When the deviation exceeds the stability threshold for more than 50 sampling periods, it is determined that the air flow cannot reach a stable state, and data is re-collected and the control parameters are optimized.

[0047] It should be noted that the three-dimensional vector (wind speed, wind direction, injection volume) of the optimized control parameters generates a specific control instruction set through the following logic: First, the wind speed parameter (0.5 - 3.0 m / s) is mapped and converted into a fan speed instruction. For example, the speed of the fan motor is controlled through a PWM (pulse width modulation) signal to ensure that the actual wind speed matches the optimized value; the wind direction parameter (0 - 360°) drives the air deflector motor to adjust the angle, converting the abstract angle value into a mechanical displacement (such as 0° corresponding to the horizontal direction and 90° corresponding to vertically upward), and dynamically adjusting the opening of the intake / exhaust valves in each region in combination with the air flow distribution ratio (such as when the wind direction in a certain region needs to be biased downward, the angle of the lower air deflector and the opening of the intake valve are linked). The injection volume parameter (50 - 200 ml / min) controls the opening and closing timing and pulse width of the solenoid valve of the injector, converting the flow demand into a duty cycle instruction for the valve. On this basis, the system further generates a timing table for valve actions according to the disinfection task stage (such as pre-disinfection, main cycle, homogenization stage). For example, in the homogenization stage, the exhaust valve in the high-concentration region is preferentially opened to balance the distribution, and at the same time, the execution order and duration of each instruction are recorded. Finally, it is integrated into a complete control instruction set including valve opening ratio, motor angle value, solenoid valve pulse parameters, and execution timing, ensuring that the hardware device can accurately execute the optimized air flow and concentration control strategy.

[0048] In step S15, the air flow balance degree is calculated according to the gas circulation state to obtain the air flow balance degree.

[0049] In one implementation, the air flow balance degree is calculated by the following formula: Where, represents the air flow balance degree, represents the total number of micro-regions, Indicates the micro-region number, Indicates the disinfectant concentration value of the th micro-region,

[0050] It should be noted that the calculation process of the air flow balance is based on the real-time monitoring data of each micro-region: First, the measured value of the current disinfectant concentration is obtained through the concentration sensors distributed in each micro-region, and the arithmetic mean of the concentration values of all micro-regions (the total number is determined by the number of disinfection partitions preset by the system) is calculated to obtain the reference concentration level; Subsequently, the square of the deviation between the concentration value of each partition and the average value is calculated and summed, and the sum is divided by the total number of partitions to obtain the variance characterization quantity of the concentration dispersion degree, and the standard deviation of the concentration distribution is obtained by taking the square root of the variance; Finally, the standard deviation is divided by the average concentration value to obtain a dimensionless air flow balance index. The smaller this value is, the more uniform the concentration distribution in each region is. When it is lower than 0.15, it is considered to reach an excellent balance state.

[0051] In step S16, when the air flow balance is lower than the preset balance threshold, secondary optimization is triggered to generate an optimized allocation plan.

[0052] In one implementation, the population is initialized according to the preset population size, where each individual gene sequence is a three-dimensional vector, including three dimensions: fan speed, air valve opening, and disinfectant cycle; Calculate the fitness according to the following formula: where, represents the fitness, represents the air flow balance, represents the normalization factor, represents the anti-zero constant; Arrange in descending order according to the magnitude of the fitness, and retain individuals with a preset proportion from largest to smallest to obtain the remaining individual population; Perform arithmetic crossover and non-uniform mutation on the remaining individual population to obtain a new generation of population, repeat the iterative process, and record the number of iterations; When the number of iterations reaches the preset upper limit of the number of iterations, the final population is obtained; Extract the individual gene sequence with the highest fitness in the final population, and after boundary verification, output the allocation plan.

[0053] It should be noted that during population initialization, the system randomly generates three-dimensional parameter combinations according to the preset feasible range: for example, the value range of the fan speed dimension is 30% - 100% (such as individual A taking 75% corresponding to medium-speed operation), the air valve opening dimension is randomly generated between 20% - 80% (such as individual B setting 60% to reserve adjustment margin), and the disinfectant spraying period dimension is limited to the range of 120 - 600 seconds (such as individual C configuring 300 seconds to balance continuous disinfection and energy consumption). This design avoids the risk of shutdown by restricting the minimum fan speed, limits the upper limit of the air valve opening to prevent airflow overshoot, and at the same time ensures that the disinfection interval meets the air treatment requirements specified in GB / T 18883 - 2022, ensuring both the diversity of the initial population and meeting the boundary conditions for the safe operation of the equipment. The population size is set to 100.

[0054] It should be noted that the anti-removal zero constant takes the value of 10^-5, a dimensionless constant. The normalization factor takes the value of 1.5, dimensionless. The fitness is dimensionless. The air flow balance degree is dimensionless. The upper limit of the number of iterations is set to 100 times.

[0055] It should be noted that during the population iteration optimization process, the coordinated execution of arithmetic crossover and non-uniform mutation follows the following engineering-oriented operation chain: First, use the tournament selection method to select the top 30% of the high-quality individuals with the highest fitness from the parent population as the crossover base group, and perform a convex combination operation on the three-dimensional gene vectors (fan speed, air valve opening, disinfection period) of each pair of parent individuals (such as the offspring individual rc = αrf+(1-α)rm, where α is a dynamic weight coefficient between [0.3, 0.7], rc represents the offspring gene vector, rf represents the father gene vector, and rm represents the mother gene vector); then perform non-uniform Gaussian mutation on the offspring generated by the crossover - apply a normal perturbation with a standard deviation of σ = 2% in the fan speed dimension (ensuring that the speed fluctuation does not exceed the safety margin of ±5%), and the air valve opening accepts a jump-type adjustment of ±8% - ±15% according to the probability inversely proportional to the fitness (the mutation amplitude of individuals with low fitness is larger); after each round of iteration, eliminate illegal individuals that violate the lower limit of the disinfection period (<120s) or fan overclocking (>100%) through the Pareto front analysis method, and finally generate a new generation of population with both inheritance and diversity until the number of iterations reaches the preset 50 generations or the optimal fitness fluctuates continuously for 5 generations <0.1% to terminate the optimization.

[0056] It should be noted that during the implementation of the boundary check, static range truncation is performed on the three-dimensional parameters (such as the fan speed is forcibly limited to the range of 30% - 100%, if 115% exists in the gene sequence, it is corrected to 100% and an overlimit alarm is triggered).

[0057] In summary, the present invention discloses an intelligent micro-region gas circulation control method for a disinfection cabinet. By monitoring local abnormal regions of the disinfection cabinet, the particle swarm optimization algorithm is first used to optimize control parameters, and software simulation is used to monitor the results. For the second optimization, actual measurement data is combined, and the genetic algorithm is used to optimize control parameters to improve the control efficiency of gas circulation.

[0058] Referring Figure 2 , the second embodiment of the present invention provides an intelligent micro-region gas circulation control system for a disinfection cabinet, including: A data acquisition module for acquiring disinfectant concentration data and air flow velocity data; A region status module for performing region status analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal region distribution map; An air flow optimization module for optimizing the air flow based on the abnormal region distribution map to obtain optimized control parameters; A primary control module for generating a control instruction based on the optimized control parameters, and after the air flow circulation is stable, acquiring the gas circulation status; An equilibrium calculation module for calculating the equilibrium degree based on the gas circulation status to obtain the air flow equilibrium degree; A secondary optimization module for triggering secondary optimization and generating a secondary optimized distribution plan when the air flow equilibrium degree is lower than a preset equilibrium degree threshold.

[0059] It should be noted that the intelligent micro-region gas circulation control system for a disinfection cabinet provided in the embodiment of the present invention is used to execute all the process steps of the intelligent micro-region gas circulation control method for a disinfection cabinet in the above embodiment. The working principles and beneficial effects of the two correspond one by one, and thus will not be elaborated here.

[0060] The embodiment of the present invention also provides an electronic device. The electronic device includes: a processor, a memory, and a computer program stored in the memory and executable on the processor, such as a data acquisition program. When the processor executes the computer program, the steps in the above embodiments of the intelligent micro-region gas circulation control method for a disinfection cabinet are implemented, such as Figure 1 The step S11 shown. Or, when the processor executes the computer program, the functions of each module / unit in the above device embodiments are implemented, such as the data acquisition module.

[0061] Exemplarily, the computer program may be divided into one or more modules / units, and the one or more modules / units are stored in the memory and executed by the processor to implement the present invention. The one or more modules / units may be a series of computer program instruction segments capable of performing specific functions, and these instruction segments are used to describe the execution process of the computer program in the electronic device.

[0062] The electronic device may be a computing device such as a desktop computer, a notebook, a palm computer, and a smart tablet. The electronic device may include, but is not limited to, a processor and a memory. Those skilled in the art can understand that the above components are only examples of the electronic device and do not constitute a limitation on the electronic device. It may include more or fewer components than the above, or combine some components, or different components. For example, the electronic device may further include input / output devices, network access devices, a bus, etc.

[0063] The so-called processor may be a central processing unit (CPU), or may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor may be a microprocessor or the processor may also be any conventional processor, etc. The processor is the control center of the electronic device and connects various parts of the entire electronic device through various interfaces and circuits.

[0064] The memory can be used to store the computer program and / or modules. By running or executing the computer program and / or modules stored in the memory, and invoking the data stored in the memory, the processor realizes various functions of the electronic device. The memory mainly includes a program storage area and a data storage area. Among them, the program storage area can store an operating system, application programs required for at least one function (such as a sound playback function, an image playback function, etc.); the data storage area can store data created according to the use of the mobile phone (such as audio data, phone book, etc.). In addition, the memory can include high-speed random access memory, and can also include non-volatile memory, such as a hard disk, memory, plug-in hard disk, smart media card (SMC), secure digital (SD) card, flash card, at least one magnetic disk storage device, flash memory device, or other volatile solid-state storage devices.

[0065] Among them, if the modules / units integrated in the electronic device are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on such an understanding, to implement all or part of the processes in the above-mentioned embodiment methods of the present invention, it can also be completed by instructing relevant hardware through a computer program. The computer program can be stored in a computer-readable storage medium. When the computer program is executed by the processor, the steps of the above-mentioned various method embodiments can be realized. Among them, the computer program includes computer program code, and the computer program code can be in the form of source code, object code, executable file, or some intermediate form, etc. The computer-readable medium can include: any entity or device capable of carrying the computer program code, recording medium, USB flash drive, mobile hard disk, magnetic disk, optical disc, computer memory, read-only memory (ROM), random access memory (RAM), electrical carrier signal, telecommunication signal, and software distribution medium, etc. It should be noted that the content included in the computer-readable medium can be appropriately increased or decreased according to the requirements of legislation and patent practice in the jurisdiction. For example, in some jurisdictions, according to legislation and patent practice, the computer-readable medium does not include electrical carrier signals and telecommunication signals.

[0066] It should be noted that the device embodiments described above are merely illustrative. The units described as separate components may or may not be physically separated, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed to multiple network units. Some or all of the modules can be selected according to actual needs to achieve the purpose of the solution of this embodiment. In addition, in the attached drawings of the device embodiments provided by the present invention, the connection relationship between the modules indicates that there is a communication connection between them, which can be specifically implemented as one or more communication buses or signal lines. Those of ordinary skill in the art can understand and implement it without creative efforts.

[0067] The specific embodiments described above further elaborate on the purpose, technical solution, and beneficial effects of the present invention. It should be understood that the above are only specific embodiments of the present invention and are not used to limit the protection scope of the present invention. In particular, for those skilled in the art, any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.

Claims

1. An intelligent micro-area gas circulation control method for a disinfection cabinet, characterized in that: include: Obtain disinfectant concentration data and air flow velocity data; Performing regional status analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal area distribution map; Perform airflow optimization according to the abnormal area distribution map to obtain optimized control parameters; Generate a control instruction according to the optimized control parameter, and obtain the gas circulation state after the gas circulation is stable; Calculate the balance degree according to the gas circulation state to obtain the airflow balance degree; When the airflow balance is lower than a preset balance threshold, secondary optimization is triggered to generate a secondary optimized distribution plan.

2. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 1 is characterized in that Performing regional status analysis based on the disinfectant concentration data and the airflow velocity data to obtain an abnormal region distribution map includes: Perform data cleaning according to the disinfectant concentration data to obtain cleaning concentration data after cleaning; Perform variance calculation based on the cleaning concentration data to obtain concentration deviation; Marking the area where the concentration deviation is higher than a preset concentration deviation threshold as an abnormal concentration area; Performing data cleaning according to the air flow velocity data to obtain cleaning flow velocity data after cleaning; Calculate the standard deviation based on the cleaning flow rate data to obtain air flow uniformity; Marking the area where the airflow uniformity is lower than a preset uniform standard value as an abnormal flow rate area; marking the overlapping area of ​​the abnormal concentration area and the abnormal velocity area as a comprehensive abnormal area; A thermal map is generated according to the abnormal concentration area, the abnormal velocity area and the comprehensive abnormal area to obtain an abnormal area distribution map.

3. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 1, characterized in that: The step of optimizing the airflow according to the abnormal area distribution map to obtain optimized control parameters includes: Divide the abnormal area distribution map into grids according to a preset grid size to obtain micro-areas; Performing particle swarm optimization on each of the micro-regions, including: Initialize the particle swarm size, inertia weight, learning factor and maximum number of iterations, where each particle represents a set of disinfection cabinet control parameters; Randomly initialize the current position and current velocity of each particle; Calculating the fitness according to the current position and the current speed to obtain the current fitness of the current particle; Update the individual optimal solution and the global optimal solution according to the current fitness; Obtain an inertia term by multiplying the inertia weight by the current speed; Calculate a cognitive item according to the learning factor, the individual optimal solution and the current position; Calculate a social item according to the learning factor, the global optimal solution and the current position; Performing weighted summation according to the inertia item, the cognitive item and the social item to obtain an update speed; Calculate the update position according to the update speed; Using the updated speed and the updated position as the current speed and the current position respectively, repeating the iterative process, and recording the number of iterations; When the number of iterations reaches the maximum number of iterations, or fitness convergence is detected, the iteration process is completed and the global optimal solution is output as the optimization control parameter.

4. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 1, characterized in that: The step of generating a control instruction according to the optimized control parameter and obtaining a gas circulation state after the gas circulation is stable comprises: Extracting the airflow distribution ratio of each micro-area from the optimized control parameters; generating a control instruction set for a gas circulation device according to the airflow distribution ratio; driving the air inlet valve and the air outlet valve of each micro-area according to the control instruction set, and recording the gas deviation; When the gas deviation is less than a preset stability threshold, it is determined that the airflow cycle is stable, and the gas circulation state is obtained; The gas circulation state includes a disinfectant concentration value and an air flow velocity value.

5. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 1, characterized in that: The calculating the balance degree according to the gas circulation state to obtain the airflow balance degree includes: The air flow balance is calculated using the following formula: in, Indicates the air flow balance. represents the total number of micro-regions, Indicates the micro-area number, Indicates The disinfectant concentration value of each micro-area, Represents the average of the disinfectant concentration values ​​of all micro-areas.

6. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 1, characterized in that: When the airflow balance is lower than a preset balance threshold, secondary optimization is triggered to generate a secondary optimized distribution plan, including: Initialize the population according to the preset population size, where each individual gene sequence is a three-dimensional vector, including three dimensions: fan speed, valve opening, and disinfectant cycle; The fitness is calculated according to the following formula: in, represents fitness, Indicates the air flow balance. represents the normalization factor, Indicates the zero-proof constant; Arrange the individuals in descending order according to the fitness, retain a preset proportion of individuals from large to small, and obtain the remaining individual population; Performing arithmetic crossover and non-uniform mutation on the remaining individual population to obtain a new generation of population, repeating the iterative process, and recording the number of iterations; When the number of iterations reaches a preset upper limit of the number of iterations, a final population is obtained; The individual gene sequence with the highest fitness in the final population is extracted, and after boundary verification, an allocation plan is output.

7. The intelligent micro-region gas circulation control method of the disinfection cabinet according to claim 3, characterized in that: The step of calculating the fitness according to the current position and the current speed to obtain the current fitness of the current particle includes: The current fitness is calculated according to the following formula: in, represents the current fitness, Indicates the maximum flow rate, Indicates the minimum flow rate, represents the average flow velocity, represents the target concentration, represents the measured concentration, represents the normalized concentration, represents the comprehensive energy consumption, represents the normalized energy consumption, Represents the indicator weight coefficient.

8. An intelligent micro-area gas circulation control system for a disinfection cabinet, characterized in that: include: A data acquisition module, used to acquire disinfectant concentration data and air flow velocity data; A regional status module, used to perform regional status analysis based on the disinfectant concentration data and the air flow velocity data to obtain an abnormal region distribution map; An airflow optimization module, used to optimize the airflow according to the abnormal area distribution map to obtain optimized control parameters; A primary control module, used to generate a control instruction according to the optimized control parameters, and obtain a gas circulation state after the gas circulation is stable; A balance calculation module, used to calculate the balance degree according to the gas circulation state to obtain the airflow balance degree; The secondary optimization module is used to trigger the secondary optimization when the airflow balance is lower than a preset balance threshold, and generate a secondary optimized distribution plan.

9. An electronic device, characterized in that: It comprises a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, and when the processor executes the computer program, it implements the intelligent micro-area gas circulation control method of the disinfection cabinet as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: The computer-readable storage medium includes a stored computer program, wherein when the computer program is running, the device where the computer-readable storage medium is located is controlled to execute the intelligent micro-area gas circulation control method of the disinfection cabinet as described in any one of claims 1 to 7.

Citation Information

Patent Citations

  • Positioning method for hazardous gas leakage diffusion accident source

    CN104834017A

  • Gas circulation control method and system of disinfection cabinet, electronic equipment and medium

    CN116099023A

  • Automatic management and safety early warning method and system for gas purification type medicine storage cabinet

    CN117744934A

  • Intelligent control method for disinfection cabinet

    CN118436833A

  • System of sterilizing and cleanjng pollution air at kitchen

    KR102548992B1