Condensation reaction process production safety risk grading early warning method

Abstract

The application discloses a kind of condensed reaction process production with security risk grading early warning method, it is related to the field of chemical production safety, comprising: from the same type condensed reaction normal production data reading reaction input, detection position thermal response and phase change discharge response, according to the thermal response transmission order caused by reaction input, the detection location with continuous receiving relationship is merged to form detection area, and the receiving proportion of detection area to reaction input under each reaction stage is formed, area thermal response benchmark and phase change discharge response benchmark;The reaction input of current condensed reaction, area thermal response and phase change discharge response are obtained, match current reaction stage, and distribute current reaction input according to receiving proportion.The application can stage processing reaction input, detection position thermal response and phase change discharge response in the same type condensed reaction normal production data, so that the monitoring data of current condensed reaction process can be compared under the same reaction stage and the same detection area.

Description

Technical Field

[0001] This invention relates to the field of chemical production safety, specifically a safety risk classification and early warning method for condensation reaction processes. Background Technology

[0002] Condensation reactions are a common type of reaction in the production of fine chemicals, pharmaceutical intermediates, resin materials, and functional chemicals. They are usually accompanied by processes such as feeding, dropping, catalysis, heating, cooling, stirring, and condensation recovery. In actual production, the thermal response inside the reactor is not only affected by the exothermic reaction, but also by the viscosity of the material, the mixing state, the local concentration distribution, the heat exchange conditions, and the phase change discharge process. Therefore, the safety monitoring of the condensation reaction process usually requires simultaneous attention to multiple data such as reactor temperature, heat exchange state, condensation load, steam discharge rate, and reflux rate. Existing production processes typically monitor safety through methods such as temperature, pressure, cooling water flow rate, condensate recovery rate, alarm thresholds, and interlock conditions. These methods can identify significant temperature increases, pressure anomalies, cooling anomalies, or discharge anomalies, and have a good operational foundation under normal operating conditions. However, in partial condensation reactions, localized heat release may initially manifest as changes in phase change discharge responses such as condensation load, steam discharge rate, or reflux rate, while the localized thermal response within the reactor may not immediately show a significant increase. The localized thermal response may only gradually become apparent when the phase change discharge's capacity to absorb localized heat release changes. Therefore, for condensation reaction production processes, relying solely on the thermal response at a single detection location, average reactor temperature, or independent phase change discharge data may not accurately distinguish between fluctuations in normal reaction stages, temporary absence of local thermal response, changes in phase change discharge, and subsequent thermal response switching. To improve the targeted nature of safety risk identification in condensation reaction production processes, a method is needed that can combine reaction input, thermal response at detection location, phase change discharge response, and benchmark relationships between different reaction stages for graded early warning. Therefore, this invention proposes a safety risk graded early warning method for condensation reaction production processes. Summary of the Invention

[0003] The purpose of this invention is to provide a safety risk classification and early warning method for condensation reaction process production, so as to solve the problems mentioned in the background art.

[0004] This invention can be achieved through the following technical solution: a safety risk classification and early warning method for condensation reaction process production, comprising: Step 1: Read the reaction input, detection location thermal response, and phase change discharge response from the normal production data of the isomorphic condensation reaction. According to the thermal response transmission order caused by the reaction input, merge the detection locations with continuous inheritance relationship to form detection areas, and form the inheritance ratio of the detection area to the reaction input, the regional thermal response benchmark, and the phase change discharge response benchmark for each reaction stage. Step 2: Obtain the reaction input, regional thermal response, and phase change discharge response of the current condensation reaction; match the current reaction stage; allocate the current reaction input according to the acceptance ratio; and form the regional input response of each detection area. Step 3: Develop a regional reference thermal response based on the regional input response and the regional thermal response benchmark. When the regional reference thermal response is higher than the current regional thermal response, develop a local thermal response missing amount based on the difference between the regional reference thermal response and the current regional thermal response. Step 4: Based on the portion of the current phase change discharge response that exceeds the phase change discharge response benchmark, a phase change discharge deviation is formed. The phase change discharge deviation is matched with the local thermal response missing amount according to the thermal response transmission order to form the phase change acceptance state of the corresponding detection area. Step 5: When the regional baseline thermal response does not decrease, the phase change discharge deviation decreases during the continuous detection period, and the current regional thermal response increases during the corresponding detection period, a thermal response switching state is formed, and the safety risk level is determined according to the phase change acceptance state, the thermal response switching state, and the progressive changes in the current regional thermal response.

[0005] A further technical improvement of the present invention is that: in step one, the method for forming the detection area and the bearing ratio includes: A detection position response sequence is formed according to the occurrence time of the thermal response at each detection position after the same reaction input; Based on the response interval, response direction relationship, and response change ratio of adjacent detection positions in the detection position response sequence, the continuous connection result between adjacent detection positions is determined; Adjacent detection positions whose continuous acceptance results meet the acceptance benchmark of the corresponding reaction stage are grouped into the same continuous acceptance position group, and adjacent detection positions whose continuous acceptance results do not meet the acceptance benchmark of the corresponding reaction stage are used as the boundary between different continuous acceptance position groups. Each group of consecutive receiving positions is grouped into a detection area, and the receiving ratio of the corresponding detection area to the reaction input is formed according to the proportion of the response change of each detection position relative to the same reaction input within the same group of consecutive receiving positions.

[0006] A further technical improvement of the present invention is that: in step one, the method for forming the regional thermal response benchmark and the phase change discharge response benchmark includes: Read the reaction input, regional thermal response, and phase change discharge response from the normal production data of the isomorphic condensation reaction corresponding to each reaction stage; Based on the proportion of the reaction input received by the detection area, the reaction input is converted into the regional input response of each detection area; Candidate benchmarks for regional thermal response are formed based on the correspondence between the changes in regional input response and regional thermal response, and candidate benchmarks for phase change discharge response are formed based on the correspondence between the changes in regional input response and phase change discharge response. The production section that is a candidate phase change receiving section is defined as the production section whose regional thermal response decreases relative to the regional thermal response candidate benchmark, whose phase change discharge response increases relative to the phase change discharge response candidate benchmark, and whose period of decrease in regional thermal response and period of increase in phase change discharge response satisfy the time correspondence benchmark of the corresponding reaction stage. Candidate phase change receiving sections are eliminated, and a regional thermal response benchmark is formed based on the correspondence between the changes in regional input response and regional thermal response within the eliminated production sections. A phase change discharge response benchmark is also formed based on the correspondence between the changes in regional input response and phase change discharge response within the eliminated production sections.

[0007] A further technical improvement of the present invention is that: in step two, the matching method for the current reaction stage includes: Based on the moment of change in the reaction input of the current condensation reaction, extract the current region's thermal response and the current phase transition discharge response; The current regional thermal response is time-series registered with the regional thermal response benchmarks corresponding to each reaction stage to form regional response registration results. The current phase change discharge response is time-series registered with the phase change discharge response reference corresponding to each reaction stage to form the phase change discharge registration result. When the regional response registration result and the phase change discharge registration result correspond to the same reaction stage, the corresponding reaction stage is determined as the current reaction stage; When the regional response registration result and the phase change discharge registration result correspond to different reaction stages, the current reaction stage is determined based on the period of decrease in the current regional thermal response relative to the regional thermal response benchmark, the period of increase in the current phase change discharge response relative to the phase change discharge response benchmark, and the time sequence corresponding benchmark of the corresponding reaction stage.

[0008] A further technical improvement of the present invention is that: in step two, the method for forming the regional input response includes: Based on the acceptance ratio corresponding to the current reaction stage, the current reaction input is allocated to each detection area to form the candidate area input response for each detection area; According to the thermal response transfer sequence corresponding to the current reaction stage, determine the regional transition period of each detection area relative to the current reaction input; Based on the input response of the candidate region and the regional thermal response benchmark, the regional thermal response that each detection region should generate during the regional acceptance period is determined. The current regional thermal response during the regional takeover period is compared with the expected regional thermal response to form the deviation of the current regional thermal response from the expected regional thermal response. The input response of the candidate region is corrected based on the response deviation to form the regional input response of each detection region.

[0009] A further technical improvement of the present invention is that: in step three, the method for forming the local thermal response deficiency includes: Based on the regional input response and the regional thermal response benchmark, the regional benchmark thermal response of each detection area is formed; The regional baseline thermal response and the current regional thermal response are matched according to the regional transition period. When the regional baseline thermal response is higher than the current regional thermal response, the candidate thermal response missing amount for the corresponding detection area is formed based on the difference between the regional baseline thermal response and the current regional thermal response. The candidate thermal response missing amount in the same detection area is read according to adjacent detection time periods, and the candidate thermal response missing amount in the next detection time period is retained when the candidate thermal response missing amount in the next detection time period is not less than the difference between the candidate thermal response missing amount in the previous detection time period and the thermal response comparison resolution. According to the thermal response transmission order, the thermal responses of adjacent detection areas are correlated with the missing amounts of retained candidate thermal responses in different time periods, and the thermal responses of adjacent detection areas and the missing amounts of retained candidate thermal responses are compared to see if they meet the transition criteria of the corresponding reaction stages. When the thermal response of adjacent detection areas and the amount of missing candidate thermal response do not meet the transition criteria for the corresponding reaction stage, the amount of missing candidate thermal response is determined as the amount of missing local thermal response in the corresponding detection area.

