CN111308286B - Method for estimating service life of traction transformer under multi-parameter empowerment combination - Google Patents

Method for estimating service life of traction transformer under multi-parameter empowerment combination Download PDF

Info

Publication number
CN111308286B
CN111308286B CN202010150333.8A CN202010150333A CN111308286B CN 111308286 B CN111308286 B CN 111308286B CN 202010150333 A CN202010150333 A CN 202010150333A CN 111308286 B CN111308286 B CN 111308286B
Authority
CN
China
Prior art keywords
oil
tensile strength
impact time
service life
test
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010150333.8A
Other languages
Chinese (zh)
Other versions
CN111308286A (en
Inventor
郭蕾
代禄喻
张玉琳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN202010150333.8A priority Critical patent/CN111308286B/en
Publication of CN111308286A publication Critical patent/CN111308286A/en
Application granted granted Critical
Publication of CN111308286B publication Critical patent/CN111308286B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
    • G01R31/1227—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials
    • G01R31/1263—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing of components, parts or materials of solid or fluid materials, e.g. insulation films, bulk material; of semiconductors or LV electronic components or parts; of cable, line or wire insulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/12—Testing dielectric strength or breakdown voltage ; Testing or monitoring effectiveness or level of insulation, e.g. of a cable or of an apparatus, for example using partial discharge measurements; Electrostatic testing
    • G01R31/16—Construction of testing vessels; Electrodes therefor
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N2030/022—Column chromatography characterised by the kind of separation mechanism
    • G01N2030/027—Liquid chromatography

Abstract

The invention discloses a method for estimating service life of a traction transformer under multi-parameter empowerment combination. Firstly, an impact load test platform is set up, impact time is set, power frequency dielectric loss factors, tensile strength and furfural content in oil of the oil paper insulation of the traction transformer at different impact time are obtained through a frequency domain dielectric spectrum tester, a tensile strength tester and a liquid chromatograph, then abnormal values of the power frequency dielectric loss factors, the tensile strength and the furfural content in the oil are removed and preprocessed through an abnormal sample detection method, then the weight of each parameter at different impact time is determined through a variation coefficient method, then characteristic parameters of the oil paper insulation service life of the traction transformer are obtained, and finally service life calculation is carried out. According to the method, abnormal data are effectively reduced through an abnormal sample detection method, the reliability of the data is improved, and the influence of a single weight on evaluation is reduced through a method for determining the weight through a coefficient of variation method.

Description

Method for estimating service life of traction transformer under multi-parameter empowerment combination
Technical Field
The invention relates to the field of service life assessment of a traction transformer, in particular to a method for estimating the service life of the traction transformer under multi-parameter empowerment combination.
Background
The traction power supply system is used as an important component of a high-speed railway and an electrified railway, the stable operation of the traction power supply system is an important basis for the development of the electrified railway, in the practical engineering, the operation condition of a traction transformer is complex, and frequent switching caused by severe traction load fluctuation, start-stop and split-phase is often faced, so that the operation condition of the traction transformer is worse than that of a common power transformer, and meanwhile, the internal thermal field of the traction transformer is more concentrated due to the limitation of the size and the weight. Under the long-term action of an impact traction load, the insulation paper in the oil paper insulation system of the traction transformer bears impact mechanical stress, the degradation of the insulation mechanical property of a winding is accelerated, and then a safety accident of the traction transformer is induced, so that the safe operation of a traction power supply system is seriously threatened.
The service life of the traction transformer mainly depends on the state of oil-paper insulation inside the traction transformer, the service life of the traction transformer under impact load is rarely researched currently, evaluation calculation is mainly carried out on the oil-paper insulation state in the traction transformer under single factor, the test result is high in randomness and is not representative, the service life of the traction transformer under impact load cannot be truly reflected, and in order to effectively calculate the service life of the traction transformer under impact load, the invention provides a method for estimating the service life of the traction transformer under multi-parameter empowerment combination.
Disclosure of Invention
In order to effectively calculate the service life of the traction transformer under impact load, the invention provides a method for estimating the service life of the traction transformer under multi-parameter empowerment combination, which comprises the following steps:
the first step is as follows: building platform
The method comprises the steps of building an impact load comprehensive test platform, wherein the test platform mainly comprises a test box body (1), a first partition plate (2), a second partition plate (3), insulating oil (4), insulating paper (5), an electric heating plate (6), an impact load simulation generator (7), a frequency domain dielectric spectrum tester (8), a liquid chromatograph (9), a tensile strength tester (10), a computer (11), a constant temperature box (12), a three-electrode system (13), a valve (14), an oil flow port (15) and a flow pipe (16), wherein the test box body (1) is divided into a first test area (17), a second test area (18) and a third test area (19) by the first partition plate (2) and the second partition plate (3), the insulating paper (5) is clamped in the electric heating plate (6), the insulating oil (4) is arranged in the first test area (17), and the insulating oil (4) is higher than the electric heating plate (6), the incubator (12) and the three-electrode system (13) are placed in the second test area (18), and the tensile strength tester (10) is placed in the third test area (19);
the second step is that: parameter measurement
Starting the impact load simulation generator (7), setting the total impact time to be T days by taking days as units, and taking out the insulating oil (4) and the insulating paper (5) after each day of impact to perform the following tests:
starting the frequency domain dielectric spectrum tester (8), placing the insulating paper (5) impacted every day in a three-electrode system (13), setting the temperature of the constant temperature box (12) to be 25 ℃, testing the power frequency dielectric loss factor P times of the insulating paper (5) after the impact time of each day, marking as tan delta (i, t), and finally storing the test data in the computer (11);
opening the oil flow port (15), opening the valve (14), taking the insulating oil (4) into the liquid chromatograph (9), testing the furfural content in the oil P times after each day of impact, marking as WC (i, t), and finally storing the test data in the computer (11);
starting the tensile strength tester (10), placing the insulating paper (5) impacted every day in the tensile strength tester (10), testing the tensile strength of the insulating paper (5) P times after each day of impact time, marking as TS (i, t), and finally storing the test data in the computer (11);
tan delta (i, t), WC (i, t), TS (i, t), wherein i represents the ith test, and t represents the impact time of the tth day;
the third step: outlier rejection processing
The method comprises the following steps of carrying out abnormal value elimination pretreatment on power frequency dielectric loss factors tan delta (i, t), tensile strength TS (i, t) and furfural content WC (i, t) in oil, which are tested at the same impact time, and comprises the following steps:
forming a data set by using the power frequency dielectric loss factor tan delta (i, t), the tensile strength TS (i, t) and the furfural content WC (i, t) in the oil at the same impact time: x is the number of1,n,t，x2,n,t，x3,n,t，…xi,n,tI denotes the ith test, n denotes the test characteristicsN-1 represents the power frequency dielectric loss factor tan delta (i, t) characteristic, n-2 represents the tensile strength TS (i, t) characteristic, and n-3 represents the furfural content WC (i, t) characteristic in the oil;
calculating the mean value u of each feature at the same impact timen,tSum variance σ2 n,t：
Figure GDA0002788135010000021
Figure GDA0002788135010000022
Removing the part not in un,t±3σ2 n,tData points in the range, and the number of the removed data points of each feature at the same impact time is recorded as k1,t、k2,t、k3,tTake k as max k1,t，k2,t，k3,tF, number k of feature n cullingn,tIf < k, then k-k is randomly eliminatedn,tAnd (3) removing abnormal points to form a new data set according to the values under the characteristic: x'1,n,t，x′2,n,t，x′3,n,t，…x′P-k,n,tRecalculating mean value and variance u of work frequency dielectric loss factor, tensile strength and furfural content in oil at the same impact time of the remaining P-k points of each characteristicn,t′,σ2 n,t', the expression is as follows:
Figure GDA0002788135010000023
Figure GDA0002788135010000024
the fourth step: determining respective parameter weights
After eliminating abnormal values, the variation coefficient method is utilized to carry out power frequency dielectric loss factor and tensile strength on each characteristic at the same impact timeWeight Q of degree and furfural content in oiln,tThe calculation method is as follows:
calculating the coefficient of variation Vn,t：
Figure GDA0002788135010000025
Calculating the weight Qn,t：
Figure GDA0002788135010000026
The fifth step: service life estimation of traction transformer
Firstly, each feature is processed under the same impact time un,t' normalization to obtain normalized un,tAcquiring a characteristic parameter chi of the service life of the traction transformert：
χt＝Q1,t×u1,t″+Q2,t×u2,t″+Q3,t×u3,t″ (7)
The characteristic parameter chi of the service life of the traction transformer under different impact timetFitting with the impact time t, wherein the fitting method comprises the following steps:
setting a fitting polynomial:
H(t)＝a0+a1t+a2t2+a3t3+a4t4 (8)
wherein H (t) is a fitting function, a0、a1、a2、a3、a4Is a fitting coefficient;
calculating the deviation square sum Loss:
Figure GDA0002788135010000031
to fitting coefficient a0、a1、a2、a3、a4Respectively solve forThe partial derivative is 0:
Figure GDA0002788135010000032
Figure GDA0002788135010000033
Figure GDA0002788135010000034
Figure GDA0002788135010000035
Figure GDA0002788135010000036
simplified formulae (10) to (14):
Figure GDA0002788135010000037
Figure GDA0002788135010000038
Figure GDA0002788135010000039
Figure GDA00027881350100000310
Figure GDA00027881350100000311
writing all coefficients containing only t into a first matrix B:
Figure GDA00027881350100000312
write all fitting coefficients alone to the second matrix a:
Figure GDA00027881350100000313
all contain chitIs written as a third matrix Y:
Figure GDA0002788135010000041
solving the matrix equation BA as Y to obtain a fitting coefficient a0、a1、a2、a3、a4And substituting the fitting coefficient into the formula (8) to obtain an expression of calculation of the service Life of the traction transformer on the impact load under the multi-parameter weighting combination, wherein the service Life is expressed as:
Life＝Γ(H(t)) (23)
where Γ (h (t)) represents the inverse function of the fitting polynomial h (t).
According to the method, abnormal data are effectively reduced through an abnormal sample detection method, the reliability of the data is improved, and the influence of a single weight on evaluation is reduced through a method for determining the weight through a coefficient of variation method. Based on the new estimation method for the service life of the traction transformer under the multi-parameter empowerment combination, the service life of the traction transformer can be effectively judged.
Drawings
FIG. 1 is a flow chart of a method for estimating service life of a traction transformer under multi-parameter weighted combination.
FIG. 2 is a schematic diagram of a comprehensive test platform.
The device comprises a test box body (1), a first partition plate (2), a second partition plate (3), insulating oil (4), insulating paper (5), an electric hot plate (6), an impact load simulation generator (7), a frequency domain dielectric spectrum tester (8), a liquid chromatograph (9), a tensile strength tester (10), a computer (11), a constant temperature box (12), a three-electrode system (13), a valve (14), an oil flow port (15) and a flow conveying pipe (16).
Detailed Description
The first step is as follows: building platform
The method comprises the steps of building an impact load comprehensive test platform, wherein the test platform mainly comprises a test box body (1), a first partition plate (2), a second partition plate (3), insulating oil (4), insulating paper (5), an electric heating plate (6), an impact load simulation generator (7), a frequency domain dielectric spectrum tester (8), a liquid chromatograph (9), a tensile strength tester (10), a computer (11), a constant temperature box (12), a three-electrode system (13), a valve (14), an oil flow port (15) and a flow pipe (16), wherein the test box body (1) is divided into a first test area (17), a second test area (18) and a third test area (19) by the first partition plate (2) and the second partition plate (3), the insulating paper (5) is clamped in the electric heating plate (6), the insulating oil (4) is arranged in the first test area (17), and the insulating oil (4) is higher than the electric heating plate (6), the incubator (12) and the three-electrode system (13) are placed in the second test area (18), and the tensile strength tester (10) is placed in the third test area (19);
the second step is that: parameter measurement
Starting the impact load simulation generator (7), setting the total impact time to be T days by taking days as units, and taking out the insulating oil (4) and the insulating paper (5) after each day of impact to perform the following tests:
starting the frequency domain dielectric spectrum tester (8), placing the insulating paper (5) impacted every day in a three-electrode system (13), setting the temperature of the constant temperature box (12) to be 25 ℃, testing the power frequency dielectric loss factor P times of the insulating paper (5) after the impact time of each day, marking as tan delta (i, t), and finally storing the test data in the computer (11);
opening the oil flow port (15), opening the valve (14), taking the insulating oil (4) into the liquid chromatograph (9), testing the furfural content in the oil P times after each day of impact, marking as WC (i, t), and finally storing the test data in the computer (11);
starting the tensile strength tester (10), placing the insulating paper (5) impacted every day in the tensile strength tester (10), testing the tensile strength of the insulating paper (5) P times after each day of impact time, marking as TS (i, t), and finally storing the test data in the computer (11);
tan delta (i, t), WC (i, t), TS (i, t), wherein i represents the ith test, and t represents the impact time of the tth day;
the third step: outlier rejection processing
The method comprises the following steps of carrying out abnormal value elimination pretreatment on power frequency dielectric loss factors tan delta (i, t), tensile strength TS (i, t) and furfural content WC (i, t) in oil, which are tested at the same impact time, and comprises the following steps:
forming a data set by using the power frequency dielectric loss factor tan delta (i, t), the tensile strength TS (i, t) and the furfural content WC (i, t) in the oil at the same impact time: x is the number of1,n,t，x2,n,t，x3,n,t，…xi,n,tI represents the ith test, n represents the test characteristic, n-1 represents the power frequency dielectric loss factor tan delta (i, t) characteristic, n-2 represents the tensile strength TS (i, t) characteristic, and n-3 represents the furfural content WC (i, t) characteristic in the oil;
calculating the mean value u of each feature at the same impact timen,tSum variance σ2 n,t：
Figure GDA0002788135010000051
Figure GDA0002788135010000052
Removing the part not in un,t±3σ2 n,tData points in the range, and the number of the removed data points of each feature at the same impact time is recorded as k1,t、k2,t、k3,tTake k as max k1,t，k2,t，k3,tF, number k of feature n cullingn,tIf < k, then k-k is randomly eliminatedn,tAnd (3) removing abnormal points to form a new data set according to the values under the characteristic: x'1,n,t，x′2,n,t，x′3,n,t，…x′P-k,n,tRecalculating mean value and variance u of work frequency dielectric loss factor, tensile strength and furfural content in oil at the same impact time of the remaining P-k points of each characteristicn,t′,σ2 n,t', the expression is as follows:
Figure GDA0002788135010000053
Figure GDA0002788135010000054
the fourth step: determining respective parameter weights
After removing abnormal numerical values, weighting Q for power frequency dielectric loss factor, tensile strength and furfural content in oil of each characteristic at the same impact time by using a coefficient of variation methodn,tThe calculation method is as follows:
calculating the coefficient of variation Vn,t：
Figure GDA0002788135010000055
Calculating the weight Qn,t：
Figure GDA0002788135010000056
The fifth step: service life estimation of traction transformer
Firstly, each feature is processed under the same impact time un,t' normalization to obtain normalized un,tAcquiring a characteristic parameter chi of the service life of the traction transformert：
χt＝Q1,t×u1,t″+Q2,t×u2,t″+Q3,t×u3,t″ (7)
Pulling down different impact timeCharacteristic parameter chi of service life of transformertFitting with the impact time t, wherein the fitting method comprises the following steps:
setting a fitting polynomial:
H(t)＝a0+a1t+a2t2+a3t3+a4t4 (8)
wherein H (t) is a fitting function, a0、a1、a2、a3、a4Is a fitting coefficient;
calculating the deviation square sum Loss:
Figure GDA0002788135010000061
to fitting coefficient a0、a1、a2、a3、a4Respectively solving the partial derivatives of the two to be 0:
Figure GDA0002788135010000062
Figure GDA0002788135010000063
Figure GDA0002788135010000064
Figure GDA0002788135010000065
Figure GDA0002788135010000066
simplified formulae (10) to (14):
Figure GDA0002788135010000067
Figure GDA0002788135010000068
Figure GDA0002788135010000069
Figure GDA00027881350100000610
Figure GDA00027881350100000611
writing all coefficients containing only t into a first matrix B:
Figure GDA00027881350100000612
write all fitting coefficients alone to the second matrix a:
Figure GDA0002788135010000071
all contain chitIs written as a third matrix Y:
Figure GDA0002788135010000072
solving the matrix equation BA as Y to obtain a fitting coefficient a0、a1、a2、a3、a4And substituting the fitting coefficient into the formula (8) to obtain an expression of calculation of the service Life of the traction transformer on the impact load under the multi-parameter weighting combination, wherein the service Life is expressed as:
Life＝Γ(H(t)) (23)
where Γ (h (t)) represents the inverse function of the fitting polynomial h (t).

Claims (1)

1. A method for estimating service life of a traction transformer under multi-parameter empowerment combination is characterized by comprising the following steps:
the first step is as follows: building platform
The method comprises the steps of building an impact load comprehensive test platform, wherein the test platform mainly comprises a test box body (1), a first partition plate (2), a second partition plate (3), insulating oil (4), insulating paper (5), an electric heating plate (6), an impact load simulation generator (7), a frequency domain dielectric spectrum tester (8), a liquid chromatograph (9), a tensile strength tester (10), a computer (11), a constant temperature box (12), a three-electrode system (13), a valve (14), an oil flow port (15) and a flow pipe (16), wherein the test box body (1) is divided into a first test area (17), a second test area (18) and a third test area (19) by the first partition plate (2) and the second partition plate (3), the insulating paper (5) is clamped in the electric heating plate (6), the insulating oil (4) is arranged in the first test area (17), and the insulating oil (4) is higher than the electric heating plate (6), the incubator (12) and the three-electrode system (13) are placed in the second test area (18), and the tensile strength tester (10) is placed in the third test area (19);
the second step is that: parameter measurement
Starting the impact load simulation generator (7), setting the total impact time to be T days by taking days as units, and taking out the insulating oil (4) and the insulating paper (5) after each day of impact to perform the following tests:
starting the frequency domain dielectric spectrum tester (8), placing the insulating paper (5) impacted every day in a three-electrode system (13), setting the temperature of the constant temperature box (12) to be 25 ℃, testing the power frequency dielectric loss factor P times of the insulating paper (5) after the impact time of each day, marking as tan delta (i, t), and finally storing the test data in the computer (11);
opening the oil flow port (15), opening the valve (14), taking the insulating oil (4) into the liquid chromatograph (9), testing the furfural content in the oil P times after each day of impact, marking as WC (i, t), and finally storing the test data in the computer (11);
starting the tensile strength tester (10), placing the insulating paper (5) impacted every day in the tensile strength tester (10), testing the tensile strength of the insulating paper (5) P times after each day of impact time, marking as TS (i, t), and finally storing the test data in the computer (11);
tan delta (i, t), WC (i, t), TS (i, t), wherein i represents the ith test, and t represents the impact time of the tth day;
the third step: outlier rejection processing
The method comprises the following steps of carrying out abnormal value elimination pretreatment on power frequency dielectric loss factors tan delta (i, t), tensile strength TS (i, t) and furfural content WC (i, t) in oil, which are tested at the same impact time, and comprises the following steps:
forming a data set by using the power frequency dielectric loss factor tan delta (i, t), the tensile strength TS (i, t) and the furfural content WC (i, t) in the oil at the same impact time: x is the number of1,n,t，x2,n,t，x3,n,t，…xi,n,tI represents the ith test, n represents the test characteristic, n-1 represents the power frequency dielectric loss factor tan delta (i, t) characteristic, n-2 represents the tensile strength TS (i, t) characteristic, and n-3 represents the furfural content WC (i, t) characteristic in the oil;
calculating the mean value u of each feature at the same impact timen,tSum variance σ2 n,t：
Figure FDA0002788134000000011
Figure FDA0002788134000000012
Removing the part not in un,t±3σ2 n,tData points in the range, and the number of the removed data points of each feature at the same impact time is recorded as k1,t、k2,t、k3,tTake k as max k1,t，k2,t，k3,tF, number k of feature n cullingn,tIf < k, then k-k is randomly eliminatedn,tAnd (3) removing abnormal points to form a new data set according to the values under the characteristic: x'1,n,t，x′2,n,t，x′3,n,t，…x′P-k,n,tRecalculating mean value and variance u of work frequency dielectric loss factor, tensile strength and furfural content in oil at the same impact time of the remaining P-k points of each characteristicn,t′,σ2 n,t', the expression is as follows:
Figure FDA0002788134000000013
Figure FDA0002788134000000021
the fourth step: determining respective parameter weights
After removing abnormal numerical values, weighting Q for power frequency dielectric loss factor, tensile strength and furfural content in oil of each characteristic at the same impact time by using a coefficient of variation methodn,tThe calculation method is as follows:
calculating the coefficient of variation Vn,t：
Figure FDA0002788134000000022
Calculating the weight Qn,t：
Figure FDA0002788134000000023
The fifth step: service life estimation of traction transformer
Firstly, each feature is processed under the same impact time un,t' normalization to obtain normalized un,tAcquiring a characteristic parameter chi of the service life of the traction transformert：
χt＝Q1,t×u1,t″+Q2,t×u2,t″+Q3,t×u3,t″ (7)
The characteristic parameter chi of the service life of the traction transformer under different impact timetFitting with the impact time t, wherein the fitting method comprises the following steps:
setting a fitting polynomial:
H(t)＝a0+a1t+a2t2+a3t3+a4t4 (8)
wherein H (t) is a fitting function, a0、a1、a2、a3、a4Is a fitting coefficient;
calculating the deviation square sum Loss:
Figure FDA0002788134000000024
to fitting coefficient a0、a1、a2、a3、a4Respectively solving the partial derivatives of the two to be 0:
Figure FDA0002788134000000025
Figure FDA0002788134000000026
Figure FDA0002788134000000027
Figure FDA0002788134000000028
Figure FDA0002788134000000029
simplified formulae (10) to (14):
Figure FDA00027881340000000210
Figure FDA0002788134000000031
Figure FDA0002788134000000032
Figure FDA0002788134000000033
Figure FDA0002788134000000034
writing all coefficients containing only t into a first matrix B:
Figure FDA0002788134000000035
write all fitting coefficients alone to the second matrix a:
Figure FDA0002788134000000036
all contain chitIs written as a third matrix Y:
Figure FDA0002788134000000037
solving the matrix equation BA as Y to obtain a fitting coefficient a0、a1、a2、a3、a4And substituting the fitting coefficient into the formula (8) to obtain an expression of calculation of the service Life of the traction transformer on the impact load under the multi-parameter weighting combination, wherein the service Life is expressed as:
Life＝Γ(H(t)) (23)
where Γ (h (t)) represents the inverse function of the fitting polynomial h (t).
CN202010150333.8A 2020-03-06 2020-03-06 Method for estimating service life of traction transformer under multi-parameter empowerment combination Active CN111308286B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010150333.8A CN111308286B (en) 2020-03-06 2020-03-06 Method for estimating service life of traction transformer under multi-parameter empowerment combination

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010150333.8A CN111308286B (en) 2020-03-06 2020-03-06 Method for estimating service life of traction transformer under multi-parameter empowerment combination

Publications (2)

Publication Number Publication Date
CN111308286A CN111308286A (en) 2020-06-19
CN111308286B true CN111308286B (en) 2021-01-12

Family

ID=71149623

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010150333.8A Active CN111308286B (en) 2020-03-06 2020-03-06 Method for estimating service life of traction transformer under multi-parameter empowerment combination

Country Status (1)

Country Link
CN (1) CN111308286B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117054259B (en) * 2023-09-27 2024-03-22 保定德利盛电器有限公司 Impact-resistant testing system for strength of transformer shell

Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4532043A (en) * 1983-08-15 1985-07-30 The Dow Chemical Company Method for characterizing the molecular weight and molecular weight distribution of ultra-high molecular weight water soluble polymers
JPS6391534A (en) * 1986-10-06 1988-04-22 Kenkichi Murakami Method and apparatus for simply and rapidly determining distribution of polymerization degree between crosslinking points of crosslinked amorphous polymer
KR20030004510A (en) * 2001-07-05 2003-01-15 이동영 Apparatus and method for multi-testing insulation of power cables
CN101887094A (en) * 2010-06-30 2010-11-17 河南省电力公司济源供电公司 Judging method of solid insulating and aging degrees of transformer
CN102735977A (en) * 2012-07-12 2012-10-17 四川大学 Electrical-thermal aging device of cable accessories under high voltage and high current and test method thereof
CN102981108A (en) * 2012-11-29 2013-03-20 重庆大学 Transformer internal insulation aging diagnosis method based on multi-feature information fusion technology
CN103293417A (en) * 2013-05-28 2013-09-11 北京华电云通电力技术有限公司 Detection method of transformer aging degree
JP2014143209A (en) * 2014-03-18 2014-08-07 Panasonic Corp Lighting device, and illuminating fixture and illumination system using the same
CN104408275A (en) * 2014-05-20 2015-03-11 国家电网公司 Transformer reliability evaluation method based on integral health index
KR20160107994A (en) * 2015-03-06 2016-09-19 (주)오은씨엔아이 Multiple waveshape impulse generator for performance-testing signal surge protective device
CN106199234A (en) * 2016-06-22 2016-12-07 国家电网公司 A kind of transformator residual life detection method based on health index
CN107329065A (en) * 2017-08-24 2017-11-07 西南交通大学 A kind of method for assessing tractive transformer sleeve pipe state of insulation in military service
CN107422234A (en) * 2017-08-24 2017-12-01 西南交通大学 A kind of method for assessing tractive transformer humidified insulation state in military service
CN108802584A (en) * 2018-06-29 2018-11-13 广东工业大学 A kind of transformer insulating paper ageing state appraisal procedure
CN110057941A (en) * 2019-05-15 2019-07-26 广西电网有限责任公司南宁供电局 Detection method based on furfural content in live transformer insulation oil
CN110297167A (en) * 2019-07-17 2019-10-01 国网浙江省电力有限公司电力科学研究院 A kind of transformer ageing state evaluation method based on Multi-source Information Fusion

Patent Citations (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4532043A (en) * 1983-08-15 1985-07-30 The Dow Chemical Company Method for characterizing the molecular weight and molecular weight distribution of ultra-high molecular weight water soluble polymers
JPS6391534A (en) * 1986-10-06 1988-04-22 Kenkichi Murakami Method and apparatus for simply and rapidly determining distribution of polymerization degree between crosslinking points of crosslinked amorphous polymer
KR20030004510A (en) * 2001-07-05 2003-01-15 이동영 Apparatus and method for multi-testing insulation of power cables
CN101887094A (en) * 2010-06-30 2010-11-17 河南省电力公司济源供电公司 Judging method of solid insulating and aging degrees of transformer
CN102735977A (en) * 2012-07-12 2012-10-17 四川大学 Electrical-thermal aging device of cable accessories under high voltage and high current and test method thereof
CN102981108A (en) * 2012-11-29 2013-03-20 重庆大学 Transformer internal insulation aging diagnosis method based on multi-feature information fusion technology
CN103293417A (en) * 2013-05-28 2013-09-11 北京华电云通电力技术有限公司 Detection method of transformer aging degree
JP2014143209A (en) * 2014-03-18 2014-08-07 Panasonic Corp Lighting device, and illuminating fixture and illumination system using the same
CN104408275A (en) * 2014-05-20 2015-03-11 国家电网公司 Transformer reliability evaluation method based on integral health index
KR20160107994A (en) * 2015-03-06 2016-09-19 (주)오은씨엔아이 Multiple waveshape impulse generator for performance-testing signal surge protective device
CN106199234A (en) * 2016-06-22 2016-12-07 国家电网公司 A kind of transformator residual life detection method based on health index
CN107329065A (en) * 2017-08-24 2017-11-07 西南交通大学 A kind of method for assessing tractive transformer sleeve pipe state of insulation in military service
CN107422234A (en) * 2017-08-24 2017-12-01 西南交通大学 A kind of method for assessing tractive transformer humidified insulation state in military service
CN108802584A (en) * 2018-06-29 2018-11-13 广东工业大学 A kind of transformer insulating paper ageing state appraisal procedure
CN110057941A (en) * 2019-05-15 2019-07-26 广西电网有限责任公司南宁供电局 Detection method based on furfural content in live transformer insulation oil
CN110297167A (en) * 2019-07-17 2019-10-01 国网浙江省电力有限公司电力科学研究院 A kind of transformer ageing state evaluation method based on Multi-source Information Fusion

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Data Fusion and Machine Learning Integration for Transformer Loss of Life Estimation;Mohsen Mahoor 等;《2018 IEEE/PES Transmission and Distribution Conference and Exposition (T&D)》;20180419;全文 *
基于振荡波的变压器绕组故障诊断方法研究;吴振宇 等;《中国电机工程学报》;20200131;全文 *

Also Published As

Publication number Publication date
CN111308286A (en) 2020-06-19

Similar Documents

Publication Publication Date Title
CN108680613B (en) Method for evaluating moisture content in insulating paper by using initial slope of complex dielectric constant
CN106021756A (en) Method for assessing insulation state of oil paper based on characteristic quantity of frequency domain dielectric spectroscopy
CN103149452A (en) Method for evaluating ageing state of paper oil insulation
CN106338336A (en) Transformer vibration on-line monitoring system
CN107957436B (en) Method for evaluating uneven moisture degree of oil-impregnated paper insulation assembly
CN111308286B (en) Method for estimating service life of traction transformer under multi-parameter empowerment combination
CN112179852B (en) Method and device for predicting residual life of composite insulator
CN202720302U (en) Simulation system for fault phenomena of turn-to-turn short circuit and insulating property reduction of winding
CN106093614A (en) A kind of method of return voltage initial slope assessment transformer insulation state
CN108896161A (en) A kind of Wound iron-core transformer winding failure analog platform and assessment method
CN106291432A (en) The method of inspection of CT saturation degree and system
Bjelić et al. Fault detection in a power transformer based on reverberation time
Liu et al. A novel curve database for moisture evaluation of transformer oil-immersed cellulose insulation using FDS and exponential decay model
CN116520101A (en) Multi-source fusion-based sleeve defect identification and insulation state evaluation method
Maur et al. An Approach for Aging State Estimation of Epoxy-Alumina Nano-Composites based High Voltage Dry-Type Insulation using Debye Model Parameters
CN114219016A (en) Oiled paper insulation moisture content comprehensive evaluation method based on adaboost
CN111638430A (en) Method and system for evaluating aging state of insulating paper under vibration-temperature synergistic effect
Xin et al. The vibration response of transformer windings under harmonic excitations and its applications
CN115468946A (en) Transformer oil aging diagnosis method and device and storage medium
CN117407631B (en) Method and system for evaluating high-voltage sleeve damp state
Nejedly et al. Evaluation of the extent of ageing of paper in oil-immersed power transformers
Guastavino et al. A study about partial discharge measurements performed applying to insulating systems square voltages with different rise times
Liu et al. Using Frequency Domain Dielectric Spectrum for Aging State Evaluation of Transformer Oil Paper Insulation
Duplessis et al. A smart way to minimize test time for transformer dielectric measurements
CN113341307B (en) Short-circuit closing capability test platform and method for open-type grounding switch

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant