1. Introduction
Power transformers are the most critical equipment in power systems, and the insulation condition of their internal oil–paper insulation directly affects the safe operation of transformers [
1]. With the continuous development of the power grid, the number of sudden short-circuit faults increases. The short-circuit electromagnetic force can cause winding deformation, thereby causing serious damage to the oil–paper insulation structure.
According to the classification of aging factors, the aging of insulation in transformers can be roughly divided into thermal aging, electrical aging, and mechanical aging [
2]. The thermal aging platform constructed in reference [
3] studied the thermal aging characteristics of insulating paper in mineral oil and vegetable oil, and the results showed that aging in vegetable oil is slower. The insulating paper thermal aging test platform built in reference [
4] studied the degradation law of the degree of polymerization of insulating paper, and the results showed that the rate of degree of polymerization decrease gradually slows down with aging time. Reference [
5] estimated the degree of polymerization (DP) and percentage of moisture content (PMCs) in oil-impregnated paper insulation using dielectric response sensing obtained from frequency-domain spectral tests. The insulating paper aging platform constructed in reference [
6] studied the microstructure and electrical characteristic changes in insulating paper under different AC/DC composite voltages and found consistency in the degradation mechanism and microstructure of insulating paper under different composite voltages. In reference [
7], accelerated thermal aging experiments were carried out on insulating paper at 120 °C and 130 °C, respectively, and the results showed that the crystallinity index and refractive index of cellulose fiber in insulating oil could be used as feasible indexes for the diagnosis of paper degradation. Reference [
8] studied the micro-slip wear between the insulation layers of transformers and found that abrupt axial mechanical forces can cause damage to insulating paper, thus affecting its insulation performance.
The above studies all investigated the effects of single aging conditions on the aging characteristics of insulating paper, while the aging of insulating paper is generally influenced by multiple factors. Reference [
9] established an electro-thermal coupled aging platform for aging research of insulating paper and used multivariate regression analysis to analyze the breakdown life of insulating paper, obtaining a lifetime model applicable to electro-thermal aging. Reference [
10] also established an electro-thermal coupled aging platform and found that the degree of polymerization decrease in electro-thermal aging samples is faster than that in single-factor thermal aging. Reference [
11] studied the frequency-domain dielectric properties of transformer oil–paper insulation under different aging stress conditions and aging times. Reference [
12] found that mechanical–thermal coupled aging has a more significant effect on the degradation of insulating paper mechanical properties compared to thermal aging, and the more intense the mechanical vibration, the more significant the decrease in the degree of polymerization and tensile strength of insulating paper. Reference [
13] conducted mechanical–electrical–thermal aging tests with experimentally given conditions and obtained an empirical model for the variation in tensile strength of insulating paper under the combined action of heating, electrical, and mechanical forces.
Most of the aforementioned studies considered the effects of two aging factors in the aging tests, and the experimental conditions were directly given. This paper will consider the combined aging of the electric field, temperature field, and mechanical stress on insulating paper. The values of each parameter will be obtained through simulation and theoretical calculations, which will be more representative of the operating conditions of insulating paper during circuit breaker reclosing.
In this paper, the short-circuit stress, electric field, and temperature field are calculated by setting up the reclosing current calculation model, the magnetic-circuit-force coupling model, and the simulation model of the electrostatic field and temperature field. Based on the above simulation results, the polymerization degree, tensile strength, elongation at break, dielectric loss factor, and breakdown field strength of insulating paper at different aging stages were measured and analyzed.
2. Finite Element Simulation of Oil–Paper Insulation
Figure 1 shows the insulating paper in the oil–paper insulation structure. The insulating paper is wrapped on the wire turns, the wire turns are wound to make a winding, and the winding is soaked in the insulating oil. At this time, the insulating paper and the insulating oil form an oil–paper insulation structure.
In this section, the applied forms of field strength, temperature, and mechanical stress are defined by simulation and theoretical calculation. A two-dimensional model of transformer main insulation is established to calculate the general field strength of insulating paper. A Simulink reclosing process simulation model and a 3D magnetic-channel-force coupling finite element simulation model of a transformer are built to calculate the mechanical stress of insulating paper. A two-dimensional model of the transformer temperature field is established to calculate the average winding temperature, and the maximum temperature of insulating paper is finally determined by theoretical calculation.
2.1. Electrodynamic Simulation of Oil–Paper Insulation
2.1.1. Reclosing Process Current Calculation
The Simulink reclosing process simulation model built based on SFSZ7-31500/110 transformer is shown in
Figure 2. The rated voltage is 110/38.5/10.5 kV, and the parameters of load loss P
k, no-load loss P
0, short-circuit voltage percentage U
k%, and no-load current percentage I
0% are shown in
Table 1.
The time of reclosing is as follows: The transformer starts to operate normally, and the three-phase short-circuit fault occurs at 0.1 s, and it is a permanent short-circuit fault. The 0.5 s circuit breaker operates in tripping and cuts out the fault part. The 0.6 s circuit breaker operates in closing and the closing angle is 0°. The current waveform of the process is shown in
Figure 3.
It can be seen from the simulation calculation that the secondary short-circuit impulse current is slightly larger than the primary short−circuit impulse current, which is due to the existence of the hysteresis effect in the opening process, so that there will be remanent magnetism inside the core and it is difficult to disappear in a short time. When the remanent magnetism and the direction of excitation are the same, the secondary short-circuit impulse current will be greater than the primary short-circuit impulse current. According to the current calculation model of the reclosing process established above, different closing angles are respectively set to calculate the impact current of a secondary short circuit under a three-phase short-circuit condition.
2.1.2. Simulation of Short-Circuit Stress in Magnetic-Structure Coupling
The electromagnetic force received by the wire cake element of the finite element model of the transformer in the electromagnetic field is [
14,
15]:
In the formula, Fi is the electromagnetic force received by the i unit of the winding, B is the magnetic induction intensity of the wire cake unit, S is the area of the wire cake unit, J is the current density of the wire cake unit, and R is the distance between the center line of the core and the center of gravity of the wire cake unit.
The total electromagnetic force of the cake is:
Through coupling, the electromagnetic force is added to the physical field of solid mechanics as a load. According to the elastic-plasticity of copper material, the electromagnetic force load is divided into several increments, and the balance equation of the finite element system based on increments is established in the physical field of solid mechanics:
where K
ep is the elastic-plastic stiffness matrix of the system, the incremental displacement phasor, and the unbalanced force phasor.
According to the constitutive relation of the material:
where: B is the strain matrix of the unit and is the strain increment.
where D
e is the elastic matrix, D
ep is the elastic–plastic matrix, and is the stress increment.
In order to calculate the short-circuit stress on the winding during reclosing, a three-dimensional magneto-circuit-force coupling finite element simulation model was established [
16,
17]. The following dimensions should be indicated in
Figure 4. “Among them, the center distance between the columns is 1335 mm, the height of the core window is 1670 mm, the diameter of the core is 600 mm, the inner diameter of the low-pressure winding is 331 mm, the outer diameter of the low-pressure winding is 396 mm, the inner diameter of the medium-pressure winding is 435 mm, the outer diameter of the medium-pressure winding is 500 mm, the inner diameter of the high-pressure winding is 547 mm, and the outer diameter of the high-pressure winding is 640.5 mm.”
The impulse current calculated by the reclosing analog circuit is used as the excitation and input into the three-dimensional magneto-circuit-force coupling finite element simulation model for coupling calculation, and the short-circuit stress distribution of the winding is obtained.
Figure 5 shows the distribution of short-circuit stress when the closing angle is 0°, and
Table 2 shows the short-circuit stress at different closing angles. When the closing angle is 0°, the maximum short-circuit stress of the insulating paper is 8.71 MPa. This test will take this as the basis for adding mechanical stress.
2.2. Oil–Paper Insulated Electric Field Simulation
Because the main insulation structure of the transformer can be approximated as a concentric cylinder with good symmetry, the two-dimensional model can meet the requirements when conducting an electric field simulation calculation. The two-dimensional simulation model established according to the actual parameters of the transformer is mainly composed of low-voltage winding, medium-voltage winding, high-voltage winding, insulating oil, insulated paper tube, and Angle ring, among which the low-voltage winding, medium-voltage winding and high-voltage winding are modeled in the unit of wire cake. The values of the relative dielectric constants of the insulation materials are shown in
Table 3 [
18].
The established calculation model of the main insulation of the transformer winding is shown in
Figure 6. The ground potential is set at the boundary of the model (a homogeneous boundary condition), 110 kV voltage is applied in the high-voltage winding, and the ground potential is also set in the middle-voltage and low-voltage winding. The calculated distribution results of the transformer main insulation composite electric field are shown in
Figure 7.
As can be seen from
Figure 7, the distribution of the electric field in the middle of the transformer between the high- and middle-voltage windings is close to uniform, roughly between 2 and 3 × 10
6 V/m. The average field strength between turns of the transformer is 2.26 × 10
6 V/m through simulation calculation. This test is to explore the aging characteristics of insulating paper, so the average field strength between turns is taken as the basis for pressurization.
2.3. Oil–Paper Insulation Temperature Field Simulation
Taking a three-phase three-column true-type transformer as the research object, in order to simplify the calculation, the following assumptions are made when establishing the two-dimensional model of the oil-immersed transformer [
19]:
Do not consider the clamp, oil tank, magnetic shield, and other structural parts, mainly including transformer oil, iron core, and each winding;
An axisymmetric model is adopted, and only the core and winding on one side of the axis are taken;
The thermal conductivity, density, and specific heat capacity of the transformer core and winding are constant and do not change with the change in temperature;
The ambient temperature is assumed to be 293 K. The two-dimensional calculation results of the transformer temperature field are shown in
Figure 8.
As can be seen from
Figure 8, the transformer winding temperature is high, the transformer oil temperature is stratified, and the temperature gradually decreases from top to bottom. The maximum average temperature of the winding during a short circuit is 107.67° C.
At the same time, according to GB/T 1094.5-2008 [
20], the average thermal stability temperature of transformer winding under a short circuit is calculated theoretically. GB/T 1094.5-2008 specifies that power transformers should be free from damage under the action of short-circuit overcurrent caused externally. It is specified that the duration of the current to withstand the short-circuit heat resistance is 2 s, and the calculation formula of the average thermal stability temperature of the copper winding after the transformer burst short circuit is given as follows:
where θ
1 is the average temperature after winding short circuit t, unit °C; and θ
0 is the starting temperature of the winding, which is the sum of the maximum ambient temperature and the average temperature rise limit of the winding, in unit °C. J is the short-circuit current density calculated according to the square mean root value of the symmetric short-circuit current, in A/mm
2; and t is the short circuit duration, unit s.
Taking a three-phase three-column true-type transformer as the research object, the cross-sectional area of each winding wire of the transformer is given in
Table 4.
When the high-voltage power supply is a short circuit, the short-circuit current value of the high voltage is 1542.3 A, and the short-circuit current value of the medium voltage is 4738.2 A. The high-voltage power supply, low-voltage short circuit, and high-voltage short-circuit current value is 880.9 A; the low-voltage short-circuit current value is 5727.4 A.
Taking high-voltage power supply and medium-voltage short-circuit working condition as an example, the short-circuit current calculation result is 1542.3 A, the cross-sectional area of the high-voltage winding is 61.44 mm
2, then the symmetric short-circuit current density J of the transformer high-voltage winding is the ratio of 25.10 A/mm
2, θ
0 is 105 °C, t is 2 s, and the Formula (6) is replaced:
Through the above calculation method, the average thermal stability temperature of the high- and low-voltage windings can be calculated, as shown in
Table 5.
Through simulation, the maximum average temperature of winding under dynamic stability is 107.67 °C, and through theoretical calculation, the maximum average temperature of winding under thermal stability is 113.7 °C. Considering the maximum temperature condition that insulating paper may withstand [
21], the maximum temperature of 113.7 °C is taken as the temperature condition of joint aging.
3. Electrical, Thermal, and Mechanical Aging Test of Oil–Paper Insulation
In order to study the performance aging law of insulating paper during transformer operation, the combined electrical, thermal, and mechanical aging test was carried out.
3.1. Electrical, Thermal, and Mechanical Aging Test Process of Oil–Paper Insulation
In order to reduce the moisture content in the insulating paper to a lower level, the insulating paper sample should be dried. First, the insulating paper (thickness 0.3 mm) should be cut into an electrode shape and placed in a beaker, and then the beaker should be placed in a vacuum drying oven at 105 ± 5 °C for 24 h. Taking into account the test error of a single-layer sample, five layers of insulating paper are put into the test oil tank, and 25# transformer oil is injected. This process is shown in
Figure 9. Electrodes on the upper and lower ends of the insulating paper are installed (the electrodes are connected to the power frequency high-voltage test console), five groups of samples are set and numbered respectively in the oil tank, and finally, the oil tank is put into the high- and low-temperature-alternating humid heat test tank. The test site diagram is shown in
Figure 10.
At the same time, 3.39 kV power frequency voltage is applied on both sides of the electrode to make the field strength between electrodes reach 2.26 × 106 V/m. The temperature is set to 113.7 °C; the sampling time of the test was 0 h, 24 h, 48 h, 96 h, 144 h, and 192 h. After aging, 8.71 MPa equivalent stress was applied to the insulating paper to simulate the mechanical stress subjected to the insulating paper.
The test was divided into two groups, one group applied only electric field and temperature conditions, called the electro-thermal combined aging group, and the other group applied electric field, temperature, and mechanical stress conditions, and was called the electro-thermal combined aging group.
3.2. Aging Parameter Measurement of Oiled Paper Insulation
The mechanical and electrical properties of insulating paper will deteriorate to different degrees under the condition of multi-element aging. In this section, the degree of polymerization, tensile strength, and elongation at break are selected to reflect the change in mechanical properties during aging. Dielectric loss factor and breakdown voltage were selected to reflect the change in electrical performance during aging.
3.2.1. Degree of Polymerization
The most intuitive way to reflect the aging of transformer oil–paper insulation is to detect the degree of polymerization of insulating paper. Insulating paper is the main reliance of oil-immersed transformer solid insulation. Cellulose composed of a glucose monomer connected to each other is the main component of insulating paper [
22]. The number of glucose monomers is the Degree of Polymerization. The average degree of polymerization of insulating paper is generally measured by the viscosity method. According to GB/T 1548-2004 [
23], the average degree of polymerization of insulating paper is calculated by measuring the viscosity of insulating paper [
24].
Firstly, the insulating paper is degreased to eliminate the influence of transformer oil in the insulating paper on viscosity measurement. Then, the degreased insulating paper is broken into fibers in water, and copied into quantitative pulp sheets on the hand paper machine before air drying. Finally, the insulated pulp sheet was immersed in 25 mL distilled water for 0.5 h, and then 25 mL copper ethylenediamine solution with a concentration of 1 mol/L was added, stirring to obtain the concentrated solution with the viscosity to be measured, and the viscosity of the cellulose solution was measured by using an Ulster viscometer.
Relative viscosity of cellulose solution:
where T
s is the cellulose solution flow time, and T
0 is the solvent flow time.
Water content of the insulating paper:
where m is the quality of the insulating paper before drying, and m
0 is the quality of the insulating paper after drying.
Cellulose solution concentration:
According to the relative viscosity, the value of the characteristic viscosity is found, and the average degree of polymerization of the insulating paper is obtained according to formula (11).
In the formula: K is 7.5 × 10−3.
3.2.2. Tensile Strength and Elongation at Break
A mechanical universal testing machine was used to test the tensile strength of insulating paper according to GB/T 12914-2018 standard [
25]. The sample size was 50 mm × 15 mm, and the drawing speed was 5 mm/min. The tensile strength and elongation at break were used as the evaluation basis.
Tensile strength refers to the maximum tension that the specimen can withstand before breaking per unit width under specified test conditions [
26]. The calculation formula is:
where
is the average tensile strength and
is the width of the insulating paper.
Elongation at break represents the ratio of the elongation of the specimen subjected to tension to fracture; that is, the ratio of the elongation length before and after stretching to the length before stretching. The formula is as follows:
where
is the length when pulled,
is the original length.
3.2.3. Dielectric Loss Factor
Dielectric loss factor
is one of the basic indicators to characterize the insulation performance of objects. It is a characteristic parameter to characterize the specific loss of insulators under the action of alternating voltage, and has nothing to do with the shape and size of the insulators [
27]. In this paper, a dielectric loss tester is used for measurement. The basic principle is the Schering bridge method. The diagram of the Schering bridge is shown in
Figure 11.
In
Figure 11, R
x and C
x represent the equivalent resistance and capacitance of the subject. By adjusting R
3 and C
4, the bridge reaches the balance, which is shown as the current of galvanometer P is zero. At this time:
Equivalent resistance and capacitance of the subject can be obtained:
Dielectric loss factor of the subject:
3.2.4. Breakdown Voltage
When the electric field intensity applied at both ends of the dielectric exceeds a certain critical value, a small number of freely moving carriers inside the dielectric move violently, colliding with atoms on the lattice, ionizing, and destroying the molecular structure, resulting in the final breakdown; this phenomenon is called electrical breakdown [
28].
The 20 s step-up voltage method was used to perform the voltage withstand test on the single-layer insulating paper. As shown in
Figure 12, the test electrodes were designed according to GB/T 1408.1-2016 [
29], and electrodes with different diameters were used: Apply to the sample at a short-time breakdown voltage of 40%. If the sample withstands this voltage for 20 s without breakdown, the voltage shall be increased step by step. Each increase in voltage should be immediately and continuously increased for 20 s until breakdown occurs.