Multi-working condition stable control method for liquid carbon dioxide energy storage process

Abstract

The application discloses a kind of multi-working-condition stable control methods of liquid carbon dioxide energy storage process, it is related to stable control technical field, comprising: according to the area data of entrance area, the occupancy sorting of entrance area is carried out, and the occupancy sorting result is obtained;Based on the area sorting list in occupancy sorting result, determine main habitat layer, secondary habitat layer and free layer, and according to the power instruction of grid corresponding to liquid carbon dioxide energy storage system, generate candidate load sequence;Analysis target load sequence corresponding to candidate load sequence, obtain target power load share;Target power load share is handled, and script release coefficient is obtained;Based on script release coefficient and grid power instruction, obtain grid-connected power token;The application improves the consistency of inlet flow state and energy distribution state under multi-working-condition.

Description

Technical Field

[0001] This invention relates to the field of stability control technology, and in particular to a multi-condition stability control method for a liquid carbon dioxide energy storage process. Background Technology

[0002] Liquid carbon dioxide energy storage technology utilizes electricity during off-peak hours to drive the compression, liquefaction, and storage process. During peak hours, it uses the heating and expansion of liquid carbon dioxide to drive turbine power generation, enabling cross-period storage and release of electrical energy. It has the advantages of large storage capacity, fast response speed, and no combustion emissions, and is gradually becoming an important energy storage method for new power systems.

[0003] In actual operation, liquid carbon dioxide energy storage devices need to frequently undergo charging, discharging, standby, hot standby, and switching between different operating conditions. The working fluid in the inlet distribution manifold of the compact heat exchanger constantly undergoes re-diversion and re-merging. Due to the small internal space of the inlet distribution manifold, each inlet area is affected by local flow resistance, thermal state changes, and working fluid energy distribution. Liquid carbon dioxide working fluid tends to form a continuously occupied flow state in some inlet areas, causing some inlet areas to bear a large working fluid flow and energy input for a long time, while other inlet areas bear a relatively low working fluid flow and energy input, thus forming a stable but dynamically changing occupancy hierarchy among the inlet areas. As the operating conditions change continuously, this occupancy hierarchy will shift, further affecting the energy exchange capacity of the corresponding heat exchange channels and the grid-connected charging and discharging process of the entire energy storage device.

[0004] Existing grid-connected power control methods for liquid carbon dioxide energy storage devices typically rely on grid power commands as the primary control basis. They directly adjust the operating status of compressors, expanders, or related actuators based on target charging or discharging power, with little consideration for the long-term hierarchical relationship of working fluid occupancy within the inlet distribution manifold. When operating conditions switch, existing control methods still adjust according to a uniform power distribution strategy, failing to reflect the differences in the actual working fluid transport and energy transfer capabilities of different inlet areas. This leads to inconsistencies between the power load relationship between heat exchange channels and the current actual occupancy status of the inlet areas. Consequently, some heat exchange channels bear high energy exchange loads for extended periods, while others are underutilized. This can easily cause increased fluctuations in grid-connected charging and discharging power during multi-condition switching, affecting the grid-connected power stability of the liquid carbon dioxide energy storage device under multiple operating conditions. Summary of the Invention

[0005] The purpose of this invention is to address the shortcomings of existing technologies that fail to reflect the differences in the actual working fluid transport and energy transfer capabilities of different inlet regions, and to propose a multi-condition stable control method for liquid carbon dioxide energy storage processes.

[0006] To address the problems existing in the prior art, the present invention adopts the following technical solution: A multi-condition stability control method for a liquid carbon dioxide energy storage process includes: S1. Based on the regional data of the entrance area, sort the entrance area to obtain the sorting result. S2. Based on the regional ranking list in the occupancy ranking results, determine the primary habitation layer, secondary habitation layer, and free layer, and generate a candidate load-bearing sequence according to the power grid power command corresponding to the liquid carbon dioxide energy storage system. S3. Analyze the target load sequence corresponding to the candidate load sequence to obtain the target power load share; S4. Process the target power carrying share to obtain the script release coefficient; based on the script release coefficient and the grid power command, obtain the grid-connected power token; S5. Based on the grid-connected power token and the grid-connected power interface data of the liquid carbon dioxide energy storage system, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

[0007] Preferably, based on the regional data of the entrance area, the entrance area is sorted by placeholders to obtain the placeholder sorting result, including: Determine the inlet region of the inlet distribution manifold connection in the compact heat exchanger; Acquire regional data for the inlet area; this regional data includes mass flow rate data, temperature data, and pressure data. Based on the temperature and pressure data in the regional data, determine the regional specific enthalpy of the inlet region; Calculate the inlet occupancy based on mass flow rate data and regional specific enthalpy; Based on the number of spaces occupied at the entrance, the entrance area is sorted to obtain the sorting result.

[0008] Preferably, based on the regional ranking list in the placeholder ranking results, the primary habitation layer, secondary habitation layer, and free layer are determined, including: Retrieve the list of regions sorted from the placeholder sorting results; The total number of entry areas is determined based on the regional sorting list; Based on the total number of entrance areas, the list of entrance areas is divided into the first sorting segment, the second sorting segment, and the third sorting segment in sequence. The entrance area in the first sorting section is designated as the primary residential floor; The entrance area in the second sorting segment is designated as the secondary residential layer; The entry region in the third sorting segment is identified as the free layer.

[0009] Preferably, generating a candidate bearer sequence includes: Determine the heat exchange channel corresponding to the inlet area; The heat exchange channels corresponding to the primary habitation layer are classified into the fixed channel family; The heat exchange channels corresponding to the secondary habitation layer are classified into the transition channel family; The heat exchange channels corresponding to the free layer are classified into the replacement channel family; Obtain the grid power command corresponding to the liquid carbon dioxide energy storage system; wherein, the grid power command includes the target charging power or the target discharging power; Based on the fixed channel family, transition channel family, replacement channel family and grid power command, several candidate bearer sequences are generated.

[0010] Preferably, parsing the target load sequence corresponding to the candidate load sequence to obtain the target power load share includes: The candidate carrier sequence is iteratively updated using the particle swarm optimization algorithm to obtain the target carrier sequence. The target power carrying capacity sequence is analyzed to obtain the target power carrying capacity share corresponding to the fixed channel family, the transition channel family, and the replacement channel family, respectively.

[0011] Preferably, the target power carrying share is processed to obtain the script release coefficient, including: The target power carrying capacity shares corresponding to the fixed channel family, the transition channel family, and the replacement channel family are normalized to obtain the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the replacement channel family. The total carrying capacity of a channel family is calculated based on the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the successor channel family. Based on the total share of the channel family, the normalized share of the fixed channel family, the transition channel family, and the successor channel family is folded to obtain the script release coefficient.