[0010] A further technical improvement of the present invention is that: in step four, the method for matching the phase change discharge deviation with the local thermal response loss includes: The phase change discharge deviation is calculated based on the portion of the current phase change discharge response that exceeds the phase change discharge response baseline. Read the formation time and duration of the local thermal response missing amount in each detection area, and according to the start offset and end offset in the time-corresponding reference of the corresponding reaction stage, shift the formation time and duration of the local thermal response missing amount to form the phase change discharge reading start time and phase change discharge reading end time. The phase change discharge deviation is read from the start time to the end time of phase change discharge reading, and the read phase change discharge deviation is matched with the local thermal response loss in the corresponding detection area according to the detection time period; According to the correspondence between the changes in regional thermal response and phase change discharge response, the read phase change discharge deviation is converted into a phase change discharge conversion quantity with the same dimensions as the missing amount in local thermal response. The phase change discharge conversion comparison resolution is formed based on the phase change discharge comparison resolution and the correspondence between the changes in regional thermal response and phase change discharge response. Calculate the deviation difference between the phase change discharge conversion amount and the local thermal response missing amount; When the deviation difference is not greater than the phase change discharge conversion comparison resolution, the correspondence between the corresponding detection area, the local thermal response missing amount and the phase change discharge deviation amount is retained. When there is only one retained correspondence for the same phase change discharge deviation, the retained correspondence is determined as the phase change acceptance status of the corresponding detection area; When there are more than two retained correspondences for the same phase change discharge deviation, calculate the overlap length between the reading time period of the phase change discharge deviation and the phase change discharge reading start time to phase change discharge reading end time for each detection area, and retain the correspondence with the largest overlap length. When there is only one correspondence with the largest overlap length, the correspondence with the largest overlap length is determined as the phase transition acceptance state of the corresponding detection area; When there are two or more correspondences with the same maximum overlap length, the correspondence with the smallest deviation difference is determined as the phase transition acceptance state of the corresponding detection area.

[0011] A further technical improvement of the present invention is that: in step five, the method for forming the thermal response switching state and the safety risk level includes: Read the regional reference thermal response, phase change discharge deviation, and current regional thermal response corresponding to the detection area that forms the phase change receiving state; When the regional reference thermal response of the next detection period is not less than the difference between the regional reference thermal response and the thermal response comparison resolution of the previous detection period, and the phase change discharge deviation of the next detection period is less than the difference between the phase change discharge deviation and the phase change discharge comparison resolution of the previous detection period, the next detection period is regarded as the phase change transition weakening period. Based on the time-series correspondence of the corresponding reaction stages, the phase transition weakening period is mapped to the detection period of the current region's thermal response; When the thermal response of the current area in the next detection period is greater than the sum of the thermal response of the current area in the previous detection period and the thermal response comparison resolution, a thermal response switching state is formed. Based on the formation of phase change acceptance state, phase change acceptance weakening period and thermal response switching state, the safety risk level of the corresponding detection area is determined, and the progressive change of the current area's thermal response after the formation of thermal response switching state is used as the basis for the progression of safety risk level.

[0012] Compared with the prior art, the present invention has the following beneficial effects: This invention can process the reaction input, detection location thermal response, and phase change discharge response in the normal production data of isomorphic condensation reaction in stages, so that the monitoring data of the current condensation reaction process can be compared under the same reaction stage and the same detection area. As a result, the direct comparison error caused by different reaction stages, different detection locations and different input types can be reduced, and the basic consistency of the safety risk judgment of the condensation reaction production process can be improved. Furthermore, when the local thermal response does not appear as the baseline, the present invention can make corresponding judgments by combining the changes in the phase change discharge response, thereby identifying the situation where local heat release is temporarily taken over by the phase change discharge process. Compared with observing only the single changes in the current region's thermal response or phase change discharge response, the present invention can more clearly distinguish the correlation between the absence of local thermal response and the deviation of phase change discharge, which is beneficial to discovering the development process of more hidden thermal risks in the condensation reaction process. On the other hand, the present invention can further determine the thermal response switching state and form a safety risk level based on the temporal relationship between the phase change acceptance state, the decrease in phase change discharge deviation, and the increase in the current region thermal response. Thus, the risk status in the condensation reaction can be transformed from a single alarm judgment into a graded early warning result, providing production personnel with a clearer risk basis for process observation, process adjustment, and safety handling. Attached Figure Description

[0013] To facilitate understanding by those skilled in the art, the present invention will be further described below with reference to the accompanying drawings.

[0014] Figure 1 This is a schematic diagram of the method logic of the present invention. Detailed Implementation

[0015] To further illustrate the technical means and effects of the present invention in achieving its intended purpose, the following detailed description of the specific implementation methods, structures, features, and effects of the present invention, in conjunction with the accompanying drawings and preferred embodiments, is provided.

[0016] Please see Figure 1 As shown, this invention provides a safety risk classification and early warning method for condensation reaction process production, including: Step 1: Read the reaction input, detection location thermal response, and phase change discharge response from the normal production data of the isomorphic condensation reaction. According to the thermal response transmission order caused by the reaction input, merge the detection locations with continuous inheritance relationship to form detection areas, and form the inheritance ratio of the detection area to the reaction input, the regional thermal response benchmark, and the phase change discharge response benchmark for each reaction stage. In practice, historical normal production batches with the same reactor structure, material system, feeding sequence, stirring method, condensation recovery path, and detection location arrangement as the current condensation reaction are selected as normal production data for the same type of condensation reaction. For each historical normal production batch, reaction input, detection location thermal response, and phase change discharge response are recorded at the same sampling interval. Reaction inputs are read according to input type, including at least one of the following: changes in feed flow rate, dropping rate, catalyst addition, heating input, cooling adjustment, and stirring speed. Detection location thermal responses are measured using one of the following: temperature change, heat flux change, or heat transfer response change. Phase change discharge responses are measured using one of the following: steam discharge rate, condensate recovery rate, condensation load change, return flow rate change, or discharge pipeline temperature change. When different input types have different dimensions, the changes of different input types are not directly added together; instead, a corresponding acceptance ratio, regional thermal response benchmark, and phase change discharge response benchmark are formed for each input type.

[0017] When forming the detection position response sequence, the sampling time when the same input type changes is taken as the starting point. The continuous sampling period closest to this starting point, where no input type changes, is read as the stable period. The stable period includes at least three consecutive sampling intervals. If a stable period meeting the conditions cannot be obtained before the starting point, the input change does not participate in the data formation in step one. The thermal response of each detection position within the stable period is checked. If the thermal response differences of any three consecutive adjacent sampling times are in the same direction and the absolute value of the difference is greater than the corresponding detection accuracy, this continuous sampling period is not considered a stable period, and a new continuous sampling period is read forward. After determining the stable period, the change in thermal response at each detection position relative to the average value of the stable period is calculated. When the absolute value of the change in thermal response at a certain detection position first exceeds the larger of the maximum fluctuation amplitude and the detection accuracy within the stable period, this sampling time is taken as the occurrence time of the thermal response at that detection position. The detection locations are arranged from earliest to latest according to the occurrence time of the thermal response at each detection location, forming a detection location response sequence; when the occurrence time of the thermal response at two or more detection locations is within the same sampling interval, they are arranged from largest to smallest according to the absolute value of the cumulative change of the thermal response at each detection location relative to the average value of the stable period at the end of the sampling interval.

[0018] When determining the continuity between adjacent detection positions, the response sequences of adjacent detection positions are read sequentially. The time difference between the occurrence time of the thermal response at the latter detection position and the occurrence time of the thermal response at the former detection position is calculated as the response interval. The occurrence time of the thermal response at the former detection position is taken as the observation start time, and the end time of the sampling interval following the occurrence time of the thermal response at the latter detection position is taken as the observation end time. The cumulative change of the thermal response at the two detection positions relative to their respective stable period averages is read at the end of the observation time, and the sign of the cumulative change is used to form the response direction relationship. Then, the ratio of the absolute value of the cumulative change at the latter detection position to the absolute value of the cumulative change at the former detection position is calculated to form the response change ratio. The continuity benchmark for the corresponding reaction stage is jointly formed by the range of response interval values, response direction relationships, and response change ratio values ​​of adjacent detection positions within the same reaction stage in the normal production data of the isomorphic condensation reaction. If the response interval, response direction relationship, and response change ratio of adjacent detection positions all meet the acceptance criteria of the corresponding reaction stage, then the adjacent detection positions are determined to be continuously connected; if any one of them is not met, then the adjacent detection positions are determined not to be continuously connected.

[0019] When merging detection regions, the continuous reception results between adjacent detection locations are read along the response sequence of the detection locations. Adjacent detection locations whose continuous reception results meet the reception benchmark of the corresponding reaction stage are grouped into the same continuous reception location group. Adjacent detection locations whose continuous reception results do not meet the reception benchmark of the corresponding reaction stage are used as the boundaries between different continuous reception location groups. Each continuous reception location group is merged to form a detection region, and the detection location identifier, detection location response sequence, original thermal response data of each detection location, response interval, response direction relationship, response change ratio, and reception benchmark of the corresponding reaction stage are retained for subsequent reading of the thermal response of the current region according to the same detection region.

[0020] When determining the acceptance ratio of a detection area to the reaction input, for the same input type change and the same reaction stage, the thermal response change of each detection position within each detection area is read relative to the average value of the stable period at the end of the observation. The absolute values ​​of the thermal response changes of each detection position within the same detection area are summed to form the regional thermal response change of that detection area relative to the input type change; the regional thermal response changes of all detection areas relative to the input type change are summed to form the total thermal response change of all detection areas. The ratio of the regional thermal response change of a single detection area to the total thermal response change of all detection areas is taken as the single-batch acceptance ratio of that detection area to the input type change. The single-batch acceptance ratio is repeatedly calculated for multiple historical normal production batches within the same reaction stage, and the average of the multiple single-batch acceptance ratios is taken to form the acceptance ratio of the detection area to the reaction input under that reaction stage.

[0021] When establishing regional thermal response benchmarks and phase change discharge response benchmarks, the normal production data of the isomorphic condensation reaction is divided into multiple reaction stages according to the process records of historical normal production batches. The reaction input, regional thermal response, and phase change discharge response corresponding to each reaction stage are then read. The regional thermal response is formed by the thermal responses of each detection location within the same detection area. Specifically, the regional weight of each detection location is determined by the proportion of its absolute change in thermal response to the sum of the absolute changes in thermal response at all detection locations within that area. The regional thermal response at the corresponding sampling time is then multiplied by the corresponding regional weight and summed to form the regional thermal response of that detection area at the corresponding sampling time. Based on the proportion of reaction input received by each detection area, the reaction input change of the corresponding input type is allocated to each detection area, forming the regional input response of each detection area.

[0022] For each detection area within the same reaction stage, the regional input response and regional thermal response are read according to the sampling time. Using the regional input response and regional thermal response at the initial sampling time of that reaction stage as starting values, the changes in regional input response and regional thermal response are calculated at each sampling time. At sampling times where the change in regional input response is not zero, the ratio of the change in regional thermal response to the change in regional input response is calculated, and the average of these ratios from multiple sampling times within the same reaction stage is taken. The product of the change in regional input response and this average value forms the candidate baseline change in regional thermal response, which, combined with the starting value of the regional thermal response, forms the candidate baseline for regional thermal response.

[0023] When forming candidate benchmarks for phase change emission response, for the same sampling time within the same reaction stage, the regional input response changes of all detection areas are read, and the regional input response changes of all detection areas are summed to form the total regional input response change. Simultaneously, the phase change emission response is read, and the phase change emission response at the beginning of the sampling time of that reaction stage is used as the starting value for the phase change emission response. The phase change emission response change is calculated for each sampling time. At sampling times where the total regional input response change is not zero, the ratio of the phase change emission response change to the total regional input response change is calculated. The average of the ratios formed at multiple sampling times within the same reaction stage is taken to form the correspondence between the phase change emission response and the total regional input response change. The product of the total regional input response change and this correspondence is used to form the phase change emission response candidate benchmark change, which, combined with the phase change emission response starting value, forms the total phase change emission response candidate benchmark. Then, according to the proportion of the change in regional input response of each detection area to the total change in regional input response, the total amount of phase change discharge response candidate references is allocated to each detection area, forming the phase change discharge response candidate reference component corresponding to each detection area; when the total change in regional input response at a certain sampling time is zero, no phase change discharge response candidate reference component is formed at that sampling time.

[0024] Before determining candidate phase change receiving sections, a thermal response comparison resolution and a phase change discharge comparison resolution are first established. For each detection area, the maximum fluctuation amplitude of the regional thermal response during a stable period is read, and the detection accuracy of each detection position within that area is also read. The detection accuracies of each detection position are weighted and summed according to regional weights to form a regional weighted detection accuracy. When the maximum fluctuation amplitude of the regional thermal response is greater than the regional weighted detection accuracy, the maximum fluctuation amplitude of the regional thermal response is used as the thermal response comparison resolution of that detection area; when the maximum fluctuation amplitude of the regional thermal response is not greater than the regional weighted detection accuracy, the regional weighted detection accuracy is used as the thermal response comparison resolution of that detection area. The maximum fluctuation amplitude of the phase change discharge response and the corresponding detection accuracy during a stable period are read. When the maximum fluctuation amplitude of the phase change discharge response is greater than the corresponding detection accuracy, the maximum fluctuation amplitude of the phase change discharge response is used as the phase change discharge comparison resolution; when the maximum fluctuation amplitude of the phase change discharge response is not greater than the corresponding detection accuracy, the corresponding detection accuracy is used as the phase change discharge comparison resolution.

[0025] When determining candidate phase change receiving sections, the regional thermal response is compared with the regional thermal response candidate benchmark at each sampling time. Continuous sampling periods where the regional thermal response is lower than the regional thermal response candidate benchmark and the difference is greater than the thermal response comparison resolution are designated as the decreasing period of the regional thermal response. Similarly, the phase change discharge response is compared with the total phase change discharge response candidate benchmark at each sampling time. Continuous sampling periods where the phase change discharge response is higher than the total phase change discharge response candidate benchmark and the difference is greater than the phase change discharge comparison resolution are designated as the increasing period of the phase change discharge response. Based on the range of start and end times of the difference between the decreasing and increasing periods of the regional thermal response and the increasing period of the phase change discharge response in multiple historical normal production batches within the same reaction stage, a time-series correspondence benchmark for the corresponding reaction stage is formed. If, within a certain production section, the regional thermal response decreases relative to the regional thermal response candidate benchmark, and the phase change discharge response increases relative to the total phase change discharge response candidate benchmark, and the decreasing period of the regional thermal response and the increasing period of the phase change discharge response satisfy the time-series correspondence benchmark for the corresponding reaction stage, then that production section is determined as a candidate phase change receiving section; otherwise, it is not determined as a candidate phase change receiving section.

[0026] After removing candidate phase change receiving sections, the regional input response, regional thermal response, and phase change discharge response within the retained production section are read. The calculation methods for the regional thermal response candidate benchmark, the total number of phase change discharge response candidate benchmarks, and the phase change discharge response candidate benchmark components corresponding to each detection area are reused to re-form the regional thermal response benchmark, the total number of phase change discharge response benchmarks, and the corresponding phase change discharge response benchmark components for each detection area under each reaction stage. The resulting detection areas, receiving ratios, regional thermal response benchmarks, total number of phase change discharge response benchmarks, corresponding phase change discharge response benchmark components for each detection area, receiving benchmarks, time-series correspondence benchmarks, thermal response comparison resolution, and phase change discharge comparison resolution serve as the reading objects for subsequent matching of the current condensation reaction stage, formation of regional input responses, formation of local thermal response missing amounts, matching of phase change discharge deviations, formation of thermal response switching states, and determination of safety risk levels.

[0027] For example, in a certain reaction stage, three detection locations P1, P2, and P3 are set with a sampling interval of 1 minute. The acceleration rate change of a drop occurs at minute 0. If the input types remain unchanged from minute -5 to minute -2, and P1, P2, and P3 do not show three consecutive adjacent sampling differences in the same direction that are greater than the temperature sampling accuracy, then minute -5 to minute -2 is considered a stable period. If P1, P2, and P3 first exceed the larger value between the maximum fluctuation amplitude of the stable period and the temperature sampling accuracy at minutes 2, 4, and 9, respectively, then the response sequence of the detection locations is P1, P2, and P3. If the response interval in the corresponding reaction stage's acceptance benchmark ranges from 1 to 3 minutes, the allowed response direction relationship is positive, positive, and the response change ratio ranges from 0.3 to 0.8; if the response interval from P1 to P2 is 2 minutes, and at the end of the observation, the cumulative change of P1 is +4, the cumulative change of P2 is +2, and the response change ratio is 0.5, then P1 and P2 satisfy continuous acceptance; if the response interval from P2 to P3 is 5 minutes, it does not satisfy the acceptance benchmark, then P2 and P3 serve as a boundary, thus forming a first detection area including P1 and P2 and a second detection area including P3. Under the change in droplet acceleration rate, if the change in regional thermal response of the first detection area is 6 and the change in regional thermal response of the second detection area is 1, then the single-batch acceptance ratio of the first detection area is 6 / 7, and the single-batch acceptance ratio of the second detection area is 1 / 7. If the thermal response of the first detection area is lower than the regional thermal response candidate benchmark and the difference is greater than the thermal response comparison resolution within the 12th to 16th minute, and the phase change discharge response is higher than the total phase change discharge response candidate benchmark and the difference is greater than the phase change discharge comparison resolution within the 14th to 18th minute, then the difference between the start time and the difference between the end time are both 2 minutes. The range of the difference between the start time and the range of the difference between the end time formed by multiple historical normal production batches are written into the time sequence benchmark of the corresponding reaction stage.

[0028] Step 2: Obtain the reaction input, regional thermal response, and phase change discharge response of the current condensation reaction; match the current reaction stage; allocate the current reaction input according to the acceptance ratio; and form the regional input response of each detection area. In specific implementation, after step one has established the detection area, acceptance ratio, regional thermal response benchmark, total phase change discharge response benchmark, phase change discharge response benchmark components corresponding to each detection area, time-series benchmark, thermal response comparison resolution, and phase change discharge comparison resolution, the reaction input of the current condensation reaction, the thermal response at each detection location, and the phase change discharge response are obtained. The current reaction input is read according to the same input type as in step one; the thermal response at the current detection location uses the same data type as the regional thermal response benchmark; the current phase change discharge response uses the same data type as the total phase change discharge response benchmark. According to the detection area and regional weights established in step one, the current thermal responses at each detection location within the same detection area are synthesized into the current regional thermal response, so that the current regional thermal response and the regional thermal response benchmark have the same detection area, the same sampling time, and the same thermal response dimensions.

[0029] When matching the current reaction stage, the sampling time at which the current reaction input changes is taken as the current reaction input change time. When two or more input types change at the same sampling time, that sampling time is still taken as the current reaction input change time, but the changes of different input types are not directly added together. The collected sampling length between the current reaction input change time and the current sampling time is read, and the minimum sampling length of all reaction stages in the normal production data of the isomorphic condensation reaction is read. When the collected sampling length is not less than the minimum sampling length, the minimum sampling length is taken as the unified registration sampling length of all reaction stages; when the collected sampling length is less than the minimum sampling length, the collected sampling length is taken as the unified registration sampling length of all reaction stages. The total amount of the current regional thermal response, the current phase change discharge response, the regional thermal response benchmark, and the phase change discharge response benchmark is extracted according to the unified registration sampling length.

[0030] The current regional thermal response is time-series registered with the regional thermal response benchmarks corresponding to each reaction stage. Specifically, the time of change of the current reaction input is aligned with the starting sampling time of each reaction stage. The absolute value of the difference between the current regional thermal response and the regional thermal response benchmark is calculated for each sampling time. The absolute value of the difference is then divided by the thermal response comparison resolution of the corresponding reaction stage to form the thermal response registration difference. The average value of the thermal response registration differences for each detection area and each sampling time is taken to form the average difference of the regional response for the corresponding reaction stage. The reaction stage with the smallest average difference of regional response is taken as the regional response registration result.

[0031] The current phase change emission response is time-series registered with the total phase change emission response reference for each reaction stage. Specifically, at each registration sampling time, the phase change emission response reference component corresponding to each detection area is read, and the phase change emission response reference components corresponding to each detection area are summed to form the total phase change emission response reference for registration at that sampling time. The absolute value of the difference between the current phase change emission response and the total phase change emission response reference for registration is calculated at each sampling time, and the absolute value of the difference is divided by the phase change emission comparison resolution of the corresponding reaction stage to form the phase change emission registration difference. The average value of the phase change emission registration difference at each sampling time is taken to form the average phase change emission difference for the corresponding reaction stage, and the reaction stage with the smallest average phase change emission difference is taken as the phase change emission registration result.

[0032] When the regional response registration result and the phase change emission registration result correspond to the same reaction stage, that reaction stage is determined as the current reaction stage. When the regional response registration result and the phase change emission registration result correspond to different reaction stages, each reaction stage is used as the object to be verified. The regional thermal response benchmark, the total amount of the phase change emission response benchmark, the thermal response comparison resolution, the phase change emission comparison resolution, and the time-series correspondence benchmark corresponding to the same reaction stage are read to form a decreasing period when the current regional thermal response is lower than the regional thermal response benchmark and an increasing period when the current phase change emission response is higher than the total amount of the phase change emission response benchmark. If the difference between the start time and the difference between the end time of the decreasing period and the increasing period satisfy the time-series correspondence benchmark of the corresponding reaction stage, then that reaction stage is regarded as the reaction stage that satisfies the time-series correspondence benchmark. If only one reaction stage meets the timing correspondence benchmark, that reaction stage is determined as the current reaction stage; if two or more reaction stages meet the timing correspondence benchmark, the reaction stage with the smallest sum of regional response average difference and phase change discharge average difference is determined as the current reaction stage; if no reaction stage meets the timing correspondence benchmark, the reaction stage with the smallest sum of regional response average difference and phase change discharge average difference is determined as the current reaction stage.

[0033] After determining the current reaction stage, the corresponding acceptance ratio for the current reaction stage is read. For the same input type, the continuous sampling period closest to the moment of change in the current reaction input, where no input type has changed, is taken as the current stable sampling period, and the average value of the current reaction input within the current stable sampling period is taken as the starting value of the current input. When at least three consecutive sampling intervals of the current stable sampling period cannot be obtained, the current reaction input value at the previous sampling moment is taken as the starting value of the current input. The change in the current reaction input relative to the starting value of the current input is multiplied by the acceptance ratio of the corresponding detection area under the current reaction stage to form the candidate region input response for that detection area. When two or more input types change simultaneously, candidate region input response components for the corresponding input types are formed and saved separately according to input type.

[0034] When determining the regional transition period, the response sequence and continuous transition results of the detection positions corresponding to the current reaction stage are read. The detection position with the earliest occurrence time in each detection region is taken as the transition start point of that detection region. The detection regions are arranged from earliest to latest according to the occurrence time of their transition start points, forming a thermal response transmission order at the detection region level. Based on the difference in occurrence time of each detection region's transition start point relative to the reaction input change time in the normal production data of the isomorphic condensation reaction, a transition start difference value is formed for each detection region. The current reaction input change time plus the transition start difference value is taken as the start time of the regional transition period, and the start time of the regional transition period of the next detection region is taken as the end time of the regional transition period of the current detection region. If there is no next detection region, the current reaction input change time plus the average duration of the current reaction stage is taken as the end time of the regional transition period. When the end time of the regional transition period is not later than the start time of the regional transition period, the end time of the regional transition period is adjusted to the time corresponding to one sampling interval after the start time of the regional transition period.

[0035] When generating the required regional thermal response, the candidate regional input response and the regional thermal response benchmark for the corresponding detection area under the current reaction stage are read, and the current regional thermal response at the start of the regional transition period is used as the starting value for the current regional thermal response. The ratio of the change in regional thermal response to the change in regional input response in the regional thermal response benchmark is read, and the candidate regional input response is multiplied by this ratio to generate the required change in regional thermal response. This is then added to the current starting value for regional thermal response to generate the regional thermal response that the detection area should generate during the regional transition period. When the candidate regional input response includes candidate regional input response components of two or more input types, the required change in regional thermal response for each input type is calculated separately and then summed.

[0036] The current regional thermal response and the expected regional thermal response are compared at each sampling time within the regional transition period to form a response deviation sequence. The response deviation is the current regional thermal response minus the expected regional thermal response. The response deviation with the largest absolute value in the response deviation sequence is used as the response deviation for correction. When the absolute value of the response deviation for correction is not greater than the thermal response comparison resolution, the candidate regional input response is used as the regional input response. When the response deviation for correction is greater than the thermal response comparison resolution, the response deviation for correction is converted into an input correction amount according to the ratio of the change in regional thermal response to the change in regional input response, and the input correction amount is added to the candidate regional input response to form the regional input response. When the response deviation for correction is less than the negative value of the thermal response comparison resolution, the candidate regional input response is retained as the regional input response, and in subsequent step three, the difference between the current regional thermal response and the expected regional thermal response forms the basis for comparison of the local thermal response missing amount.

[0037] For example, if the drop rate change occurs at minute 0 of the current condensation reaction, and data has been collected up to minute 6, the current regional thermal response and current phase change discharge response are time-registered with the regional thermal response benchmarks and phase change discharge response benchmarks corresponding to each reaction stage, respectively. If the average difference in regional response and the average difference in phase change discharge for reaction stage A are both less than those for reaction stage B, then reaction stage A is determined as the current reaction stage. The proportion of the drop rate change received by the first and second detection regions under reaction stage A is read. If the current drop rate change is 14, the proportion received by the first detection region is 6 / 7, and the proportion received by the second detection region is 1 / 7, then the candidate region input response for the first detection region is 12, and the candidate region input response for the second detection region is 2. If the regional transition period of the first detection area is from the 2nd minute to the 4th minute, and the ratio of the change in regional thermal response to the change in regional input response is 0.5, then the change in regional thermal response that the first detection area should generate is 6. When the response deviation of the first detection area at the corresponding sampling time is 1 and the thermal response comparison resolution is 0.5, the response deviation is converted into an input correction amount of 2, and the regional input response of the first detection area is corrected to 14.

[0038] Step 3: Develop a regional reference thermal response based on the regional input response and the regional thermal response benchmark. When the regional reference thermal response is higher than the current regional thermal response, develop a local thermal response missing amount based on the difference between the regional reference thermal response and the current regional thermal response. In practice, after establishing the current reaction stage, the regional input response of each detection area, the current regional thermal response, the regional transition period, the thermal response transmission order of the detection area hierarchy, the transition benchmark, and the thermal response comparison resolution, a regional benchmark thermal response for each detection area is formed. For each detection area, the regional input response, regional thermal response benchmark, and regional transition period are read in the current reaction stage. Within the regional transition period, the ratio of the change in regional thermal response to the change in regional input response in the regional thermal response benchmark is read one by one according to the sampling time. When the regional input response corresponds to only one input type, and the corresponding ratio at that sampling time has been formed, the regional input response is multiplied by the corresponding ratio at that sampling time to form the regional benchmark thermal response change at that sampling time, and this is added to the current regional thermal response at the start of the regional transition period to form the regional benchmark thermal response for that detection area at that sampling time. When the corresponding ratio has not been formed, the regional benchmark thermal response at that sampling time is not formed, and that sampling time is not included in the comparison of missing candidate thermal responses. When the regional input response includes regional input response components of two or more input types, the regional reference thermal response change of each input type at the same sampling time is formed according to the ratio corresponding to each input type. If the ratio corresponding to any non-zero regional input response component is not formed, the regional reference thermal response at that sampling time is not formed. If the ratio corresponding to each non-zero regional input response component has been formed, the regional reference thermal response change at the same sampling time is summed and added to the current regional thermal response at the start of the regional transition period to form the regional reference thermal response of the detection area at that sampling time.

[0039] After establishing the regional baseline thermal response, the regional baseline thermal response and the current regional thermal response are correlated according to the regional transition period. Using the start and end times of the regional transition period for the same detection area as the comparison range, the established regional baseline thermal response and the current regional thermal response are read at each sampling time. When the regional baseline thermal response at a certain sampling time is higher than the current regional thermal response, and the difference is greater than the thermal response comparison resolution, this difference is taken as the missing candidate thermal response for that detection area at that sampling time, and the positive sign of the difference between the regional baseline thermal response and the current regional thermal response is used as the missing direction. When the regional baseline thermal response is not higher than the current regional thermal response, or the difference is not greater than the thermal response comparison resolution, no missing candidate thermal response is generated for that sampling time. If no missing candidate thermal response is generated for the same detection area within the regional transition period, that detection area will not proceed to subsequent continuous screening.

[0040] For candidate thermal response missing values ​​formed within the same detection area, continuous screening is performed according to the sampling time from earliest to latest. The earliest formed candidate thermal response missing value is used as the starting candidate thermal response missing value for the current continuous candidate time period and is retained. The next detection time period is then read. If a candidate thermal response missing value is formed in the next detection time period, and the candidate thermal response missing value in the next detection time period is not less than the difference between the candidate thermal response missing value in the previous retained detection time period and the thermal response comparison resolution, the candidate thermal response missing value in the next detection time period is retained. If no candidate thermal response missing value is formed in the next detection time period, or the candidate thermal response missing value in the next detection time period is less than the difference between the candidate thermal response missing value in the previous retained detection time period and the thermal response comparison resolution, the current continuous candidate time period ends, and the earliest formed candidate thermal response missing value is reread from subsequent sampling times as the starting candidate thermal response missing value for the new continuous candidate time period. Candidate thermal response missing values ​​that are not retained are not included in the determination of local thermal response missing values.

[0041] After completing the continuity screening, the retained candidate thermal response missing quantities are matched with the thermal responses of adjacent detection areas according to the thermal response transmission order of the detection area hierarchy. For each retained candidate thermal response missing quantity in the current detection area, the adjacent detection area located one position after the current detection area in the thermal response transmission order is read; the response interval value range in the corresponding reference for the current reaction stage is read; the sampling time of the retained candidate thermal response missing quantity is added to the lower limit of the response interval value range as the reading start time of the adjacent detection area; the sampling time of the retained candidate thermal response missing quantity is added to the upper limit of the response interval value range as the reading end time of the adjacent detection area. When the reading end time of the adjacent detection area is later than the current acquisition time, the reference comparison is not performed at the current sampling time, and the retained candidate thermal response missing quantity, reading start time, and reading end time are saved; when the current acquisition time is not earlier than the reading end time, the thermal response of the current area of ​​the adjacent detection area is read and the reference comparison is performed.

[0042] When performing the baseline comparison, the sampling time of the retained candidate thermal response missing quantity is used as the starting time for the thermal response change of the adjacent detection area. From the start time to the end time of reading, the change in the current region's thermal response relative to the starting time is calculated for each sampling time. If the change at a certain sampling time is greater than the thermal response comparison resolution for the first time, that sampling time is taken as the response time of the adjacent detection area, and the response direction of the adjacent detection area is recorded as the positive direction. If the change at a certain sampling time is less than the negative value of the thermal response comparison resolution for the first time, that sampling time is taken as the response time of the adjacent detection area, and the response direction of the adjacent detection area is recorded as the negative direction. If there is no change greater than the thermal response comparison resolution or less than the negative value of the thermal response comparison resolution for the first time from the start time to the end time of reading, it is determined that the thermal response of the adjacent detection area and the retained candidate thermal response missing quantity do not meet the baseline for the corresponding reaction stage.

[0043] When adjacent detection regions form response times, the time difference between the response time of the adjacent detection region and the sampling time of the retained candidate thermal response missing quantity is used as the response interval; the missing direction and the response direction of the adjacent detection region together form the response direction relationship; the ratio of the absolute value of the thermal response change in the current region of the adjacent detection region to the absolute value of the retained candidate thermal response missing quantity is used as the response change ratio. Since the retained candidate thermal response missing quantity is formed by a difference greater than the thermal response comparison resolution, the absolute value of the retained candidate thermal response missing quantity is greater than zero. The response interval, response direction relationship, and response change ratio are compared with the acceptance benchmark corresponding to the current reaction stage; if the response interval, response direction relationship, and response change ratio all meet the acceptance benchmark, the retained candidate thermal response missing quantity is not determined as a local thermal response missing quantity; if any one of them does not meet the acceptance benchmark, the retained candidate thermal response missing quantity is determined as a local thermal response missing quantity of the corresponding detection region. When the current detection region does not have a subsequent adjacent detection region in the thermal response transmission order, the retained candidate thermal response missing quantity is determined as a local thermal response missing quantity of the corresponding detection region.

[0044] After local thermal response deficiencies are identified, they are grouped according to the same detection area and continuous sampling time. The earliest sampling time within the continuous sampling time to form a local thermal response deficiencies is taken as the formation sampling time, and the latest sampling time to form a local thermal response deficiencies is taken as the end sampling time. The period between the formation sampling time and the end sampling time is taken as the continuous sampling period for local thermal response deficiencies. Simultaneously, the local thermal response deficiencies corresponding to each sampling time within the continuous sampling period are saved. The corresponding detection area, formation sampling time, continuous sampling period, local thermal response deficiencies corresponding to each sampling time, and the corresponding current reaction stage are saved for subsequent matching with phase change discharge deviations based on time-series benchmarks.

[0045] For example, using two detection areas P1, P2, and P3, the first detection area includes P1 and P2, and the second detection area includes P3. If the current reaction phase's regional transition period for the first detection area is from minute 2 to minute 4, the regional input response is 14, and the ratio of the change in regional thermal response to the change in regional input response is 0.4, 0.5, and 0.5 in minutes 2, 3, and 4, respectively, then the corresponding regional baseline thermal response changes for the first detection area are 5.6, 7, and 7, respectively. If the difference between the current regional thermal response and the regional baseline thermal response in minute 2 is 0.2, and the thermal response comparison resolution is 0.5, then no candidate thermal response missing quantity is formed; if the difference is 1.2 in minute 3 and 1.1 in minute 4, then candidate thermal response missing quantities of 1.2 and 1.1 are formed in minutes 3 and 4, respectively, and the missing direction is recorded as the positive direction. Since the missing candidate thermal response value of 1.1 at the 4th minute is not less than the difference between the missing candidate thermal response value of 1.2 at the 3rd minute and the thermal response comparison resolution of 0.5, the missing candidate thermal response value at the 4th minute is retained. If the response interval in the benchmark ranges from 1 minute to 3 minutes, then the 4th to 6th minutes are taken as the reading period for the second detection area. If the change in thermal response in the current area of ​​the second detection area is less than the negative value of the thermal response comparison resolution for the first time during the reading period, then the response direction of the adjacent detection area is negative. If the response direction relationship formed by the missing direction does not meet the benchmark, the missing candidate thermal response values ​​retained in the 3rd to 4th minutes of the first detection area are merged to form a local missing thermal response value, and the sampling time, continuous sampling period, and the local missing thermal response value corresponding to each sampling time are saved.

[0046] Step 4: Based on the portion of the current phase change discharge response that exceeds the phase change discharge response benchmark, a phase change discharge deviation is formed. The phase change discharge deviation is matched with the local thermal response missing amount according to the thermal response transmission order to form the phase change acceptance state of the corresponding detection area. Specifically, after establishing the current reaction stage, current phase change emission response, total phase change emission response baseline, phase change emission response baseline components corresponding to each detection area, phase change emission comparison resolution, time-series corresponding baseline, local thermal response missing value, sampling time, and continuous sampling period, a phase change emission deviation is formed. For the current reaction stage, the current phase change emission response is read at each sampling time, and the phase change emission response baseline components corresponding to each detection area at that sampling time are read. The phase change emission response baseline components are the phase change emission response baseline values ​​formed in step one by allocating the total phase change emission response baseline to the corresponding detection areas according to the proportion of the regional input response change in each detection area to the total regional input response change. The phase change emission response baseline components corresponding to each detection area within the same sampling time are summed to form the total phase change emission response baseline for that sampling time. When the current phase change emission response is higher than the total phase change emission response baseline, and the difference is greater than the phase change emission comparison resolution, the difference is taken as the phase change emission deviation for that sampling time; when the difference is not greater than the phase change emission comparison resolution, no phase change emission deviation is formed for that sampling time. The phase change discharge deviations formed within consecutive sampling times are merged to form a phase change discharge deviation sequence, and the start sampling time, end sampling time, and corresponding phase change discharge deviation values ​​of the phase change discharge deviation sequence are saved.

[0047] When matching phase change discharge deviation with local thermal response missing amount, the detection areas that form local thermal response missing amounts are read sequentially according to the thermal response transmission order of the detection area level. For each detection area, the sampling time of the formation of local thermal response missing amount, the continuous sampling period, and the local thermal response missing amount value corresponding to each sampling time within the continuous sampling period are read; then, the start offset and end offset in the time-corresponding reference corresponding to the current reaction stage are read. The sampling time of the formation of local thermal response missing amount is added to the start offset to form the phase change discharge reading start time; the end sampling time of the continuous sampling period of local thermal response missing amount is added to the end offset to form the phase change discharge reading end time. When the phase change discharge reading end time is later than the currently acquired time, the phase change discharge deviation matching corresponding to that detection area is not performed temporarily, and the corresponding detection area, local thermal response missing amount value, phase change discharge reading start time, and phase change discharge reading end time are saved; when the currently acquired time is not earlier than the phase change discharge reading end time, the phase change discharge deviation reading and matching are performed.

[0048] During the period from the start to the end of phase change emission reading, the phase change emission deviation sequence is read. Using the sampling sequence within the continuous sampling period of the local thermal response missing quantity as the corresponding benchmark, the first sampling moment within this period is mapped to the start of the phase change emission reading, and subsequent sampling moments are mapped sequentially at the same sampling interval. If no phase change emission deviation is generated for a phase change emission sampling moment corresponding to a certain local thermal response missing quantity sampling moment, no deviation difference is generated for that sampling moment; if the number of sampling moments generating deviation differences is zero, no correspondence is established between the detection area and the phase change emission deviation sequence. Phase change emission deviations exceeding the sampling sequence within the continuous sampling period of the local thermal response missing quantity during the phase change emission reading period are only used for subsequent overlap length calculations.

[0049] When a sampling time is available for calculation, the correspondence between the changes in regional thermal response and phase change discharge response of the corresponding detection area in step one is read under the current reaction stage. This correspondence is formed by the ratio between the change in regional thermal response of the corresponding detection area and the change in the reference component of the phase change discharge response of the corresponding detection area within the production section after removing candidate phase change receiving sections in step one. When the change in the reference component of the phase change discharge response of the corresponding detection area is zero, no correspondence is formed for that sampling time. For each sampling time involved in the calculation, if a corresponding correspondence has already been formed, it is directly read; if no corresponding correspondence has been formed, the closest formed correspondence to that sampling time in the current reaction stage for the same detection area is used; if no formed correspondence exists for the detection area in the current reaction stage, no correspondence is formed between the detection area, the missing amount of local thermal response, and the phase change discharge deviation sequence.

[0050] When a readable correspondence exists between the changes, the phase change discharge deviation at the corresponding sampling time is multiplied by the correspondence to form a phase change discharge conversion quantity in the same thermal response dimension as the missing local thermal response quantity. The phase change discharge comparison resolution is multiplied by the absolute value of the correspondence to form a phase change discharge conversion comparison resolution in the thermal response dimension. The absolute value of the difference between the phase change discharge conversion quantity and the corresponding missing local thermal response quantity is calculated for each sampling time, and this is taken as the deviation difference. When the deviation difference at each sampling time involved in the calculation is not greater than the phase change discharge conversion comparison resolution at the corresponding sampling time, the correspondence between the corresponding detection area, the missing local thermal response quantity, and the phase change discharge deviation quantity sequence is retained, and the average of the deviation differences at each sampling time involved in the calculation is taken to form the average deviation difference of the retained correspondence. When the deviation difference at any sampling time involved in the calculation is greater than the phase change discharge conversion comparison resolution at the corresponding sampling time, the correspondence is not retained.

[0051] For the same phase change emission deviation sequence, all retained correspondences are read. When only one retained correspondence exists for the same phase change emission deviation sequence, that retained correspondence is determined as the phase change acceptance state of the corresponding detection area. When two or more retained correspondences exist for the same phase change emission deviation sequence, the overlap length between the time interval from the start sampling time to the end sampling time of the phase change emission deviation sequence and the time interval from the start to the end of the phase change emission reading for each detection area is calculated. The overlap length is calculated based on the number of overlapping sampling times, and both the start and end sampling times of the overlapping time interval are included in the overlap length. The correspondence with the largest overlap length is retained. When only one correspondence with the largest overlap length exists, the correspondence with the largest overlap length is determined as the phase change acceptance state of the corresponding detection area. When two or more correspondences with the same largest overlap length exist, the correspondence with the smallest average deviation difference is determined as the phase change acceptance state of the corresponding detection area. When the average deviation difference is also the same, the correspondence corresponding to the detection area with the highest thermal response transmission order is retained.

[0052] After the phase change acceptance state is formed, the detection area, local thermal response missing value, phase change discharge deviation sequence, phase change discharge reading start time, phase change discharge reading end time, average deviation difference and corresponding current reaction stage are saved for subsequent formation of thermal response switching state and determination of safety risk level.

[0053] For example, using two detection areas formed by P1, P2, and P3, the first detection area includes P1 and P2, and the second detection area includes P3. If a local thermal response gap is formed in the first detection area between the 3rd and 4th minute, and the start and end offsets in the time-series reference are both 2 minutes, then the phase change discharge reading start time is the 5th minute, and the phase change discharge reading end time is the 6th minute. If phase change discharge deviations are formed in both the 5th and 6th minutes, then the local thermal response gaps in the 3rd and 4th minutes are mapped to the phase change discharge deviations in the 5th and 6th minutes, respectively, according to the sampling sequence. If the correspondence between the change in regional thermal response and phase change emission response in the first detection area at the current reaction stage is 0.5, then the phase change emission deviations at the 5th and 6th minutes are multiplied by 0.5 to form the phase change emission conversion. When the deviation difference between the phase change emission conversion and the corresponding missing local thermal response is not greater than the phase change emission conversion comparison resolution, the correspondence between the first detection area, the missing local thermal response, and the phase change emission deviation sequence is retained. If the same phase change emission deviation sequence also forms a retained correspondence with the second detection area, then the overlap length between the phase change emission deviation sequence time period and the corresponding reading time periods of the two detection areas is compared, and the phase change acceptance status is determined accordingly.

[0054] Step 5: When the regional baseline thermal response does not decrease, the phase change discharge deviation decreases during the continuous detection period, and the current regional thermal response increases during the corresponding detection period, a thermal response switching state is formed, and the safety risk level is determined according to the phase change acceptance state, thermal response switching state, and the progressive change of the current regional thermal response. In practice, after establishing the phase change acceptance state, corresponding detection area, local thermal response missing value, continuous sampling period of local thermal response missing value, phase change discharge deviation sequence, phase change discharge reading start time, phase change discharge reading end time, current reaction stage, regional reference thermal response, current regional thermal response, thermal response comparison resolution, and phase change discharge comparison resolution, a thermal response switching state and safety risk level are established. For each detection area that has established a phase change acceptance state, the corresponding phase change discharge deviation sequence, continuous sampling period of local thermal response missing value, phase change discharge reading start time, and sampling interval are read.

[0055] Using the sampling sequence within the continuous sampling period of the local thermal response missing quantity as the corresponding benchmark, the first sampling moment within this period is mapped to the start time of phase change emission reading. Subsequent sampling moments are then mapped sequentially at the same sampling interval, forming a sampling sequence correspondence between the phase change emission deviation sequence and the local thermal response missing quantity. For each set of sampling sequence correspondences, the phase change emission sampling moment, the local thermal response missing quantity sampling moment, the regional benchmark thermal response of the same detection area at that local thermal response missing quantity sampling moment, and the current regional thermal response are saved. If a phase change emission deviation is not formed at a certain phase change emission sampling moment, that phase change emission sampling moment is not included in the phase change acceptance weakening judgment.

