A collaborative scheduling method based on non-ferrous metal smelting working condition data index driving

Abstract

The present application relates to the technical field of metallurgical automation, and particularly relates to a collaborative scheduling method based on non-ferrous metal smelting working condition data index driving, comprising the following steps: determining the source working condition event of downstream working condition change through the consistency of working procedure material transfer relationship and multiple target element change, so as to reduce the data mismatch of adjacent batches and cross-process data; establishing the working condition propagation phase index through the arrival order of material, heat and flue gas response, determining the change node and working condition change range by comparing the continuous working condition propagation phase index, and combining the current working condition and operation boundary of the equipment to gradually cut the scheduling branch which does not need to continue to propagate, so as to reduce the collaborative scheduling range; determining the scheduling quantity according to the working condition quantity exceeding the remaining adjustment range and the equipment adjustment response relationship, so that the equipment adjustment is matched with the actual working condition change, thereby reducing the cross-process response lag and excessive scheduling, and improving the scheduling accuracy, response timeliness and equipment operation stability of the smelting continuous production process.

Description

Technical Field

[0001] This invention relates to the technical field of metallurgical automation, and in particular to a collaborative scheduling method driven by non-ferrous metal smelting working condition data index. Background Technology

[0002] In a continuous production line for copper concentrate flash smelting—matte blowing—waste heat recovery—flue gas purification—acid production, copper concentrate enters the flash smelting furnace after batching and feeding. The resulting matte continues into the blowing process. The high-temperature flue gas generated during smelting and blowing passes through the waste heat boiler and dust collection equipment before entering the acid production system. Different batches of concentrate have differences in the composition and processing volume of copper, iron, and sulfur. Changes in operating conditions caused by batch switching are transmitted to subsequent processes through material migration, changes in furnace thermal state, and flue gas transport. Since matte material residence, furnace thermal response, and flue gas pipeline transport have different response times, the timing of responses to changes in the same batch at different process nodes is not consistent. Existing operating condition management methods usually associate different process data according to a unified timestamp, furnace number, or fixed time window, which easily leads to incorrect merging of upstream and downstream operating condition changes caused by different batch events, making it difficult to accurately determine the upstream source operating condition events that actually correspond to changes in downstream elemental composition, heat load, or flue gas state.

[0003] Existing smelting production scheduling typically involves overall optimization or periodic recalculation of smelting furnaces, blowing furnaces, waste heat equipment, and flue gas treatment equipment based on current furnace conditions, equipment load, and product plans. When changes in concentrate batches only propagate to some process nodes, this scheduling method struggles to distinguish between equipment not yet affected by that batch and equipment already in the response phase. Furthermore, different equipment have varying remaining space from their operating boundaries, meaning that changes in the same operating condition may not necessarily trigger adjustments in subsequent equipment if upstream equipment can still handle them. Including all related equipment in the scheduling simultaneously can easily expand the solution scope and cause unnecessary changes in control variables. Adjustments based solely on single parameter exceeding limits may lead to delayed adjustment timing due to lag in heat and flue gas propagation. Moreover, existing methods rarely consider the actual arrival locations of operating condition changes along material, heat, and flue gas propagation paths to propagate and trim the impact range node by node, making it difficult to form a localized coordinated scheduling range that matches the current propagation stage.

[0004] Given the problems of unclear attribution of cross-process operating condition events, difficulty in distinguishing different physical response propagation stages, and over-scheduling caused by operating condition changes in existing non-ferrous metal smelting operating condition scheduling technologies, this invention proposes a collaborative scheduling method driven by non-ferrous metal smelting operating condition data index. This method determines the source operating condition events of downstream operating condition changes through process material transfer relationships and changes in multiple target elements. It establishes an operating condition propagation phase index based on the arrival order of material response, heat response, and flue gas response, and determines the operating condition change range through changes between consecutive operating condition propagation phase indices. The operating condition change range propagates to subsequent equipment along process constraints, and terminates propagation branches that do not require further scheduling by combining the remaining adjustment range formed by the current operating condition and operating boundary of the equipment. Scheduling amounts are determined only for equipment corresponding to branches that have not terminated, thus solving the problems of cross-process response mismatch and difficulty in accurately determining local scheduling ranges under conditions of batch changes in concentrate batching. Summary of the Invention

[0005] The purpose of this section is to outline some aspects of the embodiments of the present invention and to briefly introduce some preferred embodiments. Some simplifications or omissions may be made in this section, as well as in the abstract and title of the present application, to avoid obscuring the purpose of this section, the abstract and title of the invention. Such simplifications or omissions shall not be used to limit the scope of the present invention.

[0006] In view of the aforementioned existing problems, the present invention is proposed.

[0007] To solve the above-mentioned technical problems, the present invention provides the following technical solution: In a first aspect, the present invention provides a collaborative scheduling method driven by non-ferrous metal smelting working condition data index, comprising: limiting the range of source working condition events according to the material transfer relationship of the process, and determining the source working condition events of downstream working condition changes according to the consistency of changes of multiple target elements; Based on the arrival order of the material, heat, and flue gas responses of the source operating condition event, the arrival node and subsequent node are determined to form an operating condition propagation phase index. Compare the propagation phase index of the current operating condition with the propagation phase index of the previous operating condition to determine the change node and the range of operating condition change, and propagate along the process constraint relationship; determine the remaining adjustment range according to the current operating condition and the operating boundary of the equipment, and terminate the branch when the range of operating condition change does not exceed the remaining adjustment range. Using the equipment corresponding to the unterminated branch as the scheduling object, the scheduling amount is determined based on the portion of the operating condition change range that exceeds the remaining adjustment range and the equipment adjustment response relationship, and the scheduling object is adjusted accordingly.

[0008] In a second aspect, the present invention provides a computer device, comprising: One or more processors; The memory stores operable instructions that, when executed by the one or more processors, cause the one or more processors to perform operations, including the flow of the cooperative scheduling method driven by the non-ferrous metal smelting working condition data index as described above.

[0009] Thirdly, the present invention provides a computer-readable medium for storing software, the software including instructions executable by one or more computers, the instructions causing the one or more computers to perform operations, the operations including the flow of the aforementioned cooperative scheduling method driven by non-ferrous metal smelting working condition data index.

[0010] The beneficial effects of this invention are as follows: This invention determines the source operating condition events of downstream operating condition changes by identifying the consistency of material transfer relationships and changes in multiple target elements, thereby reducing data mismatches between adjacent batches and across processes; it establishes an operating condition propagation phase index by establishing the arrival order of material, heat, and flue gas responses, clarifying the current arrival position of the same operating condition event in different propagation paths; it determines the change nodes and operating condition change range by comparing continuous operating condition propagation phase indices, and, in conjunction with the current operating condition and operating boundary of the equipment, it gradually trims scheduling branches that do not need to continue propagating, narrowing the scope of collaborative scheduling; furthermore, it determines the scheduling quantity based on the operating condition quantity exceeding the remaining adjustment range and the equipment adjustment response relationship, so that the equipment adjustment matches the actual operating condition changes, thereby reducing cross-process response lag and over-scheduling, and improving the scheduling accuracy, response timeliness, and equipment operation stability of the continuous smelting production process. Attached Figure Description