[0012] Preferably, the grid-connected power token is obtained based on the script release coefficient and the grid power command, including: Extract the power value and power direction identifier from the power grid command; Generate a token power value based on the script release coefficient and power value; By associating the token power value with the power direction identifier, a grid-connected power token is obtained.

[0013] Preferably, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled, including: Acquire grid-connected power interface data of the liquid carbon dioxide energy storage system; wherein, the grid-connected power interface data includes charging power interface data and discharging power interface data; When the power direction identifier in the grid-connected power token is the charging direction identifier, the token power value in the grid-connected power token is written into the charging power interface data. When the power direction identifier in the grid-connected power token is the discharge direction identifier, the token power value in the grid-connected power token is written into the discharge power interface data. Based on the grid-connected power interface data after the power value of the token is written, grid-connected power execution data is generated; Based on the grid-connected power execution data, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

[0014] Preferably, the candidate carrier sequence is iteratively updated using a particle swarm optimization algorithm to obtain the target carrier sequence, including: Each candidate bearer sequence is treated as a particle in the particle swarm algorithm. The initial velocity of the particle is determined based on the power carrying share in the candidate carrying sequence; The initial position of the particle is determined based on the channel family arrangement order and power carrying share in the candidate carrying sequence; Based on the initial velocity and initial position of the particles, construct an initial particle swarm; The fitness of particles is calculated based on the candidate carrier sequences; Based on the particle's fitness, determine the individual optimal particle and the global optimal particle; Based on the individual optimal particle and the global optimal particle, the initial particle swarm is iteratively updated to obtain the particle swarm after the iteration. The target carrying sequence is obtained from the particle swarm after the iteration ends.

[0015] Compared with the prior art, the beneficial effects of the present invention are: 1. This invention obtains mass flow rate, temperature and pressure data of the inlet region of a compact heat exchanger, calculates the inlet occupancy and establishes an occupancy ranking relationship, and divides the inlet region into a primary habitable layer, a secondary habitable layer and a free layer. This enables a quantitative description of the actual occupancy state of the working fluid inside the inlet distribution manifold, allowing subsequent power distribution to be adjusted according to the actual carrying capacity of the inlet region, and improving the consistency between the inlet flow state and energy distribution state under multiple operating conditions.

[0016] 2. This invention establishes a fixed channel family, a transition channel family, and a successor channel family based on the primary habitable layer, the secondary habitable layer, and the free layer. It also generates a candidate load-bearing sequence in conjunction with the power grid command and obtains the target load-bearing sequence using the particle swarm optimization algorithm. This allows the power load-bearing relationship corresponding to different operating conditions to be adjusted synchronously with the inlet occupancy status, avoiding different heat exchange channels from bearing unbalanced energy exchange tasks for a long time and improving the rationality of power allocation.

[0017] 3. This invention generates a script release coefficient based on the target bearing sequence, and further generates a grid-connected power token and grid-connected power execution data. The power execution process is based on the change of the occupancy level in the inlet area, so that the liquid carbon dioxide energy storage device can maintain stable grid-connected charging and discharging power during multiple operating conditions such as charging, discharging, standby and switching, reduce power fluctuations caused by operating condition switching and improve grid-connected operation stability. Attached Figure Description

[0018] The accompanying drawings, which are included to provide a further understanding of the invention and form part of this application, illustrate exemplary embodiments of the invention and, together with their description, serve to explain the invention and do not constitute an undue limitation thereof. In the drawings: Figure 1 This is a flowchart illustrating a multi-condition stable control method for a liquid carbon dioxide energy storage process, provided in an embodiment of the present invention. Detailed Implementation

[0019] The technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments.

[0020] This embodiment provides a multi-condition stability control method for a liquid carbon dioxide energy storage process. (See also...) Figure 1 Specifically, including: S1. Based on the regional data of the entrance area, sort the entrance area to obtain the sorting result. In an embodiment of the present invention, the entrance area is sorted based on the area data of the entrance area to obtain a placeholder sorting result, including: Determine the inlet region of the inlet distribution manifold connection in the compact heat exchanger; A compact heat exchanger is a heat exchange device composed of a large number of micro heat exchange channels, capable of completing efficient heat exchange within a limited space. Its internal flow channel size is small and the heat exchange area density is high, enabling liquid carbon dioxide to quickly complete heat exchange as it flows through each heat exchange channel. An inlet distribution manifold is a fluid distribution cavity located on the inlet side of the compact heat exchanger, used to distribute the incoming liquid carbon dioxide working fluid to each inlet region. It plays a role in guiding and distributing the working fluid. An inlet region refers to the independent working fluid entry position formed between the inlet distribution manifold and each heat exchange channel. Each inlet region corresponds to the inlet position of one or a group of heat exchange channels. Specifically, based on the inlet end structure diagram of the compact heat exchanger, the flow channel layout diagram of the inlet distribution manifold, and the inlet numbering table of the heat exchange channels, the location of each branch outlet of the inlet distribution manifold is confirmed. Along the flow direction of the liquid carbon dioxide working fluid from the inlet distribution manifold into the heat exchange channel, the inlet range of the heat exchange channel covered by each branch outlet is determined as an inlet area, and a corresponding record table is established for each inlet area with area number, branch outlet number, corresponding heat exchange channel number, and measuring point number.

[0021] Acquire regional data for the inlet area; this regional data includes mass flow rate data, temperature data, and pressure data. It should be noted that mass flow rate data refers to the mass of liquid carbon dioxide working fluid flowing through the corresponding inlet area per unit time. It is data that characterizes the working fluid transport capacity and working fluid supply level of the inlet area. It is used to reflect the actual mass of working fluid entering the heat exchanger under the current operating conditions and provides basic data for calculating the energy input level and inlet occupancy of the inlet area.

[0022] Specifically, mass flow rate measurement points, temperature measurement points, and pressure measurement points are arranged in each determined inlet area. Mass flow rate measurement points are set on the working medium inlet channel corresponding to the inlet area, temperature measurement points are set at the working medium inlet section of the inlet area, and pressure measurement points are set at the working medium inlet section of the inlet area. The mass flow rate data, temperature data, and pressure data collected in the same inlet area are collected according to the area number to form the area data of that inlet area.

[0023] Based on the temperature and pressure data in the regional data, determine the regional specific enthalpy of the inlet region; The specific enthalpy of a region refers to the total energy state of a unit mass of liquid carbon dioxide working fluid in the corresponding inlet region. It is used to characterize the energy level at which the working fluid can participate in subsequent energy exchange when it enters the heat exchanger. Specifically, the region number of the inlet area is used as a data index. The temperature and pressure data corresponding to the same inlet area are combined to form the carbon dioxide working fluid state data. The temperature data is used as the state temperature when the liquid carbon dioxide working fluid enters the compact heat exchanger, and the pressure data is used as the state pressure when the liquid carbon dioxide working fluid enters the compact heat exchanger. The state temperature and state pressure are input into the carbon dioxide working fluid state equation to obtain the unit mass enthalpy value corresponding to the state temperature and state pressure. The unit mass enthalpy value is used as the region specific enthalpy of the inlet area.