[0056] When forming the phase transition attenuation sampling time, the effective phase transition discharge sampling times are read from earliest to latest along the sampling time sequence of phase transition discharge deviation. An effective phase transition discharge sampling time is the sampling time within the phase transition discharge deviation sequence where a phase transition discharge deviation has already occurred and can be mapped to the sampling time of the local thermal response missing amount according to the sampling sequence. If there are fewer than two effective phase transition discharge sampling times, a phase transition attenuation sampling time is not formed. If there are more than two effective phase transition discharge sampling times, two adjacent effective phase transition discharge sampling times with a time difference equal to one sampling interval are read; when the time difference between two effective phase transition discharge sampling times is greater than one sampling interval, these two effective phase transition discharge sampling times are not compared as adjacent detection periods.

[0057] For two adjacent effective phase change emission sampling times, the phase change emission deviation of the preceding and following effective phase change emission sampling times are read respectively, and the regional reference thermal responses corresponding to the preceding and following effective phase change emission sampling times are read according to the sampling sequence correspondence. When the regional reference thermal response corresponding to the following effective phase change emission sampling time is not less than the difference between the regional reference thermal response and the thermal response comparison resolution corresponding to the preceding effective phase change emission sampling time, and the phase change emission deviation of the following effective phase change emission sampling time is less than the difference between the phase change emission deviation and the phase change emission comparison resolution of the preceding effective phase change emission sampling time, the following effective phase change emission sampling time is determined as the phase change transition weakening sampling time, and the preceding effective phase change emission sampling time compared with the following effective phase change emission sampling time is saved as the previous sampling time of phase change transition weakening; if any condition is not met, the following effective phase change emission sampling time is not determined as the phase change transition weakening sampling time.

[0058] After establishing the phase transition attenuation sampling times, these times are merged from earliest to latest. The first phase transition attenuation sampling time is taken as the starting sampling time of the current phase transition attenuation period. Continuing to read subsequent phase transition attenuation sampling times, if the time difference between the next and previous sampling times is equal to one sampling interval, the next sampling time is added to the current phase transition attenuation period. If the time difference is greater than one sampling interval, the current phase transition attenuation period ends, and the next sampling time is taken as the starting sampling time of a new phase transition attenuation period. The earliest sampling time within each phase transition attenuation period is taken as the start time, and the latest sampling time is taken as the end time. The phase transition discharge deviation values ​​corresponding to each sampling time within the phase transition attenuation period are saved. If a phase transition weakening sampling moment is not formed, then a phase transition weakening period will not be formed.

[0059] After the phase change attenuation period is established, the sampling time before and after the phase change attenuation is mapped to the detection period of the current region's thermal response according to the sampling sequence. Specifically, for each phase change attenuation sampling time, the sampling time before the phase change attenuation corresponding to that sampling time is read, and the sampling times of the local thermal response missing amount corresponding to the sampling time before and after the phase change attenuation are read according to the sampling sequence. The two sampling times of the local thermal response missing amount are used as the detection periods corresponding to the current region's thermal response. The detection periods corresponding to the current region's thermal response are arranged from earliest to latest according to the sampling time, and duplicate sampling times are removed to form the sequence of detection periods corresponding to the current region's thermal response.

[0060] When a thermal response switching state is established, the current thermal response of the same detection area is read according to the detection time sequence corresponding to the current regional thermal response. If the detection time sequence corresponding to the current regional thermal response has fewer than two sampling times, no thermal response switching state judgment is performed; if the detection time sequence corresponding to the current regional thermal response includes more than two sampling times, the detection time sequences corresponding to two adjacent current regional thermal responses are compared. If the current regional thermal response of the later detection time is greater than the sum of the current regional thermal response of the previous detection time and the thermal response comparison resolution, a thermal response switching state is established in the later detection time; if the current regional thermal response of the later detection time is not greater than the sum of the current regional thermal response of the previous detection time and the thermal response comparison resolution, a thermal response switching state is not established in the later detection time.

[0061] After a thermal response switching state is established, the next detection period after the thermal response switching state is established is used as the starting sampling time for progressive change. Sampling times are then read sequentially along the corresponding detection period sequence for the current region's thermal response. If the current region's thermal response in the next detection period is not less than the sum of the current region's thermal response and the thermal response comparison resolution in the previous detection period, the next detection period is merged into the current progressive change period. If the current region's thermal response in the next detection period is less than the sum of the current region's thermal response and the thermal response comparison resolution in the previous detection period, the current progressive change period ends. If there are sampling times between two or more thermal response switching states that do not meet the progressive change conditions, two or more progressive change periods for the current region's thermal response are formed respectively. The starting sampling time, ending sampling time, and current region thermal response values ​​corresponding to each sampling time within each progressive change period are saved. If a thermal response switching state is not established, a progressive change period for the current region's thermal response is not formed.

[0062] When determining the safety risk level, the phase change acceptance state, phase change weakening period, thermal response switching state, and the progressive change period of the current area's thermal response are read according to the same detection area. If a phase change acceptance state is formed but a phase change weakening period is not formed, the corresponding detection area is determined to be at the first safety risk level; if a phase change acceptance state and a phase change weakening period are formed but a thermal response switching state is not formed, the corresponding detection area is determined to be at the second safety risk level; if a phase change acceptance state, a phase change weakening period, and a thermal response switching state are formed, the corresponding detection area is determined to be at the third safety risk level, and the progressive change period of the current area's thermal response is used as the basis for the progression of the third safety risk level. After the safety risk level is determined, the corresponding detection area, safety risk level, phase change acceptance state, phase change weakening period, thermal response switching state, progressive change period of the current area's thermal response, the current area's thermal response value corresponding to each sampling time within the progressive change period, and the corresponding current reaction stage are saved as the current condensation reaction safety risk classification and early warning result.

[0063] For example, if the first detection area has already formed a phase change receiving state, and the corresponding phase change discharge deviation sequence is located between the 5th and 7th minutes, and the continuous sampling period for the local thermal response missing amount is from the 3rd to the 5th minute, then the phase change discharge sampling times of the 5th, 6th, and 7th minutes are respectively mapped to the sampling times of the local thermal response missing amount of the 3rd, 4th, and 5th minutes according to the sampling sequence. When it is determined that the phase change receiving is weakening, the phase change discharge deviation from the 5th to the 7th minute is read, and the regional reference thermal response from the 3rd to the 5th minute is read. If the phase change discharge deviation at the 6th minute is less than the difference between the phase change discharge deviation at the 5th minute and the phase change discharge comparison resolution, and the regional reference thermal response at the 4th minute is not less than the difference between the regional reference thermal response at the 3rd minute and the thermal response comparison resolution, then the 6th minute becomes the sampling time for weakening phase change receiving, and the 5th minute is saved as the sampling time before the weakening of phase change receiving. Subsequently, the sampling time of phase change emission at the 5th minute is mapped to the sampling time of missing local thermal response at the 3rd minute according to the sampling sequence, and the current thermal response of the first detection area is read at the 3rd minute; the sampling time of phase change emission at the 6th minute is mapped to the sampling time of missing local thermal response at the 4th minute according to the sampling sequence, and the current thermal response of the first detection area is read at the 4th minute. If the current thermal response of the first detection area at the 4th minute is greater than the sum of the current thermal response of the first detection area at the 3rd minute and the thermal response comparison resolution, then the first detection area enters a thermal response switching state, and the safety risk level is determined accordingly.

[0064] The above description is merely a preferred embodiment of the present invention and is not intended to limit the present invention in any way. Although the present invention has been disclosed above with reference to preferred embodiments, it is not intended to limit the present invention. Any person skilled in the art can make some modifications or alterations to the above-disclosed technical content to create equivalent embodiments without departing from the scope of the present invention. Any simple modifications, equivalent changes and alterations made to the above embodiments based on the technical essence of the present invention without departing from the scope of the present invention shall still fall within the scope of the present invention.

Claims

1. A method for safety risk classification and early warning in condensation reaction process production, characterized in that, include: Step 1: Read the reaction input, detection location thermal response, and phase change discharge response from the normal production data of the isomorphic condensation reaction. According to the thermal response transmission order caused by the reaction input, merge the detection locations with continuous inheritance relationship to form detection areas, and form the inheritance ratio of the detection area to the reaction input, the regional thermal response benchmark, and the phase change discharge response benchmark for each reaction stage. Step 2: Obtain the reaction input, regional thermal response, and phase change discharge response of the current condensation reaction; match the current reaction stage; allocate the current reaction input according to the acceptance ratio; and form the regional input response of each detection area. Step 3: Develop a regional reference thermal response based on the regional input response and the regional thermal response benchmark. When the regional reference thermal response is higher than the current regional thermal response, develop a local thermal response missing amount based on the difference between the regional reference thermal response and the current regional thermal response. Step 4: Based on the portion of the current phase change discharge response that exceeds the phase change discharge response benchmark, a phase change discharge deviation is formed. The phase change discharge deviation is matched with the local thermal response missing amount according to the thermal response transmission order to form the phase change acceptance state of the corresponding detection area. Step 5: When the regional baseline thermal response does not decrease, the phase change discharge deviation decreases during the continuous detection period, and the current regional thermal response increases during the corresponding detection period, a thermal response switching state is formed, and the safety risk level is determined according to the phase change acceptance state, the thermal response switching state, and the progressive changes in the current regional thermal response.

2. The safety risk classification and early warning method for condensation reaction process production according to claim 1, characterized in that, Step one, the method for forming the detection area and the acceptance ratio, includes: A detection position response sequence is formed according to the occurrence time of the thermal response at each detection position after the same reaction input; Based on the response interval, response direction relationship, and response change ratio of adjacent detection positions in the detection position response sequence, the continuous connection result between adjacent detection positions is determined; Adjacent detection positions whose continuous acceptance results meet the acceptance benchmark of the corresponding reaction stage are grouped into the same continuous acceptance position group, and adjacent detection positions whose continuous acceptance results do not meet the acceptance benchmark of the corresponding reaction stage are used as the boundary between different continuous acceptance position groups. Each group of consecutive receiving positions is grouped into a detection area, and the receiving ratio of the corresponding detection area to the reaction input is formed according to the proportion of the response change of each detection position relative to the same reaction input within the same group of consecutive receiving positions.

3. The safety risk classification and early warning method for condensation reaction process production according to claim 1, characterized in that, Step one, the method for forming the regional thermal response benchmark and the phase change discharge response benchmark, includes: Read the reaction input, regional thermal response, and phase change discharge response from the normal production data of the isomorphic condensation reaction corresponding to each reaction stage; Based on the proportion of the reaction input received by the detection area, the reaction input is converted into the regional input response of each detection area; Candidate benchmarks for regional thermal response are formed based on the correspondence between the changes in regional input response and regional thermal response, and candidate benchmarks for phase change discharge response are formed based on the correspondence between the changes in regional input response and phase change discharge response. The production section that is a candidate phase change receiving section is defined as the production section whose regional thermal response decreases relative to the regional thermal response candidate benchmark, whose phase change discharge response increases relative to the phase change discharge response candidate benchmark, and whose period of decrease in regional thermal response and period of increase in phase change discharge response satisfy the time correspondence benchmark of the corresponding reaction stage. Candidate phase change receiving sections are eliminated, and a regional thermal response benchmark is formed based on the correspondence between the changes in regional input response and regional thermal response within the eliminated production sections. A phase change discharge response benchmark is also formed based on the correspondence between the changes in regional input response and phase change discharge response within the eliminated production sections.

4. The safety risk classification and early warning method for condensation reaction process production according to claim 3, characterized in that, In step two, the matching method for the current reaction stage includes: Based on the moment of change in the reaction input of the current condensation reaction, extract the current region's thermal response and the current phase transition discharge response; The current regional thermal response is time-series registered with the regional thermal response benchmarks corresponding to each reaction stage to form regional response registration results. The current phase change discharge response is time-series registered with the phase change discharge response reference corresponding to each reaction stage to form the phase change discharge registration result. When the regional response registration result and the phase change discharge registration result correspond to the same reaction stage, the corresponding reaction stage is determined as the current reaction stage; When the regional response registration result and the phase change discharge registration result correspond to different reaction stages, the current reaction stage is determined based on the period of decrease in the current regional thermal response relative to the regional thermal response benchmark, the period of increase in the current phase change discharge response relative to the phase change discharge response benchmark, and the time sequence corresponding benchmark of the corresponding reaction stage.

5. The safety risk classification and early warning method for condensation reaction process production according to claim 1, characterized in that, Step two, the method for forming the regional input response, includes: Based on the acceptance ratio corresponding to the current reaction stage, the current reaction input is allocated to each detection area to form the candidate area input response for each detection area; According to the thermal response transfer sequence corresponding to the current reaction stage, determine the regional transition period of each detection area relative to the current reaction input; Based on the input response of the candidate region and the regional thermal response benchmark, the regional thermal response that each detection region should generate during the regional acceptance period is determined. The current regional thermal response during the regional takeover period is compared with the expected regional thermal response to form the deviation of the current regional thermal response from the expected regional thermal response. The input response of the candidate region is corrected based on the response deviation to form the regional input response of each detection region.

6. The safety risk classification and early warning method for condensation reaction process production according to claim 5, characterized in that, Step three, the method for forming the local thermal response deficiency includes: Based on the regional input response and the regional thermal response benchmark, the regional benchmark thermal response of each detection area is formed; The regional baseline thermal response and the current regional thermal response are matched according to the regional transition period. When the regional baseline thermal response is higher than the current regional thermal response, the candidate thermal response missing amount for the corresponding detection area is formed based on the difference between the regional baseline thermal response and the current regional thermal response. The candidate thermal response missing amount in the same detection area is read according to adjacent detection time periods, and the candidate thermal response missing amount in the next detection time period is retained when the candidate thermal response missing amount in the next detection time period is not less than the difference between the candidate thermal response missing amount in the previous detection time period and the thermal response comparison resolution. According to the thermal response transmission order, the thermal responses of adjacent detection areas are correlated with the missing amounts of retained candidate thermal responses in different time periods, and the thermal responses of adjacent detection areas and the missing amounts of retained candidate thermal responses are compared to see if they meet the transition criteria of the corresponding reaction stages. When the thermal response of adjacent detection areas and the amount of missing candidate thermal response do not meet the transition criteria for the corresponding reaction stage, the amount of missing candidate thermal response is determined as the amount of missing local thermal response in the corresponding detection area.

7. The safety risk classification and early warning method for condensation reaction process production according to claim 6, characterized in that, Step four, the matching method between the phase change discharge deviation and the local thermal response missing amount, includes: The phase change discharge deviation is calculated based on the portion of the current phase change discharge response that exceeds the phase change discharge response baseline. Read the formation time and duration of the local thermal response missing amount in each detection area, and according to the start offset and end offset in the time-corresponding reference of the corresponding reaction stage, shift the formation time and duration of the local thermal response missing amount to form the phase change discharge reading start time and phase change discharge reading end time. The phase change discharge deviation is read from the start time to the end time of phase change discharge reading, and the read phase change discharge deviation is matched with the local thermal response loss in the corresponding detection area according to the detection time period; According to the correspondence between the changes in regional thermal response and phase change discharge response, the read phase change discharge deviation is converted into a phase change discharge conversion quantity with the same dimensions as the missing amount in local thermal response. The phase change discharge conversion comparison resolution is formed based on the phase change discharge comparison resolution and the correspondence between the changes in regional thermal response and phase change discharge response. Calculate the deviation difference between the phase change discharge conversion amount and the local thermal response missing amount; When the deviation difference is not greater than the phase change discharge conversion comparison resolution, the correspondence between the corresponding detection area, the local thermal response missing amount and the phase change discharge deviation amount is retained. When there is only one retained correspondence for the same phase change discharge deviation, the retained correspondence is determined as the phase change acceptance status of the corresponding detection area; When there are more than two retained correspondences for the same phase change discharge deviation, calculate the overlap length between the reading time period of the phase change discharge deviation and the phase change discharge reading start time to phase change discharge reading end time for each detection area, and retain the correspondence with the largest overlap length. When there is only one correspondence with the largest overlap length, the correspondence with the largest overlap length is determined as the phase transition acceptance state of the corresponding detection area; When there are two or more correspondences with the same maximum overlap length, the correspondence with the smallest deviation difference is determined as the phase transition acceptance state of the corresponding detection area.

8. The safety risk classification and early warning method for condensation reaction process production according to claim 7, characterized in that, Step five, the method for forming the thermal response switching state and safety risk level, includes: Read the regional reference thermal response, phase change discharge deviation, and current regional thermal response corresponding to the detection area that forms the phase change receiving state; When the regional reference thermal response of the next detection period is not less than the difference between the regional reference thermal response and the thermal response comparison resolution of the previous detection period, and the phase change discharge deviation of the next detection period is less than the difference between the phase change discharge deviation and the phase change discharge comparison resolution of the previous detection period, the next detection period is regarded as the phase change transition weakening period. Based on the time-series correspondence of the corresponding reaction stages, the phase transition weakening period is mapped to the detection period of the current region's thermal response; When the thermal response of the current area in the next detection period is greater than the sum of the thermal response of the current area in the previous detection period and the thermal response comparison resolution, a thermal response switching state is formed. Based on the formation of phase change acceptance state, phase change acceptance weakening period and thermal response switching state, the safety risk level of the corresponding detection area is determined, and the progressive change of the current area's thermal response after the formation of thermal response switching state is used as the basis for the progression of safety risk level.

Images