Fine modeling method and system for energy storage battery

Abstract

The application discloses a fine modeling method of an energy storage battery, which comprises the following steps: obtaining data information of a target energy storage battery; constructing a first-order equivalent circuit model of the target energy storage battery; performing discretization processing on the constructed first-order equivalent circuit base model, and separating a charging process and a discharging process; calculating a state of charge value of the target energy storage battery; updating ohmic internal resistance, polarization resistance and polarization capacitance parameters in the first-order equivalent circuit model; and completing fine modeling of the target energy storage battery according to the obtained parameter updating result. The application further discloses a system for implementing the fine modeling method of the energy storage battery. The application not only realizes fine modeling of a terminal voltage and a state of charge value of the target energy storage battery, but also has higher reliability and better accuracy.

Description

Technical Field

[0001] This invention belongs to the field of electrical automation, and specifically relates to a refined modeling method and system for energy storage batteries. Background Technology

[0002] With economic and technological development and the improvement of people's living standards, electricity has become an indispensable secondary energy source in people's production and daily life, bringing endless convenience. Therefore, ensuring a stable and reliable supply of electricity has become one of the most important tasks of the power system.

[0003] Currently, with more and more new energy power generation systems being integrated into the power grid, the randomness and intermittency of their output pose significant challenges to the safe and stable operation of the power system. Energy storage systems, with their rapid response capabilities, play a "stabilizing" role in peak shaving, frequency regulation, output smoothing, voltage support, and energy time shifting. Therefore, the safe and stable operation of energy storage systems is of great importance to the power system.

[0004] Energy storage batteries are the fundamental components of energy storage systems, and their modeling plays a crucial role in the modeling and analysis of energy storage systems. Current modeling schemes for energy storage batteries often employ Kalman filtering, particle filtering, data-driven machine learning, and equivalent circuit models. However, Kalman filtering-based modeling relies heavily on the accuracy of the battery model, is sensitive to dynamic parameter changes, and has a high computational load; particle filtering-based schemes are complex to implement, computationally intensive, and difficult to apply in real time; data-driven machine learning schemes require large amounts of training data and have high computational resource demands; and equivalent circuit models often only consider the discharge process, failing to fully characterize the differences in dynamic charging and discharging characteristics, resulting in poor model accuracy. Summary of the Invention

[0005] One of the objectives of this invention is to provide a refined modeling method for energy storage batteries that is highly reliable and accurate.

[0006] The second objective of this invention is to provide a system for implementing a refined modeling method for the energy storage battery.

[0007] The refined modeling method for energy storage batteries provided by this invention includes the following steps:

[0008] S1. Obtain data information of the target energy storage battery;

[0009] S2. Based on the data obtained in step S1, construct a first-order equivalent circuit model of the target energy storage battery;

[0010] S3. Discretize the basic model of the first-order equivalent circuit constructed in step S2, and separate the charging process and the discharging process;

[0011] S4. Calculate the state of charge (SOC) value of the target energy storage battery based on the discretized charge-discharge model obtained in step S3.

[0012] S5. Based on the state of charge value obtained in step S4, update the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model.

[0013] S6. Based on the parameter update results obtained in step S5, complete the refined modeling of the target energy storage battery.

[0014] Step S1, which involves acquiring data information of the target energy storage battery, specifically includes the following steps:

[0015] Acquire the real-time output current, real-time output voltage, state of charge value, sampling time, charging test data, discharging test data, and open-circuit voltage test data of the target energy storage battery.

[0016] Step S2, which involves constructing a first-order equivalent circuit model of the target energy storage battery based on the data obtained in step S1, specifically includes the following steps:

[0017] The following formula is used as the first-order equivalent circuit model of the target energy storage battery:

[0018] In the formula The polarization voltage across the capacitor of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; for The first derivative with respect to time; The polarization resistance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The polarization capacitance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The ohmic internal resistance of the first-order equivalent circuit model of the target energy storage battery; The real-time output current of the target energy storage battery; The open-circuit voltage of the target energy storage battery; The real-time output voltage of the target energy storage battery.

[0019] Step S3 involves discretizing the first-order equivalent circuit model constructed in step S2, specifically including the following steps:

[0020] To facilitate calculation of the model in the discrete time domain, the first-order equivalent circuit model is discretized to obtain the recursive relationship of the polarization voltage:

[0021] In the formula Let be the polarization voltage of the polarization branch at the k-th sampling time; The polarization voltage of the polarization branch at the (k-1)th sampling time; This refers to the model sampling time, i.e., the discretization step size; The polarization resistor in the polarization branch. Let be the polarization capacitance in the polarization branch, and e be the base of the natural exponential function. Let be the battery current at the k-th sampling time, with the battery discharge direction being the positive direction of the current;

[0022] The corresponding battery terminal voltage can be expressed as:

[0023] In the formula Let be the battery terminal voltage at the k-th sampling time; Let be the battery open-circuit voltage at the k-th sampling time; This is the ohmic internal resistance of the battery.

[0024] Step S3, which involves separating the charging and discharging processes, specifically includes the following steps:

[0025] Considering the different dynamic characteristics of the charging and discharging processes of energy storage batteries, the charging and discharging process parameters are further represented separately, thus separating the charging and discharging processes and obtaining the terminal voltage model after charge-discharge separation:

[0026] In the formula This refers to the battery terminal voltage. This is the open-circuit voltage. In practical applications, the ohmic internal resistance, polarization resistance, and polarization capacitance during the charging and discharging processes will use different values, distinguished by the subscripts c and d, to accurately reflect the differences in the dynamic characteristics of charging and discharging.

[0027] Step S4, which involves calculating the state of charge (SOC) value of the target energy storage battery based on the discretized charge-discharge model obtained in step S3, specifically includes the following steps:

[0028] The open-circuit voltage of the target energy storage battery at the k-th sampling time is calculated using the following formula:

[0029] In the formula, Let be the open-circuit voltage of the target energy storage battery at the k-th sampling time; Let be the terminal voltage of the target energy storage battery at the k-th sampling time; Let be the ohmic internal resistance corresponding to the charging process at the k-th sampling time. This represents the battery current corresponding to the charging process at the k-th sampling time. The polarization voltage corresponding to the charging process at the k-th sampling time; Let be the ohmic internal resistance corresponding to the discharge process at the k-th sampling time. This represents the battery current corresponding to the discharge process at the k-th sampling time. The polarization voltage corresponding to the discharge process at the k-th sampling time;

[0030] The state of charge of the target energy storage battery at the k-th sampling time is calculated using the following formula:

[0031] In the formula The state of charge (SOC) value of the target energy storage battery at the k-th sampling time; and These are the state of charge values ​​corresponding to adjacent calibration open-circuit voltages; and These are the adjacent rated open-circuit voltages, respectively.

[0032] The fitting relationship between the open-circuit voltage and the state-of-charge value is expressed as follows:

[0033] In the formula The open-circuit voltage of the target energy storage battery; The state of charge (SOC) value of the target energy storage battery;

[0034] Based on the discretized charge-discharge model obtained in step S3, the open-circuit voltage of the target energy storage battery is first calculated, and then the state of charge value of the target energy storage battery is calculated based on the correspondence between the open-circuit voltage and the state of charge value.

[0035] Step S5, which updates the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model based on the state of charge value obtained in step S4, specifically includes the following steps:

[0036] During the charging process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula:

[0037] In the formula The ohmic internal resistance of the target energy storage battery during the charging process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the charging process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the charging process; The state of charge of the target energy storage battery.

[0038] During the discharge process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula:

[0039] In the formula The ohmic internal resistance of the target energy storage battery during the discharge process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the discharge process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the discharge process.

[0040] This invention also provides a system for implementing the refined modeling method for energy storage batteries, comprising a data acquisition module, a first construction module, a model processing module, a charge calculation module, a second construction module, and a refined modeling module; the data acquisition module, the first construction module, the model processing module, the charge calculation module, the second construction module, and the refined modeling module are connected in series; the data acquisition module is used to acquire data information of the target energy storage battery and upload the data information to the first construction module; the first construction module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received data information and the acquired data information, and upload the data information to the model processing module; the model processing module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received .... The constructed first-order equivalent circuit model is discretized, and the charging and discharging processes are separated. The data is then uploaded to the charge calculation module. The charge calculation module calculates the state of charge (SOC) of the target energy storage battery based on the received data and the resulting discretized charge / discharge model, and uploads this data to the second construction module. The second construction module updates the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model based on the received data and the obtained SOC, and uploads this data to the fine modeling module. The fine modeling module completes the fine modeling of the target energy storage battery based on the received data and the updated parameters.

[0041] The refined modeling method and system for energy storage batteries provided by this invention, by constructing a first-order equivalent circuit model and combining discretization processing, separation of charging and discharging processes, and updating of state-of-charge related parameters, not only achieves refined modeling of the terminal voltage and state-of-charge values ​​of the target energy storage battery, but also has higher reliability and better accuracy. Attached Figure Description

[0042] Figure 1 This is a schematic diagram of the method flow of the present invention.