[0024] Specifically, the regional state temperature and regional state pressure are substituted into the Punger-Robinson equation of state. The Punger-Robinson equation of state states states that pressure equals the product of the gas constant and the regional state temperature, divided by the difference between the molar volume and the co-volume parameter, and then subtracted from the combined term formed by the attraction parameter divided by the molar volume, the sum of the molar volume and the co-volume parameter, and the difference between the co-volume parameter and the molar volume. The attraction parameter and co-volume parameter are calculated based on the carbon dioxide critical temperature, carbon dioxide critical pressure, and carbon dioxide eccentricity factor. The compressibility factor of the carbon dioxide working fluid under the corresponding state in the inlet region is solved based on the regional state temperature, regional state pressure, attraction parameter, and co-volume parameter. The residual enthalpy of the carbon dioxide working fluid is calculated based on the compressibility factor, regional state temperature, regional state pressure, attraction parameter, and co-volume parameter. The ideal enthalpy of the carbon dioxide working fluid is obtained by integrating the isobaric specific heat of the ideal gas with respect to temperature based on the regional state temperature. The ideal enthalpy and the residual enthalpy are added together to obtain the molar enthalpy of the carbon dioxide working fluid under the corresponding state in the inlet region. The molar enthalpy is divided by the molar mass of carbon dioxide to obtain the regional specific enthalpy of the inlet region.

[0025] Calculate the inlet occupancy based on mass flow rate data and regional specific enthalpy; Inlet occupancy refers to a comprehensive quantity calculated based on inlet area mass flow rate data and area specific enthalpy. It is used to reflect the actual occupancy degree of the corresponding inlet area currently undertaking working fluid input and energy input. The larger the occupancy, the more concentrated the working fluid transport and energy input undertaken by the corresponding inlet area. Mass flow rate data and regional enthalpy are extracted from the same inlet region according to the region number of the inlet region. The mass flow rate data is used as the mass of liquid carbon dioxide entering the compact heat exchanger per unit time in the inlet region, and the regional enthalpy is used as the energy data carried by the liquid carbon dioxide per unit mass in the inlet region. The mass flow rate data and the regional enthalpy are multiplied to obtain the inlet occupancy of the inlet region, and the inlet occupancy is bound to the region number of the corresponding inlet region to form the inlet region occupancy data.

[0026] It should be noted that mass flow rate represents the mass of liquid carbon dioxide working fluid entering the inlet region per unit time, and specific enthalpy represents the total energy carried by a unit mass of liquid carbon dioxide working fluid. According to the thermodynamic law of conservation of energy, the product of mass flow rate and specific enthalpy represents the total energy input through the inlet region per unit time, i.e., the energy flow rate. Therefore, the larger the mass flow rate, the more working fluid enters the inlet region, and the larger the specific enthalpy, the higher the energy carried by a unit mass of working fluid. Both together determine the actual energy input capacity of the inlet region. Using only mass flow rate cannot reflect the difference in energy carried by the working fluid under different thermodynamic states, and using only specific enthalpy cannot reflect the mass of working fluid entering the inlet region per unit time. Therefore, the product of mass flow rate and specific enthalpy is used as the inlet occupancy, which can simultaneously characterize the comprehensive level of working fluid transport and energy input undertaken by the inlet region. The larger the inlet occupancy, the higher the working fluid input and energy input undertaken by the corresponding inlet region per unit time, and the greater the degree of occupation of the inlet energy exchange resources. Therefore, it can be used as a physical quantity to measure the degree of occupancy of the inlet region.

[0027] Based on the number of spaces occupied at the entrance, the entrance area is sorted to obtain the sorting result.

[0028] The occupancy sorting result refers to the ordered result formed by arranging the occupancy of each inlet area from largest to smallest, which is used to reflect the hierarchical relationship of the current working medium input and energy input undertaken by different inlet areas.

[0029] Specifically, all entrance area placeholder data are arranged in descending order of entrance placeholder quantity. Entrance areas with the same entrance placeholder quantity are arranged in order of their area numbers. After sorting, a placeholder sequence number is generated for each entrance area. The area number, entrance placeholder quantity, and placeholder sequence number are combined to form a placeholder sorting record. All placeholder sorting records are then summarized to form the placeholder sorting result.

[0030] S2. Based on the regional ranking list in the occupancy ranking results, determine the primary habitation layer, secondary habitation layer, and free layer, and generate a candidate load-bearing sequence according to the power grid power command corresponding to the liquid carbon dioxide energy storage system. In embodiments of the present invention, the primary habitation layer, secondary habitation layer, and free layer are determined based on the regional ranking list in the occupancy ranking results, including: Retrieve the list of regions sorted from the placeholder sorting results; The regional sorting list refers to an ordered record arranged from largest to smallest according to the number of entry slots corresponding to each entry region. It is used to reflect the relative order in which each entry region undertakes working fluid input and energy input. Each sorting position corresponds to an entry region and its order of occupation. Specifically, the generated placeholder sorting results are used as the data source. All placeholder sorting records contained in the placeholder sorting results are extracted. Each placeholder sorting record includes an entry area identifier, an entry space quantity, and a placeholder sequence number. All placeholder sorting records are arranged in order from front to back according to the placeholder sequence number. The arranged entry area identifiers are summarized into an area sorting list, and each entry area identifier in the area sorting list maintains the corresponding placeholder sequence number and entry space quantity.

[0031] The total number of entry areas is determined based on the regional sorting list; Specifically, each entry region identifier in the region sorting list is counted, and the total number of entry regions is determined after the count is completed. The total number of entry regions is then bound to the region sorting list, so that the total number of entry regions is used to limit the scope of the region sorting list and avoid missing or reusing entry regions.

[0032] Based on the total number of entrance areas, the list of entrance areas is divided into the first sorting segment, the second sorting segment, and the third sorting segment in sequence. The first sorting segment refers to the set of consecutive regions at the beginning of the region sorting list. The entrance regions within this segment have a high entrance occupancy rate, indicating that they bear a high level of working fluid and energy input. The second sorting segment refers to the set of consecutive regions in the middle of the region sorting list. The entrance regions within this segment bear working fluid and energy input between the preceding and following regions, indicating that they are in a medium energy-carrying state. The third sorting segment refers to the set of consecutive regions at the end of the region sorting list. The entrance regions within this segment have a low entrance occupancy rate, indicating that they bear a low level of working fluid and energy input. Divide the total number of entry areas by three to obtain the number of basic segments. Remove the decimal part of the number of basic segments to obtain the number of integer segments. Determine the entry area from the first sorting position in the area sorting list to the sorting position corresponding to the integer segment number as the first sorting segment. Determine the entry area from the sorting position corresponding to the integer segment number plus one to the sorting position corresponding to the twice integer segment number as the second sorting segment. Determine the entry area from the sorting position corresponding to the twice integer segment number plus one to the sorting position corresponding to the total number of entry areas as the third sorting segment.

[0033] The entrance area in the first sorting section is designated as the primary residential floor; The entrance area in the second sorting segment is designated as the secondary residential layer; The entry region in the third sorting segment is identified as the free layer.

[0034] The primary habitation layer refers to the set of regions formed by the inlet regions corresponding to the first sorting segment, which represents the inlet regions where the liquid carbon dioxide working fluid preferentially undertakes energy input under the current operating state; the secondary habitation layer refers to the set of regions formed by the inlet regions corresponding to the second sorting segment, which represents the inlet regions where the liquid carbon dioxide working fluid undertakes transitional energy input; the free layer refers to the set of regions formed by the inlet regions corresponding to the third sorting segment, which represents the inlet regions where the liquid carbon dioxide working fluid currently undertakes less energy input and has room for adjustment.

[0035] Specifically, extract all entry area identifiers contained in the first sorting segment, bind the segment identifier of the first sorting segment with all entry area identifiers contained in the first sorting segment to generate the primary habitation layer, and record the primary habitation layer as the set of entry areas with the highest occupancy. Extract all entry area identifiers contained in the second sorting segment, bind the segment identifier of the second sorting segment with all entry area identifiers contained in the second sorting segment to generate the secondary habitation layer, and record the secondary habitation layer as the set of entry areas with the middle occupancy. Extract all entry area identifiers contained in the third sorting segment, bind the segment identifier of the third sorting segment with all entry area identifiers contained in the third sorting segment to generate the free layer, and record the free layer as the set of entry areas with the lowest occupancy.

[0036] It should be noted that the inlet occupancy of the inlet region reflects the comprehensive capacity of the corresponding inlet region to handle working fluid transport and energy input per unit time. Different occupancy amounts mean that the attraction capacity of each inlet region for liquid carbon dioxide working fluid varies. If subsequent control is based solely on the inlet occupancy amount, there is a lack of hierarchical relationship between the inlet regions, making it difficult to reflect the preferential and continuous accumulation positions of the working fluid within the inlet distribution manifold. Therefore, based on the inlet occupancy ranking results, the inlet regions are divided into the primary habitable layer, the secondary habitable layer, and the free layer, so that the inlet regions handling high energy input, medium energy input, and low energy input form a stable hierarchical structure. This transforms the continuously changing inlet occupancy amount into a set of inlet regions with a clear hierarchical relationship, providing a unified data foundation for subsequent channel family division and grid-connected power allocation under multiple operating conditions, and ensuring that the spatial distribution relationship of manifold habitability can be continuously characterized under different operating conditions.

[0037] In an embodiment of the present invention, generating a candidate bearer sequence includes: Determine the heat exchange channel corresponding to the inlet area; A heat exchange channel refers to an independent flow channel inside a compact heat exchanger that supplies liquid carbon dioxide working fluid and completes heat exchange. Each heat exchange channel corresponds to an inlet area in the inlet distribution manifold. As the liquid carbon dioxide working fluid flows along the heat exchange channel, it exchanges heat with the heat exchange wall, thereby completing the energy storage or energy release process. Specifically, based on the inlet distribution manifold flow channel arrangement, inlet area number, and heat exchange channel inlet number of the compact heat exchanger, a correspondence between inlet areas and heat exchange channels is established. All heat exchange channels connected to the same inlet area are recorded as the heat exchange channels corresponding to that inlet area. The inlet area number, the branch outlet number of the inlet area, and the corresponding heat exchange channel number are written into the inlet area channel correspondence table.

[0038] The heat exchange channels corresponding to the primary habitation layer are classified into the fixed channel family; The fixed channel family refers to the set of all heat exchange channels corresponding to the main habitation layer. This set of channels corresponds to the inlet region with a high inlet occupancy, which represents the main flow path through which the liquid carbon dioxide working fluid maintains a high flow rate and high energy input for a long time. Specifically, based on the entrance area numbers contained in the primary habitation layer, all heat exchange channel numbers corresponding to these entrance areas are extracted from the entrance area channel correspondence table. The extracted heat exchange channel numbers are then grouped into the same channel set, and this channel set is marked as the fixed channel family.

[0039] The heat exchange channels corresponding to the secondary habitation layer are classified into the transition channel family; The transition channel family refers to the set of all heat exchange channels corresponding to the secondary habitable layer. This channel family corresponds to the inlet region with an inlet occupancy of intermediate level, which is used to undertake the flow of working fluid and energy transfer between the primary habitable layer and the free layer. Specifically, when classifying the heat exchange channels corresponding to the secondary habitation layer into the transition channel family, based on the inlet region numbers contained in the secondary habitation layer, all heat exchange channel numbers corresponding to these inlet regions are extracted from the inlet region channel correspondence table, and the extracted heat exchange channel numbers are collected into the same channel set, which is then marked as the transition channel family.

[0040] The heat exchange channels corresponding to the free layer are classified into the replacement channel family; The successor channel family refers to the channel set consisting of all heat exchange channels corresponding to the free layer. This channel set corresponds to the inlet area with a low inlet occupancy and is used to receive the transfer of working fluid flow and energy input when the operating conditions change. Specifically, based on the inlet region numbers contained in the free layer, all heat exchange channel numbers corresponding to these inlet regions are extracted from the inlet region channel correspondence table. The extracted heat exchange channel numbers are then grouped into the same channel set, and this channel set is marked as the successor channel family.

[0041] It should be noted that after liquid carbon dioxide forms manifold habitability within the inlet distribution manifold, the heat exchange channels corresponding to different inlet regions already have different working fluid carrying capacity and energy exchange capacity. If all heat exchange channels are still treated as the same object for subsequent power distribution, the already formed working fluid flow hierarchy will be disrupted, causing high-occupancy regions and low-occupancy regions to undertake the same regulation task, resulting in the working fluid migration direction being inconsistent with the actual energy input path. Therefore, based on the primary habitable layer, secondary habitable layer, and free layer, fixed channel families, transition channel families, and successor channel families are established respectively, so that heat exchange channels with different energy carrying capacities form mutually independent channel sets. The spatial occupancy relationship of the inlet region is mapped to the carrying capacity relationship of the heat exchange channels, thereby ensuring that the subsequent power distribution of the power grid and the generation of the carrying capacity sequence are based on the actual distribution state of manifold habitability, and achieving consistency between the working fluid flow path and the energy transfer path under multiple operating conditions.