[0011] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. Wherein: Figure 1 This is a flowchart illustrating the collaborative scheduling method driven by non-ferrous metal smelting working condition data index as shown in this invention. Detailed Implementation

[0012] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments.

[0013] Based on the embodiments of this invention, all other embodiments obtained by those skilled in the art without inventive effort should fall within the scope of protection of this invention.

[0014] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0015] This embodiment uses a continuous production line for copper concentrate flash smelting—matte blowing—waste heat recovery—flue gas purification—acid production as a specific application scenario. The production line includes a copper concentrate batching and feeding system, a flash smelting furnace, a matte transfer device, a blowing furnace, a waste heat boiler, a dust collection device, and an acid production system. After batching and feeding, the copper concentrate enters the flash smelting furnace, and the matte formed by smelting enters the blowing furnace via the matte transfer device. The high-temperature flue gas generated by flash smelting passes through the waste heat boiler and dust collection device in sequence and enters the acid production system. After the batch of concentrate is switched, the composition of copper, iron, and sulfur elements and the material processing volume change accordingly. These changes are propagated to subsequent process nodes through matte material migration, heat transfer, and flue gas transportation.

[0016] In this embodiment, the copper concentrate batching and feeding system, flash smelting furnace, matte transfer device, blowing furnace, waste heat boiler, dust collection equipment, and acid production system adopt the unified clock of the production control system; as a preferred implementation, the operating conditions involved in the calculation of this method are arranged according to a unified sampling interval of 1 second; 1 second is an illustrative value, and the actual production line can select other sampling intervals according to the sampling capacity of the measuring device and the response speed of the equipment; the material processing volume is obtained by the corresponding metering device, the mass fraction of the target element is obtained by the corresponding material analysis results, and the current operating condition of the equipment is obtained by the production control system.

[0017] In this embodiment, the equipment operation boundary and the equipment control boundary correspond to different technical objects. The equipment operation boundary is the upper and lower bounds of the corresponding operating condition quantity that is allowed to operate continuously under the current production process conditions. The equipment control boundary is the upper and lower bounds of the feed rate, fan speed, valve opening, or other execution control quantity that is allowed to be adjusted. The remaining adjustment range represents the range of operating condition quantities between the current operating condition distance and the equipment operation boundary corresponding to the direction of operating condition change.

[0018] The production system also stores node propagation delay benchmarks obtained from the calibration of similar historical source operating conditions. The node propagation delay benchmarks correspond to the response type, propagation path, and propagation nodes, reflecting the existing propagation delay of changes in similar concentrate batches from the source operating condition event occurrence node to the corresponding propagation node. The node propagation delay benchmarks serve as the initial basis for the arrival time of propagation nodes that have not yet formed an actual response to the current source operating condition event. Once the current source operating condition event has formed an actual response at the corresponding propagation node, time-shift cross-correlation calculations are performed based on the actual response of the current source operating condition event, and the resulting propagation delay updates the arrival time of the corresponding propagation node.

[0019] Under the above conditions of equipment connection relationship, working condition source, unified time reference and node propagation delay reference, the working condition event formed by the batch change of concentrate batching is taken as the processing object. Based on the material batch connection and the change of multiple target elements, the source working condition event that causes the downstream working condition change is determined, and the following processing is performed in sequence.

[0020] According to an embodiment of the present invention, in combination Figure 1 The flowchart shown illustrates a collaborative scheduling method driven by non-ferrous metal smelting operating condition data indexing, which specifically includes the following steps: S1. Define the scope of the source operating condition event based on the material transfer relationship of the process, and determine the source operating condition event of the downstream operating condition change based on the consistency of changes in multiple target elements. Note that the following should be noted in this step: As a preferred implementation, this embodiment uses the batch switching of concentrate feed as a candidate source operating condition event, takes the change in matte composition generated at the outlet of the flash smelting furnace as a specific example of downstream operating condition change, and selects copper, iron and sulfur as target elements; copper reflects the change of target metal entering the matte, iron has a material relationship with smelting slag formation and changes in matte composition, and sulfur is related to both matte composition and the formation process of smelting flue gas.

[0021] S1.1. Based on the material conveying path and intermediate material batch connection relationship between adjacent smelting processes, trace back the upstream operations that have material connection relationship with the downstream operating condition changes to determine the candidate source operating condition events.

[0022] The moment when a new batch of concentrate feed begins to enter the flash smelting furnace via the feeding system is taken as the starting point of the candidate source condition event; for changes in the target element composition at the matte outlet, the upstream concentrate feed batch that may correspond to the current matte material range is determined based on the concentrate feeding sequence, the effective material retention in the flash smelting furnace, the material handling rate range, and the matte discharge sequence.

[0023] In a preferred embodiment, when the effective material retention in the flash smelting furnace is H and the material processing rate is between low and high, the lower limit of the material residence time is obtained by dividing the effective material retention by the high value of the material processing rate, and the upper limit of the material residence time is obtained by dividing the effective material retention by the low value of the material processing rate. The material residence time range is shifted forward from the moment of change in matte composition, and the concentrate batches that have a material connection relationship with the current matte material range are determined in combination with the matte discharge sequence. The changes in the corresponding batches are determined as candidate source operating conditions.

[0024] S1.2 Based on the target element composition of the material corresponding to the candidate source working condition event and the material allocation relationship between processes, calculate the direction of increase or decrease and the relative change of the target elements after the candidate source working condition event is transmitted to the downstream process, and obtain the combination of changes in the candidate target elements.

[0025] For a target element, within the same material statistics interval, the product of the material processing volume after the occurrence of the candidate source condition event and the mass fraction of the target element is subtracted from the product of the material processing volume before the occurrence of the candidate source condition event and the mass fraction of the target element. The difference is determined as the source-end change of the target element.

[0026] In this embodiment, 60 seconds is used as a set of illustrative statistical intervals. For example, the material processing volume in the 60 seconds before the switching of candidate concentrate batches is 2000 kg, and the material processing volume in the 60 seconds after the switching is 2100 kg. The copper mass fraction changes from 0.280 to 0.300, the iron mass fraction changes from 0.280 to 0.270, and the sulfur mass fraction changes from 0.300 to 0.310. According to the above relationships, the change at the copper source end is 70 kg, the change at the iron source end is 7 kg, and the change at the sulfur source end is 51 kg.

[0027] Based on the material diversion ratio in the material transfer relationship and the retention ratio of the target element in the intermediate processing stage, the change at the source end is converted step by step along the material transfer path from the candidate source operating event to the matte outlet; the material diversion ratio is determined according to the ratio of the amount of material entering the corresponding material branch within the same statistical interval to the amount of material before diversion; the target element retention ratio is determined according to the ratio of the amount of target element in the outlet material of the corresponding intermediate processing stage to the amount of target element in the inlet material.

