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 PDFInfo
- 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
Links
- 238000009413 insulation Methods 0.000 title claims abstract description 47
- 230000004044 response Effects 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005259 measurement Methods 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 16
- 241000256844 Apis mellifera Species 0.000 claims description 36
- 235000012907 honey Nutrition 0.000 claims description 24
- 241000257303 Hymenoptera Species 0.000 claims description 16
- 239000003990 capacitor Substances 0.000 claims description 11
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 6
- 238000001035 drying Methods 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 6
- 238000001291 vacuum drying Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 241000256837 Apidae Species 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 claims description 3
- 238000005315 distribution function Methods 0.000 claims description 3
- 238000001453 impedance spectrum Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 229910052757 nitrogen Inorganic materials 0.000 claims description 3
- 238000002360 preparation method Methods 0.000 claims description 3
- 238000007654 immersion Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 3
- 230000008859 change Effects 0.000 abstract description 2
- 238000004870 electrical engineering Methods 0.000 abstract description 2
- 238000003745 diagnosis Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000032683 aging Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000009659 non-destructive testing Methods 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 238000005316 response function Methods 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical 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
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:
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;
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
These feasible solutions correspond to the honey sources one-to-one, producing feasible solution XiThe formula:
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:
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:
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:
in the formula, phi, wa,wbTaking values according to the formulas (10), (11) and (12) respectively;
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:
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.
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
These solutions correspond one-to-one with the honey sources, producing a feasible solution XiThe formula:
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:
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:
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:
in the formula, phi, wa,wbThe values are respectively expressed by the formulas (10), (11) and (12).
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:
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;
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
These feasible solutions correspond to the honey sources one-to-one, producing feasible solution XiThe formula:
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:
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:
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:
in the formula, phi, wa,wbTaking values according to the formulas (10), (11) and (12) respectively;
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.
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)
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 |
-
2021
- 2021-10-15 CN CN202111204060.1A patent/CN113935171A/en active Pending
Cited By (2)
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 |