[0042] Obtain the grid power command corresponding to the liquid carbon dioxide energy storage system; wherein, the grid power command includes the target charging power or the target discharging power; The grid power command refers to the power regulation demand data sent by the grid side to the liquid carbon dioxide energy storage device, which is used to characterize the charging power or discharging power that needs to be executed under the current operating conditions; the target charging power refers to the electrical energy absorption power that the liquid carbon dioxide energy storage device needs to achieve during the charging process; the target discharging power refers to the electrical energy output power that the liquid carbon dioxide energy storage device needs to achieve during the discharging process. Specifically, power command messages are received from the power grid dispatch terminal according to the current operating conditions. The power command messages are parsed to obtain the power type field and the power value field. Data with the power type field of charging is determined as the target charging power, and data with the power type field of discharging is determined as the target discharging power. The target charging power or target discharging power is bound with the current operating conditions, the time of command reception, and the device number to form a power grid power command record.

[0043] Based on the fixed channel family, transition channel family, replacement channel family and grid power command, several candidate bearer sequences are generated.

[0044] The candidate load-bearing sequence refers to multiple load-bearing schemes for channel families generated based on the fixed channel family, transition channel family, replacement channel family, and power grid command. Each candidate load-bearing sequence includes the arrangement order of each channel family and the corresponding power load-bearing share, which is used to represent the distribution method of working fluid energy input among each channel family under different operating conditions.

[0045] Specifically, the fixed channel family, transition channel family, and succession channel family are regarded as three power carrying objects. The carrying object arrangement results are generated according to the different arrangement order of the three power carrying objects. According to the target charging power or target discharging power in the power grid power instruction, the power value is allocated to the fixed channel family, transition channel family, and succession channel family in each carrying object arrangement result, forming a power carrying share combination corresponding to each carrying object arrangement result. The carrying object arrangement result, the corresponding power carrying share combination, the power type field, and the power value field are bound together to generate a candidate carrying sequence. All candidate carrying sequences are summarized into a candidate carrying sequence set.

[0046] S3. Analyze the target load sequence corresponding to the candidate load sequence to obtain the target power load share; In an embodiment of the present invention, parsing the target bearer sequence corresponding to the candidate bearer sequence to obtain the target power bearer share includes: The candidate carrier sequence is iteratively updated using the particle swarm optimization algorithm to obtain the target carrier sequence. In an embodiment of the present invention, the candidate carrier sequence is iteratively updated using a particle swarm optimization algorithm to obtain the target carrier sequence, including: Each candidate bearer sequence is treated as a particle in the particle swarm algorithm. A particle refers to an independent computational object formed by mapping a candidate bearer sequence. It is used to represent a channel family power bearer mode and the corresponding bearer sequence. Each particle corresponds to an executable power allocation state. Specifically, a particle number is established for each candidate carrier sequence in the candidate carrier sequence set. The channel family arrangement order, the power carrier share corresponding to the fixed channel family, the power carrier share corresponding to the transition channel family, and the power carrier share corresponding to the successor channel family in the candidate carrier sequence are written into the same particle record. This particle record is used as a particle in the particle swarm algorithm.

[0047] The initial velocity of the particle is determined based on the power carrying share in the candidate carrying sequence; Specifically, the power carrying share corresponding to the fixed channel family, the power carrying share corresponding to the transition channel family, and the power carrying share corresponding to the successor channel family are extracted from the same candidate carrying sequence. The power carrying share corresponding to the fixed channel family is subtracted from the power carrying share corresponding to the transition channel family to obtain the first share change. The power carrying share corresponding to the transition channel family is subtracted from the power carrying share corresponding to the successor channel family to obtain the second share change. The power carrying share corresponding to the successor channel family is subtracted from the power carrying share corresponding to the fixed channel family to obtain the third share change. The first share change, the second share change, and the third share change are combined in the order of fixed channel family, transition channel family, and successor channel family to obtain the initial velocity of the particle.

[0048] It should be noted that the power carrying share refers to the proportion of power allocated to a certain channel family from the target charging power or target discharging power corresponding to the power command of the power grid. It is used to indicate the size of the power conversion task undertaken by the fixed channel family, the transition channel family, or the replacement channel family. The larger the power carrying share, the heavier the working fluid flow, heat exchange, and energy conversion task undertaken by the corresponding channel family. The smaller the power carrying share, the lighter the energy conversion task undertaken by the corresponding channel family. The power carrying shares corresponding to each channel family together constitute the power allocation result of the power grid among the channel families, and serve as the basic data for generating candidate carrying sequences, determining the initial velocity of particles, and optimizing the target carrying sequence.

[0049] The initial position of the particle is determined based on the channel family arrangement order and power carrying share in the candidate carrying sequence; Specifically, the channel family arrangement order in the candidate carrying sequence is extracted, and the positions of the fixed channel family, transition channel family and successor channel family in the channel family arrangement order are recorded as channel family position data. The power carrying share corresponding to the fixed channel family, the power carrying share corresponding to the transition channel family and the power carrying share corresponding to the successor channel family in the candidate carrying sequence are extracted. The channel family position data and the power carrying shares corresponding to the three channel families are combined to form the initial position of the particle.

[0050] Based on the initial velocity and initial position of the particles, construct an initial particle swarm; Specifically, the particle number, initial position, and initial velocity of each particle are written into the same particle record, and all particle records are summarized according to particle number to form an initial particle swarm.

[0051] The fitness of particles is calculated based on the candidate carrier sequences; Fitness refers to a comprehensive evaluation value used to evaluate the quality of the candidate load-bearing sequence corresponding to a particle. The higher the fitness, the more the corresponding candidate load-bearing sequence meets the power grid demand and the more it conforms to the power load-bearing relationship between the fixed channel family, the transition channel family, and the replacement channel family. Specifically, the power carrying share of the fixed channel family, the power carrying share of the transition channel family, the power carrying share of the successor channel family, and the power value in the power grid command are extracted from the candidate carrying sequence corresponding to the particle. The power carrying shares of the fixed channel family, the power carrying share of the transition channel family, and the power carrying share of the successor channel family are added together to obtain the total carrying power. The power matching difference is calculated by subtracting the power carrying share of the fixed channel family from the power carrying share of the successor channel family to obtain the selective release difference. The fitness of the particle is obtained by combining the power matching difference and the selective release difference.