[0028] When there are multiple effective material transfer paths, the change in the target element is calculated for each material transfer path according to its complete path material diversion ratio and the comprehensive retention ratio of the target element. Then, the change in the target element of each material transfer path is added together to obtain the net change in the target element of the downstream process. For multiple material transfer paths with a shared front end, the converted complete path effective diversion ratio is used to avoid double counting of shared material quantities.

[0029] For example, if the effective material diversion ratio from concentrate to matte outlet is 0.95, and the comprehensive retention ratios of copper, iron, and sulfur entering matte are 0.95, 0.85, and 0.75, respectively, then the net change in the target elements transferred to matte outlet by copper, iron, and sulfur are approximately 63.18 kg, 5.65 kg, and 36.34 kg, respectively.

[0030] The zero-value judgment threshold is determined based on the measurement error of material processing volume and the detection error of target element mass fraction. The allowable error for the product of material processing volume and target element mass fraction is determined based on the measurement error of material processing volume and the detection error of target element mass fraction. The allowable errors obtained before and after the occurrence of the candidate source working condition event are combined, and then converted downstream according to the same material diversion ratio and target element retention ratio to obtain the zero-value judgment threshold corresponding to the net change of target element.

[0031] When the net change of the target element is greater than the zero value judgment threshold, the direction of increase or decrease of the target element is determined to be increase; when the net change of the target element is less than the negative value of the zero value judgment threshold, the direction of increase or decrease of the target element is determined to be decrease; when the net change of the target element is between the positive and negative zero value judgment thresholds, the direction of increase or decrease of the target element is determined to be basically unchanged.

[0032] The relative change range of the target element is determined by the ratio of the absolute value of the net change of the target element to the reference element quantity of the target element. The reference element quantity of the target element is obtained by multiplying the downstream receiving material processing volume of the same statistical interval as the candidate source operating condition event before the occurrence of the event and the mass fraction of the same target element. When the reference element quantity of the target element is lower than the lower limit of the reference element quantity, the lower limit of the reference element quantity is used in the calculation of the relative change range. The lower limit of the reference element quantity can be determined by multiplying the minimum effective material processing volume with the corresponding quantitative detection lower limit of the target element.

[0033] For example, if the amount of matte material discharged within 60 seconds before the occurrence of a candidate source event is 800 kg, with mass fractions of copper, iron, and sulfur of 0.600, 0.180, and 0.200, respectively, then the target element baseline amounts for copper, iron, and sulfur are 480 kg, 144 kg, and 160 kg, respectively. Based on the aforementioned net change in target elements, the relative change ranges for copper, iron, and sulfur are approximately 0.132, 0.039, and 0.227, respectively. Candidate target element change combinations are formed based on the direction of increase or decrease in target elements and their relative change ranges.

[0034] S1.3. Compare the increase / decrease direction of the same target element in the downstream operating condition changes with the combination of changes of candidate target elements element by element, and eliminate candidate source operating condition events where the increase / decrease direction of target elements is inconsistent.

[0035] The actual test results of matte exports are calculated using the same material statistical range and relative change magnitude as the candidate target element change combination. For example, if copper, iron, and sulfur all increase in the actual downstream operating conditions, with relative change magnitudes of 0.129, 0.041, and 0.222 respectively, then these are compared with the increase / decrease directions of copper, iron, and sulfur for each candidate source operating condition event. If the increase / decrease direction of any target element is inconsistent, the corresponding candidate source operating condition event is removed. If the increase / decrease directions of all three target elements are consistent, the corresponding candidate source operating condition event is retained.

[0036] S1.4 Calculate the absolute deviation between the relative change amplitude of each target element in the remaining candidate source operating condition events and the relative change amplitude of the same target element in the downstream operating condition changes. The candidate source operating condition event with the smallest sum of absolute deviations and all absolute deviations within the corresponding allowable deviation range is determined as the source operating condition event.

[0037] The allowable deviation range is determined after propagating the results of the target element mass fraction detection error, material processing volume measurement error, material diversion ratio error, and target element retention ratio error to the relative change amplitude. The zero value judgment threshold determines whether there is an effective increase or decrease in a single target element, and the allowable deviation range determines whether the candidate target element change combination is consistent with the downstream actual target element change. The two play different judgment roles.

[0038] Taking the aforementioned example, the candidate relative variation ranges for copper, iron, and sulfur are 0.132, 0.039, and 0.227, respectively, while the actual relative variation ranges downstream are 0.129, 0.041, and 0.222, respectively. The corresponding absolute deviations are approximately 0.003, 0.002, and 0.005, respectively, and the sum of the absolute deviations is approximately 0.010.

[0039] When multiple candidate source conditions all meet the allowable deviation range, the one with the smallest sum of absolute deviations is determined as the source condition event; when the sum of absolute deviations is the same, as a preferred implementation, the candidate source condition event that is directly continuous with the current copper matte material interval in the material batch connection relationship is determined as the source condition event; when there is no candidate source condition event that meets the allowable deviation range of all target elements, the current downstream condition change does not enter the propagation phase index processing of the corresponding source condition event.

[0040] It should be noted that this step restricts candidate upstream batches by actual material transfer relationships, and then uses the change direction and relative change magnitude of multiple target elements after smelting to determine the event attribution, so that the subsequent propagation status judgment revolves around the changes of the same concentrate batch.

[0041] S2. Based on the arrival order of materials, heat, and flue gas responses of the source operating condition event, determine the arrival nodes and subsequent nodes to form an operating condition propagation phase index. It should be noted that in this step:

[0042] In this embodiment, the material response uses the sulfur mass flow rate as an example. The sulfur mass flow rate is obtained by multiplying the material mass flow rate of the corresponding node by the sulfur mass fraction. The heat response uses the flue gas sensible heat flow rate, which is determined based on the flue gas mass flow rate and the enthalpy difference of the flue gas relative to a unified reference state. The flue gas response uses the standard state flue gas volumetric flow rate, and the actual flue gas volumetric flow rate is converted to a unified standard state based on the absolute temperature and absolute pressure of the corresponding node.

[0043] S2.1. Based on the material conveying connection, heat exchange connection and flue gas pipeline connection between the source working condition event occurrence node and the downstream process, determine the material response propagation path, heat response propagation path and flue gas response propagation path respectively, and arrange the propagation nodes according to the process connection sequence.

[0044] In this embodiment, the material response propagation path sequentially includes the concentrate feeding node, the flash smelting furnace matte outlet node, the matte transfer node, and the blowing furnace inlet node; the heat response propagation path sequentially includes the flash smelting furnace flue gas outlet node, the waste heat boiler inlet node, and the waste heat boiler outlet node; the flue gas response propagation path sequentially includes the flash smelting furnace flue gas outlet node, the waste heat boiler flue gas outlet node, the dust collection equipment flue gas outlet node, and the acid production system inlet node.

[0045] When the same device is located in different response type propagation paths, its propagation status is determined by the corresponding operating conditions of each response type.

[0046] S2.2. Based on the material residence time range, thermal response time range, and flue gas transport time range between adjacent propagation nodes, the candidate arrival time periods for the source operating condition event to reach each propagation node are determined by accumulating the data along the material response propagation path, thermal response propagation path, and flue gas response propagation path.

[0047] For material response, the range of material residence time between adjacent propagation nodes is determined based on the effective amount of material retained by the equipment and the actual material processing rate range; the lower limit of material residence time is obtained by dividing the effective amount of material retained by the higher value of the material processing rate, and the upper limit of material residence time is obtained by dividing the effective amount of material retained by the lower value of the material processing rate.

[0048] For example, the candidate arrival time from the concentrate feeding node to the flash smelting furnace matte outlet node is 360s to 480s after the source event. The time from the matte outlet node to the matte transfer node is increased by 100s to 180s, and the time from the matte transfer node to the blowing furnace inlet node is increased by another 240s to 360s. Therefore, the cumulative candidate arrival time relative to the source event is 360s to 480s, 460s to 660s, and 700s to 1020s, respectively.

[0049] For flue gas response, the flue gas transport time between adjacent propagation nodes is determined based on the effective volume of the corresponding flue gas pipe section and the actual flue gas volumetric flow rate under the operating conditions of the pipe section; the standard state flue gas volumetric flow rate is converted into the actual flue gas volumetric flow rate under the operating conditions based on the absolute pressure and absolute temperature of the corresponding pipe section and then used in the flue gas transport time calculation.

[0050] For heat response, the heat response time range of adjacent thermal nodes can be determined based on existing operating results from the generation of identifiable heat load changes at the source thermal node to the occurrence of heat load changes in the same direction at the subsequent thermal node under stable production conditions.

[0051] For propagation nodes in the production system that have been calibrated with a propagation delay reference, the occurrence time of the source condition event is added to the corresponding node propagation delay reference to obtain the initial value of the arrival time of the propagation node. At the same time, the effective propagation range of the initial value of the arrival time is limited by the candidate arrival time period formed based on the material residence time range, thermal response time range, or flue gas transport time range. Thus, when the current source condition event has not yet formed an actual response in the subsequent node, the node propagation delay reference provides the time location basis for the subsequent node.

[0052] S2.3 During the candidate arrival period, perform time-shift cross-correlation calculations on the response changes at the source condition event occurrence node and the condition changes of the same response type at the propagation node, and determine the arrival time of the propagation node by the time shift when the cross-correlation coefficient is maximized.

[0053] Once the current source condition event has generated an actual response at the corresponding propagation node, the same response type condition quantities covering the candidate arrival time period are extracted from the source condition event occurrence node and the propagation node, respectively. The average value of the stable interval before the occurrence of the source condition event is subtracted from the two sets of condition quantities, and the amplitude is normalized according to the standard deviation of the corresponding data segment. Then, the normalized cross-correlation coefficient is calculated only within the time shift range corresponding to the candidate arrival time period.

[0054] The normalized cross-correlation coefficient is determined by the ratio between the sum of the product of the source node response change and the time-shifted propagation node response change and the square root of the product of the sum of the squares of the two sets of average differences; the time shift when the cross-correlation coefficient is at its maximum during the candidate arrival period is determined as the propagation delay of the current source condition event to the propagation node, and the sum of the occurrence time of the source condition event and the propagation delay is determined as the confirmed arrival time.

[0055] When a propagation node has not yet formed an actual response to the current source event, the propagation node retains the initial value of the arrival time based on the node propagation delay reference. After the actual response to the current source event is formed and the confirmed arrival time is obtained, the original initial value of the arrival time is updated with the confirmed arrival time. Thus, the same propagation node has a reference arrival time before the current event response occurs, and is corrected by the actual cross-correlation result after the current event response occurs.

[0056] When two candidate time-shift positions achieve the same maximum cross-correlation coefficient, as a preferred implementation, the earlier time-shift position within the candidate arrival time period is selected.

[0057] S2.4. Along the corresponding propagation path, the last propagation node whose arrival time is no later than the current time is determined as the arrival node, and the propagation node that is adjacent to the arrival node and whose arrival time is later than the current time is determined as the successor node.

[0058] In this embodiment, the initial value of the arrival time is used when the propagation node has not yet obtained the confirmation arrival time, and the confirmation arrival time is used when the propagation node has obtained the confirmation arrival time.

[0059] For example, if the source event occurs at time t0, the arrival times currently saved for the matte outlet node, matte transfer node, and smelting furnace inlet node are t0+430s, t0+560s, and t0+840s, respectively. When the current time is after t0+560s but before t0+840s, the matte transfer node is determined as the arrival node, and the smelting furnace inlet node is determined as the successor node.

[0060] After propagating to the end of the path, the end propagation node is determined as the arrival node, and no new successor node is set.

[0061] S2.5. Classify response types according to material response, heat response, and flue gas response, and arrange the propagation nodes of the same response type according to the process connection sequence of the corresponding propagation path.

[0062] Preferably, the sequence of propagation nodes is determined according to the actual material conveying connection, heat exchange connection, or flue gas pipeline connection.

[0063] S2.6. Configure the arrival time of each propagation node to the corresponding propagation node. Propagation nodes whose arrival time is no later than the current time are determined to be in the arrived state, and propagation nodes whose arrival time is later than the current time are determined to be in the unarrived state.

[0064] When the propagation node has not yet formed an actual response, the currently saved arrival time is the initial value of the arrival time formed by the node propagation delay reference; after the actual response is formed and the time shift cross-correlation calculation is completed, the currently saved arrival time is updated to the confirmed arrival time.

[0065] When the state of any propagation node changes within the same source operating condition event and the same response type, a new operating condition propagation phase index is formed; the most recently saved operating condition propagation phase index of the same source operating condition event and the same response type before the current operating condition propagation phase index is formed is determined as the previous operating condition propagation phase index.

[0066] S2.7 Determine the current propagation position of the source event under the corresponding response type based on the arriving node, the successor node, and the propagation segment between the two nodes.

[0067] For example, if the material response has reached the matte transfer node, and the arrival time of the furnace inlet node is later than the current time, the propagation segment between the matte transfer node and the furnace inlet node is determined as the current propagation position of the material response.

[0068] S2.8 Establish an index hierarchy according to the source operating condition event, response type and propagation node order, and configure the arrival time and current propagation position under the corresponding propagation node to construct the operating condition propagation phase index.

[0069] The operating condition propagation phase index must at least maintain the correspondence between the source operating condition event, response type, propagation node sequence, arrival time, and current propagation position; establish a time positioning relationship between the propagation node and the formation time of the corresponding index state and the actual operating condition quantity in the production control system; when obtaining the operating condition quantity under the corresponding index state, obtain the actual operating condition quantity based on the propagation node, response type, and index state formation time.

