Icing signal analysis method and system based on circulation index

Abstract

The invention relates to the technical field of icing precursor signal prediction, in particular to an icing signal analysis method and system based on a circulation index, integrates a modern meteorological analysis tool and a data processing technology, and aims to improve the prediction accuracy and response efficiency of an icing event. The method comprises the following steps: analyzing and calculating a key atmospheric circulation index by comprehensively utilizing fifth-generation global climate atmosphere reanalysis data provided by an European mid-term weather forecast center in combination with an actual icing record of a power transmission line in a specific region; comprise Siberia high-pressure strength indexes, East Asia groove strength indexes and the like, and the possibility and severity of icing can be accurately tracked and predicted through the support of the data; a lead-lag correlation analysis method is applied to further reveal the dynamic relationship between the icing and the atmospheric circulation index, and a scientific basis is provided for making effective icing response measures.

Description

Technical Field

[0001] The present invention relates to the technical field of icing precursor signal forecasting, and in particular to an icing signal analysis method and system based on a circulation index. Background Art

[0002] The study of quantitative analysis of large-scale circulation indices to determine icing precursor signals is an important topic in modern meteorological science, with significant practical application value and scientific significance. Icing poses a serious threat to infrastructure such as power facilities, transportation and buildings, especially in cold seasons, often leading to serious consequences such as power outages, traffic accidents and building damage. By studying large-scale circulation indices, we can gain a deeper understanding of the changes in atmospheric circulation patterns and their association with surface weather conditions, thereby identifying meteorological precursor signals that may lead to icing. This quantitative analysis technique can provide early warning of potential icing risks, provide decision-making support for relevant departments, take preventive measures, and reduce economic losses and casualties. In addition, the use of large-scale circulation indices to identify icing precursor signals can also enrich and improve meteorological prediction models and improve the accuracy and reliability of meteorological forecasts.

[0003] In the background, large-scale circulation indices have been widely used to study the formation mechanism of extreme weather events. Through in-depth analysis of these indices, the key factors affecting icing phenomena can be captured. Therefore, studying this technology can not only enhance the scientific understanding of icing phenomena, but also provide effective early warning means in practical applications, which has important theoretical research and practical application significance. Summary of the invention

[0004] In view of the above problems existing in the prior art, the present invention is proposed.

[0005] Therefore, the technical problem to be solved by the present invention is: to overcome the deficiencies of the above-mentioned prior art and to provide a technical method for determining ice precursor signals based on quantitative analysis of large-scale circulation indices, so as to accurately improve the accuracy of ice prediction and identification in plateau areas.

[0006] To solve the above technical problems, the present invention provides the following technical solutions: an icing signal analysis method based on the circulation index, comprising: collecting and analyzing subsequent data, extracting key meteorological parameters as a data set; obtaining icing data records of power transmission lines in a specific area; classifying the processed data set; calculating the atmospheric circulation index; and using lead-lag correlation analysis to determine the dynamic relationship between icing and the atmospheric circulation index.

[0007] As a preferred solution of the ice cover signal analysis method based on circulation index described in the present invention, the key meteorological parameters include potential, relative humidity, specific humidity, temperature, wind u component, wind v component and pressure level data set of vertical velocity.

[0008] As a preferred embodiment of the ice signal analysis method based on circulation index described in the present invention, the acquisition of ice data records of transmission lines in a specific area includes collecting ice data of typical transmission line towers in mountainous areas of a specific province, and the geographic coordinate range is 103.33° to 104.34° east longitude and 25.74° to 26.01° north latitude.

[0009] As a preferred solution of the ice cover signal analysis method based on the circulation index described in the present invention, the classification of the processed data set includes selecting data in a specific time period as a training set and a test set.

[0010] As a preferred scheme of the ice cover signal analysis method based on circulation index described in the present invention, wherein: the atmospheric circulation index includes the Siberian high pressure intensity index, the East Asian trough intensity index, the subtropical high pressure ridge point, the subtropical high pressure area index, and the 850hPa subtropical high pressure index.

[0011] As a preferred embodiment of the ice cover signal analysis method based on the circulation index described in the present invention, the lead-lag correlation analysis includes a statistical method for analyzing the dynamic relationship between two time series data, and determines the leading and lagging relationship by calculating the correlation coefficient of one time series with respect to another time series at different time lags; specifically, it includes selecting the time series data to be analyzed, calculating the correlation coefficient of the lag period, and determining the lag period; if the lag period is a positive value, it indicates that one sequence leads the other sequence; if it is a negative value, it indicates that one sequence lags behind the other sequence.

[0012] As a preferred solution of the ice signal analysis method based on circulation index described in the present invention, the dynamic relationship between the ice and the atmospheric circulation index includes studying the dynamic interaction between the ice covering the transmission line tower and the index in the atmospheric environment; the dynamic relationship includes the change of moisture state with meteorological conditions, and the influence of the thickness and density of the ice.

[0013] Another object of the present invention is to provide an ice signal analysis system based on the circulation index, aiming to improve the accuracy and efficiency of ice prediction. By integrating advanced data collection, processing and analysis technologies, the system can effectively collect key meteorological parameters and ice data records of transmission lines in specific areas. By classifying and processing these data and combining them with the calculation of the atmospheric circulation index, the system can make a more accurate assessment of the ice risk. In addition, the system uses the lead-lag correlation analysis method to explore the complex dynamic relationship between ice and the atmospheric circulation index, so as to predict and identify potential ice events. The design of the system includes several key modules: data collection and analysis module, ice data acquisition module, data processing and classification module, circulation index calculation module and dynamic relationship analysis module, each of which focuses on performing specific functions and works together to provide a comprehensive solution. Through this comprehensive analysis, the system not only improves the accuracy of ice prediction, but also provides a powerful decision support tool for the power industry and related departments, which helps to reduce economic losses and safety risks under extreme weather conditions.

[0014] In order to solve the above technical problems, the present invention provides the following technical solutions: an ice signal analysis system based on circulation index, comprising: a data collection and analysis module, an ice data acquisition module, a data processing and classification module, a circulation index calculation module and a dynamic relationship analysis module;

[0015] The data collection and analysis module collects and analyzes the subsequent data and extracts key meteorological parameters as a data set;

[0016] The ice coverage data acquisition module acquires ice coverage data records of power transmission lines in a specific area;

[0017] The data processing classification module classifies the processed data set;

[0018] The circulation index calculation module calculates the atmospheric circulation index;

[0019] The dynamic relationship analysis module determines the dynamic relationship between ice cover and atmospheric circulation index using lead-lag correlation analysis.

[0020] A computer device comprises a memory and a processor, wherein the memory stores a computer program, and wherein the processor implements the steps of the ice cover signal analysis method based on the circulation index as described above when executing the computer program.

[0021] A computer-readable storage medium having a computer program stored thereon, characterized in that when the computer program is executed by a processor, the steps of the ice cover signal analysis method based on the circulation index as described above are implemented.

[0022] The beneficial effects of the present invention are as follows: by obtaining the actual European Centre for Medium-Range Weather Forecasts (ECMWF) fifth-generation global climate atmospheric reanalysis (ERA5) data and transmission line tower icing data in Qujing City, Yunnan Province, and through a series of data processing and lead-lag correlation analysis, the accuracy of icing prediction and identification in plateau areas can be more accurately obtained, providing scientific and technological support for icing occurrence forecasting, which has important scientific significance and application value. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the accompanying drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other accompanying drawings can be obtained based on these accompanying drawings without paying creative work.

[0024] Figure 1 A flow chart of a wind turbine icing inversion and prediction technology for an icing signal analysis method based on a circulation index provided in one embodiment of the present invention.

[0025] Figure 2 A probability distribution diagram of maximum ice coverage according to an ice coverage signal analysis method based on circulation index provided in one embodiment of the present invention.

[0026] Figure 3 A standardized schematic diagram of the daily average value of the atmospheric circulation index selected during the icing process of an icing signal analysis method based on the circulation index provided by one embodiment of the present invention.

[0027] Figure 4 A schematic diagram of the leading and lagging correlation coefficients of the average daily variation of each circulation index and the maximum daily coverage thickness in an ice cover signal analysis method based on circulation index provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0028] In order to make the above-mentioned purposes, features and advantages of the present invention more obvious and easy to understand, the specific implementation methods of the present invention are described in detail below in conjunction with the drawings of the specification. Obviously, the described embodiments are part of the embodiments of the present invention, but not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary persons in the art without creative work should fall within the scope of protection of the present invention.

[0029] Example 1, reference Figure 1-Figure 2 , is an embodiment of the present invention, which provides an ice cover signal analysis method based on the circulation index, including:

[0030] S1: Collect and analyze subsequent data and extract key meteorological parameters as a data set.

[0031] It should be noted that if Figure 1 As shown in S101, the key meteorological parameters include potential, relative humidity, specific humidity, temperature, wind u component, wind v component and pressure level data set of vertical velocity.

[0032] Furthermore, data from the fifth generation of the European Centre for Medium-Range Weather Forecasts (ECMWF) Global Climate Atmospheric Reanalysis (ERA5) from January 5 to January 14, 2021 were collected;

[0033] Furthermore, by obtaining the actual European Center for Medium-Range Weather Forecasts (ECMWF) fifth-generation global climate atmospheric reanalysis (ERA5) data and transmission line tower icing data in Qujing City, Yunnan Province, and through a series of data processing and lead-lag correlation analysis, the accuracy of icing prediction and identification in plateau areas can be more accurately determined, providing scientific and technological support for icing occurrence forecasting, which has important scientific significance and application value; a single-layer data set containing the mean sea level pressure was used, and the pressure level data set ranged from 1000hPa to 400hPa, with a total of 18 layers; these data sets cover the range of 40-160°E and 10-70°N, with a spatial resolution of 0.25°×0.25° and a temporal resolution of 1 hour.

[0034] S2: Obtain ice coverage data records of transmission lines in a specific area.

[0035] It should be noted that if Figure 1 As shown in S102, obtaining ice coverage data records of power transmission lines in a specific region includes collecting typical transmission line tower ice coverage data in a mountainous area of ​​a specific province, with a geographic coordinate range of 103.33° to 104.34° east longitude and 25.74° to 26.01° north latitude.

[0036] Furthermore, ice coverage data of typical transmission line towers located in the mountainous areas of Fuyuan, Zhanyi and Huize in the east of Qujing City, Yunnan Province were collected.

[0037] Furthermore, we collected ice coverage data of typical transmission line towers in the Fuyuan, Zhanyi and Huize mountainous areas in the east of Qujing City, Yunnan Province. These towers have experienced severe and rapid ice coverage for many years; their geographical coordinates range from approximately 103.33°E to 104.34°E, and 25.74°N to 26.01°N; from January 7 to 13, 2021, affected by the southern trough and cold air, Qujing experienced a significant drop in temperature. This caused continuous icing of UHV transmission lines in mountainous areas such as Fuyuan, Zhanyi and Huize; the severe icing areas are mainly concentrated on the Kunliulong line, and the probability distribution of maximum icing is as follows: Figure 2 As shown, the horizontal axis Maximum Ice Thickness represents the maximum ice thickness, and the vertical axis Probability represents the probability.

[0038] S3: Classify the processed data set.

[0039] It should be noted that if Figure 1 As shown in S103 , classifying the processed data set includes selecting data in a specific time period as a training set and a test set.

[0040] Furthermore, the processed data set is classified, and the data from 16:30 on January 21, 2024 to 23:50 on February 1, 2024 are selected as the training set, and the data from 00:10 on February 2, 2024 to 20:00 on February 10, 2024 are selected as the test set.

[0041] Example 2, reference Figure 1 , Figure 3-Figure 4 , is an embodiment of the present invention, which provides an ice cover signal analysis method based on the circulation index, including:

[0042] S4: Calculate the atmospheric circulation index.

[0043] It should be noted that if Figure 1 As shown in S104, the atmospheric circulation index includes the Siberian high pressure intensity index, the East Asian trough intensity index, the subtropical high pressure ridge point, the subtropical high pressure area index, and the 850hPa subtropical high pressure index.

[0044] Furthermore, six atmospheric circulation indices were calculated, including the Siberian High Pressure Intensity Index (Siberian High Pressure System), the East Asian Trough Intensity Index, the subtropical high pressure ridge point, the subtropical high pressure area (intensity) index, and the 850hPa subtropical high pressure index (50-110°E).

[0045] Furthermore, the Siberian high intensity index is calculated as the normalized sea level pressure in the region of 80-120°E and 40-65°N; the East Asian trough intensity index is calculated based on the normalized 500hPa height field in the region of 110-145°E and 25-45°N; the subtropical high pressure index in the range of 90-160°E north of 10°N is calculated; the subtropical area index is defined as the number of grid points with 500hPa potential height greater than 5880gpm in the region; the intensity index is the cumulative difference between the grids with potential height >5880gpm and 5870gpm; the subtropical ridge point is determined as the westernmost longitude of the contour line of potential height 5880gpm; the 850hPa subtropical high pressure index represents the average deviation of 850hPa potential height in the range of 50-110°E and 10-30°N.

[0046] S5: Determine the dynamic relationship between ice cover and atmospheric circulation indices using lead-lag correlation analysis.

[0047] It should be noted that if Figure 1As shown in S105, the lead-lag correlation analysis includes a statistical method for analyzing the dynamic relationship between two time series data, and determines the leading and lagging relationship by calculating the correlation coefficient of one time series with respect to another time series at different time lags; specifically, it includes selecting the time series data to be analyzed, calculating the correlation coefficient of the lag period, and determining the lag period; if the lag period is a positive value, it indicates that one sequence leads the other sequence; if it is a negative value, it indicates that one sequence lags the other sequence.

[0048] Furthermore, the dynamic relationship between ice cover and atmospheric circulation indices includes the study of the dynamic interaction between ice cover on transmission line towers and indices in the atmospheric environment; the dynamic relationship includes changes in moisture status with meteorological conditions, and the impact on ice thickness and density.

[0049] Specifically, the lead-lag correlation is used to analyze the relationship between the Siberian High, the East Asian Trough, the 50-110°E subtropical high and the ice cover.

[0050] Furthermore, the lead-lag correlation was used to analyze the relationship between the Siberian High, East Asian Trough, 50-110°E subtropical high and ice cover. The results showed that when the East Asian trough intensity and subtropical high area were 1 day ahead, their correlation coefficients with the maximum ice thickness were the lowest, which were -0.86 and -0.83 respectively; the 850hPa subtropical high intensity index was the lowest when it was 2 days ahead, with a correlation coefficient of -0.82; on the contrary, the subtropical high, Siberian high and subtropical high ridge point had the highest correlation with ice thickness, with correlation coefficients of 0.46, 0.92 and 0.86 respectively; in summary, each indicator can provide a clearer precursor signal for the daily maximum ice thickness when it is 0 to 2 days ahead.

[0051] Specifically, the lead-lag correlation analysis method is a statistical method used to analyze the dynamic relationship between two time series data; by calculating the correlation coefficient of one time series relative to another time series at different time lags, the leading and lagging relationship between them can be determined; the specific process includes selecting the time series data to be analyzed, calculating the correlation coefficient of each possible lag period, and determining the lag period corresponding to the maximum correlation coefficient; if the lag period is a positive value, it means that one sequence leads the other sequence; if it is a negative value, it means that one sequence lags behind the other sequence.

[0052] This method is widely used in economics, meteorology and engineering. For example, in economics, it is used to analyze the relationship between indicators such as consumption and income, unemployment rate and inflation; in meteorology, it is used to analyze the relationship between meteorological variables, such as temperature and precipitation; in engineering, it is used to optimize system design and control strategies. Its advantage is that it can reveal the dynamic relationship between variables and help prediction and decision-making, but it also has limitations, such as being only applicable to linear relationships and relying on the stationarity of time series; therefore, in practical applications, data needs to be properly preprocessed.

[0053] Specific as Figure 3 , Figure 4 As shown:

