CN113935171A - Broadband equivalent modeling method for dielectric response of oiled paper insulation - Google Patents

Broadband equivalent modeling method for dielectric response of oiled paper insulation Download PDF

Info

Publication number
CN113935171A
CN113935171A CN202111204060.1A CN202111204060A CN113935171A CN 113935171 A CN113935171 A CN 113935171A CN 202111204060 A CN202111204060 A CN 202111204060A CN 113935171 A CN113935171 A CN 113935171A
Authority
CN
China
Prior art keywords
dielectric response
bee
honey
temperature
broadband
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.)
Pending
Application number
CN202111204060.1A
Other languages
Chinese (zh)
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.)
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Original Assignee
Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
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 Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd filed Critical Electric Power Research Institute of State Grid Shandong Electric Power Co Ltd
Priority to CN202111204060.1A priority Critical patent/CN113935171A/en
Publication of CN113935171A publication Critical patent/CN113935171A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R31/00Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
    • G01R31/12Testing 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/1227Testing 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/1263Testing 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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

The invention relates to a broadband equivalent modeling method for dielectric response of oiled paper insulation, and belongs to the field of electrical engineering. The method comprises the following steps: s1: preparing an oiled paper insulation sample; s2: simulating artificial moisture of a test object; s3: measurement of frequency domain dielectric response (FDS) data; s4: modeling dielectric response data; s5: and (5) verifying results, and summarizing model data. The method can realize accurate modeling of frequency domain dielectric response in a wide frequency range, provides a digital and quantitative tool for further adopting the model to read the FDS information of the oilpaper insulation and obtain the insulation state of the FDS information, and provides a preliminary foundation for next research on the change rule of the model parameters.

Description

