A method and system for visualizing three-dimensional diffusion of an organic acid

By combining distributed sensors and edge computing nodes with an improved multi-parameter lookup table method and trapezoidal numerical integration method, the problems of user perception difficulty, high cost, and hardware deployment difficulties in monitoring organic acids in museums have been solved, realizing low-cost display and dynamic early warning of cultural relic corrosion risks.

CN120995928BActive Publication Date: 2026-08-04NANJING YRD ECO DEV RI CO LTD
View PDF 2 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Patents(China)
Current Assignee / Owner
NANJING YRD ECO DEV RI CO LTD
Filing Date
2025-08-07
Publication Date
2026-08-04

AI Technical Summary

Technical Problem

In existing technologies, the monitoring of organic acids in museum air faces problems such as difficulty in user perception, high monitoring costs, and difficulties in hardware deployment. It cannot effectively reflect the cumulative effect of pollutant concentration changes over time, lacks intuitive display and dynamic early warning of the corrosion rate of cultural relics, and the professional equipment and software are expensive.

Method used

Data is collected using distributed sensors, and moving average filtering and data compression are performed through edge computing nodes. The corrosion rate is calculated using an improved multi-parameter lookup table method, and the cumulative corrosion exposure is calculated using the trapezoidal numerical integration method. A cultural relic health index is generated and mapped to a corrosion risk level. A heat map is generated through a 3D simulation engine and corresponding protection measures are triggered.

Benefits of technology

It achieves low-cost, lightweight 3D diffusion simulation of organic acids, provides an intuitive display and dynamic early warning of the corrosion risk of cultural relics, reduces hardware resource consumption, is suitable for small and medium-sized museums, and improves users' ability to perceive the health status of cultural relics.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN120995928B_ABST
    Figure CN120995928B_ABST
Patent Text Reader

Abstract

The present application relates to the technical field of air quality monitoring, and aims at the problems of difficult air organic acid gas perception, high monitoring cost and difficult hardware deployment for museum users, and provides a kind of organic acid visual three-dimensional diffusion simulation method and system, current environmental parameters are collected and filtered, improved multi-parameter lookup table method and trapezoidal numerical integration method are innovatively used to calculate the cumulative corrosion exposure of the week, and the cultural relic health index is calculated;Finally, the cultural relic health index and the cumulative corrosion exposure of the week are jointly mapped into the corrosion risk level, the organic acid corrosion risk level heat map is generated through the three-dimensional simulation engine, and the adjustment mechanism is triggered.Compared with the prior art, the present application reduces the cost of organic acid gas monitoring, while improving the visualization degree of organic acid monitoring and cultural relic health degree evaluation, and is suitable for museum whole exhibition hall, bronze exhibit cabinet, painting warehouse and other scenes.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This invention relates to the field of air quality monitoring technology, and in particular to a visualization three-dimensional diffusion simulation method and system for organic acids. Background Technology

[0002] Museum indoor air pollutants are diverse, including organic acids, inorganic acids, volatile organic compounds, ozone, and ammonia. Among these pollutants, gaseous organic acids are particularly concerning, including formic acid, acetic acid, and oxalic acid, which mainly originate from old decoration materials, visitors' exhaled breath, sebum evaporation, and chemical cleaning agents. When organic acids come into contact with the surface of cultural relics, they can cause serious corrosion problems: copper, silver, and other metallic artifacts will form carboxylates, leading to rust; paper and textiles will become brittle and discolored due to cellulose hydrolysis; and carbonate components in stone and murals will dissolve and peel off. This long-term cumulative corrosion effect often causes irreversible damage to cultural relics.

[0003] Currently, there are three main types of technologies for monitoring organic acids: First, computational fluid dynamics simulation technology, or CFD simulation technology, which simulates the three-dimensional diffusion of organic acids through computational fluid dynamics to generate concentration distribution cloud maps. Common commercial tools include ANSYS Fluent fluid simulation software and COMSOL multiphysics simulation. Second, some museums have adopted equipment such as gas chromatography-mass spectrometry and proton transfer reaction mass spectrometry to monitor organic acid concentrations and generate risk heat maps. Third, in recent years, foreign scholars have been trying to integrate CFD simulation with artificial intelligence for simulation monitoring, correcting simulation errors through algorithms such as LSTM and convolutional neural networks, or accelerating calculations using reduced-order model ROM.

[0004] However, existing technologies still have two major problems: First, they are difficult for users to perceive. Current heat maps are based on single or periodic sampling data, providing only transient or steady-state concentration distributions. They cannot reflect the cumulative effect of pollutant concentration changes over time, and lack intuitive display and dynamic early warning of the corrosion rate of cultural relics, making it difficult for users to perceive. Moreover, there are few mature solutions on the market that can achieve dynamic monitoring of cultural relic corrosion. Second, they are costly and difficult to deploy hardware. Currently, professional CFD simulation software such as ANSYS costs more than $100,000 per year. The purchase and maintenance costs of equipment such as gas chromatography-mass spectrometry are expensive, and frequent calibration by professionals is required. CFD transient simulation and AI training require high-performance computing clusters and servers equipped with professional graphics cards, and require the construction of a large database, which poses a high barrier to entry for small and medium-sized museums.

[0005] Therefore, developing a low-cost, lightweight, and visualized three-dimensional diffusion simulation method for organic acids has become an urgent need in the field of cultural relic protection. Summary of the Invention

[0006] The purpose of this invention is to provide a visual three-dimensional diffusion simulation method and system for organic acids, so as to solve the technical problems of difficulty in users to perceive organic acid gases in the air, high monitoring costs, and difficulties in hardware deployment in the prior art.

[0007] To achieve the above objectives, the present invention proposes the following technical solution:

[0008] Firstly, this technical solution provides a visualized three-dimensional diffusion simulation method for organic acids, used to monitor and simulate organic acid gases in indoor air of exhibition halls, including the following steps:

[0009] S1: Deploy distributed sensors indoors in the exhibition hall to collect data on the concentration of organic acids, temperature, and humidity in the current environment, and transmit the collected sensor data to edge computing nodes;

[0010] S2: The sensor data is filtered by moving average through edge computing nodes, and then uploaded to the server after data compression and encapsulation.

[0011] S3: The organic acid concentration, temperature and humidity at the current moment are mapped to the instantaneous corrosion rate of the cultural relic by an improved multi-parameter lookup table method, and the instantaneous corrosion rate is cached into an instantaneous corrosion rate sequence over time.

[0012] S4: Read the instantaneous corrosion rate sequence, take the current moment as the end point and back 7 days as the starting point, and use the trapezoidal numerical integration method to calculate the cumulative corrosion exposure for the week;

