Regional power grid risk and power transformer depreciation age correlation analysis method

By constructing an evolutionary model of the change in the service age of transformer failure rate and combining with the sequential Monte Carlo method to simulate the operating status of transformers in the regional power grid, the problem of the failure of the existing technology to effectively evaluate the impact of transformer depreciation life on regional power grid risks is solved. The quantitative assessment of the transformer depreciation period for regional power grid risks has been achieved, providing power grid companies with scientific basis to scientifically formulate transformer depreciation periods.

CN119397785AInactive Publication Date: 2025-02-07ECONOMIC TECH RES INST OF STATE GRID HENAN ELECTRIC POWER +1
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202411500019.2
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2024-10-25
Publication Date
2025-02-07
Estimated Expiration
Not applicable · inactive patent

AI Technical Summary

Technical Problem

The failure of the prior art to effectively evaluate the impact of transformer depreciation years on regional power grid risks has led to the lack of scientific basis for power grid companies when formulating transformer depreciation years.

Method used

Through bathtub curve fitting based on transformer failure data, an evolutionary model of the change in transformer failure rate and service age was constructed, and the operating status of transformer in the regional power grid was simulated by combining the sequential Monte Carlo method to calculate the impact of transformer depreciation years on regional power grid risks. .

Benefits of technology

This method can effectively capture the impact of transformer depreciation years on the reliability of transformers in regional power grids, quantify its impact on regional power grid risks, and provide scientific basis for power grid companies to scientifically formulate transformer depreciation years.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN119397785A_ABST
    Figure CN119397785A_ABST
Patent Text Reader

Abstract

The invention belongs to the technical field of power system risk assessment, and particularly relates to a regional power grid risk and power transformer depreciation age limit correlation analysis method. The influence of the transformer depreciation age limit and the regional power grid operation time on the regional power grid risk is considered at the same time, the influence of the transformer depreciation age limit on the reliability of the transformer in the regional power grid can be effectively captured, and then the influence of the transformer depreciation age limit on the regional power grid risk is quantified; the method fills the blank of the current regional power grid risk assessment method based on the transformer depreciation age limit, can be used as an important basis for a power grid company to scientifically formulate the transformer depreciation age limit, and promotes the power grid enterprise to realize cost reduction and benefit improvement.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention belongs to the technical field of power system risk assessment, and in particular relates to a method for analyzing the correlation between regional power grid risk and the depreciation period of power transformers. Background Art

[0002] The depreciation period of transformers not only affects the reliability level of transformers in regional power grids, but also affects the proportion of over-aged assets in the power grid and the cost accounting of transmission and distribution prices. It is an important factor affecting the scientific management and healthy development of power grid companies. As a key equipment of the power grid, the selection of the depreciation period of transformers directly affects the risk level of the regional power grid. When the depreciation period is small, the reliability of transformers in the regional power grid is high. At this time, the risk of the regional power grid is low, but it will lead to premature replacement of transformers, causing the power grid company to invest in purchase costs in advance. When the depreciation period is high, although it will reduce the purchase cost of the power grid company, it will also reduce the reliability of transformers in the regional power grid, increasing the risk of the regional power grid.

[0003] In summary, there is an urgent need for a correlation analysis method between regional power grid risk and power transformer depreciation life to quantify the impact of transformer depreciation life on regional power grid risk value, and to provide a scientific basis for power grid companies to scientifically formulate transformer depreciation life. Summary of the invention

[0004] The purpose of the present invention is to overcome the shortcomings of the prior art and to provide a method for analyzing the correlation between regional power grid risk and the depreciation period of power transformers, so as to solve the problem that the current power grid risk assessment method does not involve the change of regional power grid risk with the depreciation period of transformers, quantify the impact of the depreciation period of transformers on the risk value of regional power grids, and provide a scientific basis for power grid companies to scientifically formulate the depreciation period of transformers.

[0005] The object of the present invention is achieved by: a method for analyzing the correlation between regional power grid risk and depreciation period of power transformer, which comprises the following steps:

[0006] S1. Based on the transformer fault related information, the bathtub curve is used to fit the transformer fault data, and an evolution model of the transformer failure rate changing with the transformer service age is constructed;

[0007] S2, based on the evolution model of transformer failure rate changing with transformer service age in step S1, determine the transformer failure rate and establish a transformer repairable forced outage model;

[0008] S3, based on the transformer outage model in S2, the sequential Monte Carlo method is used to simulate the operating status of the transformer in the regional power grid, and the regional power grid risk sampling value I is calculated according to the operating status of the transformer in the regional power grid;

[0009] S4. If the number of simulations of the sequential Monte Carlo method does not reach the set value, step S3 is repeated. Otherwise, the average of the regional power grid risk sampling values ​​I obtained by multiple simulations is used as the regional power grid risk sampling value II under the current regional power grid operation time and the depreciation period of the transformer, and step S5 is executed;

[0010] S5. Change the regional power grid operation time, re-execute steps S2 to S4 to calculate the regional power grid risk sampling value II, until all possible values ​​within the regional power grid operation time definition domain are traversed, and the average value of the regional power grid risk sampling values ​​II corresponding to different regional power grid operation times is used as the regional power grid risk value under the current transformer depreciation period;

[0011] S6, changing the depreciation period of the transformer, repeating steps S2 to S5, and calculating the regional power grid risk value corresponding to different depreciation periods within the definition domain of the transformer depreciation period;

[0012] S7. Taking the transformer depreciation period as an independent variable and the regional power grid risk values ​​corresponding to different depreciation periods as a function, a mathematical fitting method is used to fit the data obtained in step S6 to obtain a correlation curve between the regional power grid risk and the depreciation period of the power transformer.

[0013] Furthermore, the step S1 includes:

[0014] determining parameters of a bathtub curve based on transformer fault information;

[0015] Determine the demarcation point T according to the design life of the transformer d , when the transformer service life is less than T d The transformer stable period failure rate model is a constant, and its value is the bathtub curve at T d When the transformer age is greater than T d The failure rate model of the transformer during the wear period is a bathtub curve;

[0016] The four-parameter bathtub curve is fitted based on the transformer fault related information. The four-parameter bathtub curve function is shown as follows:

[0017]

[0018] Where, α, β, η, λ>0; parameter λ is the multiplication factor of the failure rate function; η is the interval parameter; the values ​​of shape parameters α and β jointly determine the shape of the failure rate curve, α∈(0,1), β>0, t∈(0,η); T s The service life of the transformer;

[0019] When the equipment is in service life T s Greater than T d The failure rate model of the transformer wear period is h(T s), the transformer failure rate evolution model is shown as follows:

[0020]

[0021] Where h(t) is the four-parameter bathtub curve function, T s is the transformer service life, T d It is the dividing point between the stable period failure rate curve and the wear period failure rate curve.

[0022] Further, the S2 includes:

[0023] The regional power grid operation time and transformer depreciation life are taken as input variables, and their respective definition domains are determined. The regional power grid risk is taken as the output variable. The regional power grid operation time T y The domain of is shown below:

[0024]

[0025] Among them, T y Indicates the operating life of the regional power grid and the depreciation life of the transformer T z The domain of is shown below:

[0026]

[0027] Among them, T z Indicates the depreciation period of the transformer;

[0028] The transformer age is determined based on the regional power grid operation time and transformer depreciation period, as follows:

[0029] T=T mod T

[0030] syz

[0031] Among them, T s Indicates the transformer service life, T y Indicates the regional power grid operation time, T z It represents the depreciation period of the transformer, and the symbol mod represents the remainder operation;

[0032] The outage model of the transformer is a repairable forced outage model, where μ is the repair rate, which indicates the probability that the transformer switches from a fault state to an operating state. The calculation of the repair rate μ is as follows:

[0033]

[0034] Among them, M TTR is the mean repair time of the transformer.

[0035] Furthermore, the step S3 includes:

[0036] Generate a random number R between (0,1) for each transformer in the regional power grid i , the specific formula for simulating the operating status of each transformer in the regional power grid is shown as follows:

[0037]

[0038] In the formula, s i Indicates that transformer i is in working or failure state, Q i is the transformer failure rate, R i is a uniformly distributed random number;

[0039] Calculate the operating time and repair time of the transformer as shown below,

[0040]

[0041] Among them, f(T s ) is the transformer failure rate function, μ is the probability that the transformer changes from the fault state to the operating state, ζ1, ζ2 are random numbers that obey the uniform distribution between [0,1];

[0042] When a device fails in the power supply path of a load node, the node load loss is calculated, and the regional power grid risk sampling value I is calculated according to the following formula:

[0043]

[0044] Among them, EENS I The regional power grid risk sampling value I, S is the state set where the regional power grid cannot meet the load demand, P i is the power reduction in state i, t i is the power reduction time in state i.

[0045] Furthermore, the step S5 comprises:

[0046] Increase the regional power grid operation time, re-execute steps S2 to S4 and calculate the regional power grid risk sampling value II according to the following formula until all possible values ​​within the regional power grid operation time definition domain are traversed;

[0047]

[0048] Among them, EENS II is the regional power grid risk sampling value II, T step is the number of sequential Monte Carlo simulations;

[0049] The average value of the regional power grid risk value samples II corresponding to the operation time of different regional power grids is used as the regional power grid risk value under the current transformer depreciation period;

[0050]

[0051] Among them, EENS is the regional power grid risk sampling value II, T y,max It is the maximum operating time of the regional power grid.

[0052] Further, the step S7 includes:

[0053] The regional power grid risk under each depreciation period is normalized according to the following formula:

[0054]

[0055] in, represents the normalized regional power grid risk, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period, EENS max It represents the maximum value of the regional power grid risk value corresponding to all depreciation years, EENS ,min It represents the minimum value of the regional power grid risk value corresponding to all depreciation years;

[0056] The normalized regional power grid risk and its corresponding depreciation period are fitted with the risk function based on the two-parameter Weibull distribution using the least squares method.

[0057]

[0058] Among them, EENS represents the regional power grid risk, S(m,n) is the residual sum of squares between the actual value and the fitted value, m, n, c are the risk function coefficients, t s Indicates the depreciation period of the transformer, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period of the transformer, eens i It represents the fitting value of regional power grid risk corresponding to the i-th depreciation period of the transformer.

[0059] Beneficial effects of the invention: The method for analyzing the correlation between regional power grid risk and power transformer depreciation life of the invention simultaneously considers the influence of transformer depreciation life and regional power grid operation time on regional power grid risk, and can effectively capture the influence of transformer depreciation life on transformer reliability in regional power grid, and then quantify the influence of transformer depreciation life on regional power grid risk, filling the gap in the current method for evaluating regional power grid risk based on transformer depreciation life, and can serve as an important basis for power grid companies to scientifically formulate transformer depreciation life, and promote power grid enterprises to reduce costs and increase efficiency. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0061] Figure 1 It is a flow chart of the model for associating transformer depreciation period with regional power grid risk in the present invention;

[0062] Figure 2 It is a schematic diagram of the regional power grid structure in the present invention;

[0063] Figure 3 It is a schematic diagram of depreciation period-regional power grid risk in the present invention. DETAILED DESCRIPTION

[0064] The present invention will be further described below in conjunction with the accompanying drawings.

[0065] It should be understood that the specific embodiments described herein are only used to explain the present invention, and are not used to limit the present invention.

[0066] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0067] like Figure 1-3 As shown, a method for analyzing the correlation between regional power grid risk and depreciation period of power transformers of the present invention comprises the following steps:

[0068] S1. Based on the transformer fault related information, the bathtub curve is used to fit the transformer fault data, and an evolution model of the transformer failure rate changing with the transformer service age is constructed;

[0069] S2, based on the evolution model of transformer failure rate changing with transformer service age in step S1, determine the transformer failure rate and establish a transformer repairable forced outage model;

[0070] S3, based on the transformer outage model in S2, the sequential Monte Carlo method is used to simulate the operating status of the transformer in the regional power grid, and the regional power grid risk sampling value I is calculated according to the operating status of the transformer in the regional power grid;

[0071] S4. If the number of simulations of the sequential Monte Carlo method does not reach the set value, step S3 is repeated. Otherwise, the average of the regional power grid risk sampling values ​​I obtained by multiple simulations is used as the regional power grid risk sampling value II under the current regional power grid operation time and the depreciation period of the transformer, and step S5 is executed;

[0072] S5. Change the regional power grid operation time, re-execute steps S2 to S4 to calculate the regional power grid risk sampling value II, until all possible values ​​within the regional power grid operation time definition domain are traversed, and the average value of the regional power grid risk sampling values ​​II corresponding to different regional power grid operation times is used as the regional power grid risk value under the current transformer depreciation period;

[0073] S6, changing the depreciation period of the transformer, repeating steps S2 to S5, and calculating the regional power grid risk value corresponding to different depreciation periods within the definition domain of the transformer depreciation period;

[0074] S7. Taking the transformer depreciation period as an independent variable and the regional power grid risk values ​​corresponding to different depreciation periods as a function, a mathematical fitting method is used to fit the data obtained in step S6 to obtain a correlation curve between the regional power grid risk and the depreciation period of the power transformer.

[0075] Furthermore, the step S1 includes:

[0076] determining parameters of a bathtub curve based on transformer fault information;

[0077] Determine the demarcation point T according to the design life of the transformer d , when the transformer service life is less than T d The transformer stable period failure rate model is a constant, and its value is the bathtub curve at T d When the transformer age is greater than T d The failure rate model of the transformer during the wear period is a bathtub curve;

[0078] The four-parameter bathtub curve is fitted based on the transformer fault related information. The four-parameter bathtub curve function is shown as follows:

[0079]

[0080] Where, α, β, η, λ>0; parameter λ is the multiplication factor of the failure rate function. Changing the value of λ will keep the shape of the failure rate curve unchanged while the scale of the failure rate axis will change; η is an interval parameter. Changing the value of η will not change the shape of each curve, but only the scale of the coordinate axis will change; the values ​​of shape parameters α and β jointly determine the shape of the failure rate curve, α∈(0,1), β>0, t∈(0,η); T s The service life of the transformer;

[0081] The design life of the transformer is set at 30 years, and the dividing point between the stable period failure rate curve and the wear period failure rate curve is T d When the transformer service age is less than 18 years, the transformer stable period failure rate model is a four-parameter bathtub curve at T d The value of h(18) at the equipment age T s Greater than T d The failure rate model of the transformer wear period is h(T s ), the transformer failure rate evolution model is shown as follows:

[0082]

[0083] Where h(t) is the four-parameter bathtub curve function, T s is the transformer service life, T d It is the dividing point between the stable period failure rate curve and the wear period failure rate curve.

[0084] Further, the S2 includes:

[0085] The regional power grid operation time and transformer depreciation life are taken as input variables, and their respective definition domains are determined. The regional power grid risk is taken as the output variable. The regional power grid operation time T y The domain of is shown below,

[0086]

[0087] Among them, T y Indicates the operating life of the regional power grid, with a value of 1 to 100. The depreciation life of the transformer is T z The domain of is shown below,

[0088]

[0089] Among them, T z Indicates the depreciation period of the transformer, ranging from 15 to 30 years;

[0090] The transformer age is determined based on the regional power grid operation time and transformer depreciation period, as follows:

[0091] T=T mod T

[0092] syz

[0093] Among them, T s Indicates the transformer age, T y Indicates the regional power grid operation time, T z It represents the depreciation period of the transformer, and the symbol mod represents the remainder operation;

[0094] The outage model of the transformer is a repairable forced outage model, where μ is the repair rate, which indicates the probability that the transformer switches from a fault state to an operating state. The calculation of the repair rate μ is as follows:

[0095]

[0096] Among them, M TTR is the average repair time of the transformer, and in this embodiment, the value is 8.

[0097] Furthermore, the step S3 includes:

[0098] Generate a random number R between (0,1) for each transformer in the regional power grid i , the specific formula for simulating the operating status of each transformer in the regional power grid is shown as follows:

[0099]

[0100] In the formula, s i Indicates that transformer i is in working or failure state, Q i is the transformer failure rate, R i is a uniformly distributed random number;

[0101] Calculate the operating time and repair time of the transformer as shown below,

[0102]

[0103] Among them, f(T s ) is the transformer failure rate function, μ is the probability that the transformer changes from the fault state to the operating state, ζ1, ζ2 are random numbers that obey the uniform distribution between [0,1];

[0104] The load node power supply path is determined as follows:

[0105]

[0106] Among them, lp1 indicates that the power supply path that meets the node 1 is transformer 1, transformer 2, line 10, lp2 indicates that the power supply path that meets the node 1 is transformer 1, transformer 3, line 10, lp3 indicates that the power supply path that meets the node 1 is transformer 1, transformer 4, line 11, lp4 indicates that the power supply path that meets the node 1 is transformer 1, transformer 5, line 10, line 12; lp5 indicates that the power supply path that meets the node 1 is transformer 1, transformer 7, line 11, line 13; lp6 indicates that the power supply path that meets the node 1 is transformer 1, transformer 6, line 10, line 12, line 14, lp7 indicates that the power supply path that meets the node 1 is transformer 1, transformer 8, line 11, line 13, line 15;

[0107] When a device fails in the power supply path of a load node, the node load loss is calculated, and the regional power grid risk sampling value I is calculated according to the following formula:

[0108]

[0109] Among them, EENS I The regional power grid risk sampling value I, S is the state set where the regional power grid cannot meet the load demand, P i is the power reduction in state i, t i is the power reduction time in state i.

[0110] Furthermore, the step S5 comprises:

[0111] Increase the regional power grid operation time, re-execute steps S2 to S4 and calculate the regional power grid risk sampling value II according to the following formula until all possible values ​​within the regional power grid operation time definition domain are traversed;

[0112]

[0113] Among them, EENS II is the regional power grid risk sampling value II, T step is the number of sequential Monte Carlo simulations;

[0114] The average value of the regional power grid risk value samples II corresponding to the operation time of different regional power grids is used as the regional power grid risk value under the current transformer depreciation period;

[0115]

[0116] Among them, EENS is the regional power grid risk sampling value II, T y,max It is the maximum operating time of the regional power grid.

[0117] Further, the step S7 includes:

[0118] The regional power grid risk under each depreciation period is normalized according to the following formula:

[0119]

[0120] in, represents the normalized regional power grid risk, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period, EENS max It represents the maximum value of the regional power grid risk value corresponding to all depreciation years, EENS ,min It represents the minimum value of the regional power grid risk value corresponding to all depreciation years;

[0121] The normalized regional power grid risk and its corresponding depreciation period are fitted with the risk function based on the two-parameter Weibull distribution using the least squares method.

[0122]

[0123] Among them, EENS represents the regional power grid risk, S(m,n) is the residual sum of squares between the actual value and the fitted value, m, n, c are the risk function coefficients, t s Indicates the depreciation period of the transformer, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period of the transformer, eens i It represents the fitting value of regional power grid risk corresponding to the i-th depreciation period of the transformer.

[0124] Final depreciation period-regional power grid risk diagram Figure 3 The transformer failure rate evolution model based on the bathtub curve can provide more accurate transformer failure rate information, thereby improving the accuracy and reliability of risk assessment results. The transformer failure rate is calculated based on the depreciation period and equipment age of the transformer, which can better reflect the actual status of the transformer and make the regional power grid risk assessment results more refined.