[0054] Figure 3 Standardization of the daily average of the atmospheric circulation index selected during the icing process (the horizontal axis time points represent 08:00 today to 07:00 the next day); where the vertical axis Index represents the index; the horizontal axis Month / Day represents the month / day; East Asia trough intensity index represents the East Asian trough intensity index; 850hPa subtropical high index represents the 850hPa subtropical high index; Subtropical high area represents the subtropical high area; Subtropical high intensity represents the subtropical high intensity; Siberian High Pressure System represents the Siberian high pressure system; Subtropical ridge point represents the subtropical high ridge point;

[0055] Figure 4 are the leading and lagging correlation coefficients of the average daily variation of each circulation index and the maximum daily coverage thickness; the negative (positive) horizontal axis represents the number of leading (lagging) days, and the vertical axis represents the correlation coefficient; the negative (positive) delay on the horizontal axis represents the leading (lagging) correlation between the daily average circulation index and the daily maximum coverage thickness; among them, (a) East Asia trough intensity index represents the East Asian trough intensity index; (b) 850hPa subtropical high index represents the 850hPa subtropical high index; (c) subtropical high area represents the subtropical high area; (d) subtropical high intensity represents the subtropical high intensity; (e) siberian high pressure system represents the Siberian high pressure system; (f) subtropical ridge point represents the subtropical high pressure ridge point.

[0056] To sum up, the method provided in the embodiment of the present invention can obtain the actual European Center for Medium-Term Weather Forecasts (ECMWF) fifth-generation global climate atmospheric reanalysis (ERA5) data and transmission line tower icing data in Qujing City, Yunnan Province, and through a series of data processing and lead-lag correlation analysis, it can more accurately predict and identify icing in plateau areas, provide scientific and technological support for icing occurrence forecasting, and has important scientific significance and application value.

[0057] The above is a schematic scheme of an ice signal analysis method based on circulation index in this embodiment. It should be noted that the technical scheme of the system of the ice signal analysis method based on circulation index and the technical scheme of the ice signal analysis method based on circulation index belong to the same concept. For the details not described in detail in the technical scheme of the ice signal analysis system based on circulation index in this embodiment, please refer to the description of the technical scheme of the ice signal analysis method based on circulation index.

[0058] Embodiment 3 is an embodiment of the present invention, which provides an ice signal analysis system based on the circulation index, including: a data collection and analysis module, an ice data acquisition module, a data processing and classification module, a circulation index calculation module and a dynamic relationship analysis module;

[0059] The data collection and analysis module collects and analyzes the subsequent data and extracts key meteorological parameters as a data set;

[0060] The ice coverage data acquisition module acquires ice coverage data records of power transmission lines in a specific area;

[0061] The data processing classification module classifies the processed data set;

[0062] The circulation index calculation module calculates the atmospheric circulation index;

[0063] The dynamic relationship analysis module determines the dynamic relationship between ice cover and atmospheric circulation index using lead-lag correlation analysis.

[0064] This embodiment further provides a computing device, which is applicable to the case of an ice cover signal analysis method based on a circulation index, and includes:

[0065] Memory and processor; the memory is used to store computer executable instructions, and the processor is used to execute computer executable instructions to implement the ice cover signal analysis method based on circulation index as proposed in the above embodiment.

[0066] This embodiment also provides a storage medium on which a computer program is stored. When the program is executed by a processor, the ice cover signal analysis method based on the circulation index as proposed in the above embodiment is implemented.

[0067] The storage medium proposed in this embodiment and the ice cover signal analysis method based on circulation index proposed in the above embodiment belong to the same inventive concept. The technical details not described in detail in this embodiment can be referred to the above embodiment, and this embodiment has the same beneficial effects as the above embodiment.

[0068] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium, including several instructions for a computer device (which can be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes: U disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), random access memory (RAM, Random Access Memory), disk or optical disk, etc., which can store program code.

[0069] The logic and / or steps represented in the flowchart or otherwise described herein, for example, can be considered as an ordered list of executable instructions for implementing logical functions, and can be embodied in any computer-readable medium for use by an instruction execution system, device or apparatus (such as a computer-based system, a system including a processor, or other system that can fetch instructions from an instruction execution system, device or apparatus and execute instructions), or in conjunction with such instruction execution systems, devices or apparatuses. For the purposes of this specification, "computer-readable medium" can be any device that can contain, store, communicate, propagate or transmit a program for use by an instruction execution system, device or apparatus, or in conjunction with such instruction execution systems, devices or apparatuses.

[0070] It should be understood that the various parts of the present invention can be implemented by hardware, software, firmware or a combination thereof. In the above-mentioned embodiments, a plurality of steps or methods can be implemented by software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented by hardware, as in another embodiment, it can be implemented by any one of the following technologies known in the art or their combination: a discrete logic circuit having a logic gate circuit for implementing a logic function for a data signal, a dedicated integrated circuit having a suitable combination of logic gate circuits, a programmable gate array (PGA), a field programmable gate array (FPGA), etc.

[0071] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention rather than to limit it. Although the present invention has been described in detail with reference to the preferred embodiments, those skilled in the art should understand that the technical solutions of the present invention may be modified or replaced by equivalents without departing from the spirit and scope of the technical solutions of the present invention, which should all be included in the scope of the claims of the present invention.

Claims

1. An ice cover signal analysis method based on circulation index, characterized by: include: Collect and analyze subsequent data and extract key meteorological parameters as data sets; Obtain ice coverage data records for transmission lines in a specific area; Classify the processed data set; Calculate the atmospheric circulation index; The dynamic relationship between ice cover and atmospheric circulation indices was determined using lead-lag correlation analysis.

2. The ice cover signal analysis method based on circulation index according to claim 1, characterized in that: The key meteorological parameters include geopotential, relative humidity, specific humidity, temperature, wind u component, wind v component and pressure level data set of vertical velocity.

3. The ice cover signal analysis method based on circulation index according to claim 2, characterized in that: The obtaining of ice coverage data records of power transmission lines in a specific region includes collecting ice coverage data of typical power transmission line towers in mountainous areas of a specific province, with geographic coordinates ranging from 103.33° east longitude to 104.34° east longitude and 25.74° to 26.01° north latitude.

4. The ice cover signal analysis method based on circulation index according to claim 3, characterized in that: The classifying of the processed data set includes selecting data in a specific time period as a training set and a test set.

5. The ice cover signal analysis method based on circulation index according to claim 4, characterized in that: The atmospheric circulation index includes the Siberian high pressure intensity index, the East Asian trough intensity index, the subtropical high pressure ridge point, the subtropical high pressure area index, and the 850hPa subtropical high pressure index.

6. The ice cover signal analysis method based on circulation index according to claim 5, characterized in that: The lead-lag correlation analysis includes a statistical method for analyzing the dynamic relationship between two time series data, and determines the leading and lagging relationship by calculating the correlation coefficient of one time series with respect to another time series at different time lags; specifically, it includes selecting the time series data to be analyzed, calculating the correlation coefficient of the lag period, and determining the lag period; if the lag period is a positive value, it means that one sequence leads the other sequence; if it is a negative value, it means that one sequence lags the other sequence.

7. The ice cover signal analysis method based on circulation index according to claim 6, characterized in that: The dynamic relationship between ice cover and atmospheric circulation index includes studying the dynamic interaction between ice cover on transmission line towers and indexes in the atmospheric environment; the dynamic relationship includes changes in moisture state with meteorological conditions, and the thickness and density of ice cover.

8. A system for analyzing ice cover signals based on circulation index according to any one of claims 1 to 7, characterized in that: include: Data collection and analysis module, ice cover data acquisition module, data processing and classification module, circulation index calculation module and dynamic relationship analysis module; The data collection and analysis module collects and analyzes the subsequent data and extracts key meteorological parameters as a data set; The ice coverage data acquisition module acquires ice coverage data records of power transmission lines in a specific area; The data processing classification module classifies the processed data set; The circulation index calculation module calculates the atmospheric circulation index; The dynamic relationship analysis module determines the dynamic relationship between ice cover and atmospheric circulation index using lead-lag correlation analysis.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the ice cover signal analysis method based on the circulation index described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the ice cover signal analysis method based on the circulation index described in any one of claims 1 to 7 are implemented.

Images