[0013] S5: Calculate the health index of cultural relics by calling the instantaneous corrosion rate and the cumulative corrosion exposure of the week at the same time, and combining the temperature, humidity and sensitivity coefficient of the material of the cultural relics;

[0014] S6: Map the health index of cultural relics together with the cumulative corrosion exposure of the week to a three-level corrosion risk level of red, yellow and green, and generate a heat map of organic acid corrosion risk level through a three-dimensional simulation engine; trigger sensor audible and visual alarms, SMS reminders or automatically start the display case microenvironment adjustment system according to the corrosion risk level.

[0015] Specifically, in step S1, the distributed sensor includes an organic acid sensor and an environmental parameter sensor. The organic acid sensor is used to collect the concentrations of formic acid and acetic acid in the air and filters out large molecular interferences through a filtration device. The environmental parameters include temperature, humidity, and airflow velocity. This step collects organic acid and environmental data through sensors, providing a comprehensive data source for subsequent corrosion rate calculation and risk assessment.

[0016] Furthermore, the filtration device is a filter membrane and activated carbon device installed at or inside the air inlet of the organic acid sensor; the macromolecular interfering substances include formaldehyde, ethanol, and toluene. Since macromolecular interfering substances such as formaldehyde and ethanol can cross-interfere with organic acids in the air, affecting the sensitivity and selectivity of the sensor, the filtration device can reduce their impact on the organic acid measurement results and enhance the reliability of the system.

[0017] Specifically, in S2, the moving average filtering refers to weighted averaging of data within a sliding time window to suppress random noise, making the data smoother and improving its quality and stability; the sliding time window is adjustable, ranging from 1 to 30 seconds. This increases the flexibility and adaptability of the invention, allowing it to be adjusted according to different monitoring needs.

[0018] Furthermore, in step S2, the data compression and encapsulation steps are as follows:

[0019] S21: Fix each data frame to 14 bytes;

[0020] S22: Perform secondary compression on each data frame using an offline-trained static Huffman code table;

[0021] S23: Compressed data is uploaded in batches via MQTT-SN over UDP packets, with each edge node triggering a transmission every 30 seconds or when the buffer is full (1kB).

[0022] S24: Data compression and encapsulation are complete; upload the data to the server.

[0023] Data compression and encapsulation can reduce the amount of data transmitted, improve data transmission efficiency, and also facilitate subsequent data processing by the server.

[0024] Specifically, in step S3, the step of calculating the instantaneous corrosion rate using the improved multi-parameter lookup table method is as follows:

[0025] S31: Receive sensor data, including organic acid concentration, temperature, and humidity; receive the material type of cultural relics, including bronze, paper, stone, and textiles, and find the corresponding material corrosion rate pre-stored reference table.

[0026] S32: Check whether the current environmental organic acid concentration, temperature, and humidity are all within the range of the pre-stored reference table.

[0027] S321: If within the range of the reference table, then proceed to the dynamic interpolation lookup process, using trilinear interpolation. The formula for calculating the instantaneous corrosion rate R is as follows:

[0028] Weight normalization: ,

[0029] Where R represents the instantaneous corrosion rate of the artifact, in μm / y; Ri represents the corrosion rate at the i-th vertex in the baseline table, in μm / y; and the weight w i The position of the current environmental parameters in the pre-stored baseline table is determined by:

[0030]

[0031] Among them, w i represents the weights, unitless; di represents the standardized distance from the current environment parameters to each vertex, unitless; e represents the smoothing factor, unitless; HA i ,T i ,RH i ΔHA, ΔT, and ΔRH represent the organic acid concentration, temperature, and humidity values ​​at the i-th vertex in the baseline table, respectively, in ppm, °C, and %; ΔHA, ΔT, and ΔRH represent the step sizes of the organic acid concentration, temperature, and humidity in the baseline table, respectively, in ppm, °C, and %; The advantage of this approach is that it allows for the calculation of corrosion rates based on different parameter risk levels and material characteristics, thus improving the accuracy and adaptability of corrosion rate calculations.

[0032] S322: If the value exceeds the baseline range, model compensation extrapolation is triggered, and a simplified physical model is invoked to calculate the corrosion rate. The formula for calculating the instantaneous corrosion rate R using the simplified physical model is as follows:

[0033]

[0034] Where R represents the instantaneous corrosion rate of the cultural relic, in μm / y; Rmax represents the corrosion rate corresponding to the maximum boundary of the parameter in the benchmark table, in μm / y; HAmax and RHmax represent the upper limits of the organic acid concentration and humidity parameters in the benchmark table, in ppm and % respectively; n represents the material-related organic acid concentration sensitivity index; k represents the humidity influence coefficient; and α represents the material compensation factor, all in dimensionless units.

[0035] S33: Output corrosion rate value;

[0036] The improved multi-parameter lookup table method is designed with reference to the standards for controlling organic acids and humidity in cultural relics collections, as specified in "JGJ 66-2015 Museum Building Design Code", "WH / T 24—2006 Basic Requirements for Ancient Books and Special Collections in Libraries", and "DB4403_T 565—2024 Technical Specification for Indoor Ambient Air Pollutant Control in Exhibitions". Compared with the complex models used in existing CFD simulation tools, the dynamic interpolation lookup table presets the corrosion rates corresponding to different materials, resulting in higher computational efficiency. Model compensation extrapolation calculates the corrosion rate by simplifying the physical model, ensuring the integrity and continuity of the calculation and avoiding calculation failures due to data exceeding the range. The combination of these two methods can reduce hardware resource consumption, deployment costs, and complexity while ensuring the efficiency and reliability of corrosion rate calculation.

[0037] Specifically, in S4, the trapezoidal numerical integration method refers to sampling organic acid concentration data at sampling intervals, starting from the current time and going back 7 days, and calculating the cumulative corrosion exposure for the week. The formula for calculating the cumulative corrosion exposure D(t) for the week is as follows:

[0038]

[0039] Where i is an integer representing the i-th sampling point in the time series, Ci represents the organic acid concentration of the i-th sampling, Δt represents the sampling time interval of 5 minutes, and when the concentration change per minute is greater than 0.1 ppm, the sampling interval is shortened to 1 minute.

[0040] Furthermore, the materials of the S3 cultural relics include bronzes, paper artifacts, stone artifacts, and textile artifacts. According to GB / T 30688-2014 "Diseases and Illustrations of Bronze Ware in Museum Collections," WW / T 0067-2015 "Preventive Protection Standard for Paper Cultural Relics," WW / T 0002-2007 "Classification and Illustrations of Diseases of Stone Cultural Relics," and GB / T 30236-2013 "Technical Specifications for the Protection of Textile Cultural Relics," the annual threshold for the corrosion rate of cultural relics is: 10 ppm·h / year for bronzes, 5 ppm·h / year for paper, 20 ppm·h / year for stone, and 3 ppm·h / year for textiles. This clarifies the types of cultural relics made of specific materials and their annual corrosion rate thresholds, providing specific reference standards for subsequent corrosion risk assessment.