[0125] In summary, a method for analyzing the correlation between regional power grid risk and power transformer depreciation life of the present invention simultaneously considers the influence of transformer depreciation life and regional power grid operation time on regional power grid risk, and can effectively capture the influence of transformer depreciation life on transformer reliability in regional power grid, and then quantify the influence of transformer depreciation life on regional power grid risk, which fills the gap in the current method of assessing regional power grid risk based on transformer depreciation life, and can serve as an important basis for power grid companies to scientifically formulate transformer depreciation life, and promote power grid enterprises to reduce costs and increase efficiency.

[0126] In the description of this specification, the description with reference to the terms "one embodiment", "example", "specific example", etc. means that the specific features, structures, materials or characteristics described in conjunction with the embodiment or example are included in at least one embodiment or example of the present invention. In this specification, the schematic representation of the above terms does not necessarily refer to the same embodiment or example. Moreover, the specific features, structures, materials or characteristics described can be combined in any one or more embodiments or examples in a suitable manner.

[0127] The above shows and describes the basic principles, main features and advantages of the present invention. Those skilled in the art should understand that the present invention is not limited to the above embodiments, and the above embodiments and descriptions are only for explaining the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention may have various changes and improvements, and these changes and improvements all fall within the scope of the present invention to be protected.

Claims

1. A method for analyzing the correlation between regional power grid risk and depreciation period of power transformers, characterized in that: The following steps are involved: S1. Based on the transformer fault related information, the bathtub curve is used to fit the transformer fault data, and an evolution model of the transformer failure rate changing with the transformer service age is constructed; S2, based on the evolution model of transformer failure rate changing with transformer service age in step S1, determine the transformer failure rate and establish a transformer repairable forced outage model; S3, based on the transformer outage model in S2, the sequential Monte Carlo method is used to simulate the operating status of the transformer in the regional power grid, and the regional power grid risk sampling value I is calculated according to the operating status of the transformer in the regional power grid; S4. If the number of simulations of the sequential Monte Carlo method does not reach the set value, step S3 is repeated. Otherwise, the average of the regional power grid risk sampling values ​​I obtained by multiple simulations is used as the regional power grid risk sampling value II under the current regional power grid operation time and the depreciation period of the transformer, and step S5 is executed; S5. Change the regional power grid operation time, re-execute steps S2 to S4 to calculate the regional power grid risk sampling value II, until all possible values ​​within the regional power grid operation time definition domain are traversed, and the average value of the regional power grid risk sampling values ​​II corresponding to different regional power grid operation times is used as the regional power grid risk value under the current transformer depreciation period; S6, changing the depreciation period of the transformer, repeating steps S2 to S5, and calculating the regional power grid risk value corresponding to different depreciation periods within the definition domain of the transformer depreciation period; S7. Taking the transformer depreciation period as an independent variable and the regional power grid risk values ​​corresponding to different depreciation periods as a function, a mathematical fitting method is used to fit the data obtained in step S6 to obtain a correlation curve between the regional power grid risk and the depreciation period of the power transformer.

2. A method for analyzing the correlation between regional power grid risk and depreciation period of power transformers as claimed in claim 1, characterized in that: The step S1 comprises: determining parameters of a bathtub curve based on transformer fault information; Determine the demarcation point T according to the design life of the transformer d , when the transformer service life is less than T d The transformer stable period failure rate model is a constant, and its value is the bathtub curve at T d When the transformer age is greater than T d The failure rate model of the transformer during the wear period is a bathtub curve; The four-parameter bathtub curve is fitted based on the transformer fault related information. The four-parameter bathtub curve function is shown as follows: Where, α, β, η, λ>0; parameter λ is the multiplication factor of the failure rate function; η is the interval parameter; the values ​​of shape parameters α and β jointly determine the shape of the failure rate curve, α∈(0,1), β>0, t∈(0,η); T s The service life of the transformer; When the equipment is in service life T s Greater than T d The failure rate model of the transformer wear period is h(T s ), the transformer failure rate evolution model is shown as follows: Where h(t) is the four-parameter bathtub curve function, T s is the transformer service life, T d It is the dividing point between the stable period failure rate curve and the wear period failure rate curve.