[0070] It should be noted that this step organizes the material, heat, and flue gas responses of the same source operating condition event into the same operating condition propagation phase index according to the actual propagation path and propagation time. The historical node propagation delay benchmark provides the time position before the current actual response node is generated, and the current event cross-correlation result corrects the arrival time of the actual response node, thereby forming a continuous and updatable propagation position relationship.

[0071] S3. Compare the propagation phase index of the operating condition with the propagation phase index of the previous operating condition to determine the changing node and the range of operating condition changes, and propagate along the process constraint relationship; determine the remaining adjustment range according to the current operating condition and operating boundary of the equipment, and terminate the branch when the operating condition change range does not exceed the remaining adjustment range. Note that the following points should be noted in this step: S3.1. Compare the propagation phase index of the operating condition with the propagation phase index of the previous operating condition, and determine the propagation node status of the same source operating condition event and the same response type. The propagation node that changes from the unarrived state to the arrived state is identified as the changed node.

[0072] For example, if the copper matte transfer node is in an unreached state in the previous propagation phase index, and the copper matte transfer node is in a reached state in the current propagation phase index, then the copper matte transfer node is determined as the change node of the material response.

[0073] S3.2. Compare the operating condition quantity corresponding to the response type of the changing node at the current arrival time with the operating condition quantity corresponding to the response type of the changing node in the previous operating condition propagation phase index, and determine the operating condition change range of the changing node based on the difference between the two and the measurement error of the operating condition quantity.

[0074] When the current operating condition quantity measurement value is the current value and the upper limit of the absolute value of the measurement error is the current error, the current operating condition quantity range is composed of the current value minus the current error to the current value plus the current error; the previous operating condition quantity range is formed in the same way; after subtracting the previous operating condition quantity range from the current operating condition quantity range, the operating condition change range is obtained.

[0075] Specifically, the lower limit of the operating condition variation range is equal to the current measured value minus the previous measured value, and then minus the sum of the current error and the previous error; the upper limit of the operating condition variation range is equal to the current measured value minus the previous measured value, plus the sum of the current error and the previous error.

[0076] For example, the mass flow rate of sulfur at the copper matte transfer node was 2.30 kg / s in the previous index state and 2.85 kg / s in the current state. The upper limit of the absolute value of the error between the two measurements is 0.025 kg / s. Therefore, the operating condition variation range is 0.50 kg / s to 0.60 kg / s, and the direction of increase is increasing.

[0077] S3.3 Select the material balance relationship, heat transfer relationship or flue gas transport relationship according to the response type of the changing node, convert the operating condition change range into the operating condition change range of the subsequent node, and propagate it level by level along the corresponding propagation path.

[0078] For material response, the range of change in the mass flow rate of the target element at the current propagation node is multiplied by the material diversion ratio and the target element retention ratio from that node to the subsequent node, respectively, to obtain the range of change in the mass flow rate of the target element at the subsequent node.

[0079] For example, if the product of the material diversion ratio from the matte transfer node to the inlet node of the smelting furnace and the sulfur retention ratio is 0.90, then the range of sulfur mass flow rate variation from 0.50 kg / s to 0.60 kg / s is converted to 0.45 kg / s to 0.54 kg / s.

[0080] For heat response, based on the heat balance relationship between adjacent thermal nodes, the range of changes in the sensible heat flow rate of the upstream flue gas is subtracted from the range of changes in the heat absorption and heat loss of adjacent equipment to obtain the range of changes in the sensible heat flow rate of the flue gas at subsequent nodes; when the stable operating range has an effective heat transfer ratio, the range conversion can also be completed based on the effective heat transfer ratio.

[0081] For flue gas response, based on the standard state flue gas volume balance and flue gas diversion ratio between adjacent propagation nodes, the standard state flue gas volume flow rate variation range of the current propagation node is converted into the standard state flue gas volume flow rate variation range of the subsequent node.

[0082] In this embodiment, the process constraints correspond to the material balance relationship, heat transfer relationship, or flue gas transport relationship, and the corresponding equipment operating boundaries.

[0083] S3.4. Based on the current operating condition of the device corresponding to the current propagation node and the direction of increase or decrease in the range of operating condition changes, select the corresponding upper or lower operating boundary from the device operating boundary, and determine the remaining adjustment range based on the difference between the current operating condition and the selected operating boundary.

[0084] For the increasing direction, the remaining adjustment range consists of zero to the difference between the upper limit of equipment operation and the current operating condition of the equipment; for the decreasing direction, the remaining adjustment range consists of the difference between the lower limit of equipment operation and the current operating condition of the equipment to zero.

[0085] Before comparison, the equipment operating boundary uses the same physical quantity and unit as the corresponding response type operating condition quantity; therefore, the operating condition change range and the remaining adjustment range correspond to the same propagation node, the same operating condition quantity, and the same dimension.

[0086] For example, if the current sulfur mass flow rate at the copper matte transfer node is 2.85 kg / s, and the corresponding upper limit of equipment operation is 3.25 kg / s, then the remaining adjustment range for the increased direction is 0 to 0.40 kg / s.

[0087] S3.5. Compare the range of operating condition changes with the remaining adjustment range; if the range of operating condition changes is included in the remaining adjustment range, terminate the branch corresponding to the current propagation node; if the range of operating condition changes is not included in the remaining adjustment range, determine the device corresponding to the current propagation node as the device corresponding to the branch that has not been terminated, and transfer the range of operating condition changes to the subsequent node after conversion along the process constraint relationship.

[0088] For the increasing direction, if the upper limit of the operating condition change range is not higher than the upper limit of the remaining adjustment range, the branch corresponding to the current propagation node terminates; if the upper limit of the operating condition change range is higher than the upper limit of the remaining adjustment range, the device corresponding to the current propagation node is determined to be the device corresponding to the branch that has not terminated.

[0089] For the reduction direction, the lower limit of the operating condition change range is compared with the lower limit of the remaining adjustment range.

[0090] Branch termination indicates that the device corresponding to the current propagation node does not need to enter the further impact propagation calculation of this round of cross-process collaborative scheduling, but does not mean that the actual physical propagation of materials, heat or flue gas terminates at this node.

[0091] For example, the operating condition variation range of the copper matte transfer node is 0.50 kg / s to 0.60 kg / s, and the remaining adjustment range is 0 to 0.40 kg / s. Therefore, the equipment corresponding to the copper matte transfer node is determined to be the equipment corresponding to the non-terminating branch. The complete operating condition variation range is further converted to the inlet node of the blowing furnace according to the material transfer relationship of 0.90, resulting in the subsequent node operating condition variation range of 0.45 kg / s to 0.54 kg / s.

[0092] S3.6 Repeat the process of determining and comparing the remaining adjustment range for subsequent nodes until the corresponding branch terminates or reaches the end device, and retain the operating condition change range and remaining adjustment range of the corresponding device for each non-terminated branch.