[0041] Specifically, the formula for calculating the health index H of cultural relics in S5 is as follows:

[0042]

[0043] Among them, Rc is the normalized instantaneous corrosion rate, ranging from 0 to 100; Dc is the proportion of cumulative exposure; Ee is the environmental humidity fluctuation coefficient; Mm is the material sensitivity coefficient; α, β, and γ are the weight coefficients corresponding to the normalized instantaneous corrosion rate, the proportion of cumulative exposure, and the environmental humidity fluctuation coefficient, respectively. In this step, referring to the index requirements for organic acids, temperature, and humidity in the preservation environment of cultural relics in "WW / T 0016.1-2023 Quality of the Preservation Environment for Cultural Relics in Collections - Part 1: Index Requirements", and combining the requirements for the preservation environment of collections in Chapter 6 of "JGJ 66-2015 Guide for Environmental Quality Control in Museums" of the National Cultural Heritage Administration, the weight distribution principle of the health index H is designed. By comprehensively considering the instantaneous corrosion rate, the proportion of cumulative exposure, and the environmental humidity fluctuation coefficient, it provides a scientific and quantitative basis for the evaluation of the health of cultural relics, making the evaluation of the health status of cultural relics more objective and accurate, and facilitating users to take corresponding protection measures according to the health of cultural relics.

[0044] Specifically, in S6, the corrosion risk levels of the cultural relics are divided into three types: red, yellow, and green. Among them, red indicates a high-risk level, that is, when the cultural relic health index H ≤ 70 or the weekly cumulative exposure > 30% of the annual threshold; yellow indicates a general risk level, that is, when the cultural relic health index 70 < H ≤ 85 and the weekly cumulative exposure ≤ 30% of the annual threshold; green indicates no risk, that is, when the cultural relic health index H > 85 and the weekly cumulative exposure ≤ 5% of the annual threshold.

[0045] Specifically, the triggering measures for the corrosion risk levels specifically include:

[0046] High-risk level, trigger the audible and visual alarm of the sensor and immediately send a signal to start the microenvironment regulation system of the display cabinet;

[0047] General risk level, send a text message reminder to the administrator;

[0048] No risk, prompt normal.

[0049] The corrosion risk levels of the cultural relics visually distinguish the corrosion risk levels of the cultural relics through colors, enabling users to quickly and intuitively understand the corrosion risk status of the cultural relics, facilitating timely adoption of corresponding protection measures, and making up for the technical deficiencies of the traditional difficulty in monitoring and perceiving organic acids and cultural relic corrosion in the exhibition hall by users.

[0050] On the second aspect, this technical solution also provides a visual three-dimensional diffusion simulation system for organic acids, including distributed sensors, edge computing nodes, and servers:

[0051] The distributed sensors are deployed indoors in the exhibition hall to collect data on the concentration of organic acids, temperature, and humidity in the current environment. The distributed sensors include organic acid sensors and environmental parameter sensors. The organic acid sensors have a range of 0-5 ppm, an accuracy of ≤ ±0.03 ppm, a cross-sensitivity to ethanol / formaldehyde of ≤3%, and a response time of ≤30 seconds. This is to ensure that the sensors can meet the high-precision and high-sensitivity measurement requirements of the aforementioned method for organic acids and environmental parameters. The environmental parameter sensors include temperature and humidity sensors.

[0052] The edge computing node is used to perform moving average filtering on sensor data, compress and encapsulate the data, and upload the sensor data to the server.

[0053] The server includes a 3D simulation engine and a user interaction module, wherein the user interaction module includes a cumulative exposure query module, a cultural relic corrosion early warning module, and a cultural relic health query module.

[0054] The cumulative exposure query module is used to query the cumulative corrosion exposure of cultural relics in the current week;

[0055] The cultural relic health query module is used to query the cultural relic health index.

[0056] The cultural relic corrosion early warning module is used to trigger sensor audible and visual alarms, send SMS reminders to administrators, or automatically start the display case microenvironment adjustment system based on the risk level of cultural relic corrosion.

[0057] Beneficial effects

[0058] As can be seen from the above technical solutions, the technical solution of the present invention provides a method and system for visualizing the three-dimensional diffusion simulation of organic acid gases in indoor air of exhibition halls, so as to solve the technical defects of existing technologies such as difficulty in user perception, high monitoring cost and difficulty in hardware deployment.

[0059] 1. In terms of user perception, this invention employs an improved multi-parameter lookup table method to calculate the corrosion rate, utilizes the trapezoidal numerical integration method to calculate the cumulative exposure of organic acids, and employs a cultural relic health index to quantitatively perceive the health status of cultural relics. Addressing the past difficulty in perceiving organic acids, this invention provides users with a solution for intuitively displaying and dynamically warning about the health status of cultural relics, thus overcoming the technical shortcomings of existing technologies that hinder user perception.

[0060] 2. Regarding monitoring costs, this invention innovatively employs an improved multi-parameter lookup table method, which increases the calculation speed of corrosion rates by pre-setting the corrosion rates corresponding to different materials. Furthermore, it utilizes a 3D simulation engine to display heatmaps of organic acid corrosion risk levels, further reducing software costs. Compared to the complex models, high annual fees, and equipment costs of current mainstream CFD simulation software and gas chromatography-mass spectrometry (GC-MS) systems, this invention consumes fewer hardware resources, has lower costs, and is more suitable for small and medium-sized museums.

[0061] 3. In terms of hardware deployment, traditional methods of CFD transient simulation and AI training require high-performance computing clusters and servers equipped with professional graphics cards, as well as the construction of a large database, which poses a high barrier to entry for small and medium-sized museums. This invention improves upon the traditional corrosion rate calculation method by developing an improved multi-parameter lookup table method. Through dynamic interpolation lookup and model compensation extrapolation, it reduces the cost of instantaneous calculations while ensuring the integrity and continuity of the calculations. This calculation method can reduce hardware resource consumption, deployment costs, and complexity while maintaining the efficiency and reliability of corrosion rate calculations.

[0062] It should be understood that all combinations of the foregoing concepts and the additional concepts described in more detail below can be considered part of the inventive subject matter of this disclosure, provided that such concepts do not contradict each other.

[0063] The foregoing and other aspects, embodiments, and features of the teachings of the present invention will be more fully understood from the following description in conjunction with the accompanying drawings. Other additional aspects of the invention, such as features and / or beneficial effects of exemplary embodiments, will become apparent from the following description or may be learned through practice of specific embodiments according to the teachings of the present invention. Attached Figure Description