Broadband equivalent modeling method for dielectric response of oiled paper insulation
Technical Field
The invention belongs to the field of electrical engineering, and relates to a broadband equivalent modeling method for oil paper insulation dielectric response.
Background
The problem currently faced by the grid operations sector is how to maximize the service life of electrical equipment while meeting the system load requirements. And the service life of the electrical equipment can be prolonged through proper monitoring technology and field diagnosis. In the existing various insulation state monitoring technologies, modern dielectric diagnosis is supplemented with advanced modeling and analyzing tools, and an accurate and reliable method for evaluating the insulation state of high-voltage electrical equipment is expected to be provided. On the other hand, the oil paper insulation equipment is the most important high-voltage electrical equipment in a power grid, and monitoring and analyzing the equipment is always one of key operation and maintenance work of power supply departments at all levels. Due to the combined action of a plurality of factors such as mechanical pressure, chemical corrosion, surface dirt, operating voltage and humidity, the performance of the oil paper combined insulation is gradually reduced, and the insulation breakdown can be finally caused to cause equipment failure. Frequency Domain Spectroscopy (FDS) is a new type of nondestructive testing method that is being developed in power departments at various levels in recent years, and the basic idea is that each dielectric has its own dielectric response function, which is affected by the insulation condition. Therefore, the dielectric response measurement of the insulating medium can be used as a main basis for deducing the moisture and aging degree of the insulating medium.
At present, the FDS method has accumulated a large amount of abundant experimental data from early laboratory theoretical research and experimental exploration to field application of the power department, and also discovers the rough rule that the measuring result changes along with the state of a measuring object. However, in view of the accumulated data of these stocks, how to quantitatively solve the insulation state information of moisture, aging, etc. corresponding to the FDS result and finally give a diagnosis result with high reliability has been a main technical problem that restricts further field application of the technology.
Disclosure of Invention
In view of the above, the present invention provides a broadband equivalent modeling method for dielectric response of oiled paper insulation. A broadband equivalent circuit model is correspondingly established according to frequency domain dielectric response (FDS) data of the oil paper insulation of the tested equipment, and model parameters are obtained, so that a digital and quantitative tool is provided for further adopting the model to read the FDS information of the oil paper insulation and obtaining the insulation state of the FDS information.
In order to achieve the purpose, the invention provides the following technical scheme:
a broadband equivalent modeling method for dielectric response of oiled paper insulation comprises the following steps:
s1: preparing an oiled paper insulation sample;
s2: simulating artificial moisture of a test object;
s3: measurement of frequency domain dielectric response (FDS) data;
s4: modeling dielectric response data;
s5: and (5) verifying results, and summarizing model data.
Optionally, the S1 specifically includes: aiming at the oiled paper capacitive test casing, respectively placing the coiled capacitor core and insulating oil for infiltrating the core in special vacuum tanks with different sizes for drying pretreatment, vacuumizing the vacuum tanks, and then placing the vacuum tanks into a temperature box, wherein the drying temperature is set to be 90 ℃, and performing vacuum drying for 48 hours; after the vacuum drying is finished, the capacitor core is immersed in oil for 48 hours in a vacuum environment at the temperature of 60 ℃; and finally, filling nitrogen into the vacuum tank for oil immersion for 24 hours so as to fully immerse the oil in the gaps between the sleeve and the oil paper.
Optionally, the S2 specifically includes: presetting different water content percentages within the range of 0.7% -2.5% according to different moisture degrees of the oiled paper insulating sleeve, and carrying out moisture absorption treatment on the sample prepared by S1; placing the capacitor core on a precision electronic balance, recording initial mass, calculating mass required for absorbing moisture to a target moisture, opening a humidifier, increasing air humidity, observing indication of the balance, stopping absorbing moisture once the target mass is reached, and standing for 1-2 hours.
Optionally, the S3 specifically includes: after the preparation of the damped sleeve is finished, putting the sample sleeve into a temperature box for dielectric response measurement at different temperatures, wherein the broadband measurement range is 1mHz-5 kHz;
the measurement steps are as follows:
s31: taking 30 ℃ as an initial temperature measuring point, taking 10 ℃ as a variable temperature step length, raising the temperature and keeping the temperature to the next temperature measuring point for 6 hours each time until the temperature reaches the last measuring point by 80 ℃; carrying out frequency domain dielectric spectrum measurement of the test sample sleeve and real-time monitoring of the temperature and the micro-water content of the insulating oil at each test temperature;
s32: after the frequency domain measurement is completed, the insulation resistance of the casing at the temperature is measured by a high resistance meter.
Optionally, the S4 specifically includes:
s41: an extended Debye model is used as a basic model, and based on the model, the FDS expression is listed as follows:
Figure BDA0003306222580000021
Figure BDA0003306222580000022
Figure BDA0003306222580000023
in the above formula, C' (ω) and C "(ω) represent the real part value and the imaginary part value of the complex capacitance, respectively, and tan δ is the dielectric loss factor;
s42: constructing an objective function for a parameter identification target of the model; the purpose of parameter identification is to calculate the element parameter values of all branches in the model, so that the error between a calculated spectrogram obtained according to the formulas (1) to (3) and an actually measured spectrogram is minimized, and a real part C '(omega) and an imaginary part C' (omega) of a complex capacitor are selected as reference data for parameter estimation;
Figure BDA0003306222580000031
in the formula, the subscript is "measure" to represent the related measured value, and correspondingly, the subscript "fit" to represent the calculated value of the corresponding parameter obtained in real time according to the parameter identification result at a certain time;
s43: an improved artificial bee colony algorithm is adopted, FDS data measured by S3 is used, the broadband range is 1mHz-5kHz, the objective function (4) is solved under each group of data, and the purpose is to obtain equivalent model parameter values of tested objects with different moisture degrees at different operating temperatures; the method comprises the following steps:
s431: for the time when n is 0, randomly generating XiA feasible solution
Figure BDA0003306222580000032
These feasible solutions correspond to the honey sources one-to-one, producing feasible solution XiThe formula:
Figure BDA0003306222580000033
Xminand XmaxIs the upper and lower bounds of the solution space, j takes on {1,2,3 }; calculating a feasible solution XiThen sorting the fitness function according to the size of the fitness function to be used as an initial honey bee colony X (0);
s432: and (3) updating the positions of the honey bees in the nth step:
Figure BDA0003306222580000034
wherein k is {1,2,3, …, Ne }, and k is not generated randomly from i, k and j, and rand () is [ -1,1 ];
s433: after the honey collection bees are updated to a new position, selecting a more optimal position according to greedy selection, reserving the more optimal position for the next generation of population, and after the information of the honey collection bees is shared by the following bees, following the honey collection bees through a probability distribution function:
Figure BDA0003306222580000035
ts represents the random mapping from individual space to individual space, and is expressed as S2→S;
S434: in the following bee searching stage, a bee source signal is transmitted to the following bee according to the honey bee, the probability of selecting the following bee is calculated by the formula (8), the following bee selects a bee source and searches nearby, and the following bee updates:
Figure BDA0003306222580000036
in the formula, phi, wa,wbTaking values according to the formulas (10), (11) and (12) respectively;
Figure BDA0003306222580000037
wa=w2((iter-maxcycle)/max cycle)a×(w2-w1) (11)
wb=w3-((iter-max cycle)/max cycle)b×(w4-w3) (12)
Xiis the current position of the following bee, the inertial weight w1,w2,w3,w4The value range is [0.1, 1.5 ]],w1,w3Less than w2,w4Iter is the current iteration number, maxcycle is the maximum iteration number, in is dynamically changed, when the new position fitness value is better than the original position, lambda is larger than 1, otherwise, is smaller than 1, a is in [0.8, 1]]B is in [1, 1].2]Phi is an adaptive factor;
s435: during the searching process, if honey source XiIf the honey source position information reaching the threshold limit is not changed better through the iterative search for three times, the honey source X isiWill be abandoned, meanwhile, the honey bee collecting is changed into the observation bee; observing bees randomly generating a new honey source in a search space to replace X according to the formula (4)i
S436: and if the stopping condition is met, outputting the optimal fitness value and the corresponding feasible solution, otherwise, turning to S432.
Optionally, the S5 specifically includes:
and reconstructing an FDS spectrogram by adopting the data completed by modeling, comparing the FDS spectrogram with the actually measured spectrogram, testing the modeling accuracy, completing the broadband equivalent modeling of the insulation dielectric response of the oil paper, and realizing the accurate equivalence of the insulation dielectric response of the oil paper in a broadband range of 1mHz-5kHz based on a Debye equivalent model.
The invention has the beneficial effects that: the method can realize accurate modeling of frequency domain dielectric response in a wide frequency range, provides a digital and quantitative tool for further adopting the model to read the FDS information of the oilpaper insulation and obtain the insulation state of the FDS information, and provides a preliminary foundation for next research on the change rule of the model parameters.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention. The objectives and other advantages of the invention may be realized and attained by the means of the instrumentalities and combinations particularly pointed out hereinafter.
Drawings
For the purposes of promoting a better understanding of the objects, aspects and advantages of the invention, reference will now be made to the following detailed description taken in conjunction with the accompanying drawings in which:
FIG. 1 is a circuit diagram of an extended Debye broadband equivalent model employed in the present invention;
FIG. 2 is a real part and an imaginary part of complex capacitance of a dielectric response broadband modeling result of a tested oiled paper insulation sample with the moisture content percentage of 0.71% by adopting the method;
FIG. 3 is a dielectric loss tangent of a dielectric response broadband modeling of a tested oiled paper insulation sample with a moisture content percentage of 0.71% according to the present invention;
FIG. 4 is a real part and an imaginary part of complex capacitance of a dielectric response broadband modeling result of a tested oiled paper insulation sample with the moisture content percentage of 1.98% by adopting the method;
FIG. 5 shows the dielectric loss tangent of the dielectric response of the insulation sample of the tested oiled paper with a moisture content of 1.98% according to the present invention.
Detailed Description
The embodiments of the present invention are described below with reference to specific embodiments, and other advantages and effects of the present invention will be easily understood by those skilled in the art from the disclosure of the present specification. The invention is capable of other and different embodiments and of being practiced or of being carried out in various ways, and its several details are capable of modification in various respects, all without departing from the spirit and scope of the present invention. It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present invention in a schematic way, and the features in the following embodiments and examples may be combined with each other without conflict.
Wherein the showings are for the purpose of illustrating the invention only and not for the purpose of limiting the same, and in which there is shown by way of illustration only and not in the drawings in which there is no intention to limit the invention thereto; to better illustrate the embodiments of the present invention, some parts of the drawings may be omitted, enlarged or reduced, and do not represent the size of an actual product; it will be understood by those skilled in the art that certain well-known structures in the drawings and descriptions thereof may be omitted.
The same or similar reference numerals in the drawings of the embodiments of the present invention correspond to the same or similar components; in the description of the present invention, it should be understood that if there is an orientation or positional relationship indicated by terms such as "upper", "lower", "left", "right", "front", "rear", etc., based on the orientation or positional relationship shown in the drawings, it is only for convenience of description and simplification of description, but it is not an indication or suggestion that the referred device or element must have a specific orientation, be constructed in a specific orientation, and be operated, and therefore, the terms describing the positional relationship in the drawings are only used for illustrative purposes, and are not to be construed as limiting the present invention, and the specific meaning of the terms may be understood by those skilled in the art according to specific situations.
Please refer to fig. 1 to 5, which illustrate a wideband equivalent modeling method for dielectric response of oiled paper insulation.
S1: preparing an oiled paper insulation sample;
the invention takes an oiled paper capacitive test sleeve as an example, the coiled capacitor core and the insulating oil for infiltrating the core are respectively placed in special vacuum tanks with different sizes for drying pretreatment, and the vacuum tanks are vacuumized and then placed in a temperature box (the drying temperature is set to be 90 ℃) for vacuum drying for 48 hours. And after the vacuum drying is finished, the capacitor core is immersed in oil for 48 hours in a vacuum environment at the temperature of 60 ℃. And finally, filling nitrogen into the vacuum tank for 24 hours to fully soak the oil in the oil paper gaps of the sleeve, wherein the pretreatment process is consistent with the production process flow of the real sleeve.
S2: simulating artificial moisture of a test object;
the invention takes different moisture degrees of the oiled paper insulating sleeve as an example, different water content percentages within the range of 0.7% -2.5% are preset, and the sample prepared in the step S1 is subjected to moisture absorption treatment. Placing the capacitor core on a precision electronic balance, recording initial mass, calculating mass required for absorbing moisture to a target moisture, opening a humidifier, increasing air humidity, observing indication of the balance, stopping absorbing moisture once the target mass is reached, and standing for 1-2 hours.
S3: measurement of frequency domain dielectric response (FDS) data;
after the preparation of the humidified sleeve is finished, putting the sample sleeve into a temperature box to perform dielectric response measurement (the broadband measurement range is 1mHz-5kHz) at different temperatures, wherein the measurement steps are as follows: firstly, taking 30 ℃ as an initial temperature measuring point, taking 10 ℃ as a variable temperature step length, raising the temperature and keeping the temperature to the next temperature measuring point for 6 hours each time until the temperature reaches the last measuring point, namely 80 ℃; carrying out frequency domain dielectric spectrum measurement of the test sample sleeve and real-time monitoring of the temperature and the micro-water content of the insulating oil at each test temperature; and secondly, after the frequency domain measurement is finished, measuring the insulation resistance of the sleeve at the temperature by a high-resistance meter.
S4: modeling dielectric response data;
s41: the extended debye model shown in fig. 1 is used as a basic model, and based on the model, the expression of FDS is listed as follows:
Figure BDA0003306222580000061
Figure BDA0003306222580000062
Figure BDA0003306222580000063
in the above equation, C' (ω) and C "(ω) represent the real and imaginary values of the complex capacitance, respectively, and tan δ is the dielectric loss factor.
S42: the following objective function is constructed for the parameter recognition goal of the model. The purpose of parameter identification is to calculate the element parameter values of all branches in the model, so that the error between the calculated spectrogram obtained according to the formulas (1) to (3) and the actually measured spectrogram is minimized, and the real part C '(omega) and the imaginary part C' (omega) of the complex capacitance are selected as reference data for parameter estimation.
Figure BDA0003306222580000064
In the formula, the subscript "measure" represents the measured value, and correspondingly, the subscript "fit" represents the calculated value of the corresponding parameter obtained in real time according to the parameter identification result at a certain time.
S43: and (3) respectively solving the objective function (4) under each group of data by adopting an improved artificial bee colony algorithm and using the FDS data (the broadband range is 1mHz-5kHz) measured in the step S3, so as to obtain the equivalent model parameter values of the tested object with different moisture degrees at different operating temperatures. The steps of this section are as follows:
s431: for the time when n is 0, randomly generating XiA feasible solution
Figure BDA0003306222580000065
These solutions correspond one-to-one with the honey sources, producing a feasible solution XiThe formula:
Figure BDA0003306222580000066
Xminand XmaxIs the upper and lower bounds of the solution space, j takes on {1,2,3 }. Calculating a feasible solution XiThe fitness function of (2) is then ranked according to the fitness function size as the initial bee-collecting population X (0).
S432: and (3) updating the positions of the honey bees in the nth step:
Figure BDA0003306222580000071
wherein k is {1,2, 3.,. Ne }, and k is not generated randomly from i, k and j, and rand () is [ -1,1 ].
S433: after the honey collection bees are updated to a new position, selecting a more optimal position according to greedy selection, reserving the more optimal position for the next generation of population, and after the information of the honey collection bees is shared by the following bees, following the honey collection bees through a probability distribution function:
Figure BDA0003306222580000072
ts represents the random mapping from individual space to individual space, and is expressed as S2→S。
S434: in the following bee searching stage, a bee source signal is transmitted to the following bee according to the honey bee, the probability of selecting the following bee is calculated by the formula (8), the following bee selects a bee source and searches nearby, and the following bee updates:
Figure BDA0003306222580000073
in the formula, phi, wa,wbThe values are respectively expressed by the formulas (10), (11) and (12).
Figure BDA0003306222580000074
wa=w2((iter-maxcycle)/max cycle)a×(w2-w1) (11)
wb=w3-((iter-max cycle)/max cycle)b×(w4-w3) (12)
XiIs the current position of the following bee, the inertial weight w1,w2,w3,w4The value range is [0.1, 1.5 ]],w1,w3Less than w2,w4Iter is the current iteration number, maxcycle is the maximum iteration number, in is dynamically changed, when the new position fitness value is better than the original position, lambda is larger than 1, otherwise, is smaller than 1, a is in [0.8, 1]]B is in the range of [1, 1.2 ]]Phi is the adaptive factor.
S435: during the searching process, if honey source XiIf the honey source position information reaching the threshold limit is not changed better through the iterative search for three times, the honey source X isiWill be abandoned and at the same time the honey bees will be changed into observation bees. Observing bees randomly generating a new honey source in a search space to replace X according to the formula (4)i
S436: and if the stopping condition is met, outputting the optimal fitness value and the corresponding feasible solution, otherwise, turning to the step S432.
S5: verifying results, and summarizing model data;
and reconstructing an FDS spectrogram by adopting the data completed by modeling, comparing the FDS spectrogram with the actually measured spectrogram, testing the modeling accuracy, completing the broadband equivalent modeling of the insulation dielectric response of the oil paper, and realizing the accurate equivalence of the insulation dielectric response of the oil paper in a broadband range of 1mHz-5kHz based on a Debye equivalent model.
Table 1 shows the results of the broadband dielectric response modeling (parameter values and goodness of fit of each element of the model) performed on the two tested oil paper insulation samples with the moisture contents of 0.71% and 1.98% by using the method of the present invention.
TABLE 1 oil paper insulation parameter values
Moisture content of insulation to be tested 0.71% 1.98%
C0/nF 0.250143003 0.250442603
R0/GΩ 4800 200
C1/nF 0.467695730 0.618145409
R1/GΩ 419.219612034 346.627553222
C2/nF 0.001903743 0.044946622
R2/GΩ 6.991608756 169.933788168
C3/nF 0.003034954 0.002553811
R3/GΩ 0.013460255 0.324551425
C4/nF 0.035677610 0.013360038
R4/GΩ 210.968670961 19.866349930
C5/nF 0.005904552 0.003621919
R5/GΩ 43.355216864 0.012413917
C6/nF 0.001547771 0.004152651
R6/GΩ 0.511631100 3.251985468
Goodness of fit 0.996942845 0.996262229
Finally, the above embodiments are only intended to illustrate the technical solutions of the present invention and not to limit the present invention, and although the present invention has been described in detail with reference to the preferred embodiments, it will be understood by those skilled in the art that modifications or equivalent substitutions may be made on the technical solutions of the present invention without departing from the spirit and scope of the technical solutions, and all of them should be covered by the claims of the present invention.