[0043] Figure 2 This is a schematic diagram of the simulation curve of an embodiment of the method of the present invention.

[0044] Figure 3 This is a schematic diagram of the functional modules of the system of the present invention. Detailed Implementation

[0045] like Figure 1 The diagram shown is a flowchart of the method of the present invention: The refined modeling method for energy storage batteries disclosed in this invention includes...

[0046] S1. Obtain data information of the target energy storage battery;

[0047] Acquire the real-time output current, real-time output voltage, state of charge value, sampling time, charging test data, discharging test data and open circuit voltage test data of the target energy storage battery;

[0048] Among them, real-time output current and real-time output voltage are used to characterize the operating response characteristics of the target energy storage battery; static open-circuit voltage test data are used to establish the correspondence between open-circuit voltage and state of charge value; charging test data and discharging test data are used to identify the ohmic internal resistance, polarization resistance and polarization capacitance parameters in the first-order equivalent circuit model.

[0049] S2. Based on the data obtained in step S1, construct a first-order equivalent circuit model of the target energy storage battery; specifically, this includes the following steps:

[0050] The following formula is used as the first-order equivalent circuit model of the target energy storage battery:

[0051] In the formula The polarization voltage across the capacitor of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; for The first derivative with respect to time; The polarization resistance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The polarization capacitance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The ohmic internal resistance of the first-order equivalent circuit model of the target energy storage battery; The real-time output current of the target energy storage battery; The open-circuit voltage of the target energy storage battery; The real-time output voltage of the target energy storage battery;

[0052] S3. Discretize the basic model of the first-order equivalent circuit constructed in step S2, and separate the charging and discharging processes; specifically including the following steps:

[0053] To facilitate calculation of the model in the discrete time domain, the first-order equivalent circuit model is discretized to obtain the recursive relationship of the polarization voltage:

[0054] In the formula Let be the polarization voltage of the polarization branch at the k-th sampling time; The polarization voltage of the polarization branch at the (k-1)th sampling time; This refers to the model sampling time, i.e., the discretization step size; The polarization resistor in the polarization branch. Let be the polarization capacitance in the polarization branch, and e be the base of the natural exponential function. Let be the battery current at the k-th sampling time, with the battery discharge direction being the positive direction of the current;

[0055] The corresponding battery terminal voltage can be expressed as:

[0056] In the formula Let be the battery terminal voltage at the k-th sampling time; Let be the battery open-circuit voltage at the k-th sampling time; The internal resistance of the battery is ohms.

[0057] Considering the different dynamic characteristics of the charging and discharging processes of energy storage batteries, the charging and discharging process parameters are further represented separately, thus separating the charging and discharging processes and obtaining the terminal voltage model after charge-discharge separation:

[0058] In the formula This refers to the battery terminal voltage. This is the open-circuit voltage. In practical applications, the ohmic internal resistance, polarization resistance, and polarization capacitance during the charging and discharging processes will use different values, distinguished by the subscripts c and d, to accurately reflect the differences in the dynamic characteristics of charging and discharging.

[0059] Step S4, which involves calculating the state of charge (SOC) value of the target energy storage battery based on the discretized charge-discharge model obtained in step S3, specifically includes the following steps:

[0060] The open-circuit voltage of the target energy storage battery at the k-th sampling time is calculated using the following formula:

[0061] In the formula, Let be the open-circuit voltage of the target energy storage battery at the k-th sampling time; Let be the terminal voltage of the target energy storage battery at the k-th sampling time; Let be the ohmic internal resistance corresponding to the charging process at the k-th sampling time. This represents the battery current corresponding to the charging process at the k-th sampling time. The polarization voltage corresponding to the charging process at the k-th sampling time; Let be the ohmic internal resistance corresponding to the discharge process at the k-th sampling time. This represents the battery current corresponding to the discharge process at the k-th sampling time. The polarization voltage corresponding to the discharge process at the k-th sampling time;

[0062] The state of charge of the target energy storage battery at the k-th sampling time is calculated using the following formula:

[0063] In the formula The state of charge (SOC) value of the target energy storage battery at the k-th sampling time; and These are the state of charge values ​​corresponding to adjacent calibration open-circuit voltages; and These are the adjacent rated open-circuit voltages, respectively.

[0064] The fitting relationship between the open-circuit voltage and the state-of-charge value is expressed as follows:

[0065] In the formula The open-circuit voltage of the target energy storage battery; The state of charge (SOC) value of the target energy storage battery;

[0066] Based on the discretized charge-discharge model obtained in step S3, the open-circuit voltage of the target energy storage battery is first calculated, and then the state of charge value of the target energy storage battery is calculated based on the correspondence between the open-circuit voltage and the state of charge value.

[0067] S5. Based on the state-of-charge values ​​obtained in step S4, update the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model; specifically, this includes the following steps:

[0068] During the charging process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula:

[0069] In the formula The ohmic internal resistance of the target energy storage battery during the charging process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the charging process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the charging process; The state of charge of the target energy storage battery.

[0070] During the discharge process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula:

[0071] In the formula The ohmic internal resistance of the target energy storage battery during the discharge process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the discharge process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the discharge process;

[0072] S6. Based on the parameter update results obtained in step S5, complete the refined modeling of the target energy storage battery.

[0073] The method of the present invention will be described below with reference to an embodiment:

[0074] The predicted terminal voltage obtained from the energy storage battery model constructed according to the present invention is compared with the experimentally measured value. The comparison results are as follows: Figure 2 As shown. The data points include "experimental values" and "model values," and the two are highly consistent in overall trend, indicating that the modeling scheme of this invention can well reflect the voltage dynamic response of the battery in actual operation.

[0075] like Figure 3The diagram shows the functional modules of the system of the present invention: The system for implementing the refined modeling method of the energy storage battery disclosed in this invention includes a data acquisition module, a first construction module, a model processing module, a charge calculation module, a second construction module, and a refined modeling module; the data acquisition module, the first construction module, the model processing module, the charge calculation module, the second construction module, and the refined modeling module are connected in series; the data acquisition module is used to acquire data information of the target energy storage battery and upload the data information to the first construction module; the first construction module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received data information and the acquired data information, and upload the data information to the model processing module; the model processing module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received data information and the acquired data information, and upload the data information to the model processing module; the model processing module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received data information and the acquired data information. The received data is used to discretize the basic model of the first-order equivalent circuit, separating the charging and discharging processes, and then uploading the data to the charge calculation module. The charge calculation module calculates the state of charge (SOC) of the target energy storage battery based on the received data and the resulting discretized charge / discharge model, and uploads this data to the second construction module. The second construction module updates the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model based on the received data and the obtained SOC, and uploads this data to the fine modeling module. The fine modeling module completes the fine modeling of the target energy storage battery based on the received data and the updated parameters.

Claims

1. A refined modeling method for energy storage batteries, comprising the following steps: S1. Obtain data information of the target energy storage battery; S2. Based on the data obtained in step S1, construct a first-order equivalent circuit model of the target energy storage battery; S3. Discretize the basic model of the first-order equivalent circuit constructed in step S2, and separate the charging process and the discharging process; S4. Calculate the state of charge (SOC) value of the target energy storage battery based on the discretized charge-discharge model obtained in step S3. S5. Based on the state of charge value obtained in step S4, update the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model. S6. Based on the parameter update results obtained in step S5, complete the refined modeling of the target energy storage battery.

2. The refined modeling method for energy storage batteries according to claim 1, characterized in that... Step S1, which involves acquiring data information of the target energy storage battery, specifically includes the following steps: Acquire the real-time output current, real-time output voltage, state of charge value, sampling time, charging test data, discharging test data, and open-circuit voltage test data of the target energy storage battery.

3. The refined modeling method for energy storage batteries according to claim 2, characterized in that... Step S2, which involves constructing a first-order equivalent circuit model of the target energy storage battery based on the data obtained in step S1, specifically includes the following steps: The following formula is used as the first-order equivalent circuit model of the target energy storage battery: In the formula The polarization voltage across the capacitor of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; for The first derivative with respect to time; The polarization resistance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The polarization capacitance of the polarization branch in the first-order equivalent circuit model of the target energy storage battery; The ohmic internal resistance of the first-order equivalent circuit model of the target energy storage battery; The real-time output current of the target energy storage battery; The open-circuit voltage of the target energy storage battery; The real-time output voltage of the target energy storage battery.

4. The refined modeling method for energy storage batteries according to claim 3, characterized in that... Step S3 involves discretizing the first-order equivalent circuit model constructed in step S2, specifically including the following steps: To facilitate calculation of the model in the discrete time domain, the first-order equivalent circuit model is discretized to obtain the recursive relationship of the polarization voltage: In the formula Let be the polarization voltage of the polarization branch at the k-th sampling time; The polarization voltage of the polarization branch at the (k-1)th sampling time; This refers to the model sampling time, i.e., the discretization step size; The polarization resistor in the polarization branch. Let be the polarization capacitance in the polarization branch, and e be the base of the natural exponential function. Let be the battery current at the k-th sampling time, with the battery discharge direction being the positive direction of the current; The corresponding battery terminal voltage can be expressed as: In the formula Let be the battery terminal voltage at the k-th sampling time; Let be the battery open-circuit voltage at the k-th sampling time; This is the ohmic internal resistance of the battery.

5. The refined modeling method for energy storage batteries according to claim 4, characterized in that... Step S3, which involves separating the charging and discharging processes, specifically includes the following steps: Considering the different dynamic characteristics of the charging and discharging processes of energy storage batteries, the charging and discharging process parameters are further represented separately, thus separating the charging and discharging processes and obtaining the terminal voltage model after charge-discharge separation: In the formula This refers to the battery terminal voltage. This is the open-circuit voltage. In practical applications, the ohmic internal resistance, polarization resistance, and polarization capacitance during the charging and discharging processes will use different values, distinguished by the subscripts c and d, to accurately reflect the differences in the dynamic characteristics of charging and discharging.

6. The refined modeling method for energy storage batteries according to claim 5, characterized in that... Step S4, which involves calculating the state of charge (SOC) value of the target energy storage battery based on the discretized charge-discharge model obtained in step S3, specifically includes the following steps: The open-circuit voltage of the target energy storage battery at the k-th sampling time is calculated using the following formula: In the formula, Let be the open-circuit voltage of the target energy storage battery at the k-th sampling time; Let be the terminal voltage of the target energy storage battery at the k-th sampling time; Let be the ohmic internal resistance corresponding to the charging process at the k-th sampling time. This represents the battery current corresponding to the charging process at the k-th sampling time. The polarization voltage corresponding to the charging process at the k-th sampling time; Let be the ohmic internal resistance corresponding to the discharge process at the k-th sampling time. This represents the battery current corresponding to the discharge process at the k-th sampling time. The polarization voltage corresponding to the discharge process at the k-th sampling time; The state of charge of the target energy storage battery at the k-th sampling time is calculated using the following formula: In the formula The state of charge (SOC) value of the target energy storage battery at the k-th sampling time; and These are the state of charge values ​​corresponding to adjacent calibration open-circuit voltages; and These are the adjacent rated open-circuit voltages, respectively. The fitting relationship between the open-circuit voltage and the state-of-charge value is expressed as follows: In the formula The open-circuit voltage of the target energy storage battery; The state of charge (SOC) value of the target energy storage battery; Based on the discretized charge-discharge model obtained in step S3, the open-circuit voltage of the target energy storage battery is first calculated, and then the state of charge value of the target energy storage battery is calculated based on the correspondence between the open-circuit voltage and the state of charge value.

7. The refined modeling method for energy storage batteries according to claim 6, characterized in that... Step S5, which updates the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model based on the state of charge value obtained in step S4, specifically includes the following steps: During the charging process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula: In the formula The ohmic internal resistance of the target energy storage battery during the charging process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the charging process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the charging process; The state of charge of the target energy storage battery. During the discharge process, the ohmic internal resistance, polarization resistance, and polarization capacitance parameters in the first-order equivalent circuit model are updated using the following formula: In the formula The ohmic internal resistance of the target energy storage battery during the discharge process; The polarization resistance of the first-order polarization branch of the target energy storage battery during the discharge process; The polarization capacitor of the first-order polarization branch of the target energy storage battery during the discharge process.

8. A system for implementing the refined modeling method for energy storage batteries according to any one of claims 1 to 7, characterized in that... It includes a data acquisition module, a first construction module, a model processing module, a charge calculation module, a second construction module, and a fine modeling module; the data acquisition module, the first construction module, the model processing module, the charge calculation module, the second construction module, and the fine modeling module are connected in series; the data acquisition module is used to acquire data information of the target energy storage battery and upload the data information to the first construction module; the first construction module is used to construct a first-order equivalent circuit model of the target energy storage battery based on the received data information and the acquired data information, and upload the data information to the model processing module; The model processing module is used to discretize the basic model of the first-order equivalent circuit based on the received data information, separate the charging process and the discharging process, and upload the data information to the charge calculation module; the charge calculation module is used to calculate the state of charge value of the target energy storage battery based on the received data information and the obtained discretized charging and discharging model, and upload the data information to the second construction module. The second construction module is used to update the ohmic internal resistance, polarization resistance and polarization capacitance parameters in the first-order equivalent circuit model based on the received data information and the obtained state of charge value, and upload the data information to the fine modeling module. The fine modeling module is used to complete the fine modeling of the target energy storage battery by updating the results based on the received data and parameters.

Images