[0064] The accompanying drawings are not drawn to scale. In the drawings, each identical or nearly identical component shown in the various figures may be denoted by the same reference numeral. For clarity, not every component is labeled in each figure. Embodiments of various aspects of the invention will now be described by way of example and with reference to the accompanying drawings, wherein:

[0065] Figure 1 This is a flowchart of a visualization three-dimensional diffusion simulation method for organic acids as described in this embodiment;

[0066] Figure 2 This is a flowchart illustrating the calculation of corrosion rate based on the improved multi-parameter lookup table method in this embodiment;

[0067] Figure 3 This is a structural block diagram of a visualization three-dimensional diffusion simulation system for organic acids described in this embodiment. Detailed Implementation

[0068] To make the objectives, technical solutions, and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some, not all, of the embodiments of the present invention. All other embodiments obtained by those skilled in the art based on the described embodiments of the present invention without creative effort are within the scope of protection of the present invention. Unless otherwise defined, the technical or scientific terms used herein should have the ordinary meaning understood by those skilled in the art.

[0069] The terms "first," "second," and similar words used in the specification and claims of this patent application do not indicate any order, quantity, or importance, but are merely used to distinguish different components. Similarly, unless the context clearly indicates otherwise, the singular forms of "an," "a," or "the," etc., do not indicate a quantity limitation, but rather indicate the presence of at least one. Terms such as "comprising" or "including" mean that the element or object preceding "comprising" encompasses the features, wholes, steps, operations, elements, and / or components listed following "comprising" or "including," and do not exclude the presence or addition of one or more other features, wholes, steps, operations, elements, components, and / or collections thereof. Terms such as "upper," "lower," "left," and "right" are used only to indicate relative positional relationships; when the absolute position of the described object changes, the relative positional relationship may also change accordingly.

[0070] Monitoring organic acid gases such as formic acid and acetic acid in museum interiors is a crucial aspect of cultural relic preservation. Organic acid gases accelerate metal corrosion, paper acidification, and textile brittleness and discoloration; long-term cumulative corrosion causes irreversible damage to cultural relics. Existing technologies rely on CFD simulation software and gas chromatography-mass spectrometry (GC-MS), which present technical challenges such as difficulty in user perception of organic acids, high monitoring costs, and difficulties in hardware deployment. Therefore, this embodiment aims to provide a visualized three-dimensional diffusion simulation method and system for organic acids to simultaneously address the aforementioned technical problems.

[0071] The following description, in conjunction with the accompanying drawings, details a method for visualizing three-dimensional diffusion simulation of organic acids as described in this embodiment.

[0072] Combination Figure 1 As shown, the method includes:

[0073] S1: Deploy distributed sensors indoors in the exhibition hall to collect data on the concentration of organic acids, temperature, and humidity in the current environment, and transmit the collected sensor data to edge computing nodes;

[0074] S2: The sensor data is filtered by moving average through edge computing nodes, and then uploaded to the server after data compression and encapsulation.

[0075] S3: The organic acid concentration, temperature and humidity at the current moment are mapped to the instantaneous corrosion rate of the cultural relic by an improved multi-parameter lookup table method, and the instantaneous corrosion rate is cached into an instantaneous corrosion rate sequence over time.

[0076] S4: Read the instantaneous corrosion rate sequence, take the current moment as the end point and back 7 days as the starting point, and use the trapezoidal numerical integration method to calculate the cumulative corrosion exposure for the week;

[0077] S5: Calculate the health index of cultural relics by calling the instantaneous corrosion rate and the cumulative corrosion exposure of the week at the same time, and combining the temperature, humidity and sensitivity coefficient of the material of the cultural relics;

[0078] S6: Map the health index of cultural relics together with the cumulative corrosion exposure of the week to a three-level corrosion risk level of red, yellow and green, and generate a heat map of organic acid corrosion risk level through a three-dimensional simulation engine; trigger sensor audible and visual alarms, SMS reminders or automatically start the display case microenvironment adjustment system according to the corrosion risk level.

[0079] Specifically, this method divides the indoor area of ​​the exhibition hall into key areas and auxiliary areas, and deploys different numbers of distributed sensors in each area. The key areas are those that need to be monitored, while the auxiliary areas are those that need to be monitored in general.

[0080] Specifically, this embodiment uses a museum exhibition hall as an example to illustrate the options and deployment of distributed sensors:

[0081] In this embodiment, the key areas are those where the stored cultural relics are made of sensitive materials, have a high national level, or are prone to pollutants due to high visitor traffic, including key cultural relic display areas such as bronze ware display cases, ancient book and document exhibition halls, and silk fabric display areas; one node is deployed every 5 to 10 square meters in the key areas, and the data sampling frequency is once per minute.

[0082] The auxiliary area is a storage area for cultural relics made of corrosion-resistant materials, replicas of cultural relics, or temporary exhibition areas, as well as well-ventilated edge areas, including stone carving corridors, visitor rest areas, equipment rooms, etc. One node is deployed every 20-30 square meters in the auxiliary area, and the data sampling frequency is once every 5 minutes.

[0083] Furthermore, the sensor for collecting organic acids has a range of 0-5 ppm, an accuracy of ≤ ±0.03 ppm, a cross-sensitivity of ethanol / formaldehyde of ≤3%, and a response time of ≤30 seconds. This selection is to ensure that the sensor can meet the high-precision and high-sensitivity measurement requirements of this invention for organic acids and environmental parameters. In this embodiment, the GASS-3000 series electrochemical sensor is selected as the organic acid sensor for use in the key area to meet the above-mentioned high-precision acquisition requirements. The CHEMS-100 series electrochemical sensor is selected as the organic acid sensor for use in the auxiliary area, which can meet the above parameters. At the same time, compared with the GASS-3000 series, it has a lower cost to achieve economical coverage. It should be noted that the measured values ​​and sensitivity of the organic acid electrochemical sensor will drift slightly due to factors such as electrolyte evaporation, temperature, and humidity. Therefore, the above-mentioned electrochemical sensors also need to be calibrated every 3-6 months to ensure measurement accuracy and sensitivity.

[0084] Furthermore, the environmental parameter sensors include temperature sensors and humidity sensors; in this embodiment, the ENV-TH1 series temperature and humidity sensors are selected, with one set deployed in key areas for every 10㎡ and one set deployed in auxiliary areas for every 30㎡, in order to achieve full coverage.