[0093] After propagation to the end device, even if there is no successor node, as long as the range of operating condition changes is not included in the remaining adjustment range of the end device, the end device will still be retained as the device corresponding to the unterminated branch and will enter the subsequent scheduling quantity determination process.

[0094] This yields the equipment corresponding to the unterminated branches involved in this round of cross-process collaborative scheduling, the range of corresponding operating condition changes, and the remaining adjustment range.

[0095] S4. Using the equipment corresponding to the unterminated branch as the scheduling object, determine the scheduling quantity based on the portion of the operating condition change range exceeding the remaining adjustment range and the equipment adjustment response relationship, and adjust the scheduling object. Note that the following should be noted in this step:

[0096] This embodiment takes the aforementioned non-terminating branch of material response as a complete scheduling example, and takes the matte transfer feeding device and the blowing furnace feeding adjustment device as the corresponding scheduling objects.

[0097] S4.1 For each non-terminated branch corresponding equipment, the part of the operating condition change range that exceeds the remaining adjustment range is determined as the operating condition quantity to be compensated, and the adjustment direction of the scheduling object is determined according to the increase or decrease direction of the operating condition quantity to be compensated.

[0098] For the direction of increase, the portion of the operating condition change range that is higher than the upper limit of the remaining adjustment range is defined as the out-of-bounds interval; for the direction of decrease, the portion of the operating condition change range that is lower than the lower limit of the remaining adjustment range is defined as the out-of-bounds interval.

[0099] As a preferred implementation method, the larger value among the absolute amplitudes of the out-of-bounds intervals is used as the working condition quantity to be compensated in this round.

[0100] For example, if the operating condition varies from 0.50 kg / s to 0.60 kg / s and the upper limit of the remaining adjustment range is 0.40 kg / s, then the out-of-range amplitude is 0.10 kg / s to 0.20 kg / s, and the operating condition to be compensated is 0.20 kg / s.

[0101] The direction of increase or decrease in the actual control quantity of the equipment is determined by the direction of the working condition quantity to be compensated and the positive or negative sign of the equipment adjustment response relationship, so that the change in the equipment control quantity produces a working condition compensation quantity that is opposite to the direction of the working condition quantity to be compensated.

[0102] S4.2 Select adjacent operating points in the equipment regulation characteristic curve of the scheduling object, where the operating conditions are located on both sides of the current operating condition of the equipment and the control quantities are different. Determine the equipment regulation response relationship of the current operating range based on the difference in operating conditions and the difference in control quantities of the adjacent operating points.

[0103] The equipment adjustment characteristic curve can be obtained from the performance data generated during the equipment commissioning phase or from the corresponding data of control quantity and operating condition quantity under stable production conditions; the control quantity is the actual set quantity of the equipment, and the operating condition quantity is the actual physical quantity corresponding to the current propagation branch response type.

[0104] The equipment regulation response relationship is uniformly defined as the ratio of the change in operating condition quantity to the change in control quantity, that is, the change in operating condition quantity corresponding to a unit change in control quantity.

[0105] For example, if the feeding rates of two adjacent stable operating points of the copper matte transfer and feeding device are 14.0 kg / s and 15.0 kg / s, respectively, and the corresponding sulfur mass flow rates are 2.80 kg / s and 3.00 kg / s, then the equipment adjustment response relationship is 0.20 kg / s sulfur mass flow rate per 1 kg / s feeding rate.

[0106] If the difference in operating conditions between adjacent operating points is lower than the corresponding measurement resolution, or if the current operating interval does not have a definite monotonic relationship, the group of operating points will not be used to determine the equipment adjustment response relationship; instead, the nearest stable operating point with a definite monotonic relationship will be selected on both sides of the current operating condition of the equipment.

[0107] S4.3. Based on the equipment adjustment response relationship, the operating condition quantity to be compensated is converted into a candidate scheduling quantity. The upper or lower control boundary of the equipment is selected according to the adjustment direction. The candidate scheduling quantity is limited according to the difference between the current control quantity of the scheduling object and the control boundary of the selected equipment to obtain the executable scheduling quantity.

[0108] The scheduling quantity represents the control increment relative to the current control quantity of the equipment; the absolute value of the candidate scheduling quantity is obtained by dividing the operating condition quantity to be compensated by the absolute value of the equipment adjustment response relationship; the sign of the candidate scheduling quantity is determined according to the direction of operating condition change that can produce the opposite direction to the operating condition quantity to be compensated.

[0109] For example, if the mass flow rate of sulfur to be compensated is 0.20 kg / s, and the equipment adjustment response relationship is 0.20 kg / s sulfur mass flow rate per 1 kg / s feed rate, then the absolute value of the candidate scheduling quantity is 1.00 kg / s; an increase in the feed rate corresponds to an increase in the sulfur mass flow rate, therefore, the current positive sulfur disturbance corresponds to a decrease in the feed rate, and the candidate scheduling quantity is -1.00 kg / s.

[0110] If the current control quantity of the copper matte transfer and feeding device is 14.75 kg / s and the lower limit of equipment control is 14.25 kg / s, then the current allowable reduction in control quantity is 0.50 kg / s, and therefore the executable scheduling quantity is limited to -0.50 kg / s.

[0111] S4.4 When the same non-terminating branch has multiple scheduling objects, according to the process connection order of the scheduling objects in the corresponding propagation path, the executable scheduling quantity of the current scheduling object is converted into the working condition compensation quantity based on the equipment adjustment response relationship, the working condition compensation quantity is deducted from the working condition quantity to be compensated, and the uncompensated part is converted into the working condition quantity to be compensated of the next scheduling object along the process constraint relationship.

[0112] The absolute value of the dispatchable quantity of the copper matte transfer and feeding device is 0.50 kg / s. The equipment adjustment response relationship is 0.20 kg / s sulfur mass flow rate per 1 kg / s feeding rate. Therefore, the equipment forms a compensation quantity of 0.10 kg / s sulfur mass flow rate.

[0113] The original compensation quantity was 0.20 kg / s. After deducting the compensation quantity of 0.10 kg / s, the current propagation node retains an uncompensated quantity of 0.10 kg / s. The uncompensated quantity is converted according to the material transfer relationship of 0.90 from the matte transfer node to the inlet node of the smelting furnace. The compensation quantity of the next scheduling object is 0.09 kg / s.

[0114] If the current operating range of the feed regulating device for the blowing furnace is a feed rate of 0.20 kg / s per 1 kg / s of sulfur mass flow rate, then the 0.09 kg / s compensation condition quantity corresponds to the absolute value of the candidate control increment of 0.45 kg / s. When the equipment control boundary allows a reduction of 0.50 kg / s, the executable scheduling quantity is -0.45 kg / s, corresponding to a compensation quantity of 0.09 kg / s, and the current uncompensated condition quantity is reduced to zero.

[0115] When all scheduled objects in the same non-terminating branch have reached the equipment control boundary and there are still uncompensated operating conditions, each scheduled object only executes the executable scheduled quantities within the equipment control boundary, and retains the remaining uncompensated operating conditions as the part not eliminated in this round of scheduling.

[0116] S4.5. Adjust the scheduling objects sequentially according to the process connection sequence and the executable scheduling quantity.

[0117] The production control system converts the executable scheduling quantities of each scheduling object into the set value increments of the corresponding equipment controllers. Taking the aforementioned example, the set value of the feeding rate of the matte transfer feeding device is reduced by 0.50 kg / s, the set value of the feeding rate of the blowing furnace feeding regulating device is reduced by 0.45 kg / s, and each equipment controller adjusts the corresponding feeding actuator according to the updated feeding rate set value.

[0118] For non-terminating branches formed by heat response or flue gas response, the candidate scheduling quantity is determined according to the same processing relationship, based on the equipment adjustment response relationship between the actual control quantity of the corresponding scheduling object and the sensible heat flow rate of flue gas or the volumetric flow rate of flue gas under standard conditions, and then the executable scheduling quantity is determined according to the corresponding equipment control boundary.

[0119] Finally, the scheduling objects corresponding to the current source operating condition event propagation location and the executable scheduling quantities of each scheduling object are obtained, and the corresponding equipment controllers perform equipment control quantity adjustments.

[0120] It should be noted that this embodiment determines the actual concentrate batch corresponding to the downstream operating condition change based on the material batch connection relationship and the combined changes of multiple target elements such as copper, iron, and sulfur after smelting. On this basis, the propagation time is given based on the existing node propagation delay benchmark when the current event has not yet generated a downstream actual response. After the current event response is formed, the arrival time of the propagation node is updated by the time-shift cross-correlation result limited by the candidate arrival time period, so that the material, heat, and flue gas response form a continuously updated operating condition propagation phase index according to the actual process path. Furthermore, the operating condition change range is determined based on the node state changes between the continuous operating condition propagation phase indices, and propagated to subsequent equipment along the material balance relationship, heat transfer relationship, or flue gas transport relationship. The corresponding collaborative scheduling branch ends at the position where the operating condition change does not touch the remaining space of the equipment operation. For branches that have not been terminated, the out-of-bounds operating condition quantity is converted into a control increment based on the control quantity and operating condition response relationship of the current operating range of the equipment, and allocated to the corresponding scheduling object according to the equipment control boundary and process connection sequence. Thus, the equipment range and control quantity of cross-process collaborative scheduling correspond to the actual propagation stage and equipment operating status of the concentrate batch change.

[0121] Other aspects disclosed in the embodiments of the present invention also provide a computer device including one or more processors and a memory.

[0122] The memory is used to store operable instructions that, when executed by the one or more processors, cause the one or more processors to perform operations, including the flow of the cooperative scheduling method based on non-ferrous metal smelting working condition data index driven by the foregoing embodiments, especially... Figure 1 The flowchart of the method is shown.

[0123] Other aspects disclosed in the embodiments of the present invention also propose a computer-readable medium for storing software including instructions executable by one or more computers, which, upon execution, cause the one or more computers to perform operations including the flow of the cooperative scheduling method based on non-ferrous metal smelting working condition data index driven by the foregoing embodiments, particularly... Figure 1 The flowchart of the method is shown.

[0124] It should be recognized that embodiments of the present invention may be implemented or carried out by computer hardware, a combination of hardware and software, or by computer instructions stored in a non-transitory computer-readable storage medium.

[0125] The method can be implemented using standard programming techniques, including a non-transitory computer-readable storage medium configured with a computer program in the computer program, wherein the storage medium is configured such that the computer operates in a specific and predefined manner.

[0126] Each program can be implemented in a high-level procedural or object-oriented programming language to communicate with the computer system; however, if required, the program can be implemented in assembly or machine language.

[0127] In any case, the language can be a compiled or interpreted language.

[0128] Furthermore, for this purpose, the program can run on a programmed application-specific integrated circuit.

[0129] The processes described herein (or variations and / or combinations thereof) can be executed under the control of one or more computer systems configured with executable instructions, and can be implemented by hardware or a combination thereof as code (e.g., executable instructions, one or more computer programs, or one or more applications) that commonly executes on one or more processors. The computer program includes a plurality of instructions executable by one or more processors.

[0130] Furthermore, the method can be implemented in any suitable computing platform, including but not limited to personal computers, minicomputers, mainframes, workstations, networked or distributed computing environments, standalone or integrated computer platforms, or in communication with charged particle tools or other imaging devices.

[0131] Various aspects of the present invention can be implemented in machine-readable code stored on a non-transitory storage medium or device, whether portable or integrated into a computing platform, such as a hard disk, optical read and / or write storage medium, RAM, ROM, etc., such that it can be read by a programmable computer, and when the storage medium or device is read by the computer, it can be used to configure and operate the computer to perform the processes described herein.

[0132] Furthermore, machine-readable code, or parts thereof, can be transmitted via wired or wireless networks.

[0133] When such media includes instructions or programs that combine with a microprocessor or other data processor to implement the steps described above, the invention described herein includes these and other different types of non-transitory computer-readable storage media.

[0134] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A collaborative scheduling method driven by non-ferrous metal smelting operating condition data index, characterized in that, include: The scope of the source condition event is limited based on the material transfer relationship of the process, and the source condition event of the downstream condition change is determined according to the consistency of changes in multiple target elements. Based on the arrival order of the material, heat, and flue gas responses of the source operating condition event, the arrival node and subsequent node are determined to form an operating condition propagation phase index. Compare the propagation phase index of the current operating condition with the propagation phase index of the previous operating condition to determine the change node and the range of operating condition change, and propagate along the process constraint relationship; determine the remaining adjustment range according to the current operating condition and the operating boundary of the equipment, and terminate the branch when the range of operating condition change does not exceed the remaining adjustment range. Using the equipment corresponding to the unterminated branch as the scheduling object, the scheduling amount is determined based on the portion of the operating condition change range that exceeds the remaining adjustment range and the equipment adjustment response relationship, and the scheduling object is adjusted accordingly.

2. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 1, characterized in that, Determining the source condition event includes: Based on the material conveying path and intermediate material batch connection relationship between adjacent smelting processes, the upstream operations that have a material connection relationship with the downstream operating condition change are traced back to determine the candidate source operating condition event. Based on the target element composition of the material corresponding to the candidate source working condition event and the material allocation relationship between processes, the direction of increase or decrease and the relative change of the target elements after the candidate source working condition event is transmitted to the downstream process are calculated to obtain the candidate target element change combination. The direction of increase or decrease of the same target element in the downstream operating condition change is compared with the combination of changes of the candidate target element element element by element, and candidate source operating condition events with inconsistent direction of increase or decrease of target element are eliminated. Calculate the absolute deviation between the relative change magnitude of each target element in the remaining candidate source operating condition events and the relative change magnitude of the same target element in the downstream operating condition changes, and determine the candidate source operating condition event whose absolute deviation is within the corresponding allowable deviation range and whose sum of absolute deviation is the smallest as the source operating condition event.

3. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 2, characterized in that, The calculation of the direction and relative change of the target element after the candidate source condition event is transmitted to the downstream process includes: Based on the material processing volume and target element mass fraction within the same statistical interval before and after the occurrence of the candidate source working condition event, the product of the material processing volume after the event and the target element mass fraction is subtracted from the corresponding product before the event to determine the source-end change of each target element. Based on the material diversion ratio in the material transfer relationship and the retention ratio of the target element in the intermediate processing stage, the element change amount of the source change amount is calculated along each material transfer path to the downstream process, and the element change amounts of each material transfer path are added together to determine the net change amount of the target element in the downstream process. The direction of increase or decrease of the target element is determined based on the comparison result between the net change of the target element and the zero value determination threshold. The zero value determination threshold is determined based on the target element mass fraction detection error and the material processing quantity measurement error. The product of the material processing volume of the downstream process receiving the material before the occurrence of the candidate source condition event and the mass fraction of the target element is determined as the reference element quantity of the target element. The relative change range of the target element is determined by the ratio of the absolute value of the net change of the target element to the larger value between the reference element quantity of the target element and the lower limit of the reference element quantity. The candidate target element variation combination is formed based on the direction of increase or decrease of the target element and the relative change range of the target element.

4. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 1, characterized in that, The determination of the arriving node and the successor node includes: Based on the material conveying connection, heat exchange connection, and flue gas pipeline connection between the source operating condition event occurrence node and the downstream process, the material response propagation path, heat response propagation path, and flue gas response propagation path are determined respectively, and the propagation nodes are arranged according to the process connection sequence. Based on the material residence time range, thermal response time range, and flue gas transport time range between adjacent propagation nodes, the candidate arrival time periods of the source condition event to each propagation node are determined by accumulating them step by step along the material response propagation path, thermal response propagation path, and flue gas response propagation path. Within the candidate arrival time period, the response change at the source condition event occurrence node and the condition change of the same response type at the propagation node are cross-correlated by time shift calculation, and the arrival time of the propagation node is determined by the time shift when the cross-correlation coefficient is maximized. Along the corresponding propagation path, the last propagation node whose arrival time is no later than the current time is determined as the arrival node, and the propagation node that is adjacent to the arrival node and whose arrival time is later than the current time is determined as the successor node.

5. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 4, characterized in that, The formation of the operating condition propagation phase index includes: The response types are classified according to material response, heat response, and flue gas response, and the propagation nodes of the same response type are arranged in the process connection order of the corresponding propagation path. The arrival time of each propagation node is configured to the corresponding propagation node. Propagation nodes whose arrival time is no later than the current time are determined to be in the arrived state, and propagation nodes whose arrival time is later than the current time are determined to be in the unarrived state. The current propagation position of the source event under the corresponding response type is determined by the arriving node, the successor node, and the propagation segment between the two nodes. An index hierarchy is established according to the source operating condition event, response type, and propagation node order, and the arrival time and current propagation position are configured under the corresponding propagation node to construct the operating condition propagation phase index.

6. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 5, characterized in that, The determination of the changing nodes and the range of operating condition changes includes: Compare the propagation phase index of the operating condition with the propagation phase index of the previous operating condition, and determine the propagation node status of the same source operating condition event and the same response type. The propagation node that changes from the non-arrived state to the arrived state is identified as the changed node. The operating condition change range of the changing node is determined by comparing the operating condition quantity corresponding to the response type of the changing node at the current arrival time with the operating condition quantity corresponding to the response type of the changing node in the previous operating condition propagation phase index, and by the difference between the two and the measurement error of the operating condition quantity. Select the material balance relationship, heat transfer relationship, or flue gas transport relationship according to the response type of the changed node, convert the operating condition change range into the operating condition change range of the subsequent node, and propagate it level by level along the corresponding propagation path.

7. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 6, characterized in that, Also includes: Based on the current operating condition of the device corresponding to the current propagation node and the direction of increase or decrease of the range of change of the operating condition, select the corresponding upper or lower operating boundary from the device operating boundary, and determine the remaining adjustment range based on the difference between the current operating condition and the selected operating boundary. Compare the range of operating condition changes with the remaining adjustment range; if the range of operating condition changes is included in the remaining adjustment range, terminate the branch corresponding to the current propagation node. When the operating condition change range is not included in the remaining adjustment range, the device corresponding to the current propagation node is determined as the device corresponding to the non-terminated branch, and the operating condition change range is converted along the process constraint relationship and then passed to the successor node. The process of determining and comparing the remaining adjustment range is repeated for subsequent nodes until the corresponding branch terminates or reaches the end device, and the operating condition change range and remaining adjustment range at the corresponding device of each non-terminated branch are retained.

8. The collaborative scheduling method based on non-ferrous metal smelting working condition data index driven according to claim 7, characterized in that, Adjusting the scheduling object includes: For each device corresponding to a non-terminated branch, the portion of the operating condition change range that exceeds the remaining adjustment range is identified as the operating condition quantity to be compensated, and the adjustment direction of the scheduling object is determined according to the increase or decrease direction of the operating condition quantity to be compensated. In the equipment regulation characteristic curve of the scheduling object, select adjacent operating points whose operating conditions are located on both sides of the current operating condition of the equipment and whose control quantities are different. Determine the equipment regulation response relationship of the current operating range based on the difference in operating conditions and the difference in control quantities of the adjacent operating points. Based on the equipment adjustment response relationship, the operating condition quantity to be compensated is converted into a candidate scheduling quantity. The upper or lower control boundary of the equipment is selected according to the adjustment direction. The candidate scheduling quantity is limited according to the difference between the current control quantity of the scheduling object and the selected equipment control boundary to obtain the executable scheduling quantity. When the same non-terminating branch has multiple scheduling objects, the executable scheduling quantity of the current scheduling object is converted into the working condition compensation quantity according to the process connection order of the scheduling object in the corresponding propagation path. The working condition compensation quantity is deducted from the working condition quantity to be compensated, and the uncompensated part is converted into the working condition quantity to be compensated of the next scheduling object along the process constraint relationship. The scheduling objects are adjusted sequentially according to the process connection sequence and the executable scheduling quantity.

9. A computer device, characterized in that, include: One or more processors; The memory stores operable instructions that, when executed by the one or more processors, cause the one or more processors to perform operations, including the flow of the cooperative scheduling method based on non-ferrous metal smelting working condition data index as described in any one of claims 1 to 8.

10. A computer-readable medium for storing software, characterized in that: The software includes instructions executable by one or more computers, which cause the one or more computers to perform operations, including the flow of the collaborative scheduling method based on non-ferrous metal smelting working condition data index as described in any one of claims 1 to 8.

Images