3. The method for analyzing the correlation between regional power grid risk and depreciation period of power transformers according to claim 1, characterized in that: The S2 includes: The regional power grid operation time and transformer depreciation life are taken as input variables, and their respective definition domains are determined. The regional power grid risk is taken as the output variable. The regional power grid operation time T y The domain of is shown below, Among them, T y Indicates the operating life of the regional power grid and the depreciation life of the transformer T z The domain of is shown below: Among them, T z Indicates the depreciation period of the transformer; The transformer age is determined based on the regional power grid operation time and transformer depreciation period, as follows: T=T mod T syz Among them, T s Indicates the transformer service life, T y Indicates the regional power grid operation time, T z It represents the depreciation period of the transformer, and the symbol mod represents the remainder operation; The outage model of the transformer is a repairable forced outage model, where μ is the repair rate, which indicates the probability that the transformer switches from a fault state to an operating state. The calculation of the repair rate μ is as follows: Among them, M TTR is the mean repair time of the transformer.

4. The method for analyzing the correlation between regional power grid risk and depreciation period of power transformers according to claim 1, characterized in that: The step S3 comprises: Generate a random number R between (0,1) for each transformer in the regional power grid i , the specific formula for simulating the operating status of each transformer in the regional power grid is shown as follows: In the formula, s i Indicates that transformer i is in working or failure state, Q i is the transformer failure rate, R i is a uniformly distributed random number; Calculate the operating time and repair time of the transformer as shown below, Among them, f(T s ) is the transformer failure rate function, μ is the probability of the transformer changing from a fault state to an operating state, ζ1, ζ2 are random numbers that obey a uniform distribution between [0,1]; When a device fails in the power supply path of a load node, the node load loss is calculated, and the regional power grid risk sampling value I is calculated according to the following formula: Among them, EENS I The regional power grid risk sampling value I, S is the state set where the regional power grid cannot meet the load demand, P i is the power reduction in state i, t i is the power reduction time in state i.

5. The method for analyzing the correlation between regional power grid risk and depreciation period of power transformers according to claim 1, characterized in that: The step S5 comprises: Increase the regional power grid operation time, re-execute steps S2 to S4 and calculate the regional power grid risk sampling value II according to the following formula until all possible values ​​within the regional power grid operation time definition domain are traversed; Among them, EENS II is the regional power grid risk sampling value II, T step is the number of sequential Monte Carlo simulations; The average value of the regional power grid risk value samples II corresponding to the operation time of different regional power grids is used as the regional power grid risk value under the current transformer depreciation period; Among them, EENS is the regional power grid risk sampling value II, T y,max It is the maximum operating time of the regional power grid.

6. A method for analyzing the correlation between regional power grid risk and power transformer depreciation period as claimed in claim 1, characterized in that: The step S7 comprises: The regional power grid risk under each depreciation period is normalized according to the following formula: Among them, EENS i * represents the normalized regional power grid risk, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period, EENS max It represents the maximum value of the regional power grid risk value corresponding to all depreciation years, EENS ,min It represents the minimum value of the regional power grid risk value corresponding to all depreciation years; The normalized regional power grid risk and its corresponding depreciation period are fitted with the risk function based on the two-parameter Weibull distribution using the least squares method. Among them, EENS represents the regional power grid risk, S(m,n) is the residual sum of squares between the actual value and the fitted value, m, n, c are the risk function coefficients, t s Indicates the depreciation period of the transformer, EENS i represents the regional power grid risk value corresponding to the i-th depreciation period of the transformer, eens i It represents the fitting value of regional power grid risk corresponding to the i-th depreciation period of the transformer.