[0085] Furthermore, the organic acid sensor is used to collect the content of formic acid and acetic acid in the air, and filters out large molecular interferences through a filtration device; the environmental parameter sensor is used to collect temperature, humidity and airflow speed; by collecting organic acid and environmental data through the sensors, a comprehensive data source is provided for subsequent corrosion rate calculation and risk assessment.

[0086] Furthermore, in this embodiment, the filtration device is a filter membrane and activated carbon device installed at or inside the air inlet of the organic acid sensor; the macromolecular interfering substances include formaldehyde, ethanol, and toluene. Since macromolecular interfering substances such as formaldehyde and ethanol can cross-interfere with organic acids in the air, affecting the sensitivity and selectivity of the sensor, the filtration device can reduce their impact on the organic acid measurement results and enhance the reliability of the system.

[0087] Specifically, in S2, the moving average filtering refers to weighted averaging of data within a sliding time window to suppress random noise, making the data smoother and improving its quality and stability; the sliding time window is adjustable, ranging from 1 to 30 seconds. This increases the flexibility and adaptability of the invention, allowing it to be adjusted according to different monitoring needs.

[0088] Furthermore, this embodiment employs edge computing with a single node computing power of ≥1 GOPS to achieve real-time filtering and compression of sensor data; power consumption is ≤2W, as some display cases lack external power supplies and rely on battery power; it features at least 2 UART channels + 1 SPI channel to support connections to organic acid and temperature / humidity sensors; and due to the museum's constant temperature but potential dust accumulation inside equipment cabinets, its operating temperature is between -10℃ and 50℃, ensuring dust protection. This embodiment utilizes Chipstar CMS32 series and Espressif ESP32 series computing nodes.

[0089] Furthermore, in step S2, the data compression and encapsulation steps are as follows:

[0090] S21: Fix each data frame to 14 bytes;

[0091] S22: Perform secondary compression on each data frame using an offline-trained static Huffman code table;

[0092] S23: Compressed data is uploaded in batches via MQTT-SN over UDP packets, with each edge node triggering a transmission every 30 seconds or when the buffer is full (1kB).

[0093] S24: Data compression and encapsulation are complete; upload the data to the server.

[0094] Data compression and encapsulation can reduce the amount of data transmitted, improve data transmission efficiency, and also facilitate subsequent data processing by the server.

[0095] Combination Figure 2 As shown, the improved multi-parameter lookup table method is used to map the current organic acid concentration, temperature, and humidity to the instantaneous corrosion rate of the cultural relic. The steps are as follows:

[0096] S31: Receive sensor data, including organic acid concentration, temperature, and humidity; receive the material type of cultural relics, including bronze, paper, stone, and textiles, and find the corresponding material corrosion rate pre-stored reference table.

[0097] S32: Check whether the current environmental organic acid concentration, temperature, and humidity are all within the range of the pre-stored reference table.

[0098] S321: If within the range of the reference table, then proceed to the dynamic interpolation lookup process, using trilinear interpolation. The formula for calculating the instantaneous corrosion rate R is as follows:

[0099] Weight normalization: ,

[0100] Where R represents the instantaneous corrosion rate of the artifact, in μm / y; Ri represents the corrosion rate at the i-th vertex in the baseline table, in μm / y; and the weight w iThe position of the current environmental parameters in the pre-stored baseline table is determined by:

[0101]

[0102] Among them, w i The weights are unitless; di represents the standardized distance from the current environment parameters to each vertex, unitless; e represents the smoothing factor, unitless; HA i ,T i ,RH i ΔHA, ΔT, and ΔRH represent the organic acid concentration, temperature, and humidity values ​​at the i-th vertex in the baseline table, respectively, in ppm, °C, and %; ΔHA, ΔT, and ΔRH represent the step sizes of the organic acid concentration, temperature, and humidity in the baseline table, respectively, in ppm, °C, and %; The advantage of this approach is that it allows for the calculation of corrosion rates based on different parameter risk levels and material characteristics, thus improving the accuracy and adaptability of corrosion rate calculations.

[0103] Furthermore, taking the calculation of the corrosion rate of bronze artifacts as an example, the environmental parameters are as follows:

[0104] The organic acid concentration HA is 0.07 ppm (no direct data is available in the table, but it is between 0.05 and 0.1 ppm).

[0105] Temperature T = 22°C (between 20 and 25°C)

[0106] Humidity RH = 60% (between 50% and 70%)

[0107] Baseline data (corrosion rate of bronze artifacts, in μm / y).

[0108] Referring to GB / T 16545-2015 Method for Determination of Corrosion Rate of Metals, GB / T 30688-2014 Diseases and Illustrations of Bronze Ware in Museum Collections, and WW / T 0066-2015 Specification for Monitoring of Preservation Environment of Cultural Relics in Museum Collections, the following table provides a pre-stored benchmark table for the corrosion rate of bronze ware in relation to organic acid concentration, temperature, and humidity:

[0109] 1 0.05 20 50 0.08 2 0.05 20 70 0.12 3 0.05 25 50 0.10 4 0.05 20 70 0.15 5 0.10 20 50 0.15 6 0.10 20 70 0.20 7 0.10 25 50 0.18 8 0.10 25 70 0.25

[0110] First, calculate the standardized distances from the environmental parameters to each vertex, i.e., step size ΔHA = 0.05ppm, ΔT = 5°C, ΔRH = 20%. Substituting these values ​​into the formula, we get:

[0111] ,

[0112] Similarly, calculate the other values ​​of di, and we get d2=0.70, d3=0.85, ..., d8=0.65;

[0113] Then take the smoothing factor e=e -6 The weights are calculated as follows:

[0114]

[0115] Calculate the weighted average corrosion rate R:

[0116]

[0117] S33: Output the instantaneous corrosion rate of bronze artifacts under this environment as 0.14 μm / y.

[0118] S322: If the value exceeds the baseline range, model compensation extrapolation is triggered, and a simplified physical model is invoked to calculate the corrosion rate. The formula for calculating the instantaneous corrosion rate R using the simplified physical model is as follows:

[0119]

[0120] Where R represents the instantaneous corrosion rate of the cultural relic, in μm / y; Rmax represents the corrosion rate corresponding to the maximum boundary of the parameter in the benchmark table, in μm / y; HAmax and RHmax represent the upper limits of the organic acid concentration and humidity parameters in the benchmark table, respectively, in ppm and %;

[0121] In addition, n represents the material-related organic acid concentration sensitivity index, k represents the humidity influence coefficient, and α represents the material compensation factor; all units are dimensionless. Referring to GB / T 16545-2015 Method for Determination of Corrosion Rate of Metals, GB / T 30688-2014 Diseases and Illustrations of Bronze Ware in Museum Collections, and WW / T 0066-2015 Specification for Monitoring of Preservation Environment of Cultural Relics in Museum Collections, the specific parameters for the corrosion rate of cultural relics in relation to organic acid concentration, temperature, and humidity are as follows:

[0122] Furthermore, taking the calculation of the corrosion rate of bronze artifacts as an example, the environmental parameters are as follows:

[0123] Organic acid concentration HA = 6 ppm

[0124] Temperature T=25°C

[0125] Humidity RH=65%

[0126] According to the table above, the sensitivity index of organic acid concentration related to the material of bronze is n=0.9, the humidity influence coefficient is k=0.05, and the material compensation factor is α=1.0.

[0127] Because the organic acid concentration exceeds the maximum value of 5 ppm in the baseline table, model compensation extrapolation is triggered:

[0128] First, calculate the acid concentration correction term:

[0129]

[0130] Then calculate the humidity correction term:

[0131]

[0132] Overall extrapolated corrosion rate:

[0133]

[0134] S33: Output the instantaneous corrosion rate of bronze in this environment as 1.83 μm / y.

[0135] The improved multi-parameter lookup table method is designed with reference to the standards for controlling organic acids and humidity in cultural relics collections, as specified in "JGJ 66-2015 Museum Building Design Code", "WH / T 24—2006 Basic Requirements for Ancient Books and Special Collections in Libraries", and "DB4403_T 565—2024 Technical Specification for Indoor Ambient Air Pollutant Control in Exhibitions". Compared with the complex models used in existing CFD simulation tools, the dynamic interpolation lookup table presets the corrosion rates corresponding to different materials, resulting in higher computational efficiency. Model compensation extrapolation calculates the corrosion rate by simplifying the physical model, ensuring the integrity and continuity of the calculation and avoiding calculation failures due to data exceeding the range. The combination of these two methods can reduce hardware resource consumption, deployment costs, and complexity while ensuring the efficiency and reliability of corrosion rate calculation.

[0136] Furthermore, the materials of cultural relics in S3 include bronzes, paper artifacts, stone artifacts, and textile artifacts. According to GB / T 30688-2014 "Diseases and Illustrations of Bronze Artifacts in Museum Collections," WW / T 0067-2015 "Preventive Protection Standard for Paper Artifacts," WW / T 0002-2007 "Classification and Illustrations of Diseases of Stone Artifacts," and GB / T 30236-2013 "Technical Specifications for the Protection of Textile Artifacts," the annual threshold for the corrosion rate of cultural relics is: 10 ppm·h / year for bronzes, 5 ppm·h / year for paper, 20 ppm·h / year for stone, and 3 ppm·h / year for textiles. This clarifies the types of cultural relics made of specific materials and their annual corrosion rate thresholds, providing specific reference standards for subsequent corrosion risk assessments.

[0137] Specifically, in S4, the trapezoidal numerical integration method refers to sampling organic acid concentration data at sampling intervals, starting from the current time and going back 7 days, and calculating the cumulative corrosion exposure D(t) for the week, in ppm·h. The formula for calculating the cumulative corrosion exposure D(t) for the week is as follows:

[0138]

[0139] Where i is an integer, representing the i-th sampling point in the time series, c i The concentration of organic acid in the i-th sample is expressed in ppm; Δt represents the sampling time interval of 5 minutes; when the concentration change per minute is greater than 0.1 ppm, the sampling interval is shortened to 1 minute.

[0140] The following example uses paper as an example to illustrate the calculation process of the cumulative exposure of cultural relics:

[0141] The current time is the 7th day of closure, and it is necessary to calculate the cumulative corrosion exposure of paper artifacts over the past 7 days (168 hours). The default sampling interval Δt = 5 minutes; the annual threshold for the corrosion rate of paper artifacts is 5 ppm·h / year.

[0142] Taking the first 20 minutes after the exhibition hall opens on the first day as an example, the concentration of organic acids in the air fluctuates once during the first 20 minutes after the opening, as visitors flock in. Specifically, the organic acid concentrations at 9:00, 9:05, 9:10, 9:15, and 9:20 are 0.08 ppm, 0.12 ppm, 0.18 ppm, 0.22 ppm, and 0.15 ppm, respectively. The cumulative exposure is calculated segment by segment using the formula described above.

[0143] 9:00-9:05 time period:

[0144] And so on, with each 5-minute interval constituting a time period, the corresponding exposure amount is:

[0145] The total accumulated amount in the current 20 minutes is: 0.5 + 0.75 + 1 + 0.923 = 3.175 ppm·min, which is converted to ppm·h = 3.175 / 60 ​​= 0.053 pm·h.

[0146] The daily records show 10 instances of the aforementioned 20-minute fluctuation per day, with a total daily fluctuation duration of 20 minutes / event × 10 times = 200 minutes = 3.33 hours; the concentration remained stable at 0.05 ppm during other periods, with a stable duration of 24 hours - 3.33 hours = 20.67 hours; therefore, the daily exposure is: 0.053 × 10 + 0.05 × 22.67 = 1.56 ppm·h

[0147] If the daily exposure is the same, then the cumulative corrosion exposure D(t) for the week is 1.56*7=10.92 ppm·h;

[0148] Compared with the annual threshold of 5 ppm·h / year for paper cultural relics, the cumulative exposure is 10.92 / 5*100%=218.4% of the annual threshold, which is seriously exceeded, providing data support for subsequent safety risk warning.

[0149] Specifically, in step S5, the formula for calculating the health index H of the cultural relics is as follows:

[0150]

[0151] Where H is the Cultural Relics Health Index, which is defined according to the standards for controlling organic acids and humidity in cultural relics collections in "JGJ 66-2015 Museum Building Design Code", "WH / T 24—2006 Basic Requirements for Ancient Books and Special Collections in Libraries", and "DB4403_T 565—2024 Technical Specification for Indoor Ambient Air Pollutant Control in Exhibitions". It is used to express the health of cultural relics, and the unit is dimensionless.

[0152] Rc is the normalized instantaneous corrosion rate = instantaneous corrosion rate R / maximum corrosion rate Rmax * 100; Rc ranges from 0 to 100, where 0 represents no corrosion and 100 represents corrosion reaching its maximum value.

[0153] Dc is the percentage of cumulative exposure = cumulative corrosion exposure D(t) in the current week / annual threshold * 100;

[0154] Ee is the environmental humidity fluctuation coefficient = standard deviation of environmental humidity in the past hour / average environmental humidity in the past hour * 100;

[0155] Mmt is a material sensitivity coefficient preset according to the type of cultural relic: 1.0 for bronze, 1.2 for paper, 0.8 for stone, and 1.1 for textiles.