[0052] It should be noted that the grid power command determines the energy exchange rate that the liquid carbon dioxide energy storage device should achieve under the current operating conditions, while the power carrying share of each channel family collectively determines the actual energy transfer rate undertaken by the heat exchange channel. According to the law of conservation of energy, the sum of the power carrying shares of each channel family should be consistent with the target power corresponding to the grid power command. Therefore, the difference between the total carrying power and the target power can reflect the degree to which the candidate carrying sequence meets the grid power demand. At the same time, the formation of manifold selectivity indicates that the fixed channel family has long undertaken a high working fluid flow and energy input task, while the replacement channel... The host family undertakes lower working fluid flow and energy input tasks. With the switching of multiple operating conditions, the power load should be gradually released from the stationary channel family to the successor channel family in an orderly manner to avoid the continuous concentration of local energy input. Therefore, the difference between the power load share of the successor channel family and the power load share of the stationary channel family is used to characterize the degree of selective release. The power matching difference and the selective release difference are used together as evaluation criteria. This can reflect both the degree to which the candidate load sequence meets the power demand of the grid and the degree to which the candidate load sequence conforms to the evolution law of manifold selective state. Therefore, the combined result can be used as the fitness of the particle.

[0053] Based on the particle's fitness, determine the individual optimal particle and the global optimal particle; The individual optimal particle refers to the particle state corresponding to the highest fitness of the same particle in the entire iteration process, which is used to represent the best carrying scheme that the particle itself can achieve; the global optimal particle refers to the particle state corresponding to the highest fitness of all particles in the entire iteration process, which is used to represent the carrying scheme with the highest comprehensive evaluation among all candidate carrying sequences. Specifically, all fitness records obtained by the same particle before the current iteration are compared with the fitness of the current iteration. The highest fitness record corresponding to the particle is determined as the individual optimal record of the particle, and the carrier sequence corresponding to the individual optimal record is determined as the individual optimal particle. The individual optimal records of all particles in the initial particle swarm are compared, and the carrier sequence corresponding to the highest fitness record is determined as the global optimal particle.

[0054] Based on the individual optimal particle and the global optimal particle, the initial particle swarm is iteratively updated to obtain the particle swarm after the iteration. The initial particle swarm refers to the set of particles formed by transforming all candidate carrying sequences. Each particle maintains an independent carrying order and power carrying share, and is used as the initial object for iterative optimization of the particle swarm algorithm. The particle swarm after iteration refers to the set of particles formed after multiple rounds of particle position updates and carrying scheme adjustments, in which the carrying scheme of each particle has gradually approached a highly adaptive allocation state. Specifically, for each particle in the initial particle swarm, the current channel family arrangement of the particle is compared with the channel family arrangement of the individual best particle, and the current channel family arrangement of the particle is compared with the channel family arrangement of the globally best particle. Based on the results of the two comparisons, the channel family arrangement of the particle is adjusted. The current power carrying share of the particle is calculated with the power carrying share of the individual best particle, and the current power carrying share of the particle is calculated with the power carrying share of the globally best particle. Based on the results of the two share difference calculations, the power carrying share of the particle is adjusted. The adjusted channel family arrangement and the adjusted power carrying share are combined to form the updated particle. After updating all particles, the particle swarm after the iteration is completed is obtained.

[0055] The target carrying sequence is obtained from the particle swarm after the iteration ends.

[0056] The target load-bearing sequence refers to the candidate load-bearing sequence formed by the particle with the highest fitness in the particle swarm after the iteration, which is used to represent the final power allocation scheme that meets the current power demand of the power grid and conforms to the load-bearing relationship of the channel family.

[0057] Specifically, the fitness of each particle in the particle swarm after the iteration is completed is extracted, the fitness of all particles is compared, the particle with the highest fitness is determined as the target particle, the channel family arrangement order, the fixed channel family power carrying share, the transition channel family power carrying share, and the successor channel family power carrying share of the target particle are extracted, and the channel family arrangement order and the power carrying share of each channel family are combined to form the target carrying sequence.

[0058] The target power carrying capacity sequence is analyzed to obtain the target power carrying capacity share corresponding to the fixed channel family, the transition channel family, and the replacement channel family, respectively.

[0059] Specifically, the target load sequence is used as the data source. The channel family identifier, channel family ranking, and power load share are extracted from the target load sequence. The power load share corresponding to the record identified as a fixed channel family is determined as the target power load share of the fixed channel family. The power load share corresponding to the record identified as a transitional channel family is determined as the target power load share of the transitional channel family. The power load share corresponding to the record identified as a successor channel family is determined as the target power load share of the successor channel family. The target power load shares of the fixed channel family, the transitional channel family, and the successor channel family are combined to form the target power load share.

[0060] S4. Process the target power carrying share to obtain the script release coefficient; based on the script release coefficient and the grid power command, obtain the grid-connected power token; In an embodiment of the present invention, the target power carrying share is processed to obtain the script release coefficient, including: The target power carrying capacity shares corresponding to the fixed channel family, the transition channel family, and the replacement channel family are normalized to obtain the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the replacement channel family. Normalized load share refers to the dimensionless ratio value obtained by proportionally converting the target power load share of each channel family according to the sum of the target power load shares of all channel families. It is used to eliminate the influence of changes in the target power on the power distribution relationship of each channel family, so that only the relative load relationship is retained between each channel family. Specifically, the target power carrying capacity share corresponding to the fixed channel family, the target power carrying capacity share corresponding to the transition channel family, and the target power carrying capacity share corresponding to the replacement channel family are added together to obtain the total target power carrying capacity share. The target power carrying capacity share corresponding to the fixed channel family is divided by the total target power carrying capacity share to obtain the normalized carrying capacity share corresponding to the fixed channel family. The target power carrying capacity share corresponding to the transition channel family is divided by the total target power carrying capacity share to obtain the normalized carrying capacity share corresponding to the transition channel family. The target power carrying capacity share corresponding to the replacement channel family is divided by the total target power carrying capacity share to obtain the normalized carrying capacity share corresponding to the replacement channel family.

[0061] The total carrying capacity of a channel family is calculated based on the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the successor channel family. The total load share of the channel family refers to the sum of the normalized load share corresponding to the fixed channel family, the normalized load share corresponding to the transition channel family, and the normalized load share corresponding to the replacement channel family. It is used to characterize the overall load level after all channel families jointly bear the target power, and serves as a unified benchmark for subsequent share folding processing. The normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the successor channel family are summed to obtain the total carrying capacity share of the channel family.

[0062] Based on the total share of the channel family, the normalized share of the fixed channel family, the transition channel family, and the successor channel family is folded to obtain the script release coefficient.

[0063] Share folding refers to the process of compressing and merging the normalized load shares of fixed channel families, transition channel families, and successor channel families into a single release amount based on their relative load relationships, using the total load share of the channel family as a benchmark. This transforms the power load relationship of multiple channel families into a single result that can characterize the overall load status. The script release coefficient is a single coefficient obtained after share folding, used to characterize the overall power release degree corresponding to the current target load sequence and reflect the comprehensive release status after the combined effect of the power load relationships among fixed channel families, transition channel families, and successor channel families. This provides a unified power adjustment basis for generating grid-connected power tokens.

[0064] Specifically, the normalized load share corresponding to the successor channel family and the normalized load share corresponding to the transition channel family are summed to obtain the release-side load share. The normalized load share corresponding to the fixed channel family is taken as the fixed-side load share. The load release difference is obtained by subtracting the fixed-side load share from the release-side load share. The load release difference is divided by the total load share of the channel family to obtain the script release coefficient.

[0065] It should be noted that after the target power carrying share is normalized, only the relative power distribution relationship between each channel family is retained, and their sum remains constant. Therefore, the normalized carrying share can reflect the relative proportion of energy transfer tasks undertaken by different channel families. According to the energy redistribution law, when the normalized carrying share undertaken by the successor channel family and the transition channel family gradually increases, it indicates that the energy transfer task originally concentrated in the fixed channel family is being transferred to other channel families. The carrying state of liquid carbon dioxide working fluid in the heat exchange channel is gradually changing from concentrated carrying to dispersed carrying. Conversely, it indicates that the energy is still concentrated in the fixed channel family. Therefore, the normalized carrying share corresponding to the successor channel family and the transition channel family is taken as the release-side carrying share, and the normalized carrying share corresponding to the fixed channel family is taken as the fixed-side carrying share. The difference between the two can characterize the degree of energy carrying from the fixed state to the release state. Then, the total carrying share of the channel family is proportionalized to eliminate the influence of the target power magnitude change on the calculation results, so that the result only reflects the release degree of the carrying relationship. Therefore, this proportional result can be used as the script release coefficient to represent the overall power release state corresponding to the target carrying sequence.

[0066] In an embodiment of the present invention, a grid-connected power token is obtained based on the script release coefficient and the grid power command, including: Extract the power value and power direction identifier from the power grid command; Power value refers to the data in the power grid command that indicates the charging power or discharging power, used to characterize the amount of electrical power that the liquid carbon dioxide energy storage process needs to absorb or output under the current operating conditions; power direction indicator refers to the data in the power grid command used to distinguish the direction of power flow, used to characterize whether the liquid carbon dioxide energy storage process is currently performing charging to absorb electrical energy or discharging to output electrical energy. Specifically, the power grid power command is used as the data source. The power field and direction field in the power grid power command are parsed. The data corresponding to the power field is determined as the power value, and the data corresponding to the direction field is determined as the power direction identifier. Corresponding records are established for the power value, the power direction identifier and the power grid power command number.

[0067] Generate a token power value based on the script release coefficient and power value; The token power value refers to the power execution value obtained after adjusting the power value according to the script release coefficient, which is used to characterize the grid-connected power after the channel family bearing relationship is corrected; Specifically, the script release coefficient is used as the power correction ratio, the power value is used as the power benchmark to be executed, the script release coefficient is multiplied by the power value to obtain the token power value, and a corresponding record is established between the token power value, the script release coefficient, and the power value.

[0068] By associating the token power value with the power direction identifier, a grid-connected power token is obtained.

[0069] Grid-connected power tokens refer to station-level power execution data formed by associating token power values ​​and power direction identifiers. They are used to simultaneously express the grid-connected power magnitude and grid-connected power flow direction, and serve as input data for subsequent grid-connected charging and discharging power execution.

[0070] Specifically, the token power value is written into the power value field of the grid-connected power token, the power direction identifier is written into the direction field of the grid-connected power token, and the grid power instruction number is written into the source field of the grid-connected power token, forming a grid-connected power token containing the power value field, the direction field, and the source field.

[0071] S5. Based on the grid-connected power token and the grid-connected power interface data of the liquid carbon dioxide energy storage system, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

[0072] In embodiments of the present invention, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled, including: Acquire grid-connected power interface data of the liquid carbon dioxide energy storage system; wherein, the grid-connected power interface data includes charging power interface data and discharging power interface data; Grid-connected power interface data refers to the interface data corresponding to the power interaction between the liquid carbon dioxide energy storage device and the power grid. It is used to carry the execution value of charging power or discharging power and serves as the direct input data for the device to perform grid-connected power regulation. Specifically, the power write interface between the liquid carbon dioxide energy storage device and the grid power control unit is determined. The charging power interface field and the discharging power interface field in the power write interface are called. The data currently recorded in the charging power interface field is determined as the charging power interface data, and the data currently recorded in the discharging power interface field is determined as the discharging power interface data. The charging power interface data, the discharging power interface data and the current operating conditions are bound together to form the grid-connected power interface data.

[0073] When the power direction identifier in the grid-connected power token is the charging direction identifier, the token power value in the grid-connected power token is written into the charging power interface data. Charging power interface data refers to the interface data used to receive charging direction power values, which represents the target power corresponding to the liquid carbon dioxide energy storage device when absorbing electrical energy from the grid. Specifically, when writing the token power value in the grid-connected power token into the charging power interface data, the power direction identifier and token power value in the grid-connected power token are parsed, the power direction identifier is confirmed to correspond to the charging direction, the token power value is written as the charging power execution value into the charging power interface field, the discharge power interface field is kept in an unwritten state, and the written charging power interface data is used as the grid-connected power interface data under the current operating conditions.

[0074] When the power direction identifier in the grid-connected power token is the discharge direction identifier, the token power value in the grid-connected power token is written into the discharge power interface data. Discharge power interface data refers to the interface data used to receive the power value in the discharge direction, which represents the target power corresponding to the liquid carbon dioxide energy storage device outputting electrical energy to the grid. Specifically, the power direction identifier and token power value in the grid-connected power token are parsed to confirm the discharge direction corresponding to the power direction identifier. The token power value is written as the discharge power execution value into the discharge power interface field, while the charging power interface field remains unwritten. The written discharge power interface data is then used as the grid-connected power interface data under the current operating conditions.

[0075] Based on the grid-connected power interface data after the power value of the token is written, grid-connected power execution data is generated; Grid-connected power execution data refers to the actual power execution data generated based on charging power interface data or discharging power interface data. It is used to characterize the grid-connected charging power or grid-connected discharging power that the liquid carbon dioxide energy storage device should currently execute, and directly drives the energy storage device to complete the corresponding energy exchange process. Specifically, the grid-connected power interface data under the current operating condition is used as the data source. The charging power interface data or discharging power interface data that has been written with the token power value is extracted. The written token power value is determined as the grid-connected power execution value under the current operating condition. The corresponding power direction identifier is determined as the grid-connected power execution direction under the current operating condition. The grid-connected power execution value, the grid-connected power execution direction and the current operating condition are bound together to generate grid-connected power execution data.

[0076] Based on the grid-connected power execution data, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

[0077] Multiple operating conditions refer to the operating conditions of a liquid carbon dioxide energy storage device under different operating states, such as charging, discharging, standby, and charge / discharge switching. Each operating state corresponds to different working fluid flow directions, energy transfer directions, and grid-connected power requirements. Grid-connected charging and discharging power refers to the electrical power actually exchanged between the liquid carbon dioxide energy storage device and the power grid. In the charging state, it represents the electrical power absorbed from the power grid, and in the discharging state, it represents the electrical power output to the power grid. It is the direct power quantity for the liquid carbon dioxide energy storage device to realize the energy storage and release process.

[0078] Specifically, the system analyzes the grid-connected power execution data to obtain the grid-connected power execution value, grid-connected power execution direction, and current operating condition. The grid-connected power execution value when the grid-connected power execution direction is charging is sent to the charging power control terminal, and the grid-connected power execution value when the grid-connected power execution direction is discharging is sent to the discharging power control terminal. This allows the liquid carbon dioxide energy storage system to absorb or output electrical energy according to the grid-connected power execution value corresponding to the current operating condition. The system continuously uses the grid-connected power execution data as the basis for power execution under charging, discharging, standby switching, and hot standby switching conditions to stably control the grid-connected charging and discharging power.

[0079] The above description is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any equivalent substitutions or modifications made by those skilled in the art within the scope of the technology disclosed in the present invention, based on the technical solution and inventive concept of the present invention, should be covered within the scope of protection of the present invention.

Claims

1. A multi-condition stability control method for a liquid carbon dioxide energy storage process, characterized in that, Includes the following steps: S1. Based on the regional data of the entrance area, sort the entrance area to obtain the sorting result. S2. Based on the regional ranking list in the occupancy ranking results, determine the primary habitation layer, secondary habitation layer, and free layer, and generate a candidate load-bearing sequence according to the power grid power command corresponding to the liquid carbon dioxide energy storage system. S3. Analyze the target load sequence corresponding to the candidate load sequence to obtain the target power load share; S4. Process the target power carrying share to obtain the script release coefficient; based on the script release coefficient and the grid power command, obtain the grid-connected power token; S5. Based on the grid-connected power token and the grid-connected power interface data of the liquid carbon dioxide energy storage system, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

2. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 1, characterized in that, Based on the regional data of the entrance area, the entrance area is sorted by placeholders, and the placeholder sorting results are obtained, including: Determine the inlet region of the inlet distribution manifold connection in the compact heat exchanger; Acquire regional data for the inlet area; this regional data includes mass flow rate data, temperature data, and pressure data. Based on the temperature and pressure data in the regional data, determine the regional specific enthalpy of the inlet region; Calculate the inlet occupancy based on mass flow rate data and regional specific enthalpy; Based on the number of spaces occupied at the entrance, the entrance area is sorted to obtain the sorting result.

3. The multi-condition stable control method for a liquid carbon dioxide energy storage process according to claim 1, characterized in that, Based on the regional ranking list in the placeholder ranking results, the primary habitation layer, secondary habitation layer, and free layer are determined, including: Retrieve the list of regions sorted from the placeholder sorting results; The total number of entry areas is determined based on the regional sorting list; Based on the total number of entrance areas, the list of entrance areas is divided into the first sorting segment, the second sorting segment, and the third sorting segment. The entrance area in the first sorting section is designated as the primary residential floor; The entrance area in the second sorting segment is designated as the secondary residential layer; The entry region in the third sorting segment is identified as the free layer.

4. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 1, characterized in that, Generate candidate bearer sequences, including: Determine the heat exchange channel corresponding to the inlet area; The heat exchange channels corresponding to the primary habitation layer are classified into the fixed channel family; The heat exchange channels corresponding to the secondary habitation layer are classified into the transition channel family; The heat exchange channels corresponding to the free layer are classified into the replacement channel family; Obtain the grid power command corresponding to the liquid carbon dioxide energy storage system; wherein, the grid power command includes the target charging power or the target discharging power; Based on the fixed channel family, transition channel family, replacement channel family and grid power command, several candidate bearer sequences are generated.

5. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 1, characterized in that, Parse the target load sequence corresponding to the candidate load sequence to obtain the target power load share, including: The candidate carrier sequence is iteratively updated using the particle swarm optimization algorithm to obtain the target carrier sequence. The target power carrying capacity sequence is analyzed to obtain the target power carrying capacity share corresponding to the fixed channel family, the transition channel family, and the replacement channel family, respectively.

6. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 5, characterized in that, The target power carrying share is processed to obtain the script release coefficient, including: The target power carrying capacity shares corresponding to the fixed channel family, the transition channel family, and the replacement channel family are normalized to obtain the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the replacement channel family. The total carrying capacity of a channel family is calculated based on the normalized carrying capacity share corresponding to the fixed channel family, the normalized carrying capacity share corresponding to the transition channel family, and the normalized carrying capacity share corresponding to the successor channel family. Based on the total share of the channel family, the normalized share of the fixed channel family, the transition channel family, and the successor channel family is folded to obtain the script release coefficient.

7. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 1, characterized in that, Based on the script release coefficient and grid power command, the grid-connected power token is obtained, including: Extract the power value and power direction identifier from the power grid command; Generate a token power value based on the script release coefficient and power value; By associating the token power value with the power direction identifier, a grid-connected power token is obtained.

8. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 7, characterized in that, Stable control of grid-connected charging and discharging power of liquid carbon dioxide energy storage systems under multiple operating conditions, including: Acquire grid-connected power interface data of the liquid carbon dioxide energy storage system; wherein, the grid-connected power interface data includes charging power interface data and discharging power interface data; When the power direction identifier in the grid-connected power token is the charging direction identifier, the token power value in the grid-connected power token is written into the charging power interface data. When the power direction identifier in the grid-connected power token is the discharge direction identifier, the token power value in the grid-connected power token is written into the discharge power interface data. Based on the grid-connected power interface data after the power value of the token is written, grid-connected power execution data is generated; Based on the grid-connected power execution data, the grid-connected charging and discharging power of the liquid carbon dioxide energy storage system under multiple operating conditions is stably controlled.

9. The multi-condition stability control method for a liquid carbon dioxide energy storage process according to claim 5, characterized in that, The candidate carrier sequence is iteratively updated using the particle swarm optimization algorithm to obtain the target carrier sequence, including: Each candidate bearer sequence is treated as a particle in the particle swarm algorithm. The initial velocity of the particle is determined based on the power carrying share in the candidate carrying sequence; The initial position of the particle is determined based on the channel family arrangement order and power carrying share in the candidate carrying sequence; Based on the initial velocity and initial position of the particles, construct an initial particle swarm; The fitness of particles is calculated based on the candidate carrier sequences; Based on the particle's fitness, determine the individual optimal particle and the global optimal particle; Based on the individual optimal particle and the global optimal particle, the initial particle swarm is iteratively updated to obtain the particle swarm after the iteration. The target carrying sequence is obtained from the particle swarm after the iteration ends.

Images