Claims (6)

1. A broadband equivalent modeling method for dielectric response of oiled paper insulation is characterized by comprising the following steps: the method comprises the following steps:
s1: preparing an oiled paper insulation sample;
s2: simulating artificial moisture of a test object;
s3: measurement of frequency domain dielectric response (FDS) data;
s4: modeling dielectric response data;
s5: and (5) verifying results, and summarizing model data.
2. The broadband equivalent modeling method for dielectric response of oiled paper insulation according to claim 1, characterized in that: the S1 specifically includes: aiming at the oiled paper capacitive test casing, respectively placing the coiled capacitor core and insulating oil for infiltrating the core in special vacuum tanks with different sizes for drying pretreatment, vacuumizing the vacuum tanks, and then placing the vacuum tanks into a temperature box, wherein the drying temperature is set to be 90 ℃, and performing vacuum drying for 48 hours; after the vacuum drying is finished, the capacitor core is immersed in oil for 48 hours in a vacuum environment at the temperature of 60 ℃; and finally, filling nitrogen into the vacuum tank for oil immersion for 24 hours so as to fully immerse the oil in the gaps between the sleeve and the oil paper.
3. The broadband equivalent modeling method for dielectric response of oiled paper insulation according to claim 2, characterized in that: the S2 specifically includes: presetting different water content percentages within the range of 0.7% -2.5% according to different moisture degrees of the oiled paper insulating sleeve, and carrying out moisture absorption treatment on the sample prepared by S1; placing the capacitor core on a precision electronic balance, recording initial mass, calculating mass required for absorbing moisture to a target moisture, opening a humidifier, increasing air humidity, observing indication of the balance, stopping absorbing moisture once the target mass is reached, and standing for 1-2 hours.
4. The broadband equivalent modeling method for dielectric response of oiled paper insulation according to claim 3, characterized in that: the S3 specifically includes: after the preparation of the damped sleeve is finished, putting the sample sleeve into a temperature box for dielectric response measurement at different temperatures, wherein the broadband measurement range is 1mHz-5 kHz;
the measurement steps are as follows:
s31: taking 30 ℃ as an initial temperature measuring point, taking 10 ℃ as a variable temperature step length, raising the temperature and keeping the temperature to the next temperature measuring point for 6 hours each time until the temperature reaches the last measuring point by 80 ℃; carrying out frequency domain dielectric spectrum measurement of the test sample sleeve and real-time monitoring of the temperature and the micro-water content of the insulating oil at each test temperature;
s32: after the frequency domain measurement is completed, the insulation resistance of the casing at the temperature is measured by a high resistance meter.
5. The broadband equivalent modeling method for dielectric response of oiled paper insulation according to claim 4, characterized in that: the S4 specifically includes:
s41: an extended Debye model is used as a basic model, and based on the model, the FDS expression is listed as follows:
Figure FDA0003306222570000011
Figure FDA0003306222570000021
Figure FDA0003306222570000022
in the above formula, C' (ω) and C "(ω) represent the real part value and the imaginary part value of the complex capacitance, respectively, and tan δ is the dielectric loss factor;
s42: constructing an objective function for a parameter identification target of the model; the purpose of parameter identification is to calculate the element parameter values of all branches in the model, so that the error between a calculated spectrogram obtained according to the formulas (1) to (3) and an actually measured spectrogram is minimized, and a real part C '(omega) and an imaginary part C' (omega) of a complex capacitor are selected as reference data for parameter estimation;
Figure FDA0003306222570000023
in the formula, the subscript is "measure" to represent the related measured value, and correspondingly, the subscript "fit" to represent the calculated value of the corresponding parameter obtained in real time according to the parameter identification result at a certain time;
s43: an improved artificial bee colony algorithm is adopted, FDS data measured by S3 is used, the broadband range is 1mHz-5kHz, the objective function (4) is solved under each group of data, and the purpose is to obtain equivalent model parameter values of tested objects with different moisture degrees at different operating temperatures; the method comprises the following steps:
s431: for the time when n is 0, randomly generating XiA feasible solution
Figure FDA0003306222570000024
These feasible solutions correspond to the honey sources one-to-one, producing feasible solution XiThe formula:
Figure FDA0003306222570000025
Xminand XmaxIs the upper and lower bounds of the solution space, j takes on {1,2,3 }; calculating a feasible solution XiThen sorting the fitness function according to the size of the fitness function to be used as an initial honey bee colony X (0);
s432: and (3) updating the positions of the honey bees in the nth step:
Figure FDA0003306222570000026
wherein k is {1,2,3, …, Ne }, and k is not generated randomly from i, k and j, and rand () is [ -1,1 ];
s433: after the honey collection bees are updated to a new position, selecting a more optimal position according to greedy selection, reserving the more optimal position for the next generation of population, and after the information of the honey collection bees is shared by the following bees, following the honey collection bees through a probability distribution function:
Figure FDA0003306222570000027
ts represents the random mapping from individual space to individual space, and is expressed as S2→S;
S434: in the following bee searching stage, a bee source signal is transmitted to the following bee according to the honey bee, the probability of selecting the following bee is calculated by the formula (8), the following bee selects a bee source and searches nearby, and the following bee updates:
Figure FDA0003306222570000031
in the formula, phi, wa,wbTaking values according to the formulas (10), (11) and (12) respectively;
Figure FDA0003306222570000032
wa=w2((iter-maxcycle)/max cycle)a×(w2-w1) (11)
wb=w3-((iter-max cycle)/max cycle)b×(w4-w3) (12)
Xiis the current position of the following bee, the inertial weight w1,w2,w3,w4The value range is [0.1, 1.5 ]],w1,w3Less than w2,w4Iter is the current iteration number, maxcycle is the maximum iteration number, in is dynamically changed, when the new position fitness value is better than the original position, lambda is larger than 1, otherwise, is smaller than 1, a is in [0.8, 1]]B is in the range of [1, 1.2 ]]Phi is an adaptive factor;
s435: during the searching process, if honey source XiIf the honey source position information reaching the threshold limit is not changed better through the iterative search for three times, the honey source X isiWill be abandoned, meanwhile, the honey bee collecting is changed into the observation bee; observing bees randomly generating a new honey source in a search space to replace X according to the formula (4)i
S436: and if the stopping condition is met, outputting the optimal fitness value and the corresponding feasible solution, otherwise, turning to S432.
6. The broadband equivalent modeling method for dielectric response of oiled paper insulation according to claim 5, characterized in that: the S5 specifically includes:
and reconstructing an FDS spectrogram by adopting the data completed by modeling, comparing the FDS spectrogram with the actually measured spectrogram, testing the modeling accuracy, completing the broadband equivalent modeling of the insulation dielectric response of the oil paper, and realizing the accurate equivalence of the insulation dielectric response of the oil paper in a broadband range of 1mHz-5kHz based on a Debye equivalent model.
CN202111204060.1A 2021-10-15 2021-10-15 Broadband equivalent modeling method for dielectric response of oiled paper insulation Pending CN113935171A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111204060.1A CN113935171A (en) 2021-10-15 2021-10-15 Broadband equivalent modeling method for dielectric response of oiled paper insulation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111204060.1A CN113935171A (en) 2021-10-15 2021-10-15 Broadband equivalent modeling method for dielectric response of oiled paper insulation

Publications (1)

Publication Number Publication Date
CN113935171A true CN113935171A (en) 2022-01-14

Family

ID=79279551

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111204060.1A Pending CN113935171A (en) 2021-10-15 2021-10-15 Broadband equivalent modeling method for dielectric response of oiled paper insulation

Country Status (1)

Country Link
CN (1) CN113935171A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116151064A (en) * 2022-12-30 2023-05-23 西南大学 Dielectric parameter multi-partition inversion method for oilpaper insulating sleeve core

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116151064A (en) * 2022-12-30 2023-05-23 西南大学 Dielectric parameter multi-partition inversion method for oilpaper insulating sleeve core
CN116151064B (en) * 2022-12-30 2023-08-01 西南大学 Dielectric parameter multi-partition inversion method for oilpaper insulating sleeve core

Similar Documents

Publication Publication Date Title
CN109031014B (en) Transformer comprehensive reliability assessment and prediction method based on operation data
CN112505494B (en) Method and device for evaluating insulation water content of oiled paper
CN110009236A (en) A kind of oil-immersed power transformer built-in electrical insulation degree of aging quantitative evaluating method
CN113049962A (en) LSTM-based energy storage device operation situation deduction method
CN108387824B (en) Oiled paper insulation state evaluation method for improving return voltage measurement mode
CN107944571B (en) Method for predicting residual service life of power transformer
CN108872820B (en) Method and system for evaluating insulation aging state of oil impregnated paper in high-voltage current transformer
CN113484813B (en) Intelligent ammeter fault rate prediction method and system under multi-environment stress
CN113935171A (en) Broadband equivalent modeling method for dielectric response of oiled paper insulation
CN114814501B (en) On-line diagnosis method for capacitor breakdown fault of capacitor voltage transformer
CN112668604A (en) Sleeve insulation moisture prediction method based on frequency domain dielectric response and support vector machine
CN104091213B (en) Floor surface temperature predicting method and system
CN111797566A (en) Method and system for determining key characteristic quantity for representing health state of transformer
CN108846239A (en) A kind of elastic epoxy resin accelerated storage test and appraisal procedure based on temperature and humidity
CN114624526A (en) Stepping stress accelerated life test method for evaluating reliability of electric meter
Ding et al. NSGA-II model-based dielectric frequency response parameters for aging and moisture analysis of transformer insulation
CN107727939A (en) Three branch roads for paper oil insulation improve equivalent model
Fan et al. Moisture evaluation of oil-immersed insulation in bushing based on frequency domain spectroscopy and grey relational analysis
CN111896849A (en) Oil paper insulation performance analysis method and device considering non-uniform medium
CN115329672A (en) Converter transformer state reliability assessment method based on improved TOPSIS gray correlation
CN112329825B (en) Transformer mechanical fault diagnosis method based on information dimension division and decision tree lifting
CN114528726A (en) Method and device for correcting dielectric spectrum curve of oiled paper insulation frequency domain at time-varying temperature
CN112419093A (en) Load model characteristic parameter extraction method and device based on clustering algorithm
CN112668145A (en) FDS and exponential decay model-based transformer oiled paper insulation moisture assessment method
CN113283176B (en) Optimization method and system for dielectric response extended debye model

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
TA01 Transfer of patent application right

Effective date of registration: 20220302

Address after: No.2000 Wangyue Road, Shizhong District, Jinan City, Shandong Province

Applicant after: ELECTRIC POWER RESEARCH INSTITUTE OF STATE GRID SHANDONG ELECTRIC POWER Co.

Address before: No.2000 Wangyue Road, Shizhong District, Jinan City, Shandong Province

Applicant before: SHANDONG ELECTRIC POWER Research Institute

TA01 Transfer of patent application right