[0156] α, β, and γ are the weighting coefficients corresponding to the normalized instantaneous corrosion rate, the proportion of cumulative exposure, and the environmental fluctuation coefficient, respectively. The coefficients can be adjusted, and the default values ​​are 0.4, 0.3, and 0.2.

[0157] The following example uses bronze to illustrate the calculation process for the health status of cultural relics:

[0158] In this embodiment, the instantaneous corrosion rate is 0.14 μm / y, the maximum corrosion rate is 10 μm / y, and the normalized corrosion rate is 0.14 / 10*100=1.4; the weekly cumulative corrosion exposure D(t) is 1.8 ppm·h, the annual threshold is 10 ppm·h / year, and the cumulative exposure percentage is weekly cumulative corrosion exposure D(t) / annual threshold*100=18; the standard deviation of the ambient humidity in the past hour is 15%, the average ambient humidity in the past hour is 60%, the environmental fluctuation coefficient is 15% / 60%*100=25, and the material sensitivity coefficient is 1.0. These values ​​are then substituted into the health formula for calculation.

[0159]

[0160] Since the system defines that a health level greater than 85 is risk-free and the improvement is normal. If the health level is lower than 85, it is a general risk, and a text message is sent to the administrator.

[0161] This step also refers to the index requirements for organic acids, temperature, and humidity in the preservation environment of cultural relics in "WW / T 0016.1-2023 Quality of the Preservation Environment of Collected Cultural Relics - Part 1: Index Requirements", and combines the requirements for the preservation environment of collections in Chapter 6 of "JGJ 66-2015 Guide for Environmental Quality Control in Museums" of the National Cultural Heritage Administration to design the weight distribution principle of the health index H. By comprehensively considering the instantaneous corrosion rate, the proportion of cumulative exposure, and the environmental humidity fluctuation coefficient, it provides a scientific and quantitative basis for the evaluation of the health of cultural relics, making the evaluation of the health status of cultural relics more objective and accurate, and facilitating users to take corresponding protection measures according to the health of cultural relics.

[0162] Specifically, in S6, the corrosion risk levels of cultural relics are divided into three types: red, yellow, and green. Among them, red indicates a high-risk, that is, when the health index H of cultural relics ≤ 70 or the weekly cumulative exposure > 30% of the annual threshold; yellow indicates a general risk, that is, when the health index of cultural relics 70 < H ≤ 85 and the weekly cumulative exposure ≤ 30% of the annual threshold; green indicates risk-free, that is, when the health index H of cultural relics > 85 and the weekly cumulative exposure ≤ 5% of the annual threshold; the risk level triggering measures include:

[0163] High-risk, trigger the audible and visual alarm of the sensor and immediately send a signal to start the display case microenvironment regulation system;

[0164] General risk, send a text message reminder to the administrator;

[0165] Risk-free, prompt normal.

[0166] This step also refers to the index requirements for organic acids, temperature, and humidity in the preservation environment of cultural relics in "WW / T 0016.1-2023 Quality of the Preservation Environment of Collected Cultural Relics - Part 1: Index Requirements", and combines the requirements for the preservation environment of collections in Chapter 6 of "JGJ 66-2015 Guide for Environmental Quality Control in Museums" of the National Cultural Heritage Administration to design the threshold of the cultural relic health index. The corrosion risk levels of cultural relics visually distinguish the corrosion risk levels of cultural relics by colors, enabling users to quickly and intuitively understand the corrosion risk status of cultural relics, facilitating timely adoption of corresponding protection measures, and making up for the technical deficiencies of the traditional difficulty in monitoring and perceiving organic acids and cultural relic corrosion in exhibition halls.

[0167] Combined Figure 3 As shown, this embodiment also provides a visual three-dimensional diffusion simulation system for organic acids, including distributed sensors, edge computing nodes, and servers:

[0168] The distributed sensors are deployed indoors in the exhibition hall to collect data on the concentration of organic acids, temperature, and humidity in the current environment. The distributed sensors include organic acid sensors and environmental parameter sensors. The organic acid sensors have a range of 0-5 ppm, an accuracy of ≤ ±0.03 ppm, a cross-sensitivity to ethanol / formaldehyde of ≤3%, and a response time of ≤30 seconds. This ensures that the sensors can meet the high-precision and high-sensitivity measurement requirements of the aforementioned method for organic acids and environmental parameters. The environmental parameter sensors include temperature and humidity sensors. Specifically, this embodiment selects the GASS-3000 series or CHEMS-100 series electrochemical sensors as the organic acid sensors and the ENV-TH1 series temperature and humidity sensors, which can meet the above requirements.

[0169] The edge computing node is used to perform moving average filtering on the sensor data, compress and encapsulate the data, and upload the sensor data to the server. Specifically, in this embodiment, the Chipstar CMS32 series or Espressif ESP32 series computing nodes are used as edge computing nodes.

[0170] The server includes a 3D simulation engine and a user interaction module, wherein the user interaction module includes a cumulative exposure query module, a cultural relic corrosion early warning module, and a cultural relic health query module.

[0171] The cumulative exposure query module is used to query the cumulative corrosion exposure of cultural relics in the current week;

[0172] The cultural relic health query module is used to query the cultural relic health index.

[0173] The cultural relic corrosion early warning module is used to trigger sensor audible and visual alarms, send SMS reminders to administrators, or automatically start the display case microenvironment adjustment system based on the risk level of cultural relic corrosion.

[0174] Specifically, this embodiment uses one Inspur NF5180 series server for a minimal deployment, supporting data from 20 organic acid sensors and 50 temperature and humidity sensors. The server configuration is as follows: CPU: Phytium FT-2000 quad-core, 2.0GHz or equivalent; 16GB DDR4 memory; storage: 1 x 480GB SATA SSD + 1 x 2TB HDD; network port: one 1Gbps port; operating system: domestic Kylin V10 operating system. Servers with other configurations not lower than the above requirements are also considered. As the number of sensor nodes increases, the corresponding number of servers can be increased to meet the core functions of the visualized three-dimensional diffusion simulation system for organic acids, including cumulative exposure query, cultural relic corrosion early warning, and cultural relic health query.

[0175] Since the system is built based on the method, it also has the technical advantage of high testing accuracy in practical applications.

[0176] While the present invention has been disclosed above with reference to preferred embodiments, it is not intended to limit the invention. Those skilled in the art can make various modifications and refinements without departing from the spirit and scope of the invention. Therefore, the scope of protection of the present invention shall be determined by the claims.

Claims

1. A visualization-based three-dimensional diffusion simulation method for organic acids, used for monitoring and simulating organic acid gases in indoor air of exhibition halls, characterized in that... Includes the following steps: S1: Deploy distributed sensors indoors in the exhibition hall to collect data on the concentration of organic acids, temperature, and humidity in the current environment, and transmit the collected sensor data to edge computing nodes; S2: The sensor data is filtered by moving average through edge computing nodes, and then uploaded to the server after data compression and encapsulation. S3: The organic acid concentration, temperature, and humidity at the current moment are mapped to the instantaneous corrosion rate of the cultural relic using an improved multi-parameter lookup table method, and the instantaneous corrosion rate is cached over time as an instantaneous corrosion rate sequence; the step of calculating the instantaneous corrosion rate using the improved multi-parameter lookup table method in S3 includes: S31: Receive sensor data, including organic acid concentration, temperature, and humidity; receive the material type of cultural relics, including bronze, paper, stone, and textiles, and find the corresponding material corrosion rate pre-stored reference table; S32: Check whether the current environmental organic acid concentration, temperature, and humidity are all within the range of the pre-stored reference table. If the value is within the baseline table range, then the dynamic interpolation lookup process begins, using trilinear interpolation. The formula for calculating the instantaneous corrosion rate R is as follows: Weight normalization: , Where R represents the instantaneous corrosion rate of the artifact, Ri represents the corrosion rate at the i-th vertex in the baseline table, and the weight wi is determined by the position of the current environmental parameter in the pre-stored baseline table: Where wi represents the weight, di represents the standardized distance from the current environmental parameters to each vertex, e represents the smoothing factor, HA, T, R represent the organic acid concentration, temperature and humidity of the current environment, HAi, Ti, RHi represent the organic acid concentration, temperature and humidity values ​​of the i-th vertex in the benchmark table, and ΔHA, ΔT, ΔRH represent the step size of the organic acid concentration, temperature and humidity in the benchmark table. If the value exceeds the baseline range, model compensation extrapolation is triggered, and a simplified physical model is invoked to calculate the corrosion rate. The formula for calculating the instantaneous corrosion rate R using the simplified physical model is as follows: Where R represents the instantaneous corrosion rate of the cultural relic, Rmax represents the corrosion rate corresponding to the maximum boundary of the parameter in the benchmark table, HAmax and RHmax represent the upper limits of the organic acid concentration and humidity parameters in the benchmark table, respectively, n represents the material-related organic acid concentration sensitivity index, k represents the humidity influence coefficient, and α represents the material compensation factor. S33: Output instantaneous corrosion rate R; S4: Read the instantaneous corrosion rate sequence, take the current moment as the end point and back 7 days as the starting point, and use the trapezoidal numerical integration method to calculate the cumulative corrosion exposure for the week; S5: Calculate the health index of cultural relics by calling the instantaneous corrosion rate and the cumulative corrosion exposure of the week at the same time, and combining the temperature, humidity and sensitivity coefficient of the material of the cultural relics; S6: Map the health index of cultural relics together with the cumulative corrosion exposure of the week to a three-level corrosion risk level of red, yellow and green, and generate a heat map of organic acid corrosion risk level through a three-dimensional simulation engine; trigger sensor audible and visual alarms, SMS reminders or automatically start the display case microenvironment adjustment system according to the corrosion risk level.

2. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S1, the distributed sensor includes an organic acid sensor and an environmental parameter sensor. The organic acid sensor is used to collect the concentrations of formic acid and acetic acid in the air and filter out macromolecular interfering substances through a filtering device. The environmental parameters include temperature, humidity, and air flow velocity.

3. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S2, the moving average filtering refers to performing weighted averaging on the data within a sliding time window to suppress random noise. The sliding time window is adjustable.

4. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S4, the trapezoidal numerical integration method means that with the current moment as the end point and seven days back as the starting point, sampling the organic acid concentration data at the sampling interval and calculating the weekly cumulative corrosion exposure D(t). The calculation formula for the weekly cumulative corrosion exposure D(t) is as follows: where i is an integer representing the i-th sampling point in the time series, Ci represents the concentration of organic acid at the i-th sampling, Δt represents the sampling time interval of 5 minutes. When the concentration change per minute is greater than 0.1 ppm, the sampling interval is shortened to 1 minute.

5. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S5, the calculation formula for the cultural relic health index H is: where Rc is the normalized instantaneous corrosion rate, ranging from 0 to 100; Dc is the proportion of cumulative exposure; Ee is the environmental humidity fluctuation coefficient; Mmt is the material sensitivity coefficient; α, β, and γ are the weight coefficients corresponding to the normalized instantaneous corrosion rate, the proportion of cumulative exposure, and the environmental humidity fluctuation coefficient, respectively.

6. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S6, the corrosion risk levels are divided into three types: red, yellow, and green. Red indicates a high-risk level, that is, when the cultural relic health index H ≤ 70 or the weekly cumulative exposure > 30% of the annual threshold; yellow indicates a normal risk level, that is, when 70 < H ≤ 85 and the weekly cumulative exposure ≤ 30% of the annual threshold; green indicates no risk, that is, when the cultural relic health index H > 85 and the weekly cumulative exposure ≤ 5% of the annual threshold.

7. The method for visualizing three-dimensional diffusion simulation of organic acids according to claim 1, characterized in that, In S6, the specific measures triggered by the corrosion risk levels include: High-risk level, triggering the audible and visual alarm of the sensor and immediately sending a signal to start the exhibition cabinet microenvironment regulation system; Normal risk level, sending a text message reminder to the administrator; No risk, indicating normal.

8. A visual three-dimensional diffusion simulation system for organic acids, including a distributed sensor, an edge computing node, and a server, for performing a visual three-dimensional diffusion simulation method for organic acids as described in any one of claims 1 to 7. It is characterized in that: The distributed sensor is deployed indoors in the exhibition hall and is used to collect the concentration of organic acids, temperature, and humidity in the current environment. The distributed sensor includes an organic acid sensor and an environmental parameter sensor. The range of the organic acid sensor is 0 - 5 ppm, the accuracy is ≤ ±0.03 ppm, the cross-sensitivity to ethanol / formaldehyde is ≤ 3%, the response time is ≤ 30 seconds, and the environmental parameter sensor includes a temperature sensor and a humidity sensor; The edge computing node is used to perform moving average filtering on the sensor data, compress and encapsulate the data, and upload the sensor data to the server; The server includes a three-dimensional simulation engine and a user interaction module. The user interaction module includes a cumulative exposure query module, a cultural relic corrosion warning module, and a cultural relic health query module: The cumulative exposure query module is used to query the cumulative corrosion exposure of cultural relics in the current week; The cultural relic health query module is used to query the cultural relic health index. The cultural relic corrosion early warning module is used to trigger sensor audible and visual alarms, send SMS reminders to administrators, or automatically start the display case microenvironment adjustment system based on the risk level of cultural relic corrosion.