WO2021088453A1 - Dielectric state analysis method and system, computer and storage medium - Google Patents

Dielectric state analysis method and system, computer and storage medium Download PDF

Info

Publication number
WO2021088453A1
WO2021088453A1 PCT/CN2020/109652 CN2020109652W WO2021088453A1 WO 2021088453 A1 WO2021088453 A1 WO 2021088453A1 CN 2020109652 W CN2020109652 W CN 2020109652W WO 2021088453 A1 WO2021088453 A1 WO 2021088453A1
Authority
WO
WIPO (PCT)
Prior art keywords
singular point
real
imaginary part
real part
frequency
Prior art date
Application number
PCT/CN2020/109652
Other languages
French (fr)
Chinese (zh)
Inventor
高岩峰
卢毅
张旭
王馨
范硕超
王书渊
薛文祥
王辉
陈原
张吉飞
苏斌
Original Assignee
国网冀北电力有限公司电力科学研究院
华北电力科学研究院有限责任公司
国家电网有限公司
国网冀北电力有限公司
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 国网冀北电力有限公司电力科学研究院, 华北电力科学研究院有限责任公司, 国家电网有限公司, 国网冀北电力有限公司 filed Critical 国网冀北电力有限公司电力科学研究院
Publication of WO2021088453A1 publication Critical patent/WO2021088453A1/en

Links

Images

Classifications

    • 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
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/26Measuring inductance or capacitance; Measuring quality factor, e.g. by using the resonance method; Measuring loss factor; Measuring dielectric constants ; Measuring impedance or related variables
    • G01R27/2617Measuring dielectric properties, e.g. constants

Definitions

  • the present invention relates to the field of dielectrics, and in particular to a method, system, computer and storage medium for analyzing the state of dielectrics.
  • the main purpose of the embodiments of the present invention is to provide a dielectric state analysis method, system, computer, and storage medium to accurately analyze and determine the dielectric state, avoid misjudgment, missed judgment, and misjudgment, and save the cost of power equipment operation and maintenance.
  • an embodiment of the present invention provides a dielectric state analysis method, including:
  • the limited frequency band is located between the lowest frequency and the highest frequency;
  • the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency
  • the data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency
  • the real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency
  • the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part
  • Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral
  • Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral
  • the conductance information the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
  • the embodiment of the present invention also provides a dielectric state analysis system, including:
  • An acquiring unit for acquiring real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies
  • the determining unit is used to determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
  • the finite frequency band fitting function unit is used to fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in the finite frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part data in the finite frequency band.
  • the imaginary part dielectric response fitting function; among them, the limited frequency band is located between the lowest frequency and the highest frequency;
  • the epitaxial frequency band fitting function unit is used to obtain the real part dielectric response fitting function in the real low epitaxial frequency band according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, and according to the real part corresponding to the highest frequency Data and the real part data corresponding to the second highest frequency obtain the real part dielectric response fitting function in the real part high epitaxial frequency band; obtain the imaginary part low epitaxial frequency band according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency.
  • the imaginary part dielectric response fitting function in the inner part, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band is obtained; the real part low epitaxial frequency band Located between the lowest frequency and the lowest epitaxial frequency of the real part, the high epitaxial frequency band of the real part is between the highest frequency and the highest epitaxial frequency of the real part, the low epi
  • the real part iterative unit is used to perform the following iterative processing:
  • Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral
  • the imaginary part iterative unit is used to perform the following iterative processing:
  • Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral
  • the dielectric state unit is used to obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
  • the embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and the processor implements the steps of the dielectric state analysis method when the computer program is executed.
  • the embodiment of the present invention also provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the steps of the dielectric state analysis method are implemented.
  • the dielectric state analysis method, system, computer, and storage medium of the embodiments of the present invention iteratively calculate the simulated imaginary part integral according to the real part interpolation step, until the simulated imaginary part integral is less than the preset resolution, and the simulated imaginary part integral is compared with the actual imaginary part
  • the integral is compared to obtain the conductance information and the real part information of the polarization process; according to the imaginary part interpolation step, iteratively calculate the simulated real part integral until the simulated real part integral is less than the preset resolution, and the simulated real part integral is compared with the actual real part integral Comparing to obtain the capacitance information and the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state and avoid misjudgment. , Missed judgments and wrong judgments, saving the cost of operation and maintenance of power equipment.
  • FIG. 1 is a flowchart of a method for analyzing the state of a dielectric in an embodiment of the present invention
  • Figure 3 is a flowchart of calculating imaginary singular point frequency band integrals in an embodiment of the present invention.
  • FIG. 5 is a schematic diagram of imaginary data in the first embodiment of the present invention.
  • FIG. 6 is a schematic diagram of imaginary data in the second embodiment of the present invention.
  • FIG. 7 is a schematic diagram of comparison of polarizability in an embodiment of the present invention.
  • Fig. 8 is a schematic diagram of comparison of polarizability in an embodiment of the present invention.
  • FIG. 9 is a schematic diagram of the simulated real part integral and the simulated imaginary part integral of a complex capacitor in an embodiment of the present invention.
  • FIG. 10 is a schematic diagram of the real part information and the imaginary part information of the polarization process of a complex capacitor in an embodiment of the present invention
  • FIG. 11 is a schematic diagram of conductance information and infinite frequency capacitance information of a complex capacitor in an embodiment of the present invention.
  • FIG. 12 is a schematic diagram of dielectric response parameters of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention.
  • FIG. 13 is a schematic diagram of polarization process information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention
  • FIG. 14 is a schematic diagram of conductance information and infinite frequency capacitance information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention.
  • 15 is a schematic diagram showing the comparison of the real part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitor in the embodiment of the present invention.
  • 16 is a schematic diagram of comparison of the imaginary part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention.
  • Fig. 17 is a structural block diagram of a dielectric state analysis system in an embodiment of the present invention.
  • the embodiments of the present invention can be implemented as a system, device, device, method, or computer program product. Therefore, the present disclosure may be specifically implemented in the following forms, namely: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.
  • the embodiment of the present invention provides a dielectric state analysis method, Accurate analysis and judgment of the dielectric state can avoid misjudgment, missed judgment and misjudgment, and save the cost of operation and maintenance of power equipment.
  • FIG. 1 is a flowchart of a method for analyzing the state of a dielectric in an embodiment of the present invention. As shown in Figure 1, the dielectric state analysis method includes:
  • S101 Acquire real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies.
  • the dielectric response parameter can be complex capacitance, dielectric constant or polarizability.
  • S102 Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency.
  • S103 Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting in the limited frequency band function.
  • the limited frequency band is located between the lowest frequency and the highest frequency.
  • the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency obtain the real part dielectric response fitting function in the real part low epitaxial frequency band, and according to the real part data corresponding to the highest frequency and the second highest frequency corresponding
  • the real part data obtains the real part dielectric response fitting function in the real part high epitaxial frequency band
  • the imaginary part dielectric response in the imaginary part low epitaxial frequency band is obtained according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency Fitting function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band
  • the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest Among the epitaxial frequencies, the real high epitaxial frequency band is located between the highest frequency and the real highest epitaxial frequency, the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of
  • S105 Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; where the singular point is located between the highest frequency and the lowest frequency.
  • the same part of the simulated imaginary part integral and the actual imaginary part integral can be used as the real part information of the polarization process, and the different part of the simulated imaginary part integral and the actual imaginary part integral can be used as the conductance information.
  • the real interpolation step size is updated by the following formula:
  • ln( ⁇ n ) is the real part interpolation step in the nth iteration
  • ln( ⁇ n+1 ) is the real part interpolation step in the n+1th iteration.
  • S112 Determine whether the analog real part integral is less than the preset resolution.
  • the same part of the simulated real part integral and the actual real part integral can be used as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral can be used as the capacitance information.
  • the imaginary part interpolation step size is updated by the following formula:
  • ln( ⁇ ' n ) is the imaginary part interpolation step in the nth iteration
  • ln( ⁇ ' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
  • S115 Obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
  • the capacitance information is infinite frequency capacitance information.
  • the execution subject of the dielectric state analysis method shown in FIG. 1 may be a computer.
  • the dielectric state analysis method of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step, until the simulated imaginary part integral is less than the preset resolution, and the simulated imaginary part integral is compared with the actual The imaginary part integral is compared to obtain the conductance information and the real part information of the polarization process; according to the imaginary part interpolation step, the simulated real part integral is iteratively calculated until the simulated real part integral is less than the preset resolution, and the simulated real part integral is compared with the actual real part.
  • Part integrals are compared to obtain the capacitance information and the imaginary part information of the polarization process.
  • the dielectric state can be obtained, which can accurately analyze and judge the dielectric state to avoid causing Misjudgment, omission and misjudgment saves the cost of operation and maintenance of power equipment.
  • S106 includes:
  • calculating the real part truncated frequency domain integral includes:
  • the frequency domain integral of the real part truncated is calculated.
  • the truncated frequency domain integral on the real part is calculated by the following formula:
  • F 1 is the truncated frequency domain integral on the real part
  • ⁇ '( ⁇ Ext-H ) is the real part data corresponding to the highest epitaxial frequency of the real part
  • ⁇ S is the singularity point
  • ⁇ Ext-H is the highest epitaxial frequency of the real part.
  • the lowest epitaxial frequency of the real part is calculated.
  • the frequency domain integral with the lower truncated real part is calculated by the following formula:
  • F 2 is the lower truncated frequency domain integral of the real part
  • ⁇ '( ⁇ Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part
  • ⁇ Ext-L is the lowest epitaxial frequency of the real part.
  • the upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
  • the real part dielectric response fitting function in the limited frequency band calculates the real part extension frequency band integral.
  • the simulation imaginary part integration is calculated.
  • Fig. 2 is a flow chart of calculating real part singular point frequency band integral in an embodiment of the present invention. As shown in Figure 2, the calculation of the real part singular point frequency band integral includes:
  • S201 Determine according to the real data corresponding to the upper bound of the real singular point, the upper bound of the real singular point, the real data corresponding to the singular point, the real data corresponding to the singular point, the lower bound of the real singular point, and the real data corresponding to the lower bound of the real singular point.
  • the real part dielectric response fitting function in the real part low epitaxial frequency band and the real part dielectric response fitting function in the real part high epitaxial frequency band are linear functions, it can be based on the upper bound of the real part singular point , Real part data corresponding to the upper bound of the real part singular point, singular point and real part data corresponding to the singular point determine the first singular point real part coefficient and the second singular point real part coefficient; according to the real part singular point lower bound, real part singularity The real part data corresponding to the lower bound of the point, the singular point and the real part data corresponding to the singular point determine the third singular point real part coefficient and the fourth singular point real part coefficient.
  • S202 Calculate the upper band integral of the real singular point according to the real coefficient of the first singular point, the real coefficient of the second singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point.
  • the frequency band integral on the singular point of the real part is calculated by the following formula:
  • E 1 is the frequency band integral on the real singular point
  • a is the real coefficient of the first singular point
  • b is the real coefficient of the second singular point
  • ⁇ S is the singular point
  • ⁇ SH is the upper bound of the real singular point
  • ⁇ SL is the lower bound of the singular point of the real part.
  • S203 Calculate the lower band integral of the real singular point according to the real part coefficient of the third singular point, the real part coefficient of the fourth singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point.
  • the real part singular point lower band integral is calculated by the following formula:
  • E 2 is the frequency band integral of the real singular point
  • c is the real coefficient of the third singular point
  • d is the real coefficient of the fourth singular point.
  • S111 includes:
  • the imaginary part dielectric response fitting function in the limited frequency band calculates the imaginary part extension frequency band integral.
  • Fig. 3 is a flowchart of calculating the imaginary part singular point frequency band integral in an embodiment of the present invention. As shown in Figure 3, calculating the imaginary part singular point frequency band integral includes:
  • the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band and the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band are both linear functions, it can be based on the upper bound of the imaginary part singular point .
  • the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point and the imaginary part data corresponding to the singular point determine the first singular point imaginary part coefficient and the second singular point imaginary part coefficient; according to the lower bound of the imaginary part singular point, the imaginary part singularity
  • the imaginary part data corresponding to the lower bound of the point, the singular point and the imaginary part data corresponding to the singular point determine the third singular point imaginary part coefficient and the fourth singular point imaginary part coefficient.
  • S302 Calculate the upper band integral of the imaginary part singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point.
  • the frequency band integral on the imaginary singular point is calculated by the following formula:
  • E '1 is the imaginary part of the integral bands singular point
  • a' is a first singular point of the imaginary part of the coefficient
  • b ' is the imaginary part of the second singular point coefficients
  • ⁇ 'SL Is the lower bound of the imaginary part of the singular point.
  • S303 Calculate the lower band integral of the imaginary part singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point.
  • the lower band integral of the imaginary part singular point is calculated by the following formula:
  • E '2 is a lower frequency band integrating the imaginary part of a singular point
  • c' is the coefficient of the imaginary part of the third singular point
  • d ' is the imaginary part of a singular point of the fourth coefficient.
  • Fig. 4 is a schematic diagram of real data in an embodiment of the present invention.
  • Fig. 5 is a schematic diagram of imaginary part data in the first embodiment of the present invention.
  • Fig. 6 is a schematic diagram of imaginary part data in the second embodiment of the present invention.
  • the abscissas of Fig. 4 to Fig. 6 are the logarithm of frequency (log (frequency)), and the ordinates are the logarithm of complex capacitance/dielectric constant/polarizability (dielectric response parameter) (log(complex capacitance/ Dielectric constant/polarizability)).
  • the specific embodiments of the present invention are as follows:
  • the actual measured dielectric response frequency band is [ ⁇ L , ⁇ H ].
  • ⁇ L is the lowest frequency
  • the actual measured dielectric response frequency band is [ ⁇ L , ⁇ H ].
  • ⁇ L is the lowest frequency
  • ⁇ H is the highest frequency.
  • the real part dielectric response fitting function in the low epitaxial band of the real part can be obtained according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second-lowest frequency, and according to the real part corresponding to the highest frequency
  • the real part data corresponding to the data and the second-highest frequency obtains the real part dielectric response fitting function in the real part high epitaxial frequency band.
  • the real part dielectric response fitting function in the real part low epitaxial frequency band and the real part dielectric response fitting function in the real part high epitaxial frequency band are both linear functions.
  • the real low epitaxial frequency band [ ⁇ Ext-L , ⁇ L ] is located between the lowest frequency ( ⁇ L ) and the lowest real epitaxial frequency ( ⁇ Ext-L ), and the real high epitaxial frequency band [ ⁇ H , ⁇ Ext-H ] Located between the highest frequency ( ⁇ H ) and the highest epitaxial frequency of the real part ( ⁇ Ext-H ).
  • the method of truncation is used for assignment.
  • the real part data are all assigned the value ⁇ '( ⁇ Ext-L ), and in the frequency band higher than ⁇ Ext-H , the real part data are all assigned the value ⁇ '( ⁇ Ext-H ).
  • Figure 5 is a schematic diagram of the imaginary part data when the characteristic frequency is not in the epitaxial frequency band (the imaginary part low epitaxial frequency band and the imaginary part high epitaxial frequency band), and Figure 6 is the characteristic frequency when the characteristic frequency is in the epitaxial frequency band (imaginary part low epitaxial frequency band and the imaginary part high epitaxial frequency band) Schematic diagram of the imaginary part of the data.
  • the imaginary part dielectric response fitting function in the imaginary low epitaxial frequency band can be obtained according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency.
  • Part data and the imaginary part data corresponding to the second highest frequency obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band.
  • the imaginary dielectric response fitting function in the imaginary low epitaxial frequency band and the imaginary dielectric response fitting function in the imaginary high epitaxial frequency band are linear functions.
  • the imaginary low-epitaxial frequency band [ ⁇ ' Ext-L , ⁇ L ] is located between the lowest frequency ( ⁇ L ) and the lowest imaginary-part epitaxial frequency ( ⁇ ' Ext-L ), and the imaginary high-epitaxial frequency band [ ⁇ H , ⁇ ' Ext -H ] is located between the highest frequency ( ⁇ H ) and the highest epitaxial frequency of the imaginary part ( ⁇ ' Ext-H ).
  • the method of truncation is used for assignment. That is, in the frequency band lower than ⁇ Ext-L and in the frequency band higher than ⁇ Ext-H , the imaginary part data is assigned a value of 0.
  • the relationship between the upper bound of the real part of the singular point and the singular point is as follows:
  • ln( ⁇ ) is the real part interpolation step.
  • the frequency band integral on the singular point of the real part is the integral of the frequency band [ ⁇ S , ⁇ SH ] on the singular point of the real part.
  • the real part singular point lower frequency band integral is the real part singular point lower frequency band [ ⁇ SL , ⁇ S ] integral.
  • the upper truncated frequency domain integral of the real part is the integral of the upper truncated frequency domain (0, ⁇ Ext-L );
  • the lower truncated frequency domain integral of the real part is the integral of the lower truncated frequency domain of the real part ( ⁇ Ext-H , ⁇ ).
  • the real part dielectric response fitting function in the limited frequency band calculates the real part of the extension frequency band integral.
  • the real part dielectric response fitting function in the limited frequency band [ ⁇ L , ⁇ H ] and the real part dielectric response fitting in the real low epitaxial frequency band [ ⁇ Ext-L , ⁇ L ] have been determined Function and real part of the high-extension frequency band [ ⁇ H , ⁇ Ext-H ] in the real part of the dielectric response fitting function, so the integral of the real part of the epitaxial frequency band [ ⁇ Ext-L , ⁇ SL) can be calculated according to the Simpson formula.
  • the relationship between the upper bound of the imaginary part of the singular point and the singular point is as follows:
  • ln( ⁇ ') is the imaginary part interpolation step.
  • the frequency band integral on the imaginary singular point is the integral of the frequency band [ ⁇ S , ⁇ ' SH ] on the imaginary singular point.
  • the frequency band integration under the imaginary part singular point is the integration of the frequency band under the imaginary part singular point [ ⁇ ' SL , ⁇ S ].
  • the imaginary part dielectric response fitting function in the limited frequency band the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band, the imaginary part The lowest extension frequency, the highest extension frequency of the imaginary part, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point, calculate the imaginary part extension frequency band integral.
  • the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral are added together to obtain the analog real part integral.
  • the imaginary part interpolation step is updated, and step 13 is executed again.
  • the resolution is less than the preset resolution, the same part of the simulated real part integral and the actual real part integral is used as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral is used as the capacitance information.
  • Fig. 7 is a schematic diagram of comparison of polarizability in an embodiment of the present invention.
  • the abscissa is the frequency and the unit is Hz; the ordinate is the polarization rate.
  • ⁇ '( ⁇ ) is the integral of the actual imaginary part of the polarizability
  • ⁇ ”( ⁇ ) is the integral of the simulated imaginary part of the polarizability.
  • replace ⁇ ” ( ⁇ ) is shifted vertically downward by two orders of magnitude. It can be seen from Fig. 7 that the actual imaginary part integral of the polarizability is consistent with the simulated imaginary part integral of the polarizability.
  • Fig. 8 is a schematic diagram of comparison of polarizability in an embodiment of the present invention.
  • the abscissa is the frequency and the unit is Hz; the ordinate is the polarization rate.
  • ⁇ '( ⁇ ) is the actual integral of the real part of the polarizability
  • ⁇ ”( ⁇ ) is the integral of the simulated real part of the susceptibility.
  • add ⁇ ” ( ⁇ ) is shifted vertically downward by two orders of magnitude. It can be seen from Fig. 8 that the actual real part integral of the polarizability is consistent with the simulated real part integral of the polarizability.
  • Fig. 9 is a schematic diagram of the simulated real part integral and the simulated imaginary part integral of a complex capacitor in an embodiment of the present invention.
  • FIG. 10 is a schematic diagram of the real part information and the imaginary part information of the polarization process of a complex capacitor in an embodiment of the present invention.
  • FIG. 11 is a schematic diagram of conductance information and infinite frequency capacitance information of a complex capacitor in an embodiment of the present invention.
  • the abscissas of Fig. 9-11 are all frequency, and the unit is Hz; the ordinates are all capacitance, and the unit is pF.
  • C'( ⁇ ) in Figure 9 is the integration of the simulated real part of the complex capacitor
  • C”( ⁇ ) is the integration of the simulated imaginary part of the complex capacitor
  • ⁇ 1'( ⁇ ) in Figure 10 is the real part of the polarization process of the complex capacitor Information
  • ⁇ 1”( ⁇ ) is the imaginary part information of the polarization process of the complex capacitor
  • ⁇ in Fig. 11 is the conductance information of the complex capacitor
  • C ⁇ is the infinite frequency capacitance information of the complex capacitor.
  • the present invention can be applied in the dielectric response analysis of various power equipment.
  • the following takes the most common high-temperature vulcanized silicone rubber in the field of external insulation of power systems as an example to illustrate the application of the present invention in detail.
  • the main components of high-temperature vulcanized silicone rubber are polydimethylsiloxane, nano-silica filler and nano-aluminum hydroxide filler.
  • the test temperature is 80°C, and the measurement frequency range is 10 -4 Hz to 10 3 Hz.
  • Fig. 12 is a schematic diagram of dielectric response parameters of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention.
  • Fig. 13 is a schematic diagram of polarization process information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention.
  • Fig. 14 is a schematic diagram of conductivity information and infinite frequency capacitance information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention.
  • FIG. 15 is a schematic diagram showing the comparison of the real part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention.
  • FIG. 16 is a schematic diagram of comparison of imaginary part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention.
  • the abscissas of Figure 12-16 are all frequency, the unit is Hz; the ordinates are all capacitance, the unit is pF.
  • C'( ⁇ ) in Figure 12 is the actual real integral of the dielectric response parameter of the complex capacitor
  • C"( ⁇ ) is the actual imaginary integral of the dielectric response parameter of the complex capacitor.
  • the actual implementation of the integral portion remains substantially unchanged in the measured frequency band, the actual presence of a significant imaginary part of the integral in the relaxation peak 10 0 Hz ⁇ 10 3 Hz band, the 10 -4 Hz ⁇ 10 0 Hz band obvious Conductance process (that is, the relationship between the imaginary part and frequency is -1 power).
  • ⁇ 1'( ⁇ ) in Figure 13 is the real part information of the polarization process of the complex capacitor
  • ⁇ 1"( ⁇ ) is the imaginary part information of the polarization process of the complex capacitor.
  • the real part information of the polarization process is at 10 -4 Hz ⁇ 10 -2 Hz frequency band and 10 1 Hz ⁇ 10 3 Hz frequency band have obvious dispersion phenomenon.
  • the imaginary part information of the polarization process also has dispersion phenomenon in these two frequency bands.
  • polarization The value of the real part information of the process and the value of the imaginary part information of the polarization process increase as the frequency decreases in the 10 -4 Hz ⁇ 10 -2 Hz frequency band, which is a characteristic of the typical low-frequency dispersion process. 1 Hz ⁇ 10 3 Hz frequency range, the imaginary part of the information occurs during polarization relaxation peak and the 10 2 Hz ⁇ 10 3 Hz band can be observed polarization value and the imaginary part of the process of the real information part of the polarization process
  • the parallel phenomenon of the numerical value of the information in the double exponential coordinate, that is, the universal exponential law of polarizability is observed.
  • ⁇ in Figure 14 is the conductance information of the complex capacitor
  • C ⁇ is the infinite frequency capacitance information of the complex capacitor.
  • C'( ⁇ ) in Figure 15 is the actual real part integral of the dielectric response parameter of the complex capacitor
  • ⁇ 1'( ⁇ ) is the real part information of the polarization process of the complex capacitor.
  • C"( ⁇ ) in Fig. 16 is the actual imaginary part integral of the dielectric response parameter of the complex capacitor
  • ⁇ 1"( ⁇ ) is the simulated imaginary part integral of the dielectric response parameter of the polarizability. It can be found that in the high frequency band (10 1 Hz ⁇ 10 3 Hz frequency band), the actual imaginary part integral is consistent with the simulated imaginary part integral (considering the shape coefficient). 10 0 Hz in frequency bands below, the actual score was significantly larger than the imaginary part analog integrated imaginary part (considering the shape factor), is the difference between the two conductivity information. Figure 16 clearly shows that the conductance information mainly affects the measurement results at the low frequency of the imaginary part of the complex capacitance.
  • the dielectric state analysis method of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral, Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain Capacitance information and imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, avoiding misjudgments, missed judgments, and Misjudgment saves the cost of operation and maintenance of power equipment.
  • the embodiment of the present invention also provides a dielectric state analysis system. Since the principle of the system to solve the problem is similar to the dielectric state analysis method, the implementation of the system can refer to the implementation of the method, and the repetition will not be repeated. .
  • Fig. 17 is a structural block diagram of a dielectric state analysis system in an embodiment of the present invention. As shown in Figure 17, the dielectric state analysis system includes:
  • An acquiring unit for acquiring real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies
  • the determining unit is used to determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
  • the finite frequency band fitting function unit is used to fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in the finite frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part data in the finite frequency band.
  • the imaginary part dielectric response fitting function; among them, the limited frequency band is located between the lowest frequency and the highest frequency;
  • the epitaxial frequency band fitting function unit is used to obtain the real part dielectric response fitting function in the real low epitaxial frequency band according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, and according to the real part corresponding to the highest frequency Data and the real part data corresponding to the second highest frequency obtain the real part dielectric response fitting function in the real part high epitaxial frequency band; obtain the imaginary part low epitaxial frequency band according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency.
  • the imaginary part dielectric response fitting function in the inner part, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band is obtained; the real part low epitaxial frequency band Located between the lowest frequency and the lowest epitaxial frequency of the real part, the high epitaxial frequency band of the real part is between the highest frequency and the highest epitaxial frequency of the real part, the low epi
  • the real part iterative unit is used to perform the following iterative processing:
  • Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral
  • the imaginary part iterative unit is used to perform the following iterative processing:
  • Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral
  • the dielectric state unit is used to obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
  • the real part iteration unit is specifically used for:
  • the real part dielectric response fitting function in the limited frequency band calculates the real part extension frequency band integral;
  • the imaginary part iteration unit is specifically used for:
  • the imaginary part dielectric response fitting function in the limited frequency band calculates the imaginary part extension frequency band integral;
  • the real part iteration unit is specifically used for:
  • the real data corresponding to the upper bound of the real singular point determines the first Singular point real part coefficients, second singular point real part coefficients, third singular point real part coefficients and fourth singular point real part coefficients;
  • the imaginary part iteration unit is specifically used for:
  • the frequency band integral on the singular point of the real part is calculated by the following formula:
  • E 1 is the frequency band integral on the real singular point
  • a is the real coefficient of the first singular point
  • b is the real coefficient of the second singular point
  • ⁇ S is the singular point
  • ⁇ SH is the upper bound of the real singular point
  • ⁇ SL is the lower bound of the singular point of the real part
  • E 2 is the frequency band integral of the real part of the singular point
  • c is the real part coefficient of the third singular point
  • d is the real part coefficient of the fourth singular point
  • E '1 is the imaginary part of the integral bands singular point
  • a' is a first singular point of the imaginary part of the coefficient
  • b ' is the imaginary part of the second singular point coefficients
  • ⁇ 'SL Is the lower bound of the imaginary part singular point
  • E '2 is a lower frequency band integrating the imaginary part of a singular point
  • c' is the coefficient of the imaginary part of the third singular point
  • d ' is the imaginary part of a singular point of the fourth coefficient.
  • the real part iteration unit is specifically used for:
  • the upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
  • the truncated frequency domain integral on the real part is calculated by the following formula:
  • F 1 is the truncated frequency domain integral on the real part
  • ⁇ '( ⁇ Ext-H ) is the real part data corresponding to the highest extension frequency of the real part
  • ⁇ S is the singularity point
  • ⁇ Ext-H is the highest extension frequency of the real part
  • F 2 is the lower truncated frequency domain integral of the real part
  • ⁇ '( ⁇ Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part
  • ⁇ Ext-L is the lowest epitaxial frequency of the real part.
  • the real interpolation step size is updated by the following formula:
  • ln( ⁇ n ) is the real part interpolation step in the nth iteration
  • ln( ⁇ n+1 ) is the real part interpolation step in the n+1th iteration
  • ln( ⁇ ' n ) is the imaginary part interpolation step in the nth iteration
  • ln( ⁇ ' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
  • the real part iteration unit is specifically used for:
  • the imaginary part iteration unit is specifically used for:
  • the same part of the simulated real part integral and the actual real part integral is regarded as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral is regarded as the capacitance information.
  • the dielectric state analysis system of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral, Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain Capacitance information and imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, avoiding misjudgments, missed judgments, and Misjudgment saves the cost of operation and maintenance of power equipment.
  • the embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and running on the processor.
  • the processor can implement all or part of the dielectric state analysis method when the computer program is executed. For example, when a processor executes a computer program, the following can be realized:
  • the limited frequency band is located between the lowest frequency and the highest frequency;
  • the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency
  • the data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency
  • the real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency
  • the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part
  • Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral
  • Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral
  • the conductance information the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
  • the computer device of the embodiment of the present invention iteratively calculates the simulated imaginary integral according to the real interpolation step, until the simulated imaginary integral is less than the preset resolution, and compares the simulated imaginary integral with the actual imaginary integral to obtain the conductance Information and real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step, until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain capacitance information And the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, and avoid misjudgment, missed judgment and wrong judgment. , Saving the cost of power equipment operation and maintenance.
  • the embodiment of the present invention also provides a computer-readable storage medium on which a computer program is stored.
  • a computer program When the computer program is executed by a processor, all or part of the content of the dielectric state analysis method can be realized. For example, when the processor executes the computer program, To achieve the following:
  • the limited frequency band is located between the lowest frequency and the highest frequency;
  • the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency
  • the data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency
  • the real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency
  • the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part
  • Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral
  • Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral
  • the conductance information the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
  • the computer-readable storage medium of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral , Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral, Obtain the capacitance information and the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, and avoid misjudgments and missed judgments. And misjudgment, saving the cost of power equipment operation and maintenance.
  • the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.
  • computer-usable storage media including but not limited to disk storage, CD-ROM, optical storage, etc.
  • These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device.
  • the device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
  • These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment.
  • the instructions provide steps for implementing the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement Of Resistance Or Impedance (AREA)
  • Design And Manufacture Of Integrated Circuits (AREA)

Abstract

Disclosed are a dielectric state analysis method and system, a computer and a storage medium. The method comprises: iteratively calculating a simulated imaginary part integral according to a real part interpolation step size until the simulated imaginary part integral is less than a preset resolution, and comparing the simulated imaginary part integral with an actual imaginary part integral to obtain conductance information and polarization process real part information; iteratively calculating a simulated real part integral according to an imaginary part interpolation step size until the simulated real part integral is less than the preset resolution, and comparing the simulated real part integral with an actual real part integral to obtain capacitance information and polarization process imaginary part information; and finally obtaining a dielectric state according to the conductance information, the polarization process real part information, the capacitance information and the polarization process imaginary part information. In this way, a dielectric state can be accurately analyzed and determined, incorrect identification, omitted identification and erroneous identification are avoided, and the operation and maintenance costs of a power device are saved.

Description

电介质状态分析方法、系统、计算机及存储介质Method, system, computer and storage medium for analyzing state of dielectric 技术领域Technical field
本发明涉及电介质领域,具体地,涉及一种电介质状态分析方法、系统、计算机及存储介质。The present invention relates to the field of dielectrics, and in particular to a method, system, computer and storage medium for analyzing the state of dielectrics.
背景技术Background technique
在电能的生产、传输及应用过程中,所有的电气设备中,都需要电介质材料的支持,如在发电机的匝间绝缘系统,变压器的绝缘系统,电缆的绝缘系统,绝缘子的绝缘系统等等。如何表征和研究各类电介质材料的绝缘状态,是电力设备运行维护中的重要研究课题。In the production, transmission and application of electric energy, all electrical equipment requires the support of dielectric materials, such as the turn-to-turn insulation system of generators, the insulation system of transformers, the insulation system of cables, the insulation system of insulators, etc. . How to characterize and study the insulation state of various dielectric materials is an important research topic in the operation and maintenance of power equipment.
已有的研究及工程实际应用中,尚无法实现电介质\电力设备的介电响应测量结果的解耦分析,因而就无法获得电介质\电力设备的电导过程、极化过程和无穷频率电容过程,这就导致无法准确判断电介质\电力设备的状态,给电力系统的电力设备状态诊断带来诸多不便,容易造成误判、漏判和错判。In the existing research and practical engineering applications, the decoupling analysis of the measurement results of the dielectric response of the dielectric/power equipment has not been achieved, so the conductance process, the polarization process and the infinite frequency capacitance process of the dielectric/power equipment cannot be obtained. This leads to the inability to accurately determine the status of the dielectric/power equipment, which brings a lot of inconvenience to the power equipment status diagnosis of the power system, and it is easy to cause misjudgments, missed judgments and misjudgments.
发明内容Summary of the invention
本发明实施例的主要目的在于提供一种电介质状态分析方法、系统、计算机及存储介质,以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。The main purpose of the embodiments of the present invention is to provide a dielectric state analysis method, system, computer, and storage medium to accurately analyze and determine the dielectric state, avoid misjudgment, missed judgment, and misjudgment, and save the cost of power equipment operation and maintenance.
为了实现上述目的,本发明实施例提供一种电介质状态分析方法,包括:In order to achieve the foregoing objective, an embodiment of the present invention provides a dielectric state analysis method, including:
获取多个频率对应的介电响应参数的实部数据和虚部数据;Acquire real data and imaginary data of dielectric response parameters corresponding to multiple frequencies;
根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;有限频段位于最低频率与最高频率之间;Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting function in the limited frequency band; The limited frequency band is located between the lowest frequency and the highest frequency;
根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;According to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency The data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency The real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency, the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary high epitaxial frequency band is between the highest frequency and the highest epitaxial frequency of the imaginary part. Between frequencies
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; among them, the singular point is located between the highest frequency and the lowest frequency;
根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part in the low epitaxial frequency band of the real part. Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral;
判断模拟虚部积分是否小于预设分辨率;当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain the conductance information and the real part information of the polarization process, otherwise update the real part interpolation step long;
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part in the low epitaxial frequency band of the imaginary part. Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral;
判断模拟实部积分是否小于预设分辨率;当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain the capacitance information and the imaginary part information of the polarization process, otherwise update the imaginary part interpolation step long;
根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。According to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
本发明实施例还提供一种电介质状态分析系统,包括:The embodiment of the present invention also provides a dielectric state analysis system, including:
获取单元,用于获取多个频率对应的介电响应参数的实部数据和虚部数据;An acquiring unit for acquiring real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies;
确定单元,用于根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;The determining unit is used to determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
有限频段拟合函数单元,用于拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;其中,有限频段位于最低频率与最高频率之间;The finite frequency band fitting function unit is used to fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in the finite frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part data in the finite frequency band. The imaginary part dielectric response fitting function; among them, the limited frequency band is located between the lowest frequency and the highest frequency;
外延频段拟合函数单元,用于根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;The epitaxial frequency band fitting function unit is used to obtain the real part dielectric response fitting function in the real low epitaxial frequency band according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, and according to the real part corresponding to the highest frequency Data and the real part data corresponding to the second highest frequency obtain the real part dielectric response fitting function in the real part high epitaxial frequency band; obtain the imaginary part low epitaxial frequency band according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency The imaginary part dielectric response fitting function in the inner part, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band is obtained; the real part low epitaxial frequency band Located between the lowest frequency and the lowest epitaxial frequency of the real part, the high epitaxial frequency band of the real part is between the highest frequency and the highest epitaxial frequency of the real part, the low epitaxial frequency band of the imaginary part is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary part is high epitaxial frequency The frequency band is located between the highest frequency and the highest epitaxial frequency of the imaginary part;
实部迭代单元,用于执行如下迭代处理:The real part iterative unit is used to perform the following iterative processing:
根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; among them, the singular point is located between the highest frequency and the lowest frequency;
根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part in the low epitaxial frequency band of the real part. Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral;
判断模拟虚部积分是否小于预设分辨率;当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信 息和极化过程实部信息,否则更新实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain the conductance information and the real part information of the polarization process, otherwise update the real part interpolation step long;
虚部迭代单元,用于执行如下迭代处理:The imaginary part iterative unit is used to perform the following iterative processing:
根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part in the low epitaxial frequency band of the imaginary part. Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral;
判断模拟实部积分是否小于预设分辨率;当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain the capacitance information and the imaginary part information of the polarization process, otherwise update the imaginary part interpolation step long;
电介质状态单元,用于根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。The dielectric state unit is used to obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时实现所述的电介质状态分析方法的步骤。The embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, and the processor implements the steps of the dielectric state analysis method when the computer program is executed.
本发明实施例还提供一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现所述的电介质状态分析方法的步骤。The embodiment of the present invention also provides a computer-readable storage medium on which a computer program is stored, and when the computer program is executed by a processor, the steps of the dielectric state analysis method are implemented.
本发明实施例的电介质状态分析方法、系统、计算机及存储介质根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。The dielectric state analysis method, system, computer, and storage medium of the embodiments of the present invention iteratively calculate the simulated imaginary part integral according to the real part interpolation step, until the simulated imaginary part integral is less than the preset resolution, and the simulated imaginary part integral is compared with the actual imaginary part The integral is compared to obtain the conductance information and the real part information of the polarization process; according to the imaginary part interpolation step, iteratively calculate the simulated real part integral until the simulated real part integral is less than the preset resolution, and the simulated real part integral is compared with the actual real part integral Comparing to obtain the capacitance information and the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state and avoid misjudgment. , Missed judgments and wrong judgments, saving the cost of operation and maintenance of power equipment.
附图说明Description of the drawings
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。In order to explain the embodiments of the present invention or the technical solutions in the prior art more clearly, the following will briefly introduce the accompanying drawings needed in the description of the embodiments. Obviously, the accompanying drawings in the following description are only of the present invention. For some embodiments, those of ordinary skill in the art can obtain other drawings based on these drawings without creative work.
图1是本发明实施例中电介质状态分析方法的流程图;FIG. 1 is a flowchart of a method for analyzing the state of a dielectric in an embodiment of the present invention;
图2是本发明实施例中计算实部奇异点频段积分的流程图;2 is a flowchart of calculating real part singular point frequency band integral in an embodiment of the present invention;
图3是本发明实施例中计算虚部奇异点频段积分的流程图;Figure 3 is a flowchart of calculating imaginary singular point frequency band integrals in an embodiment of the present invention;
图4是本发明实施例中实部数据示意图;4 is a schematic diagram of real data in an embodiment of the present invention;
图5是本发明第一实施例中虚部数据示意图;5 is a schematic diagram of imaginary data in the first embodiment of the present invention;
图6是本发明第二实施例中虚部数据示意图;6 is a schematic diagram of imaginary data in the second embodiment of the present invention;
图7是本发明实施例中极化率的对比示意图;FIG. 7 is a schematic diagram of comparison of polarizability in an embodiment of the present invention;
图8是本发明实施例中极化率的对比示意图;Fig. 8 is a schematic diagram of comparison of polarizability in an embodiment of the present invention;
图9是本发明实施例中复电容的模拟实部积分和模拟虚部积分的示意图;9 is a schematic diagram of the simulated real part integral and the simulated imaginary part integral of a complex capacitor in an embodiment of the present invention;
图10是本发明实施例中复电容的极化过程实部信息和极化过程虚部信息的示意图;10 is a schematic diagram of the real part information and the imaginary part information of the polarization process of a complex capacitor in an embodiment of the present invention;
图11是本发明实施例中复电容的电导信息和无穷频率电容信息的示意图;FIG. 11 is a schematic diagram of conductance information and infinite frequency capacitance information of a complex capacitor in an embodiment of the present invention;
图12是本发明实施例中高温硫化硅橡胶复电容的介电响应参数示意图;12 is a schematic diagram of dielectric response parameters of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention;
图13是本发明实施例中高温硫化硅橡胶复电容的极化过程信息示意图;13 is a schematic diagram of polarization process information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention;
图14是本发明实施例中高温硫化硅橡胶复电容的电导信息和无穷频率电容信息的示意图;FIG. 14 is a schematic diagram of conductance information and infinite frequency capacitance information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention;
图15是本发明实施例中高温硫化硅橡胶复电容的极化率的实部数据对比示意图;15 is a schematic diagram showing the comparison of the real part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitor in the embodiment of the present invention;
图16是本发明实施例中高温硫化硅橡胶复电容的极化率的虚部数据对比示意图;16 is a schematic diagram of comparison of the imaginary part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention;
图17是本发明实施例中电介质状态分析系统的结构框图。Fig. 17 is a structural block diagram of a dielectric state analysis system in an embodiment of the present invention.
具体实施方式Detailed ways
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of the present invention.
本领域技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。Those skilled in the art know that the embodiments of the present invention can be implemented as a system, device, device, method, or computer program product. Therefore, the present disclosure may be specifically implemented in the following forms, namely: complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.
鉴于现有技术无法准确判断电介质\电力设备的状态,给电力系统的电力设备状态诊断带来诸多不便,容易造成误判、漏判和错判,本发明实施例提供一种电介质状态分析方法,以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。以下结合附图对本发明进行详细说明。In view of the fact that the prior art cannot accurately determine the status of the dielectric/power equipment, which brings a lot of inconvenience to the power equipment status diagnosis of the power system, and is likely to cause misjudgment, missed judgment and misjudgment, the embodiment of the present invention provides a dielectric state analysis method, Accurate analysis and judgment of the dielectric state can avoid misjudgment, missed judgment and misjudgment, and save the cost of operation and maintenance of power equipment. The present invention will be described in detail below with reference to the accompanying drawings.
图1是本发明实施例中电介质状态分析方法的流程图。如图1所示,电介质状态分析方法包括:FIG. 1 is a flowchart of a method for analyzing the state of a dielectric in an embodiment of the present invention. As shown in Figure 1, the dielectric state analysis method includes:
S101:获取多个频率对应的介电响应参数的实部数据和虚部数据。S101: Acquire real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies.
其中,介电响应参数可以为复电容、介电常数或极化率。Among them, the dielectric response parameter can be complex capacitance, dielectric constant or polarizability.
S102:根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率。S102: Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency.
S103:拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数。其中,有限频段位于最低频率与最高频率之间。S103: Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting in the limited frequency band function. Among them, the limited frequency band is located between the lowest frequency and the highest frequency.
S104:根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高 频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;S104: According to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, obtain the real part dielectric response fitting function in the real part low epitaxial frequency band, and according to the real part data corresponding to the highest frequency and the second highest frequency corresponding The real part data obtains the real part dielectric response fitting function in the real part high epitaxial frequency band; the imaginary part dielectric response in the imaginary part low epitaxial frequency band is obtained according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency Fitting function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest Among the epitaxial frequencies, the real high epitaxial frequency band is located between the highest frequency and the real highest epitaxial frequency, the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary high epitaxial frequency band is located between the highest frequency and the imaginary part. Between the highest epitaxial frequencies;
执行如下迭代处理:Perform the following iterative processing:
S105:根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间。S105: Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; where the singular point is located between the highest frequency and the lowest frequency.
S106:根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分。S106: According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part low epitaxial frequency band The real part dielectric response fitting function and the real part dielectric response fitting function in the high epitaxial frequency band of the real part are calculated to simulate the imaginary part integral.
S107:判断模拟虚部积分是否小于预设分辨率。S107: Determine whether the analog imaginary part integral is less than the preset resolution.
S108:当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息。S108: When the resolution is less than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain conductance information and real part information of the polarization process.
具体实施时,可以将模拟虚部积分与实际虚部积分中的相同部分作为极化过程实部信息,将模拟虚部积分与实际虚部积分中的不同部分作为电导信息。In specific implementation, the same part of the simulated imaginary part integral and the actual imaginary part integral can be used as the real part information of the polarization process, and the different part of the simulated imaginary part integral and the actual imaginary part integral can be used as the conductance information.
S109:当大于或等于预设分辨率时,更新实部插值步长。S109: When the resolution is greater than or equal to the preset resolution, update the real part interpolation step.
一实施例中,通过如下公式更新实部插值步长:In one embodiment, the real interpolation step size is updated by the following formula:
Figure PCTCN2020109652-appb-000001
Figure PCTCN2020109652-appb-000001
其中,ln(Δω n)为第n次迭代中的实部插值步长,ln(Δω n+1)为第n+1次迭代中的实部插值步长。 Among them, ln(Δω n ) is the real part interpolation step in the nth iteration, and ln(Δω n+1 ) is the real part interpolation step in the n+1th iteration.
执行如下迭代处理:Perform the following iterative processing:
S110:根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界。S110: Calculate the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point according to the singular point and the interpolation step length of the imaginary part.
S111:根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分。S111: According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part singular point, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part low epitaxial frequency band The imaginary part dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band are calculated to simulate the real part integral.
S112:判断模拟实部积分是否小于预设分辨率。S112: Determine whether the analog real part integral is less than the preset resolution.
S113:当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息。S113: When the resolution is less than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain capacitance information and imaginary part information of the polarization process.
具体实施时,可以将模拟实部积分与实际实部积分中的相同部分作为极化过程虚部信息,将模拟实部积分与实际实部积分中的不同部分作为电容信息。In specific implementation, the same part of the simulated real part integral and the actual real part integral can be used as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral can be used as the capacitance information.
S114:当大于或等于预设分辨率时,更新虚部插值步长。S114: When the resolution is greater than or equal to the preset resolution, update the imaginary part interpolation step.
一实施例中,通过如下公式更新虚部插值步长:In an embodiment, the imaginary part interpolation step size is updated by the following formula:
Figure PCTCN2020109652-appb-000002
Figure PCTCN2020109652-appb-000002
其中,ln(Δω' n)为第n次迭代中的虚部插值步长,ln(Δω' n+1)为第n+1次迭代中的虚部插值步长。 Among them, ln(Δω' n ) is the imaginary part interpolation step in the nth iteration, and ln(Δω' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
S115:根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。S115: Obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
其中,电容信息为无穷频率电容信息。Among them, the capacitance information is infinite frequency capacitance information.
图1所示的电介质状态分析方法的执行主体可以为计算机。由图1所示的流程可知,本发明实施例的电介质状态分析方法根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。The execution subject of the dielectric state analysis method shown in FIG. 1 may be a computer. As can be seen from the process shown in Figure 1, the dielectric state analysis method of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step, until the simulated imaginary part integral is less than the preset resolution, and the simulated imaginary part integral is compared with the actual The imaginary part integral is compared to obtain the conductance information and the real part information of the polarization process; according to the imaginary part interpolation step, the simulated real part integral is iteratively calculated until the simulated real part integral is less than the preset resolution, and the simulated real part integral is compared with the actual real part. Part integrals are compared to obtain the capacitance information and the imaginary part information of the polarization process. Finally, according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state can be obtained, which can accurately analyze and judge the dielectric state to avoid causing Misjudgment, omission and misjudgment saves the cost of operation and maintenance of power equipment.
一实施例中,S106包括:In an embodiment, S106 includes:
根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点、奇异点对应的实部数据、实部奇异点下界和实部奇异点下界对应的实部数据计算实部奇异点频段积分。Calculate the real part based on the real data corresponding to the upper bound of the real singular point and the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point, and the real data corresponding to the lower bound of the real singular point Singularity frequency band integration.
根据实部最低外延频率对应的实部数据、实部最低外延频率、实部最高外延频率对应的实部数据、实部最高外延频率和奇异点计算实部截断频域积分。Calculate the real part truncated frequency domain integral based on the real part data corresponding to the lowest real part epitaxial frequency, the lowest real part epitaxial frequency, the real part data corresponding to the highest real part epitaxial frequency, the highest real part epitaxial frequency and the singular point.
具体实施时,计算实部截断频域积分包括:In specific implementation, calculating the real part truncated frequency domain integral includes:
根据实部最高外延频率对应的实部数据、实部最高外延频率和奇异点计算实部上截断频域积分。According to the real data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part and the singular point, the frequency domain integral of the real part truncated is calculated.
一实施例中,通过如下公式计算实部上截断频域积分:In an embodiment, the truncated frequency domain integral on the real part is calculated by the following formula:
Figure PCTCN2020109652-appb-000003
Figure PCTCN2020109652-appb-000003
其中,F 1为实部上截断频域积分,χ'(ω Ext-H)为实部最高外延频率对应的实部数据,ω S为奇异点,ω Ext-H为实部最高外延频率。 Among them, F 1 is the truncated frequency domain integral on the real part, χ'(ω Ext-H ) is the real part data corresponding to the highest epitaxial frequency of the real part, ω S is the singularity point, and ω Ext-H is the highest epitaxial frequency of the real part.
根据实部最低外延频率对应的实部数据、实部最低外延频率和奇异点计算实部下截断频域积分。According to the real data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part and the singular point, the lower truncated frequency domain integral of the real part is calculated.
一实施例中,通过如下公式计算实部下截断频域积分:In an embodiment, the frequency domain integral with the lower truncated real part is calculated by the following formula:
Figure PCTCN2020109652-appb-000004
Figure PCTCN2020109652-appb-000004
其中,F 2为实部下截断频域积分,χ'(ω Ext-L)为实部最低外延频率对应的实部数据,ω Ext-L为实部最低外延频率。 Among them, F 2 is the lower truncated frequency domain integral of the real part, χ'(ω Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part, and ω Ext-L is the lowest epitaxial frequency of the real part.
将实部上截断频域积分与实部下截断频域积分相加,得到实部截断频域积分。The upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
根据有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数、实部高外延频段内的实部介电响应拟合函数、实部最低外延频率、实部最高外延频率、实部奇异点上界和实部奇异点下界,计算实部外延频段积分。According to the real part dielectric response fitting function in the limited frequency band, the real part dielectric response fitting function in the real part low epitaxial frequency band, the real part dielectric response fitting function in the real part high epitaxial frequency band, and the real part lowest epitaxy Frequency, real part maximum extension frequency, real part singular point upper bound and real part singular point lower bound, calculate the real part extension frequency band integral.
根据实部奇异点频段积分、实部截断频域积分和实部外延频段积分计算模拟虚部积分。According to the real part singular point frequency band integration, real part truncated frequency domain integration and real part extension frequency band integration, the simulation imaginary part integration is calculated.
图2是本发明实施例中计算实部奇异点频段积分的流程图。如图2所示,计算实部奇异点频段积分包括:Fig. 2 is a flow chart of calculating real part singular point frequency band integral in an embodiment of the present invention. As shown in Figure 2, the calculation of the real part singular point frequency band integral includes:
S201:根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点、奇异点对应的实部数据、实部奇异点下界和实部奇异点下界对应的实部数据确定第一奇异点实部系数、第二奇异点实部系数、第三奇异点实部系数和第四奇异点实部系数。S201: Determine according to the real data corresponding to the upper bound of the real singular point, the upper bound of the real singular point, the real data corresponding to the singular point, the real data corresponding to the singular point, the lower bound of the real singular point, and the real data corresponding to the lower bound of the real singular point The first singular point real part coefficient, the second singular point real part coefficient, the third singular point real part coefficient, and the fourth singular point real part coefficient.
具体实施时,因为实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数均为线性函数,因此可以根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点和奇异点对应的实部数据确定第一奇异点实部系数和第二奇异点实部系数;根据实部奇异点下界、实部奇异点下界对应的实部数据、奇异点和奇异点对应的实部数据确定第三奇异点实部系数和第四奇异点实部系数。In specific implementation, because the real part dielectric response fitting function in the real part low epitaxial frequency band and the real part dielectric response fitting function in the real part high epitaxial frequency band are linear functions, it can be based on the upper bound of the real part singular point , Real part data corresponding to the upper bound of the real part singular point, singular point and real part data corresponding to the singular point determine the first singular point real part coefficient and the second singular point real part coefficient; according to the real part singular point lower bound, real part singularity The real part data corresponding to the lower bound of the point, the singular point and the real part data corresponding to the singular point determine the third singular point real part coefficient and the fourth singular point real part coefficient.
S202:根据第一奇异点实部系数、第二奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点上频段积分。S202: Calculate the upper band integral of the real singular point according to the real coefficient of the first singular point, the real coefficient of the second singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point.
一实施例中,通过如下公式计算实部奇异点上频段积分:In an embodiment, the frequency band integral on the singular point of the real part is calculated by the following formula:
Figure PCTCN2020109652-appb-000005
Figure PCTCN2020109652-appb-000005
其中,E 1为实部奇异点上频段积分,a为第一奇异点实部系数,b为第二奇异点实部系数,ω S为奇异点,ω S-H为实部奇异点上界,ω S-L为实部奇异点下界。 Among them, E 1 is the frequency band integral on the real singular point, a is the real coefficient of the first singular point, b is the real coefficient of the second singular point, ω S is the singular point, ω SH is the upper bound of the real singular point, ω SL is the lower bound of the singular point of the real part.
S203:根据第三奇异点实部系数、第四奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点下频段积分。S203: Calculate the lower band integral of the real singular point according to the real part coefficient of the third singular point, the real part coefficient of the fourth singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point.
一实施例中,通过如下公式计算实部奇异点下频段积分:In an embodiment, the real part singular point lower band integral is calculated by the following formula:
Figure PCTCN2020109652-appb-000006
Figure PCTCN2020109652-appb-000006
其中,E 2为实部奇异点下频段积分,c为第三奇异点实部系数,d为第四奇异点实部系数。 Among them, E 2 is the frequency band integral of the real singular point, c is the real coefficient of the third singular point, and d is the real coefficient of the fourth singular point.
S204:将实部奇异点上频段积分与实部奇异点下频段积分相加,得到实部奇异点频段积分。S204: Add the upper frequency band integral of the real part singular point and the lower frequency band integral of the real part singular point to obtain the real part singular point frequency band integral.
一实施例中,S111包括:In an embodiment, S111 includes:
根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点、奇异点对应的虚部数据、虚部奇异点下界和虚部奇异点下界对应的虚部数据计算虚部奇异点频段积分。Calculate the imaginary part based on the imaginary part data corresponding to the upper bound of the imaginary part singular point, the imaginary part corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the imaginary part data corresponding to the lower bound of the imaginary part singular point Singularity frequency band integration.
根据有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数、虚部高外延频段内的虚部介电响应拟合函数、虚部最低外延频率、虚部最高外延频率、虚部奇异点上界和虚部奇异点下界,计算虚部外延频段积分。According to the imaginary part dielectric response fitting function in the limited frequency band, the imaginary part dielectric response fitting function in the imaginary part low extension frequency band, the imaginary part dielectric response fitting function in the imaginary part high extension frequency band, and the imaginary part minimum extension Frequency, imaginary part maximum extension frequency, imaginary part singular point upper bound and imaginary part singular point lower bound, calculate the imaginary part extension frequency band integral.
根据虚部奇异点频段积分和虚部外延频段积分计算模拟实部积分。Calculate the analog real part integral based on the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral.
图3是本发明实施例中计算虚部奇异点频段积分的流程图。如图3所示,计算虚部奇异点频段积分包括:Fig. 3 is a flowchart of calculating the imaginary part singular point frequency band integral in an embodiment of the present invention. As shown in Figure 3, calculating the imaginary part singular point frequency band integral includes:
S301:根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点、奇异点对应的虚部数据、虚部奇异点下界和虚部奇异点下界对应的虚部数据确定第一奇异点虚部系数、第二奇异点虚部系数、第三奇异点虚部系数和第四奇异点虚部系数。S301: Determine according to the imaginary part data corresponding to the upper bound of the imaginary part singular point, the imaginary part corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point, and the imaginary part data corresponding to the lower bound of the imaginary part singular point The first singular point imaginary part coefficient, the second singular point imaginary part coefficient, the third singular point imaginary part coefficient, and the fourth singular point imaginary part coefficient.
具体实施时,因为虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数均为线性函数,因此可以根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点和奇异点对应的虚部数据确定第一奇异点虚部系数和第二奇异点虚部系数;根据虚部奇异点下界、虚部奇异点下界对应的虚部数据、奇异点和奇异点对应的虚部数据确定第三奇异点虚部系数和第四奇异点虚部系数。In specific implementation, because the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band and the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band are both linear functions, it can be based on the upper bound of the imaginary part singular point , The imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point and the imaginary part data corresponding to the singular point determine the first singular point imaginary part coefficient and the second singular point imaginary part coefficient; according to the lower bound of the imaginary part singular point, the imaginary part singularity The imaginary part data corresponding to the lower bound of the point, the singular point and the imaginary part data corresponding to the singular point determine the third singular point imaginary part coefficient and the fourth singular point imaginary part coefficient.
S302:根据第一奇异点虚部系数、第二奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点上频段积分。S302: Calculate the upper band integral of the imaginary part singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point.
一实施例中,通过如下公式计算虚部奇异点上频段积分:In an embodiment, the frequency band integral on the imaginary singular point is calculated by the following formula:
Figure PCTCN2020109652-appb-000007
Figure PCTCN2020109652-appb-000007
其中,E' 1为虚部奇异点上频段积分,a'为第一奇异点虚部系数,b'为第二奇异点虚部系数,ω' S-H为虚部奇异点上界,ω' S-L为虚部奇异点下界。 Wherein, E '1 is the imaginary part of the integral bands singular point, a' is a first singular point of the imaginary part of the coefficient, b 'is the imaginary part of the second singular point coefficients, ω' SH bounded singularity point on the imaginary part, ω 'SL Is the lower bound of the imaginary part of the singular point.
S303:根据第三奇异点虚部系数、第四奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点下频段积分。S303: Calculate the lower band integral of the imaginary part singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point.
一实施例中,通过如下公式计算虚部奇异点下频段积分:In one embodiment, the lower band integral of the imaginary part singular point is calculated by the following formula:
Figure PCTCN2020109652-appb-000008
Figure PCTCN2020109652-appb-000008
其中,E' 2为虚部奇异点下频段积分,c'为第三奇异点虚部系数,d'为第四奇异点虚部系数。 Wherein, E '2 is a lower frequency band integrating the imaginary part of a singular point, c' is the coefficient of the imaginary part of the third singular point, d 'is the imaginary part of a singular point of the fourth coefficient.
S304:将虚部奇异点上频段积分与虚部奇异点下频段积分相加,得到虚部奇异点频段积分。S304: Add the upper frequency band integral of the imaginary part singular point and the lower frequency band integral of the imaginary part singular point to obtain the frequency band integral of the imaginary part singular point.
图4是本发明实施例中实部数据示意图。图5是本发明第一实施例中虚部数据示意图。图6是本发明第二实施例中虚部数据示意图。图4-图6的横坐标均为频率的对数(log(频率)),纵坐标均为复电容/介电常数/极化率(介电响应参数)的对数(log(复电容/介电常数/极化率))。如图4-图6所示,本发明的具体实施例如下:Fig. 4 is a schematic diagram of real data in an embodiment of the present invention. Fig. 5 is a schematic diagram of imaginary part data in the first embodiment of the present invention. Fig. 6 is a schematic diagram of imaginary part data in the second embodiment of the present invention. The abscissas of Fig. 4 to Fig. 6 are the logarithm of frequency (log (frequency)), and the ordinates are the logarithm of complex capacitance/dielectric constant/polarizability (dielectric response parameter) (log(complex capacitance/ Dielectric constant/polarizability)). As shown in Figures 4 to 6, the specific embodiments of the present invention are as follows:
1、获取多个频率对应的介电响应参数的实部数据和虚部数据。根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率。1. Obtain the real part data and imaginary part data of the dielectric response parameters corresponding to multiple frequencies. Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency.
如图4所示,实际测量的介电响应频段(有限频段)为[ω LH]。ω L为最低频率,ω H为最高频率。将数据外延2~3个数量级,即将最低频率外延至实部最低外延频率ω Ext-L,即ω Ext-L=10 -2×ω L或ω Ext-L=10 -3×ω L;将最高频率外延至实部最高外延频率ω Ext-H,即ω Ext-H=10 -2×ω H或ω Ext-H=10 -3×ω HAs shown in Figure 4, the actual measured dielectric response frequency band (limited frequency band) is [ω LH ]. ω L is the lowest frequency, and ω H is the highest frequency. Extend the data by 2 to 3 orders of magnitude, that is, extend the lowest frequency to the lowest epitaxial frequency ω Ext-L of the real part, that is, ω Ext-L = 10 -2 × ω L or ω Ext-L = 10 -3 × ω L ; The highest frequency is extended to the highest extension frequency ω Ext-H of the real part, that is, ω Ext-H =10 -2 ×ω H or ω Ext-H =10 -3 ×ω H.
如图5-图6所示,实际测量的介电响应频段(有限频段)为[ω LH]。ω L为最低频率,ω H为最高频率。将数据外延2~3个数量级,即将最低频率外延至虚部最低外延频率ω' Ext-L,即ω' Ext-L=10 -2×ω L或ω' Ext-L=10 -3×ω L;将最高频率外延至虚部最高外延频率ω' Ext-H,即ω' Ext-H=10 -2×ω H或ω' Ext-H=10 -3×ω HAs shown in Figure 5-6, the actual measured dielectric response frequency band (limited frequency band) is [ω LH ]. ω L is the lowest frequency, and ω H is the highest frequency. Extend the data by 2 to 3 orders of magnitude, that is, extend the lowest frequency to the lowest epitaxial frequency ω'Ext-L of the imaginary part, that is, ω'Ext-L = 10 -2 ×ω L or ω'Ext-L = 10 -3 ×ω L ; Extend the highest frequency to the highest extension frequency ω'Ext-H of the imaginary part, that is, ω'Ext-H = 10 -2 × ω H or ω'Ext-H = 10 -3 × ω H.
2、拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数。其中,有限频段位于最低频率与最高频率之间。实部介电响应拟合函数和虚部介电响应拟合函数均为二次函数。2. Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting in the limited frequency band function. Among them, the limited frequency band is located between the lowest frequency and the highest frequency. Both the real part dielectric response fitting function and the imaginary part dielectric response fitting function are quadratic functions.
3、如图4所示,可以根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数。3. As shown in Figure 4, the real part dielectric response fitting function in the low epitaxial band of the real part can be obtained according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second-lowest frequency, and according to the real part corresponding to the highest frequency The real part data corresponding to the data and the second-highest frequency obtains the real part dielectric response fitting function in the real part high epitaxial frequency band.
其中,实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数均为线性函数。实部低外延频段[ω Ext-LL]位于最低频率(ω L)与实部最低外延频率(ω Ext-L)之间,实部高外延频段[ω HExt-H]位于最高频率(ω H)与实部最高外延频率(ω Ext-H)之间。在[ω Ext-LExt-H]以外的频段,采用截断处理的方法进行赋值。即在低于ω Ext-L的频段,实部数据均赋值为χ'(ω Ext-L),而在高于ω Ext-H的频段,实部数据均赋值为χ'(ω Ext-H)。 Among them, the real part dielectric response fitting function in the real part low epitaxial frequency band and the real part dielectric response fitting function in the real part high epitaxial frequency band are both linear functions. The real low epitaxial frequency band [ω Ext-LL ] is located between the lowest frequency (ω L ) and the lowest real epitaxial frequency (ω Ext-L ), and the real high epitaxial frequency band [ω HExt-H ] Located between the highest frequency (ω H ) and the highest epitaxial frequency of the real part (ω Ext-H ). In the frequency bands other than [ω Ext-LExt-H ], the method of truncation is used for assignment. That is, in the frequency band lower than ω Ext-L , the real part data are all assigned the value χ'(ω Ext-L ), and in the frequency band higher than ω Ext-H , the real part data are all assigned the value χ'(ω Ext-H ).
图5为特征频率不在外延频段(虚部低外延频段和虚部高外延频段)时的虚部数据示意图,图6为特征频率在外延频段(虚部低外延频段和虚部高外延频段)时的虚部数据示意图。Figure 5 is a schematic diagram of the imaginary part data when the characteristic frequency is not in the epitaxial frequency band (the imaginary part low epitaxial frequency band and the imaginary part high epitaxial frequency band), and Figure 6 is the characteristic frequency when the characteristic frequency is in the epitaxial frequency band (imaginary part low epitaxial frequency band and the imaginary part high epitaxial frequency band) Schematic diagram of the imaginary part of the data.
如图5或图6所示,可以根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数。As shown in Figure 5 or Figure 6, the imaginary part dielectric response fitting function in the imaginary low epitaxial frequency band can be obtained according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency. Part data and the imaginary part data corresponding to the second highest frequency obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band.
如图5-图6所示,虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数均为线性函数。虚部低外延频段[ω' Ext-LL]位于最低频率(ω L)与虚部最低外延频率(ω' Ext-L)之间,虚部高外延频段[ω H,ω' Ext-H]位于最高频率(ω H)与虚部最高外延频率(ω' Ext-H)之间。在[ω' Ext-L,ω' Ext-H]以外的频段,采用截断处理的方法进行赋值。即在低于ω Ext-L的频段和在高于ω Ext-H的频段,虚部数据均赋值为0。 As shown in Figure 5-6, the imaginary dielectric response fitting function in the imaginary low epitaxial frequency band and the imaginary dielectric response fitting function in the imaginary high epitaxial frequency band are linear functions. The imaginary low-epitaxial frequency band [ω' Ext-LL ] is located between the lowest frequency (ω L ) and the lowest imaginary-part epitaxial frequency (ω' Ext-L ), and the imaginary high-epitaxial frequency band [ω H ,ω' Ext -H ] is located between the highest frequency (ω H ) and the highest epitaxial frequency of the imaginary part (ω' Ext-H ). In the frequency bands other than [ω' Ext-L ,ω' Ext-H ], the method of truncation is used for assignment. That is, in the frequency band lower than ω Ext-L and in the frequency band higher than ω Ext-H , the imaginary part data is assigned a value of 0.
4、根据奇异点ω S和实部插值步长计算实部奇异点上界ω S-H和实部奇异点下界ω S-L;其中,奇异点ω S位于最高频率ω H与最低频率ω L之间。 4. Calculate the upper bound of the real singular point ω SH and the lower bound of the real singular point ω SL according to the singular point ω S and the real part interpolation step length; among them, the singular point ω S is located between the highest frequency ω H and the lowest frequency ω L.
一实施例中,实部奇异点上界与奇异点的关系如下:In an embodiment, the relationship between the upper bound of the real part of the singular point and the singular point is as follows:
ln(ω S-H)=ln(ω S)+ln(Δω); ln(ω SH )=ln(ω S )+ln(Δω);
实部奇异点下界与奇异点的关系如下:The relationship between the lower bound of the real part of the singular point and the singular point is as follows:
ln(ω S-L)=ln(ω S)-ln(Δω); ln(ω SL )=ln(ω S )-ln(Δω);
其中,ln(Δω)为实部插值步长。Among them, ln(Δω) is the real part interpolation step.
5、根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点和奇异点对应的实部数据确定第一奇异点实部系数和第二奇异点实部系数;根据实部奇异点下界、实部奇异点下界对应的实部数据、奇异点和奇异点对应的实部数据确定第三奇异点实部系数和第四奇异点实部系数。5. Determine the first singular point real part coefficient and the second singular point real part coefficient according to the upper bound of the real part singular point, the real part data corresponding to the upper bound of the real part singular point, the singular point and the real part data corresponding to the singular point; The real part singular point lower bound, the real part data corresponding to the real part singular point lower bound, the singular point and the real part data corresponding to the singular point determine the third singular point real coefficient and the fourth singular point real coefficient.
6、根据第一奇异点实部系数、第二奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点上频段积分。其中,实部奇异点上频段积分为实部奇异点上频段[ω SS-H]的积分。 6. Calculate the upper band integral of the real singular point according to the real part coefficient of the first singular point, the real part coefficient of the second singular point, the singular point, the upper bound of the real singular point and the lower bound of the real singular point. Among them, the frequency band integral on the singular point of the real part is the integral of the frequency band [ω SSH ] on the singular point of the real part.
7、根据第三奇异点实部系数、第四奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点下频段积分。其中,实部奇异点下频段积分为实部奇异点下频段[ω S-LS]的积分。 7. Calculate the lower band integral of the real singular point according to the real part coefficient of the third singular point, the real part coefficient of the fourth singular point, the singular point, the upper bound of the real singular point and the lower bound of the real singular point. Among them, the real part singular point lower frequency band integral is the real part singular point lower frequency band [ω SLS ] integral.
8、将实部奇异点上频段积分与实部奇异点下频段积分相加,得到实部奇异点频段积分。8. Add the upper frequency band integral of the real part singular point and the lower frequency band integral of the real part singular point to obtain the real part singular point frequency band integral.
9、根据实部最高外延频率对应的实部数据、实部最高外延频率和奇异点计算实部上截断频域积分;根据实部最低外延频率对应的实部数据、实部最低外延频率和奇异点计算实部下截断频域积分。将实部上截断频域积分与实部下截断频域积分相加,得到实部截断频域积分。9. Calculate the truncated frequency domain integral of the real part according to the real data corresponding to the highest extension frequency of the real part, the highest extension frequency of the real part and the singular point; according to the real part data corresponding to the lowest extension frequency of the real part, the lowest extension frequency and singularity of the real part Point to calculate the truncated frequency domain integral under the real part. The upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
其中,实部上截断频域积分为实部上截断频域(0,ω Ext-L)的积分;实部下截断频域积分为实部下截断频域(ω Ext-H,∞)的积分。 Among them, the upper truncated frequency domain integral of the real part is the integral of the upper truncated frequency domain (0,ω Ext-L ); the lower truncated frequency domain integral of the real part is the integral of the lower truncated frequency domain of the real part (ω Ext-H ,∞).
10、根据有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数、实部高外延频段内的实部介电响应拟合函数、实部最低外延频率、实部最高外延频率、实部奇异点上界和实部奇异点下界,计算实部外延频段积分。10. According to the real part dielectric response fitting function in the limited frequency band, the real part dielectric response fitting function in the real part low epitaxial frequency band, the real part dielectric response fitting function in the real part high epitaxial frequency band, the real part The lowest extension frequency, the highest extension frequency of the real part, the upper bound of the singular point of the real part and the lower bound of the singular point of the real part, calculate the real part of the extension frequency band integral.
具体实施时,因为已经确定有限频段[ω LH]内的实部介电响应拟合函数、实部低外延频段[ω Ext-LL]内的实部介电响应拟合函数和实部高外延频段[ω HExt-H]内的实部介电响应拟合函数,因此可以根据辛普森公式计算实部外延频段[ω Ext-LS-L)的积分。 In the specific implementation, because the real part dielectric response fitting function in the limited frequency band [ω LH ] and the real part dielectric response fitting in the real low epitaxial frequency band [ω Ext-LL ] have been determined Function and real part of the high-extension frequency band [ω H , ω Ext-H ] in the real part of the dielectric response fitting function, so the integral of the real part of the epitaxial frequency band [ω Ext-L , ω SL) can be calculated according to the Simpson formula.
11、根据实部奇异点频段积分、实部截断频域积分和实部外延频段积分计算模拟虚部积分。11. Calculate the analog imaginary part integral based on the real part singular point frequency band integral, the real part truncated frequency domain integral and the real part epitaxial frequency band integral.
将实部奇异点频段积分、实部截断频域积分和实部外延频段积分相加,可以得到模拟虚部积分。Add the real part singular point frequency band integral, the real part truncated frequency domain integral and the real part epitaxial frequency band integral to get the analog imaginary part integral.
12、判断模拟虚部积分是否小于预设分辨率。当大于或等于预设分辨率时,更新实部插值步长,返回重新执行步骤4;当小于预设分辨率时,将模拟虚部积分与实际虚部积分中的相同部分作为极化过程实部信息,将模拟虚部积分与实际虚部积分中的不同部分作为电导信息。12. Judge whether the integral of the simulated imaginary part is less than the preset resolution. When the resolution is greater than or equal to the preset resolution, update the real part interpolation step and return to step 4; when the resolution is less than the preset resolution, the same part of the simulated imaginary part integral and the actual imaginary part integral will be regarded as the real part of the polarization process. Part information, the difference between the simulated imaginary part integral and the actual imaginary part integral is used as the conductance information.
13、根据奇异点ω S和虚部插值步长计算虚部奇异点上界ω' S-H和虚部奇异点下界ω' S-L13 The singular points ω S and imaginary interpolation step calculating the imaginary part of the boundary singular point ω 'SH and imaginary singular points lower bound ω' SL.
一实施例中,虚部奇异点上界与奇异点的关系如下:In an embodiment, the relationship between the upper bound of the imaginary part of the singular point and the singular point is as follows:
ln(ω' S-H)=ln(ω S)+ln(Δω'); ln(ω' SH )=ln(ω S )+ln(Δω');
虚部奇异点下界与奇异点的关系如下:The relationship between the lower bound of the imaginary part of the singular point and the singular point is as follows:
ln(ω' S-L)=ln(ω S)-ln(Δω'); ln(ω' SL )=ln(ω S )-ln(Δω');
其中,ln(Δω')为虚部插值步长。Among them, ln(Δω') is the imaginary part interpolation step.
14、根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点和奇异点对应的虚部数据确定第一奇异点虚部系数和第二奇异点虚部系数;根据虚部奇异点下界、虚部奇异点下界对应的虚部数据、奇异点和奇异点对应的虚部数据确定第三奇异点虚部系数和第四奇异点虚部系数。14. Determine the first singular point imaginary part coefficient and the second singular point imaginary part coefficient according to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point upper bound, the singular point and the imaginary part data corresponding to the singular point; The lower bound of the imaginary part singular point, the imaginary part data corresponding to the lower bound of the imaginary part singular point, the singular point and the imaginary part data corresponding to the singular point determine the third singular point imaginary part coefficient and the fourth singular point imaginary part coefficient.
15、根据第一奇异点虚部系数、第二奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点上频段积分。其中,虚部奇异点上频段积分为虚部奇异点上频段[ω S,ω' S-H]的积分。 15. Calculate the upper band integral of the imaginary singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point. Among them, the frequency band integral on the imaginary singular point is the integral of the frequency band [ω S ,ω' SH ] on the imaginary singular point.
16、根据第三奇异点虚部系数、第四奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点下频段积分。其中,虚部奇异点下频段积分为虚部奇异点下频段[ω' S-LS]的积分。 16. Calculate the lower band integral of the imaginary singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point. Among them, the frequency band integration under the imaginary part singular point is the integration of the frequency band under the imaginary part singular point [ω' SLS ].
17、将虚部奇异点上频段积分与虚部奇异点下频段积分相加,得到虚部奇异点频段积分。17. Add the upper frequency band integral of the imaginary part singular point and the lower frequency band integral of the imaginary part singular point to obtain the frequency band integral of the imaginary part singular point.
18、根据有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数、虚部高外延频段内的虚部介电响应拟合函数、虚部最低外延频率、虚部最高外延频率、虚部奇异点上界和虚部奇异点下界,计算虚部外延频段积分。18.According to the imaginary part dielectric response fitting function in the limited frequency band, the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band, the imaginary part The lowest extension frequency, the highest extension frequency of the imaginary part, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point, calculate the imaginary part extension frequency band integral.
19、根据虚部奇异点频段积分和虚部外延频段积分计算模拟实部积分。19. Calculate the analog real part integral based on the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral.
将据虚部奇异点频段积分和虚部外延频段积分相加,可以得到模拟实部积分。The imaginary part singular point frequency band integral and the imaginary part extension frequency band integral are added together to obtain the analog real part integral.
20、判断模拟实部积分是否小于预设分辨率。当大于或等于预设分辨率时,更新虚部插值步长,返回重新执行步骤13。当小于预设分辨率时,将模拟实部积分与实际实部积分中的相同部分作为极化过程虚部信息,将模拟实部积分与实际实部积分中的不同部分作为电容信息。20. Judge whether the integral of the simulated real part is less than the preset resolution. When the resolution is greater than or equal to the preset resolution, the imaginary part interpolation step is updated, and step 13 is executed again. When the resolution is less than the preset resolution, the same part of the simulated real part integral and the actual real part integral is used as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral is used as the capacitance information.
21、根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。21. Obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process.
图7是本发明实施例中极化率的对比示意图。如图7所示,横坐标为频率,单位为Hz;纵坐标为极化率。χ'(ω)为极化率的实际虚部积分,χ”(ω)为极化率的模拟虚部积分。为更清晰地展示计算获得的极化率的模拟虚部积分,将χ”(ω)竖直向下平移了两个数量级。由图7可见,极化率的实际虚部积分与极化率的模拟虚部积分一致。Fig. 7 is a schematic diagram of comparison of polarizability in an embodiment of the present invention. As shown in Figure 7, the abscissa is the frequency and the unit is Hz; the ordinate is the polarization rate. χ'(ω) is the integral of the actual imaginary part of the polarizability, and χ”(ω) is the integral of the simulated imaginary part of the polarizability. In order to show the calculated integral of the simulated imaginary part of the polarizability more clearly, replace χ” (ω) is shifted vertically downward by two orders of magnitude. It can be seen from Fig. 7 that the actual imaginary part integral of the polarizability is consistent with the simulated imaginary part integral of the polarizability.
图8是本发明实施例中极化率的对比示意图。如图8所示,横坐标为频率,单位为Hz;纵坐标为极化率。χ'(ω)为极化率的实际实部积分,χ”(ω)为极化率的模拟实部积分。为更清晰地展示计算获得的极化率的模拟实部积分,将χ”(ω)竖直向下平移了两个数量级。由图8可见,极化率的实际实部积分与极化率的模拟实部积分一致。Fig. 8 is a schematic diagram of comparison of polarizability in an embodiment of the present invention. As shown in Figure 8, the abscissa is the frequency and the unit is Hz; the ordinate is the polarization rate. χ'(ω) is the actual integral of the real part of the polarizability, and χ”(ω) is the integral of the simulated real part of the susceptibility. In order to show the calculated integral of the real part of the polarizability more clearly, add χ” (ω) is shifted vertically downward by two orders of magnitude. It can be seen from Fig. 8 that the actual real part integral of the polarizability is consistent with the simulated real part integral of the polarizability.
图9是本发明实施例中复电容的模拟实部积分和模拟虚部积分的示意图。图10是本发明实施例中复电容的极化过程实部信息和极化过程虚部信息的示意图。图11是本发明实施例中复电容的电导信息和无穷频率电容信息的示意图。图9-图11的横坐标均为频率,单位为Hz;纵坐标均为电容,单位为pF。图9中的C'(ω)为复电容的模拟实部积分,C”(ω)为复电容的模拟虚部积分;图10中的χ1'(ω)为复电容的极化过程实部信息,χ1”(ω)为复电容的极化过程虚部信息;图11中的σ为复电容的电导信息,C∞为复电容的无穷频率电容信息。Fig. 9 is a schematic diagram of the simulated real part integral and the simulated imaginary part integral of a complex capacitor in an embodiment of the present invention. FIG. 10 is a schematic diagram of the real part information and the imaginary part information of the polarization process of a complex capacitor in an embodiment of the present invention. FIG. 11 is a schematic diagram of conductance information and infinite frequency capacitance information of a complex capacitor in an embodiment of the present invention. The abscissas of Fig. 9-11 are all frequency, and the unit is Hz; the ordinates are all capacitance, and the unit is pF. C'(ω) in Figure 9 is the integration of the simulated real part of the complex capacitor, C”(ω) is the integration of the simulated imaginary part of the complex capacitor; χ1'(ω) in Figure 10 is the real part of the polarization process of the complex capacitor Information, χ1”(ω) is the imaginary part information of the polarization process of the complex capacitor; σ in Fig. 11 is the conductance information of the complex capacitor, and C∞ is the infinite frequency capacitance information of the complex capacitor.
本发明可以在各种电力设备的介电响应分析中进行应用,下面以电力系统外绝缘领域最常见的高温硫化硅橡胶为例,详细说明本发明的应用方式。The present invention can be applied in the dielectric response analysis of various power equipment. The following takes the most common high-temperature vulcanized silicone rubber in the field of external insulation of power systems as an example to illustrate the application of the present invention in detail.
高温硫化硅橡胶的主要成分是聚二甲基硅氧烷、纳米白炭黑填料和纳米氢氧化铝填料,试验测试温度为80℃,测量频段为10 -4Hz~10 3Hz。 The main components of high-temperature vulcanized silicone rubber are polydimethylsiloxane, nano-silica filler and nano-aluminum hydroxide filler. The test temperature is 80°C, and the measurement frequency range is 10 -4 Hz to 10 3 Hz.
图12是本发明实施例中高温硫化硅橡胶复电容的介电响应参数示意图。图13是本发明实施例中高温硫化硅橡胶复电容的极化过程信息示意图。图14是本发明实施例中高温硫化硅橡胶复电容的电导信息和无穷频率电容信息的示意图。图15是本发明实施例中高温硫化硅橡胶复电容的极化率的实部数据对比示意图。图16是本发明实施例中高温硫化硅橡胶复电容的极化率的虚部数据对比示意图。图12-图16的横坐标均为频率,单位为Hz;纵坐标均为电容,单位为pF。Fig. 12 is a schematic diagram of dielectric response parameters of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention. Fig. 13 is a schematic diagram of polarization process information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention. Fig. 14 is a schematic diagram of conductivity information and infinite frequency capacitance information of a high-temperature vulcanized silicone rubber compound capacitor in an embodiment of the present invention. FIG. 15 is a schematic diagram showing the comparison of the real part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention. 16 is a schematic diagram of comparison of imaginary part data of the polarizability of the high-temperature vulcanized silicone rubber compound capacitance in the embodiment of the present invention. The abscissas of Figure 12-16 are all frequency, the unit is Hz; the ordinates are all capacitance, the unit is pF.
如图12所示,图12中的C'(ω)为复电容的介电响应参数的实际实部积分,C”(ω)为复电容的介电响应参数的实际虚部积分。直观地,实际实部积分在所测频段中基本保持不变,实际虚部积分在10 0Hz~10 3Hz频段中存在一明显的弛豫峰,在10 -4Hz~10 0Hz频段存在明显的电导过程(即虚部与频率呈-1次幂的关系)。 As shown in Figure 12, C'(ω) in Figure 12 is the actual real integral of the dielectric response parameter of the complex capacitor, and C"(ω) is the actual imaginary integral of the dielectric response parameter of the complex capacitor. Intuitively , the actual implementation of the integral portion remains substantially unchanged in the measured frequency band, the actual presence of a significant imaginary part of the integral in the relaxation peak 10 0 Hz ~ 10 3 Hz band, the 10 -4 Hz ~ 10 0 Hz band obvious Conductance process (that is, the relationship between the imaginary part and frequency is -1 power).
如图13所示,图13中的χ1'(ω)为复电容的极化过程实部信息,χ1”(ω)为复电容的极化过程虚部信息。极化过程实部信息在10 -4Hz~10 -2Hz频段和10 1Hz~10 3Hz频段存在明显的弥散现象,相应地,极化过程虚部信息在这两频段上也存在弥散现象。更具体而言,极化过程的实部信息的数值和极化过程的虚部信息的数值在10 -4Hz~10 -2Hz频段上都随着频率的降低而增加,这是典型的低频弥散过程的特征。在10 1Hz~10 3Hz频段,极化过程的虚部信息出现弛豫峰,且在10 2Hz~10 3Hz频段,可以观察到极化过程的实部信息的数值和极化过程的虚部信息的数值在双指数坐标中的平行现象,即观察到了普适的极化率的指数规律。 As shown in Figure 13, χ1'(ω) in Figure 13 is the real part information of the polarization process of the complex capacitor, and χ1"(ω) is the imaginary part information of the polarization process of the complex capacitor. The real part information of the polarization process is at 10 -4 Hz~10 -2 Hz frequency band and 10 1 Hz~10 3 Hz frequency band have obvious dispersion phenomenon. Correspondingly, the imaginary part information of the polarization process also has dispersion phenomenon in these two frequency bands. More specifically, polarization The value of the real part information of the process and the value of the imaginary part information of the polarization process increase as the frequency decreases in the 10 -4 Hz~10 -2 Hz frequency band, which is a characteristic of the typical low-frequency dispersion process. 1 Hz ~ 10 3 Hz frequency range, the imaginary part of the information occurs during polarization relaxation peak and the 10 2 Hz ~ 10 3 Hz band can be observed polarization value and the imaginary part of the process of the real information part of the polarization process The parallel phenomenon of the numerical value of the information in the double exponential coordinate, that is, the universal exponential law of polarizability is observed.
如图14所示,图14中的σ为复电容的电导信息,C∞为复电容的无穷频率电容信息。As shown in Figure 14, σ in Figure 14 is the conductance information of the complex capacitor, and C∞ is the infinite frequency capacitance information of the complex capacitor.
如图15所示,图15中的C'(ω)为复电容的介电响应参数的实际实部积分,χ1'(ω)为复电容的极化过程实部信息。在复 电容的介电响应参数的实际实部积分中,由于无穷频率电容信息的数值较大(如图14所示),因此极化过程的实部信息的弥散过程几乎完全被无穷频率电容所遮盖,这就导致在复电容的介电响应参数的实际实部积分中,几乎观察不到弥散现象。若仅分析复电容的介电响应参数的实际实部积分,会误以为高温硫化硅橡胶在所测频段不发生介电弥散过程,从图15中可以看出,这个判断是与实际物理实质明显不符的。As shown in Figure 15, C'(ω) in Figure 15 is the actual real part integral of the dielectric response parameter of the complex capacitor, and χ1'(ω) is the real part information of the polarization process of the complex capacitor. In the actual real part integration of the dielectric response parameters of the complex capacitor, since the value of the infinite frequency capacitance information is relatively large (as shown in Figure 14), the dispersion process of the real part information in the polarization process is almost completely affected by the infinite frequency capacitance. Covering, which leads to almost no dispersion phenomenon observed in the actual real part integral of the dielectric response parameter of the complex capacitor. If only the actual real part integral of the dielectric response parameters of the complex capacitor is analyzed, it will be misunderstood that the high-temperature vulcanized silicone rubber does not undergo a dielectric dispersion process in the measured frequency band. It can be seen from Figure 15 that this judgment is obviously in essence with the actual physical Inconsistent.
如图16所示,图16中的C”(ω)为复电容的介电响应参数的实际虚部积分,χ1”(ω)为极化率的介电响应参数的模拟虚部积分。可以发现,在高频段(10 1Hz~10 3Hz频段),实际虚部积分与模拟虚部积分(考虑形状系数)是一致的。而在10 0Hz以下频段,实际虚部积分明显大于模拟虚部积分(考虑形状系数),二者的差值即为电导信息。图16清晰地显示,电导信息主要会影响复电容虚部低频下的测量结果。若仅分析复电容的介电响应参数的实际虚部积分,会误以为复电容的虚部在低频段仅由电导过程所决定,从图16中可以看出,这个判断也是与实际物理实质明显不符的。高温硫化硅橡胶复电容的实际虚部积分不仅存在电导过程的贡献,还存在低频弥散过程的贡献,且电导过程掩盖了复电容的实际虚部积分在低频段的极化过程的贡献,只有通过极化率虚部的分析,才可以观察到硅橡胶介电响应中的低频弥散现象。 As shown in Fig. 16, C"(ω) in Fig. 16 is the actual imaginary part integral of the dielectric response parameter of the complex capacitor, and χ1"(ω) is the simulated imaginary part integral of the dielectric response parameter of the polarizability. It can be found that in the high frequency band (10 1 Hz~10 3 Hz frequency band), the actual imaginary part integral is consistent with the simulated imaginary part integral (considering the shape coefficient). 10 0 Hz in frequency bands below, the actual score was significantly larger than the imaginary part analog integrated imaginary part (considering the shape factor), is the difference between the two conductivity information. Figure 16 clearly shows that the conductance information mainly affects the measurement results at the low frequency of the imaginary part of the complex capacitance. If only the actual imaginary part integral of the dielectric response parameter of the complex capacitor is analyzed, it will be mistaken to think that the imaginary part of the complex capacitor is determined only by the conductance process in the low frequency range. It can be seen from Figure 16 that this judgment is also obvious from the actual physical essence. Inconsistent. The actual imaginary part integral of the high-temperature vulcanized silicone rubber compound capacitor has not only the contribution of the conductance process, but also the contribution of the low-frequency dispersion process, and the conductance process conceals the contribution of the actual imaginary part integral of the complex capacitance to the polarization process at the low frequency. Only by analyzing the imaginary part of the polarizability can the low-frequency dispersion phenomenon in the dielectric response of silicone rubber be observed.
综上,本发明实施例的电介质状态分析方法根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。In summary, the dielectric state analysis method of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral, Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain Capacitance information and imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, avoiding misjudgments, missed judgments, and Misjudgment saves the cost of operation and maintenance of power equipment.
基于同一发明构思,本发明实施例还提供了一种电介质状态分析系统,由于该系统解决问题的原理与电介质状态分析方法相似,因此该系统的实施可以参见方法的实施,重复之处不再赘述。Based on the same inventive concept, the embodiment of the present invention also provides a dielectric state analysis system. Since the principle of the system to solve the problem is similar to the dielectric state analysis method, the implementation of the system can refer to the implementation of the method, and the repetition will not be repeated. .
图17是本发明实施例中电介质状态分析系统的结构框图。如图17所示,电介质状态分析系统包括:Fig. 17 is a structural block diagram of a dielectric state analysis system in an embodiment of the present invention. As shown in Figure 17, the dielectric state analysis system includes:
获取单元,用于获取多个频率对应的介电响应参数的实部数据和虚部数据;An acquiring unit for acquiring real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies;
确定单元,用于根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;The determining unit is used to determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
有限频段拟合函数单元,用于拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;其中,有限频段位于最低频率与最高频率之间;The finite frequency band fitting function unit is used to fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in the finite frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part data in the finite frequency band. The imaginary part dielectric response fitting function; among them, the limited frequency band is located between the lowest frequency and the highest frequency;
外延频段拟合函数单元,用于根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;The epitaxial frequency band fitting function unit is used to obtain the real part dielectric response fitting function in the real low epitaxial frequency band according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, and according to the real part corresponding to the highest frequency Data and the real part data corresponding to the second highest frequency obtain the real part dielectric response fitting function in the real part high epitaxial frequency band; obtain the imaginary part low epitaxial frequency band according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency The imaginary part dielectric response fitting function in the inner part, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency, the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band is obtained; the real part low epitaxial frequency band Located between the lowest frequency and the lowest epitaxial frequency of the real part, the high epitaxial frequency band of the real part is between the highest frequency and the highest epitaxial frequency of the real part, the low epitaxial frequency band of the imaginary part is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary part is high epitaxial frequency The frequency band is located between the highest frequency and the highest epitaxial frequency of the imaginary part;
实部迭代单元,用于执行如下迭代处理:The real part iterative unit is used to perform the following iterative processing:
根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; among them, the singular point is located between the highest frequency and the lowest frequency;
根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part in the low epitaxial frequency band of the real part. Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral;
判断模拟虚部积分是否小于预设分辨率;当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain the conductance information and the real part information of the polarization process, otherwise update the real part interpolation step long;
虚部迭代单元,用于执行如下迭代处理:The imaginary part iterative unit is used to perform the following iterative processing:
根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part in the low epitaxial frequency band of the imaginary part. Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral;
判断模拟实部积分是否小于预设分辨率;当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain the capacitance information and the imaginary part information of the polarization process, otherwise update the imaginary part interpolation step long;
电介质状态单元,用于根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。The dielectric state unit is used to obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
在其中一种实施例中,实部迭代单元具体用于:In one of the embodiments, the real part iteration unit is specifically used for:
根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点、奇异点对应的实部数据、实部奇异点下界和实部奇异点下界对应的实部数据计算实部奇异点频段积分;Calculate the real part based on the real data corresponding to the upper bound of the real singular point and the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point, and the real data corresponding to the lower bound of the real singular point Singularity frequency band integration;
根据实部最低外延频率对应的实部数据、实部最低外延频率、实部最高外延频率对应的实部数据、实部最高外延频率和奇异点计算实部截断频域积分;Calculate the real part truncated frequency domain integral according to the real part data corresponding to the lowest real part epitaxial frequency, the lowest real part epitaxial frequency, the real part data corresponding to the highest real part epitaxial frequency, the highest real part epitaxial frequency and the singular point;
根据有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数、实部高外延频段内的实部介电响应拟合函数、实部最低外延频率、实部最高外延频率、实部奇异点上界和实部奇异点下界,计算实部外延频段积分;According to the real part dielectric response fitting function in the limited frequency band, the real part dielectric response fitting function in the real part low epitaxial frequency band, the real part dielectric response fitting function in the real part high epitaxial frequency band, and the real part lowest epitaxy Frequency, real part maximum extension frequency, real part singular point upper bound and real part singular point lower bound, calculate the real part extension frequency band integral;
根据实部奇异点频段积分、实部截断频域积分和实部外延频段积分计算模拟虚部积分;Calculate the analog imaginary part integral based on the real part singular point frequency band integral, the real part truncated frequency domain integral and the real part extension frequency band integral;
虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点、奇异点对应的虚部数据、虚部奇异点下界和虚部奇异点下界对应的虚部数据计算虚部奇异点频段积分;Calculate the imaginary part based on the imaginary data corresponding to the upper bound of the imaginary part singular point, the imaginary part corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the imaginary part data corresponding to the lower bound of the imaginary part singular point Singularity frequency band integration;
根据有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数、虚部高外延频段内的虚部介电响应拟合函数、虚部最低外延频率、虚部最高外延频率、虚部奇异点上界和虚部奇异点下界,计算虚部外延频段积分;According to the imaginary part dielectric response fitting function in the limited frequency band, the imaginary part dielectric response fitting function in the imaginary part low extension frequency band, the imaginary part dielectric response fitting function in the imaginary part high extension frequency band, and the imaginary part minimum extension Frequency, imaginary part maximum extension frequency, imaginary part singular point upper bound and imaginary part singular point lower bound, calculate the imaginary part extension frequency band integral;
根据虚部奇异点频段积分和虚部外延频段积分计算模拟实部积分。Calculate the analog real part integral based on the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral.
在其中一种实施例中,实部迭代单元具体用于:In one of the embodiments, the real part iteration unit is specifically used for:
根据实部奇异点上界、实部奇异点上界对应的实部数据、奇异点、奇异点对应的实部数据、实部奇异点下界和实部奇异点下界对应的实部数据确定第一奇异点实部系数、第二奇异点实部系数、第三奇异点实部系数和第四奇异点实部系数;According to the real data corresponding to the upper bound of the real singular point, the upper bound of the real singular point, the real data corresponding to the singular point, the real data corresponding to the singular point, the lower bound of the real singular point, and the real data corresponding to the lower bound of the real singular point, determine the first Singular point real part coefficients, second singular point real part coefficients, third singular point real part coefficients and fourth singular point real part coefficients;
根据第一奇异点实部系数、第二奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点上频段积分;Calculate the upper band integral of the real singular point according to the real coefficient of the first singular point, the real coefficient of the second singular point, the singular point, the upper bound of the real singular point and the lower bound of the real singular point;
根据第三奇异点实部系数、第四奇异点实部系数、奇异点、实部奇异点上界和实部奇异点下界计算实部奇异点下频段积分;Calculate the lower band integral of the real singular point according to the third singular point real part coefficient, the fourth singular point real part coefficient, the singular point, the upper bound of the real part singular point and the lower bound of the real part singular point;
将实部奇异点上频段积分与实部奇异点下频段积分相加,得到实部奇异点频段积分;Add the upper frequency band integral of the real singular point and the lower frequency band integral of the real singular point to obtain the real singular point frequency band integral;
虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
根据虚部奇异点上界、虚部奇异点上界对应的虚部数据、奇异点、奇异点对应的虚部数据、虚部奇异点下界和虚部奇异点下界对应的虚部数据确定第一奇异点虚部系数、第二奇异点虚部系数、第三奇异点虚部系数和第四奇异点虚部系数;Determine the first according to the imaginary part data corresponding to the upper bound of the imaginary part singular point, the imaginary part corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point, and the imaginary part data corresponding to the lower bound of the imaginary part singular point. Singular point imaginary part coefficient, second singular point imaginary part coefficient, third singular point imaginary part coefficient and fourth singular point imaginary part coefficient;
根据第一奇异点虚部系数、第二奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点上频段积分;Calculate the upper band integral of the imaginary singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary singular point and the lower bound of the imaginary singular point;
根据第三奇异点虚部系数、第四奇异点虚部系数、奇异点、虚部奇异点上界和虚部奇异点下界计算虚部奇异点下频段积分;Calculate the lower band integral of the imaginary singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point;
将虚部奇异点上频段积分与虚部奇异点下频段积分相加,得到虚部奇异点频段积分。Add the upper frequency band integral of the imaginary part singular point and the lower frequency band integral of the imaginary part singular point to obtain the frequency band integral of the imaginary part singular point.
在其中一种实施例中,通过如下公式计算实部奇异点上频段积分:In one of the embodiments, the frequency band integral on the singular point of the real part is calculated by the following formula:
Figure PCTCN2020109652-appb-000009
Figure PCTCN2020109652-appb-000009
其中,E 1为实部奇异点上频段积分,a为第一奇异点实部系数,b为第二奇异点实部系数,ω S为奇异点,ω S-H为实部奇异点上界,ω S-L为实部奇异点下界; Among them, E 1 is the frequency band integral on the real singular point, a is the real coefficient of the first singular point, b is the real coefficient of the second singular point, ω S is the singular point, ω SH is the upper bound of the real singular point, ω SL is the lower bound of the singular point of the real part;
通过如下公式计算实部奇异点下频段积分:Calculate the frequency band integral under the singular point of the real part by the following formula:
Figure PCTCN2020109652-appb-000010
Figure PCTCN2020109652-appb-000010
其中,E 2为实部奇异点下频段积分,c为第三奇异点实部系数,d为第四奇异点实部系数; Among them, E 2 is the frequency band integral of the real part of the singular point, c is the real part coefficient of the third singular point, and d is the real part coefficient of the fourth singular point;
通过如下公式计算虚部奇异点上频段积分:The frequency band integral on the singular point of the imaginary part is calculated by the following formula:
Figure PCTCN2020109652-appb-000011
Figure PCTCN2020109652-appb-000011
其中,E' 1为虚部奇异点上频段积分,a'为第一奇异点虚部系数,b'为第二奇异点虚部系数,ω' S-H为虚部奇异点上界,ω' S-L为虚部奇异点下界; Wherein, E '1 is the imaginary part of the integral bands singular point, a' is a first singular point of the imaginary part of the coefficient, b 'is the imaginary part of the second singular point coefficients, ω' SH bounded singularity point on the imaginary part, ω 'SL Is the lower bound of the imaginary part singular point;
通过如下公式计算虚部奇异点下频段积分:Calculate the lower band integral of the imaginary part singular point by the following formula:
Figure PCTCN2020109652-appb-000012
Figure PCTCN2020109652-appb-000012
其中,E' 2为虚部奇异点下频段积分,c'为第三奇异点虚部系数,d'为第四奇异点虚部系数。 Wherein, E '2 is a lower frequency band integrating the imaginary part of a singular point, c' is the coefficient of the imaginary part of the third singular point, d 'is the imaginary part of a singular point of the fourth coefficient.
在其中一种实施例中,实部迭代单元具体用于:In one of the embodiments, the real part iteration unit is specifically used for:
根据实部最高外延频率对应的实部数据、实部最高外延频率和奇异点计算实部上截断频域积分;Calculate the truncated frequency domain integral on the real part according to the real data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part and the singular point;
根据实部最低外延频率对应的实部数据、实部最低外延频率和奇异点计算实部下截断频域积分;Calculate the lower truncated frequency domain integral of the real part according to the real data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part and the singular point;
将实部上截断频域积分与实部下截断频域积分相加,得到实部截断频域积分。The upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
在其中一种实施例中,通过如下公式计算实部上截断频域积分:In one of the embodiments, the truncated frequency domain integral on the real part is calculated by the following formula:
Figure PCTCN2020109652-appb-000013
Figure PCTCN2020109652-appb-000013
其中,F 1为实部上截断频域积分,χ'(ω Ext-H)为实部最高外延频率对应的实部数据,ω S为奇异点,ω Ext-H为实部最高外延频率; Among them, F 1 is the truncated frequency domain integral on the real part, χ'(ω Ext-H ) is the real part data corresponding to the highest extension frequency of the real part, ω S is the singularity point, and ω Ext-H is the highest extension frequency of the real part;
通过如下公式计算实部下截断频域积分:Calculate the truncated frequency domain integral under the real part by the following formula:
Figure PCTCN2020109652-appb-000014
Figure PCTCN2020109652-appb-000014
其中,F 2为实部下截断频域积分,χ'(ω Ext-L)为实部最低外延频率对应的实部数据,ω Ext-L为实部最低外延频率。 Among them, F 2 is the lower truncated frequency domain integral of the real part, χ'(ω Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part, and ω Ext-L is the lowest epitaxial frequency of the real part.
在其中一种实施例中,通过如下公式更新实部插值步长:In one of the embodiments, the real interpolation step size is updated by the following formula:
Figure PCTCN2020109652-appb-000015
Figure PCTCN2020109652-appb-000015
其中,ln(Δω n)为第n次迭代中的实部插值步长,ln(Δω n+1)为第n+1次迭代中的实部插值步长; Among them, ln(Δω n ) is the real part interpolation step in the nth iteration, and ln(Δω n+1 ) is the real part interpolation step in the n+1th iteration;
通过如下公式更新虚部插值步长:Update the imaginary part interpolation step size by the following formula:
Figure PCTCN2020109652-appb-000016
Figure PCTCN2020109652-appb-000016
其中,ln(Δω' n)为第n次迭代中的虚部插值步长,ln(Δω' n+1)为第n+1次迭代中的虚部插值步长。 Among them, ln(Δω' n ) is the imaginary part interpolation step in the nth iteration, and ln(Δω' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
在其中一种实施例中,实部迭代单元具体用于:In one of the embodiments, the real part iteration unit is specifically used for:
将模拟虚部积分与实际虚部积分中的相同部分作为极化过程实部信息,将模拟虚部积分与实际虚部积分中的不同部分作为电导信息;Use the same part of the simulated imaginary part integral and the actual imaginary part integral as the real part information of the polarization process, and use the different part of the simulated imaginary part integral and the actual imaginary part integral as the conductance information;
虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
将模拟实部积分与实际实部积分中的相同部分作为极化过程虚部信息,将模拟实部积分与实际实部积分中的不同部分作为电容信息。The same part of the simulated real part integral and the actual real part integral is regarded as the imaginary part information of the polarization process, and the different part of the simulated real part integral and the actual real part integral is regarded as the capacitance information.
综上,本发明实施例的电介质状态分析系统根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。In summary, the dielectric state analysis system of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral, Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain Capacitance information and imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, avoiding misjudgments, missed judgments, and Misjudgment saves the cost of operation and maintenance of power equipment.
本发明实施例还提供了一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,处理器执行计算机程序时可以实现电介质状态分析方法的全部或部分内容,例如,处理器执行计算机程序时可以实现如下内容:The embodiment of the present invention also provides a computer device, including a memory, a processor, and a computer program stored in the memory and running on the processor. The processor can implement all or part of the dielectric state analysis method when the computer program is executed. For example, when a processor executes a computer program, the following can be realized:
获取多个频率对应的介电响应参数的实部数据和虚部数据;Acquire real data and imaginary data of dielectric response parameters corresponding to multiple frequencies;
根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;有限频段位于最低频率与最高频率之间;Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting function in the limited frequency band; The limited frequency band is located between the lowest frequency and the highest frequency;
根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;According to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency The data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency The real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency, the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary high epitaxial frequency band is between the highest frequency and the highest epitaxial frequency of the imaginary part. Between frequencies
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; among them, the singular point is located between the highest frequency and the lowest frequency;
根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part in the low epitaxial frequency band of the real part. Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral;
判断模拟虚部积分是否小于预设分辨率;当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain the conductance information and the real part information of the polarization process, otherwise update the real part interpolation step long;
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part in the low epitaxial frequency band of the imaginary part. Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral;
判断模拟实部积分是否小于预设分辨率;当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain the capacitance information and the imaginary part information of the polarization process, otherwise update the imaginary part interpolation step long;
根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。According to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
综上,本发明实施例的计算机设备根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。In summary, the computer device of the embodiment of the present invention iteratively calculates the simulated imaginary integral according to the real interpolation step, until the simulated imaginary integral is less than the preset resolution, and compares the simulated imaginary integral with the actual imaginary integral to obtain the conductance Information and real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step, until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral to obtain capacitance information And the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, and avoid misjudgment, missed judgment and wrong judgment. , Saving the cost of power equipment operation and maintenance.
本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时可以实现电介质状态分析方法的全部或部分内容,例如,处理器执行计算机程序时可以实现如下内容:The embodiment of the present invention also provides a computer-readable storage medium on which a computer program is stored. When the computer program is executed by a processor, all or part of the content of the dielectric state analysis method can be realized. For example, when the processor executes the computer program, To achieve the following:
获取多个频率对应的介电响应参数的实部数据和虚部数据;Acquire real data and imaginary data of dielectric response parameters corresponding to multiple frequencies;
根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
拟合多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;有限频段位于最低频率与最高频率之间;Fit the real part data corresponding to multiple frequencies to obtain the real part dielectric response fitting function in a limited frequency band; fit the imaginary part data corresponding to multiple frequencies to obtain the imaginary part dielectric response fitting function in the limited frequency band; The limited frequency band is located between the lowest frequency and the highest frequency;
根据最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;实部低外延频段位于最低频率与实部最低外延频率之间,实部高外延频段位于最高频率与实部最高外延频率之间,虚部低外延频段位于最低频率与虚部最低外延频率之间,虚部高外延频段位于最高频率与虚部最高外延频率之间;According to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the real part corresponding to the second highest frequency The data obtains the real part dielectric response fitting function in the real high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part dielectric response fitting in the imaginary low epitaxial frequency band is obtained Function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band; the real part low epitaxial frequency band is located at the lowest frequency and the real part lowest epitaxial frequency The real high epitaxial frequency band is between the highest frequency and the highest real epitaxial frequency, the imaginary low epitaxial frequency band is between the lowest frequency and the lowest epitaxial frequency of the imaginary part, and the imaginary high epitaxial frequency band is between the highest frequency and the highest epitaxial frequency of the imaginary part. Between frequencies
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,奇异点位于最高频率与最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; among them, the singular point is located between the highest frequency and the lowest frequency;
根据实部奇异点上界、实部奇异点下界、奇异点、实部最低外延频率、实部最高外延频率、有限频段内的实部介电响应拟合函数、实部低外延频段内的实部介电响应拟合函数和实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the singular point of the real part, the lower bound of the singular point of the real part, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the real part dielectric response fitting function in the limited frequency band, and the real part in the low epitaxial frequency band of the real part. Partial dielectric response fitting function and real part dielectric response fitting function in the high epitaxial band of the real part calculate and simulate the imaginary part integral;
判断模拟虚部积分是否小于预设分辨率;当小于预设分辨率时,将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain the conductance information and the real part information of the polarization process, otherwise update the real part interpolation step long;
执行如下迭代处理:Perform the following iterative processing:
根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
根据虚部奇异点上界、虚部奇异点下界、奇异点、虚部最低外延频率、虚部最高外延频率、有限频段内的虚部介电响应拟合函数、虚部低外延频段内的虚部介电响应拟合函数和虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the imaginary part dielectric response fitting function in the limited frequency band, and the imaginary part in the low epitaxial frequency band of the imaginary part. Partial dielectric response fitting function and the imaginary part dielectric response fitting function in the imaginary high epitaxial frequency band calculate and simulate the real part integral;
判断模拟实部积分是否小于预设分辨率;当小于预设分辨率时,将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain the capacitance information and the imaginary part information of the polarization process, otherwise update the imaginary part interpolation step long;
根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态。According to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, the dielectric state is obtained.
综上,本发明实施例的计算机可读存储介质根据实部插值步长迭代计算模拟虚部积分,直至模拟虚部积分小于预设分辨率,并将模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息;根据虚部插值步长迭代计算模拟实部积分,直至模拟实部积分小于预设分辨率,并将模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,最后根据电导信息、极化过程实部信息、电容信息和极化过程虚部信息获得电介质状态,可以准确分析判断电介质状态,避免造成误判、漏判和错判,节约了电力设备运行维护的成本。In summary, the computer-readable storage medium of the embodiment of the present invention iteratively calculates the simulated imaginary part integral according to the real part interpolation step length until the simulated imaginary part integral is less than the preset resolution, and compares the simulated imaginary part integral with the actual imaginary part integral , Obtain the conductance information and the real part information of the polarization process; iteratively calculate the simulated real part integral according to the imaginary part interpolation step length until the simulated real part integral is less than the preset resolution, and compare the simulated real part integral with the actual real part integral, Obtain the capacitance information and the imaginary part information of the polarization process, and finally obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information and the imaginary part information of the polarization process, which can accurately analyze and judge the dielectric state, and avoid misjudgments and missed judgments. And misjudgment, saving the cost of power equipment operation and maintenance.
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further describe the purpose, technical solutions and beneficial effects of the present invention in further detail. It should be understood that the above descriptions are only specific embodiments of the present invention and are not intended to limit the scope of the present invention. The protection scope, any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included in the protection scope of the present invention.
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of the present invention can be provided as a method, a system, or a computer program product. Therefore, the present invention may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present invention may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to flowcharts and/or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and/or block in the flowchart and/or block diagram, and the combination of processes and/or blocks in the flowchart and/or block diagram can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing equipment to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing equipment are generated It is a device that realizes the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions can also be stored in a computer-readable memory that can guide a computer or other programmable data processing equipment to work in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture including the instruction device. The device implements the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded on a computer or other programmable data processing equipment, so that a series of operation steps are executed on the computer or other programmable equipment to produce computer-implemented processing, so as to execute on the computer or other programmable equipment. The instructions provide steps for implementing the functions specified in one process or multiple processes in the flowchart and/or one block or multiple blocks in the block diagram.
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。In the present invention, specific examples are used to illustrate the principles and implementation of the present invention. The description of the above examples is only used to help understand the method and core idea of the present invention; at the same time, for those of ordinary skill in the art, according to this The idea of the invention will change in the specific implementation and the scope of application. In summary, the content of this specification should not be construed as limiting the invention.

Claims (18)

  1. 一种电介质状态分析方法,其特征在于,包括:A method for analyzing the state of a dielectric, which is characterized in that it includes:
    获取多个频率对应的介电响应参数的实部数据和虚部数据;Acquire real data and imaginary data of dielectric response parameters corresponding to multiple frequencies;
    根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;Determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
    拟合所述多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合所述多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;所述有限频段位于所述最低频率与所述最高频率之间;Fit the real part data corresponding to the multiple frequencies to obtain the real part dielectric response fitting function in the limited frequency band; fit the imaginary part data corresponding to the multiple frequencies to obtain the imaginary part dielectric response in the limited frequency band Fitting function; the limited frequency band is located between the lowest frequency and the highest frequency;
    根据所述最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据所述最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据所述最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据所述最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;所述实部低外延频段位于所述最低频率与所述实部最低外延频率之间,所述实部高外延频段位于所述最高频率与所述实部最高外延频率之间,所述虚部低外延频段位于所述最低频率与所述虚部最低外延频率之间,所述虚部高外延频段位于所述最高频率与所述虚部最高外延频率之间;According to the real part data corresponding to the lowest frequency and the real part data corresponding to the second lowest frequency, the real part dielectric response fitting function in the low epitaxial band of the real part is obtained, and according to the real part data corresponding to the highest frequency and the second highest frequency The corresponding real part data obtains the real part dielectric response fitting function in the real part high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency, the imaginary part in the imaginary part low epitaxial frequency band is obtained. Partial dielectric response fitting function, according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second-highest frequency to obtain the imaginary part dielectric response fitting function in the imaginary part high-extension frequency band; the real part low-extension The frequency band is located between the lowest frequency and the lowest epitaxial frequency of the real part, the high epitaxial frequency band of the real part is located between the highest frequency and the highest epitaxial frequency of the real part, and the low epitaxial frequency band of the imaginary part is located at the Between the lowest frequency and the lowest epitaxial frequency of the imaginary part, the imaginary high epitaxial frequency band is located between the highest frequency and the highest epitaxial frequency of the imaginary part;
    执行如下迭代处理:Perform the following iterative processing:
    根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,所述奇异点位于所述最高频率与所述最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; wherein the singular point is located between the highest frequency and the lowest frequency;
    根据所述实部奇异点上界、所述实部奇异点下界、所述奇异点、所述实部最低外延频率、所述实部最高外延频率、所述有限频段内的实部介电响应拟合函数、所述实部低外延频段内的实部介电响应拟合函数和所述实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the real part singular point, the lower bound of the real part singular point, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, and the real part dielectric response in the limited frequency band A fitting function, a real part dielectric response fitting function in the real part low epitaxial frequency band, and a real part dielectric response fitting function in the real part high epitaxial frequency band to calculate a simulated imaginary part integral;
    判断所述模拟虚部积分是否小于预设分辨率;当小于所述预设分辨率时,将所述模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新所述实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain conductance information and polarization process real part information, Otherwise update the real part interpolation step;
    执行如下迭代处理:Perform the following iterative processing:
    根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
    根据所述虚部奇异点上界、所述虚部奇异点下界、所述奇异点、所述虚部最低外延频率、所述虚部最高外延频率、所述有限频段内的虚部介电响应拟合函数、所述虚部低外延频段内的虚部介电响应拟合函数和所述虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part singular point, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, and the imaginary part dielectric response in the limited frequency band A fitting function, an imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, and an imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band to calculate an analog real part integral;
    判断所述模拟实部积分是否小于预设分辨率;当小于所述预设分辨率时,将所述模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新所述虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain capacitance information and imaginary part information of the polarization process, Otherwise update the imaginary part interpolation step;
    根据所述电导信息、所述极化过程实部信息、所述电容信息和所述极化过程虚部信息获得电介质状态。The dielectric state is obtained according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
  2. 根据权利要求1所述的电介质状态分析方法,其特征在于,The dielectric state analysis method according to claim 1, wherein:
    计算模拟虚部积分包括:The calculation of the simulated imaginary integral includes:
    根据所述实部奇异点上界、实部奇异点上界对应的实部数据、所述奇异点、奇异点对应的实部数据、所述实部奇异点下界和实部奇异点下界对应的实部数据计算实部奇异点频段积分;According to the real data corresponding to the upper bound of the real singular point, the real data corresponding to the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point and the lower bound of the real singular point correspond to Real part data to calculate the real part singular point frequency band integral;
    根据实部最低外延频率对应的实部数据、所述实部最低外延频率、实部最高外延频率对应的实部数据、所述实部最高外延频率和所述奇异点计算实部截断频域积分;Calculate the real part truncated frequency domain integral according to the real part data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part, the real part data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part and the singular point ;
    根据所述有限频段内的实部介电响应拟合函数、所述实部低外延频段内的实部介电响应拟合函数、所述实部高外延频段内的实部介电响应拟合函数、所述实部最低外延频率、所述实部最高外延频率、所述实部奇异点上界和所述实部奇异点下界,计算实部外延频段积分;According to the real part dielectric response fitting function in the limited frequency band, the real part dielectric response fitting function in the real part low epitaxial frequency band, and the real part dielectric response fitting function in the real part high epitaxial frequency band Function, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the upper bound of the real part singular point and the lower bound of the real part singular point, calculating the real part epitaxial frequency band integral;
    根据所述实部奇异点频段积分、所述实部截断频域积分和所述实部外延频段积分计算模拟虚部积分;Calculating an analog imaginary part integral according to the real part singular point frequency band integral, the real part truncated frequency domain integral, and the real part epitaxial frequency band integral;
    计算模拟实部积分包括:Calculating the integral of the real part of the simulation includes:
    根据所述虚部奇异点上界、虚部奇异点上界对应的虚部数据、所述奇异点、奇异点对应的虚部数据、所述虚部奇异点下界和虚部奇异点下界对应的虚部数据计算虚部奇异点频段积分;According to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the lower bound of the imaginary part singular point corresponding to The imaginary part data calculates the imaginary part singular point frequency band integral;
    根据所述有限频段内的虚部介电响应拟合函数、所述虚部低外延频段内的虚部介电响应拟合函数、所述虚部高外延频段内的虚部介电响应拟合函数、所述虚部最低外延频率、所述虚部最高外延频率、所述虚部奇异点上界和所述虚部奇异点下界,计算虚部外延频段积分;According to the imaginary part dielectric response fitting function in the limited frequency band, the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, and the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band Function, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point, calculating the integral of the imaginary part epitaxial frequency band;
    根据所述虚部奇异点频段积分和所述虚部外延频段积分计算模拟实部积分。The analog real part integral is calculated according to the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral.
  3. 根据权利要求2所述的电介质状态分析方法,其特征在于,The dielectric state analysis method according to claim 2, wherein:
    根据所述实部奇异点上界、实部奇异点上界对应的实部数据、所述奇异点、奇异点对应的实部数据、所述实部奇异点下界和实部奇异点下界对应的实部数据计算实部奇异点频段积分,包括:According to the real data corresponding to the upper bound of the real singular point, the real data corresponding to the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point and the lower bound of the real singular point correspond to The real part data calculates the real part singular point frequency band integral, including:
    根据所述实部奇异点上界、实部奇异点上界对应的实部数据、所述奇异点、奇异点对应的实部数据、所述实部奇异点下界和实部奇异点下界对应的实部数据确定第一奇异点实部系数、第二奇异点实部系数、第三奇异点实部系数和第四奇异点实部系数;According to the real data corresponding to the upper bound of the real singular point, the real data corresponding to the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point and the lower bound of the real singular point correspond to The real part data determines the first singular point real part coefficient, the second singular point real part coefficient, the third singular point real part coefficient, and the fourth singular point real part coefficient;
    根据所述第一奇异点实部系数、所述第二奇异点实部系数、所述奇异点、所述实部奇异点上界和所述实部奇异点下界计算实部奇异点上频段积分;Calculate the real part singular point upper band integral based on the first singular point real part coefficient, the second singular point real part coefficient, the singular point, the real part singular point upper bound, and the real part singular point lower bound ;
    根据所述第三奇异点实部系数、所述第四奇异点实部系数、所述奇异点、所述实部奇异点上界和所述实部奇异点下界计算实部奇异点下频段积分;Calculate the lower band integral of the real singular point according to the real part coefficient of the third singular point, the real part coefficient of the fourth singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point ;
    将所述实部奇异点上频段积分与所述实部奇异点下频段积分相加,得到实部奇异点频段积分;Adding the upper frequency band integral of the real part singular point and the lower frequency band integral of the real part singular point to obtain the real part singular point frequency band integral;
    根据所述虚部奇异点上界、虚部奇异点上界对应的虚部数据、所述奇异点、奇异点对应的虚部数据、所述虚部奇异点下 界和虚部奇异点下界对应的虚部数据计算虚部奇异点频段积分,包括:According to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the lower bound of the imaginary part singular point corresponding to The imaginary part data calculates the imaginary part singular point frequency band integral, including:
    根据所述虚部奇异点上界、虚部奇异点上界对应的虚部数据、所述奇异点、奇异点对应的虚部数据、所述虚部奇异点下界和虚部奇异点下界对应的虚部数据确定第一奇异点虚部系数、第二奇异点虚部系数、第三奇异点虚部系数和第四奇异点虚部系数;According to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the lower bound of the imaginary part singular point corresponding to The imaginary part data determines the first singular point imaginary part coefficient, the second singular point imaginary part coefficient, the third singular point imaginary part coefficient, and the fourth singular point imaginary part coefficient;
    根据所述第一奇异点虚部系数、所述第二奇异点虚部系数、所述奇异点、所述虚部奇异点上界和所述虚部奇异点下界计算虚部奇异点上频段积分;Calculate the upper band integral of the imaginary part singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point ;
    根据所述第三奇异点虚部系数、所述第四奇异点虚部系数、所述奇异点、所述虚部奇异点上界和所述虚部奇异点下界计算虚部奇异点下频段积分;Calculate the lower band integral of the imaginary part singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point ;
    将所述虚部奇异点上频段积分与所述虚部奇异点下频段积分相加,得到虚部奇异点频段积分。The upper frequency band integral of the imaginary part singular point and the lower frequency band integral of the imaginary part singular point are added to obtain the frequency band integral of the imaginary part singular point.
  4. 根据权利要求3所述的电介质状态分析方法,其特征在于,通过如下公式计算实部奇异点上频段积分:The dielectric state analysis method according to claim 3, wherein the frequency band integral on the singular point of the real part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100001
    Figure PCTCN2020109652-appb-100001
    其中,E 1为实部奇异点上频段积分,a为第一奇异点实部系数,b为第二奇异点实部系数,ω S为奇异点,ω S-H为实部奇异点上界,ω S-L为实部奇异点下界; Among them, E 1 is the frequency band integral on the real singular point, a is the real coefficient of the first singular point, b is the real coefficient of the second singular point, ω S is the singular point, ω SH is the upper bound of the real singular point, ω SL is the lower bound of the singular point of the real part;
    通过如下公式计算实部奇异点下频段积分:Calculate the frequency band integral under the singular point of the real part by the following formula:
    Figure PCTCN2020109652-appb-100002
    Figure PCTCN2020109652-appb-100002
    其中,E 2为实部奇异点下频段积分,c为第三奇异点实部系数,d为第四奇异点实部系数; Among them, E 2 is the frequency band integral of the real part of the singular point, c is the real part coefficient of the third singular point, and d is the real part coefficient of the fourth singular point;
    通过如下公式计算虚部奇异点上频段积分:The frequency band integral on the singular point of the imaginary part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100003
    Figure PCTCN2020109652-appb-100003
    其中,E' 1为虚部奇异点上频段积分,a'为第一奇异点虚部系数,b'为第二奇异点虚部系数,ω' S-H为虚部奇异点上界,ω' S-L为虚部奇异点下界; Wherein, E '1 is the imaginary part of the integral bands singular point, a' is a first singular point of the imaginary part of the coefficient, b 'is the imaginary part of the second singular point coefficients, ω' SH bounded singularity point on the imaginary part, ω 'SL Is the lower bound of the imaginary part singular point;
    通过如下公式计算虚部奇异点下频段积分:Calculate the lower band integral of the imaginary part singular point by the following formula:
    Figure PCTCN2020109652-appb-100004
    Figure PCTCN2020109652-appb-100004
    其中,E' 2为虚部奇异点下频段积分,c'为第三奇异点虚部系数,d'为第四奇异点虚部系数。 Wherein, E '2 is a lower frequency band integrating the imaginary part of a singular point, c' is the coefficient of the imaginary part of the third singular point, d 'is the imaginary part of a singular point of the fourth coefficient.
  5. 根据权利要求2所述的电介质状态分析方法,其特征在于,根据实部最低外延频率对应的实部数据、所述实部最低外延频率、实部最高外延频率对应的实部数据、所述实部最高外延频率和所述奇异点计算实部截断频域积分,包括:The dielectric state analysis method according to claim 2, wherein the real part data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part, the real part data corresponding to the highest epitaxial frequency of the real part, the real part The highest extension frequency of the part and the singular point calculation real part truncated frequency domain integral, including:
    根据所述实部最高外延频率对应的实部数据、所述实部最高外延频率和所述奇异点计算实部上截断频域积分;Calculating a truncated frequency domain integral on the real part according to the real part data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, and the singular point;
    根据所述实部最低外延频率对应的实部数据、所述实部最低外延频率和所述奇异点计算实部下截断频域积分;Calculating the lower truncated frequency domain integral of the real part according to the real part data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part, and the singular point;
    将所述实部上截断频域积分与所述实部下截断频域积分相加,得到实部截断频域积分。The upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
  6. 根据权利要求5所述的电介质状态分析方法,其特征在于,The dielectric state analysis method according to claim 5, wherein:
    通过如下公式计算实部上截断频域积分:The truncated frequency domain integral on the real part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100005
    Figure PCTCN2020109652-appb-100005
    其中,F 1为实部上截断频域积分,χ'(ω Ext-H)为实部最高外延频率对应的实部数据,ω S为奇异点,ω Ext-H为实部最高外延频率; Among them, F 1 is the truncated frequency domain integral on the real part, χ'(ω Ext-H ) is the real part data corresponding to the highest extension frequency of the real part, ω S is the singularity point, and ω Ext-H is the highest extension frequency of the real part;
    通过如下公式计算实部下截断频域积分:Calculate the truncated frequency domain integral under the real part by the following formula:
    Figure PCTCN2020109652-appb-100006
    Figure PCTCN2020109652-appb-100006
    其中,F 2为实部下截断频域积分,χ'(ω Ext-L)为实部最低外延频率对应的实部数据,ω Ext-L为实部最低外延频率。 Among them, F 2 is the lower truncated frequency domain integral of the real part, χ'(ω Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part, and ω Ext-L is the lowest epitaxial frequency of the real part.
  7. 根据权利要求1所述的电介质状态分析方法,其特征在于,The dielectric state analysis method according to claim 1, wherein:
    通过如下公式更新所述实部插值步长:Update the real part interpolation step size by the following formula:
    Figure PCTCN2020109652-appb-100007
    Figure PCTCN2020109652-appb-100007
    其中,ln(Δω n)为第n次迭代中的实部插值步长,ln(Δω n+1)为第n+1次迭代中的实部插值步长; Among them, ln(Δω n ) is the real part interpolation step in the nth iteration, and ln(Δω n+1 ) is the real part interpolation step in the n+1th iteration;
    通过如下公式更新所述虚部插值步长:Update the imaginary part interpolation step size by the following formula:
    Figure PCTCN2020109652-appb-100008
    Figure PCTCN2020109652-appb-100008
    其中,ln(Δω' n)为第n次迭代中的虚部插值步长,ln(Δω' n+1)为第n+1次迭代中的虚部插值步长。 Among them, ln(Δω' n ) is the imaginary part interpolation step in the nth iteration, and ln(Δω' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
  8. 根据权利要求1所述的电介质状态分析方法,其特征在于,The dielectric state analysis method according to claim 1, wherein:
    获得电导信息和极化过程实部信息包括:Obtaining conductance information and real part information of polarization process includes:
    将所述模拟虚部积分与所述实际虚部积分中的相同部分作为极化过程实部信息,将所述模拟虚部积分与所述实际虚部积分中的不同部分作为电导信息;Taking the same part of the simulated imaginary part integral and the actual imaginary part integral as the real part information of the polarization process, and taking the different part of the simulated imaginary part integral and the actual imaginary part integral as the conductance information;
    获得电容信息和极化过程虚部信息包括:Obtaining capacitance information and imaginary part information of the polarization process include:
    将所述模拟实部积分与所述实际实部积分中的相同部分作为极化过程虚部信息,将所述模拟实部积分与所述实际实部积分中的不同部分作为电容信息。The same part in the simulated real part integral and the actual real part integral is used as the imaginary part information of the polarization process, and the different part in the simulated real part integral and the actual real part integral is used as capacitance information.
  9. 一种电介质状态分析系统,其特征在于,包括:A dielectric state analysis system, which is characterized in that it comprises:
    获取单元,用于获取多个频率对应的介电响应参数的实部数据和虚部数据;An acquiring unit for acquiring real part data and imaginary part data of dielectric response parameters corresponding to multiple frequencies;
    确定单元,用于根据最低频率确定实部最低外延频率和虚部最低外延频率,根据最高频率确定实部最高外延频率和虚部最高外延频率;The determining unit is used to determine the lowest epitaxial frequency of the real part and the lowest epitaxial frequency of the imaginary part according to the lowest frequency, and determine the highest epitaxial frequency of the real part and the highest epitaxial frequency of the imaginary part according to the highest frequency;
    有限频段拟合函数单元,用于拟合所述多个频率对应的实部数据,获得有限频段内的实部介电响应拟合函数;拟合所述多个频率对应的虚部数据,获得有限频段内的虚部介电响应拟合函数;其中,所述有限频段位于所述最低频率与所述最高频率之间;The finite frequency band fitting function unit is used to fit the real part data corresponding to the multiple frequencies to obtain the real part dielectric response fitting function in the finite frequency band; fit the imaginary part data corresponding to the multiple frequencies to obtain Imaginary part dielectric response fitting function in a limited frequency band; wherein, the limited frequency band is located between the lowest frequency and the highest frequency;
    外延频段拟合函数单元,用于根据所述最低频率对应的实部数据和次低频率对应的实部数据获得实部低外延频段内的实部介电响应拟合函数,根据所述最高频率对应的实部数据和次高频率对应的实部数据获得实部高外延频段内的实部介电响应拟合函数;根据所述最低频率对应的虚部数据和次低频率对应的虚部数据获得虚部低外延频段内的虚部介电响应拟合函数,根据所述最高频率对应的虚部数据和次高频率对应的虚部数据获得虚部高外延频段内的虚部介电响应拟合函数;所述实部低外延频段位于所述最低频率与所述实部最低外延频率之间,所述实部高外延频段位于所述最高频率与所述实部最高外延频率之间,所述虚部低外延频段位于所述最低频率与所述虚部最低外延频率之间,所述虚部高外延频段位于所述最高频率与所述虚部最高外延频率之间;The epitaxial frequency band fitting function unit is used to obtain the real part dielectric response fitting function in the real low epitaxial frequency band according to the real part data corresponding to the lowest frequency and the real part data corresponding to the second-lowest frequency, and according to the highest frequency The corresponding real part data and the real part data corresponding to the second highest frequency obtain the real part dielectric response fitting function in the real part high epitaxial frequency band; according to the imaginary part data corresponding to the lowest frequency and the imaginary part data corresponding to the second lowest frequency The imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band is obtained, and the imaginary part dielectric response simulation function in the imaginary part high epitaxial frequency band is obtained according to the imaginary part data corresponding to the highest frequency and the imaginary part data corresponding to the second highest frequency. Resultant function; the real low epitaxial frequency band is located between the lowest frequency and the real lowest epitaxial frequency, the real high epitaxial frequency band is located between the highest frequency and the real highest epitaxial frequency, so The imaginary part low epitaxial frequency band is located between the lowest frequency and the imaginary part lowest epitaxial frequency, and the imaginary part high epitaxial frequency band is located between the highest frequency and the imaginary part highest epitaxial frequency;
    实部迭代单元,用于执行如下迭代处理:The real part iterative unit is used to perform the following iterative processing:
    根据奇异点和实部插值步长计算实部奇异点上界和实部奇异点下界;其中,所述奇异点位于所述最高频率与所述最低频率之间;Calculate the upper bound of the real singular point and the lower bound of the real singular point according to the singular point and the real part interpolation step length; wherein the singular point is located between the highest frequency and the lowest frequency;
    根据所述实部奇异点上界、所述实部奇异点下界、所述奇异点、所述实部最低外延频率、所述实部最高外延频率、所述有限频段内的实部介电响应拟合函数、所述实部低外延频段内的实部介电响应拟合函数和所述实部高外延频段内的实部介电响应拟合函数计算模拟虚部积分;According to the upper bound of the real part singular point, the lower bound of the real part singular point, the singular point, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, and the real part dielectric response in the limited frequency band A fitting function, a real part dielectric response fitting function in the real part low epitaxial frequency band, and a real part dielectric response fitting function in the real part high epitaxial frequency band to calculate a simulated imaginary part integral;
    判断所述模拟虚部积分是否小于预设分辨率;当小于所述预设分辨率时,将所述模拟虚部积分与实际虚部积分进行比较,获得电导信息和极化过程实部信息,否则更新所述实部插值步长;Determine whether the simulated imaginary part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated imaginary part integral with the actual imaginary part integral to obtain conductance information and polarization process real part information, Otherwise update the real part interpolation step;
    虚部迭代单元,用于执行如下迭代处理:The imaginary part iterative unit is used to perform the following iterative processing:
    根据奇异点和虚部插值步长计算虚部奇异点上界和虚部奇异点下界;Calculate the upper bound and lower bound of the imaginary part singular point according to the singular point and the imaginary part interpolation step;
    根据所述虚部奇异点上界、所述虚部奇异点下界、所述奇异点、所述虚部最低外延频率、所述虚部最高外延频率、所述有限频段内的虚部介电响应拟合函数、所述虚部低外延频段内的虚部介电响应拟合函数和所述虚部高外延频段内的虚部介电响应拟合函数计算模拟实部积分;According to the upper bound of the imaginary part singular point, the lower bound of the imaginary part singular point, the singular point, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, and the imaginary part dielectric response in the limited frequency band A fitting function, an imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, and an imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band to calculate an analog real part integral;
    判断所述模拟实部积分是否小于预设分辨率;当小于所述预设分辨率时,将所述模拟实部积分与实际实部积分进行比较,获得电容信息和极化过程虚部信息,否则更新所述虚部插值步长;Determine whether the simulated real part integral is less than the preset resolution; when it is smaller than the preset resolution, compare the simulated real part integral with the actual real part integral to obtain capacitance information and imaginary part information of the polarization process, Otherwise update the imaginary part interpolation step;
    电介质状态单元,用于根据所述电导信息、所述极化过程实部信息、所述电容信息和所述极化过程虚部信息获得电介质状态。The dielectric state unit is configured to obtain the dielectric state according to the conductance information, the real part information of the polarization process, the capacitance information, and the imaginary part information of the polarization process.
  10. 根据权利要求9所述的电介质状态分析系统,其特征在于,The dielectric state analysis system according to claim 9, wherein:
    所述实部迭代单元具体用于:The real part iteration unit is specifically used for:
    根据所述实部奇异点上界、实部奇异点上界对应的实部数据、所述奇异点、奇异点对应的实部数据、所述实部奇异点下界和实部奇异点下界对应的实部数据计算实部奇异点频段积分;According to the real data corresponding to the upper bound of the real singular point, the real data corresponding to the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point and the lower bound of the real singular point correspond to Real part data to calculate the real part singular point frequency band integral;
    根据实部最低外延频率对应的实部数据、所述实部最低外延频率、实部最高外延频率对应的实部数据、所述实部最高外延频率和所述奇异点计算实部截断频域积分;Calculate the real part truncated frequency domain integral according to the real part data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part, the real part data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part and the singular point ;
    根据所述有限频段内的实部介电响应拟合函数、所述实部低外延频段内的实部介电响应拟合函数、所述实部高外延频段内的实部介电响应拟合函数、所述实部最低外延频率、所述实部最高外延频率、所述实部奇异点上界和所述实部奇异点下界,计算实部外延频段积分;According to the real part dielectric response fitting function in the limited frequency band, the real part dielectric response fitting function in the real part low epitaxial frequency band, and the real part dielectric response fitting function in the real part high epitaxial frequency band Function, the lowest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, the upper bound of the real part singular point and the lower bound of the real part singular point, calculating the real part epitaxial frequency band integral;
    根据所述实部奇异点频段积分、所述实部截断频域积分和所述实部外延频段积分计算模拟虚部积分;Calculating an analog imaginary part integral according to the real part singular point frequency band integral, the real part truncated frequency domain integral, and the real part epitaxial frequency band integral;
    所述虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
    根据所述虚部奇异点上界、虚部奇异点上界对应的虚部数据、所述奇异点、奇异点对应的虚部数据、所述虚部奇异点下界和虚部奇异点下界对应的虚部数据计算虚部奇异点频段积分;According to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the lower bound of the imaginary part singular point corresponding to The imaginary part data calculates the imaginary part singular point frequency band integral;
    根据所述有限频段内的虚部介电响应拟合函数、所述虚部低外延频段内的虚部介电响应拟合函数、所述虚部高外延频段内的虚部介电响应拟合函数、所述虚部最低外延频率、所述虚部最高外延频率、所述虚部奇异点上界和所述虚部奇异点下界,计算虚部外延频段积分;According to the imaginary part dielectric response fitting function in the limited frequency band, the imaginary part dielectric response fitting function in the imaginary part low epitaxial frequency band, and the imaginary part dielectric response fitting function in the imaginary part high epitaxial frequency band Function, the lowest epitaxial frequency of the imaginary part, the highest epitaxial frequency of the imaginary part, the upper bound of the imaginary part singular point and the lower bound of the imaginary part singular point, calculating the integral of the imaginary part epitaxial frequency band;
    根据所述虚部奇异点频段积分和所述虚部外延频段积分计算模拟实部积分。The analog real part integral is calculated according to the imaginary part singular point frequency band integral and the imaginary part extension frequency band integral.
  11. 根据权利要求10所述的电介质状态分析系统,其特征在于,所述实部迭代单元具体用于:The dielectric state analysis system according to claim 10, wherein the real part iteration unit is specifically configured to:
    根据所述实部奇异点上界、实部奇异点上界对应的实部数据、所述奇异点、奇异点对应的实部数据、所述实部奇异点下界和实部奇异点下界对应的实部数据确定第一奇异点实部系数、第二奇异点实部系数、第三奇异点实部系数和第四奇异点实部系数;According to the real data corresponding to the upper bound of the real singular point, the real data corresponding to the upper bound of the real singular point, the singular point, the real data corresponding to the singular point, the lower bound of the real singular point and the lower bound of the real singular point correspond to The real part data determines the first singular point real part coefficient, the second singular point real part coefficient, the third singular point real part coefficient, and the fourth singular point real part coefficient;
    根据所述第一奇异点实部系数、所述第二奇异点实部系数、所述奇异点、所述实部奇异点上界和所述实部奇异点下界计 算实部奇异点上频段积分;Calculate the real part singular point upper band integral based on the first singular point real part coefficient, the second singular point real part coefficient, the singular point, the real part singular point upper bound, and the real part singular point lower bound ;
    根据所述第三奇异点实部系数、所述第四奇异点实部系数、所述奇异点、所述实部奇异点上界和所述实部奇异点下界计算实部奇异点下频段积分;Calculate the lower band integral of the real singular point according to the real part coefficient of the third singular point, the real part coefficient of the fourth singular point, the singular point, the upper bound of the real singular point, and the lower bound of the real singular point ;
    将所述实部奇异点上频段积分与所述实部奇异点下频段积分相加,得到实部奇异点频段积分;Adding the upper frequency band integral of the real part singular point and the lower frequency band integral of the real part singular point to obtain the real part singular point frequency band integral;
    所述虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
    根据所述虚部奇异点上界、虚部奇异点上界对应的虚部数据、所述奇异点、奇异点对应的虚部数据、所述虚部奇异点下界和虚部奇异点下界对应的虚部数据确定第一奇异点虚部系数、第二奇异点虚部系数、第三奇异点虚部系数和第四奇异点虚部系数;According to the upper bound of the imaginary part singular point, the imaginary part data corresponding to the upper bound of the imaginary part singular point, the singular point, the imaginary part data corresponding to the singular point, the lower bound of the imaginary part singular point and the lower bound of the imaginary part singular point corresponding to The imaginary part data determines the first singular point imaginary part coefficient, the second singular point imaginary part coefficient, the third singular point imaginary part coefficient, and the fourth singular point imaginary part coefficient;
    根据所述第一奇异点虚部系数、所述第二奇异点虚部系数、所述奇异点、所述虚部奇异点上界和所述虚部奇异点下界计算虚部奇异点上频段积分;Calculate the upper band integral of the imaginary part singular point according to the imaginary part coefficient of the first singular point, the imaginary part coefficient of the second singular point, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point ;
    根据所述第三奇异点虚部系数、所述第四奇异点虚部系数、所述奇异点、所述虚部奇异点上界和所述虚部奇异点下界计算虚部奇异点下频段积分;Calculate the lower band integral of the imaginary part singular point according to the third singular point imaginary part coefficient, the fourth singular point imaginary part coefficient, the singular point, the upper bound of the imaginary part singular point, and the lower bound of the imaginary part singular point ;
    将所述虚部奇异点上频段积分与所述虚部奇异点下频段积分相加,得到虚部奇异点频段积分。The upper frequency band integral of the imaginary part singular point and the lower frequency band integral of the imaginary part singular point are added to obtain the frequency band integral of the imaginary part singular point.
  12. 根据权利要求11所述的电介质状态分析系统,其特征在于,通过如下公式计算实部奇异点上频段积分:The dielectric state analysis system according to claim 11, wherein the frequency band integral on the singular point of the real part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100009
    Figure PCTCN2020109652-appb-100009
    其中,E 1为实部奇异点上频段积分,a为第一奇异点实部系数,b为第二奇异点实部系数,ω S为奇异点,ω S-H为实部奇异点上界,ω S-L为实部奇异点下界; Among them, E 1 is the frequency band integral on the real singular point, a is the real coefficient of the first singular point, b is the real coefficient of the second singular point, ω S is the singular point, ω SH is the upper bound of the real singular point, ω SL is the lower bound of the singular point of the real part;
    通过如下公式计算实部奇异点下频段积分:Calculate the frequency band integral under the singular point of the real part by the following formula:
    Figure PCTCN2020109652-appb-100010
    Figure PCTCN2020109652-appb-100010
    其中,E 2为实部奇异点下频段积分,c为第三奇异点实部系数,d为第四奇异点实部系数; Among them, E 2 is the frequency band integral of the real part of the singular point, c is the real part coefficient of the third singular point, and d is the real part coefficient of the fourth singular point;
    通过如下公式计算虚部奇异点上频段积分:The frequency band integral on the singular point of the imaginary part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100011
    Figure PCTCN2020109652-appb-100011
    其中,E' 1为虚部奇异点上频段积分,a'为第一奇异点虚部系数,b'为第二奇异点虚部系数,ω' S-H为虚部奇异点上界,ω' S-L为虚部奇异点下界; Wherein, E '1 is the imaginary part of the integral bands singular point, a' is a first singular point of the imaginary part of the coefficient, b 'is the imaginary part of the second singular point coefficients, ω' SH bounded singularity point on the imaginary part, ω 'SL Is the lower bound of the imaginary part singular point;
    通过如下公式计算虚部奇异点下频段积分:Calculate the lower band integral of the imaginary part singular point by the following formula:
    Figure PCTCN2020109652-appb-100012
    Figure PCTCN2020109652-appb-100012
    其中,E' 2为虚部奇异点下频段积分,c'为第三奇异点虚部系数,d'为第四奇异点虚部系数。 Wherein, E '2 is a lower frequency band integrating the imaginary part of a singular point, c' is the coefficient of the imaginary part of the third singular point, d 'is the imaginary part of a singular point of the fourth coefficient.
  13. 根据权利要求10所述的电介质状态分析系统,其特征在于,所述实部迭代单元具体用于:The dielectric state analysis system according to claim 10, wherein the real part iteration unit is specifically configured to:
    根据所述实部最高外延频率对应的实部数据、所述实部最高外延频率和所述奇异点计算实部上截断频域积分;Calculating a truncated frequency domain integral on the real part according to the real part data corresponding to the highest epitaxial frequency of the real part, the highest epitaxial frequency of the real part, and the singular point;
    根据所述实部最低外延频率对应的实部数据、所述实部最低外延频率和所述奇异点计算实部下截断频域积分;Calculating the lower truncated frequency domain integral of the real part according to the real part data corresponding to the lowest epitaxial frequency of the real part, the lowest epitaxial frequency of the real part, and the singular point;
    将所述实部上截断频域积分与所述实部下截断频域积分相加,得到实部截断频域积分。The upper truncated frequency domain integral of the real part and the lower truncated frequency domain integral of the real part are added to obtain the real part truncated frequency domain integral.
  14. 根据权利要求13所述的电介质状态分析系统,其特征在于,The dielectric state analysis system according to claim 13, wherein:
    通过如下公式计算实部上截断频域积分:The truncated frequency domain integral on the real part is calculated by the following formula:
    Figure PCTCN2020109652-appb-100013
    Figure PCTCN2020109652-appb-100013
    其中,F 1为实部上截断频域积分,χ'(ω Ext-H)为实部最高外延频率对应的实部数据,ω S为奇异点,ω Ext-H为实部最高外延频率; Among them, F 1 is the truncated frequency domain integral on the real part, χ'(ω Ext-H ) is the real part data corresponding to the highest extension frequency of the real part, ω S is the singularity point, and ω Ext-H is the highest extension frequency of the real part;
    通过如下公式计算实部下截断频域积分:Calculate the truncated frequency domain integral under the real part by the following formula:
    Figure PCTCN2020109652-appb-100014
    Figure PCTCN2020109652-appb-100014
    其中,F 2为实部下截断频域积分,χ'(ω Ext-L)为实部最低外延频率对应的实部数据,ω Ext-L为实部最低外延频率。 Among them, F 2 is the lower truncated frequency domain integral of the real part, χ'(ω Ext-L ) is the real part data corresponding to the lowest epitaxial frequency of the real part, and ω Ext-L is the lowest epitaxial frequency of the real part.
  15. 根据权利要求9所述的电介质状态分析系统,其特征在于,The dielectric state analysis system according to claim 9, wherein:
    通过如下公式更新所述实部插值步长:Update the real part interpolation step size by the following formula:
    Figure PCTCN2020109652-appb-100015
    Figure PCTCN2020109652-appb-100015
    其中,ln(Δω n)为第n次迭代中的实部插值步长,ln(Δω n+1)为第n+1次迭代中的实部插值步长; Among them, ln(Δω n ) is the real part interpolation step in the nth iteration, and ln(Δω n+1 ) is the real part interpolation step in the n+1th iteration;
    通过如下公式更新所述虚部插值步长:Update the imaginary part interpolation step size by the following formula:
    Figure PCTCN2020109652-appb-100016
    Figure PCTCN2020109652-appb-100016
    其中,ln(Δω' n)为第n次迭代中的虚部插值步长,ln(Δω' n+1)为第n+1次迭代中的虚部插值步长。 Among them, ln(Δω' n ) is the imaginary part interpolation step in the nth iteration, and ln(Δω' n+1 ) is the imaginary part interpolation step in the n+1th iteration.
  16. 根据权利要求9所述的电介质状态分析系统,其特征在于,所述实部迭代单元具体用于:The dielectric state analysis system according to claim 9, wherein the real part iteration unit is specifically used for:
    将所述模拟虚部积分与所述实际虚部积分中的相同部分作为极化过程实部信息,将所述模拟虚部积分与所述实际虚部积分中的不同部分作为电导信息;Taking the same part of the simulated imaginary part integral and the actual imaginary part integral as the real part information of the polarization process, and taking the different part of the simulated imaginary part integral and the actual imaginary part integral as the conductance information;
    所述虚部迭代单元具体用于:The imaginary part iteration unit is specifically used for:
    将所述模拟实部积分与所述实际实部积分中的相同部分作为极化过程虚部信息,将所述模拟实部积分与所述实际实部积分中的不同部分作为电容信息。The same part in the simulated real part integral and the actual real part integral is used as the imaginary part information of the polarization process, and the different part in the simulated real part integral and the actual real part integral is used as capacitance information.
  17. 一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至8任一项所述的电介质状态分析方法的步骤。A computer device, comprising a memory, a processor, and a computer program stored on the memory and running on the processor, wherein the processor implements the computer program described in any one of claims 1 to 8 when the processor executes the computer program. The steps of the dielectric state analysis method described.
  18. 一种计算机可读存储介质,其上存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现权利要求1至8任一项所述的电介质状态分析方法的步骤。A computer-readable storage medium with a computer program stored thereon, wherein the computer program implements the steps of the dielectric state analysis method according to any one of claims 1 to 8 when the computer program is executed by a processor.
PCT/CN2020/109652 2019-11-07 2020-08-17 Dielectric state analysis method and system, computer and storage medium WO2021088453A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201911081050.6 2019-11-07
CN201911081050.6A CN111103511B (en) 2019-11-07 2019-11-07 Dielectric state analysis method, system, computer, and storage medium

Publications (1)

Publication Number Publication Date
WO2021088453A1 true WO2021088453A1 (en) 2021-05-14

Family

ID=70420622

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2020/109652 WO2021088453A1 (en) 2019-11-07 2020-08-17 Dielectric state analysis method and system, computer and storage medium

Country Status (2)

Country Link
CN (1) CN111103511B (en)
WO (1) WO2021088453A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111103511B (en) * 2019-11-07 2021-06-15 国网冀北电力有限公司电力科学研究院 Dielectric state analysis method, system, computer, and storage medium

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4905170A (en) * 1987-11-12 1990-02-27 Forouhi Abdul R Method and apparatus of determining optical constants of amorphous semiconductors and dielectrics
CN103105566A (en) * 2013-01-14 2013-05-15 江苏省电力公司电力科学研究院 Oil paper insulation electrical equipment aging state detection method based on universal relaxation law
CN106980095A (en) * 2017-05-24 2017-07-25 中国电子科技集团公司第四十研究所 A kind of Meta Materials electromagnetic parameter inversion method based on improvement K K algorithms
CN107679327A (en) * 2017-10-10 2018-02-09 国网江苏省电力公司电力科学研究院 Paper oil insulation extension Debye model parameter identification method based on FDS
US20180045758A1 (en) * 2015-03-03 2018-02-15 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Arrangement for spatially resolved determination of the specific electrical resistance and/or the specific electrical conductivity of samples
CN111103511A (en) * 2019-11-07 2020-05-05 国网冀北电力有限公司电力科学研究院 Dielectric state analysis method, system, computer, and storage medium

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016139233A1 (en) * 2015-03-03 2016-09-09 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Arrangement for spatially resolved determination of the specific electrical resistance and/or the specific electrical conductivity of samples

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4905170A (en) * 1987-11-12 1990-02-27 Forouhi Abdul R Method and apparatus of determining optical constants of amorphous semiconductors and dielectrics
CN103105566A (en) * 2013-01-14 2013-05-15 江苏省电力公司电力科学研究院 Oil paper insulation electrical equipment aging state detection method based on universal relaxation law
US20180045758A1 (en) * 2015-03-03 2018-02-15 Fraunhofer-Gesellschaft Zur Foerderung Der Angewandten Forschung E.V. Arrangement for spatially resolved determination of the specific electrical resistance and/or the specific electrical conductivity of samples
CN106980095A (en) * 2017-05-24 2017-07-25 中国电子科技集团公司第四十研究所 A kind of Meta Materials electromagnetic parameter inversion method based on improvement K K algorithms
CN107679327A (en) * 2017-10-10 2018-02-09 国网江苏省电力公司电力科学研究院 Paper oil insulation extension Debye model parameter identification method based on FDS
CN111103511A (en) * 2019-11-07 2020-05-05 国网冀北电力有限公司电力科学研究院 Dielectric state analysis method, system, computer, and storage medium

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
刘东航 等 (LIU, DONGHANG ET AL.): "频域介电测量谱的Kramers-Kronig修正 (Kramers-Kronig Revision of Frequency-Domain Dielectric Measurement)", 《西安交通大学学报》 (《JOURNAL OF XI’AN JIAOTONG UNIVERSITY》), vol. 25, no. 5, 15 July 1991 (1991-07-15), ISSN: 0253-987X *
张正 等 (ZHANG, ZHENG ET AL.): "Kramers-Kronig关系中积分奇点的处理方法研究 (Calculation of Integral Singularities in Kramers-Kronig Relations)", 《北京师范大学学报(自然科学版)》 (JOURNAL OF BEIJING NORMAL UNIVERSITY(NATURAL SCIENCE)), vol. 53, no. 3, 15 June 2017 (2017-06-15), ISSN: 0476-0301 *
牟媛 (MU, YUAN): "电大尺寸粗糙目标太赫兹散射特性研究 (Investigation of THz Scattering Features of Electrically Large Targets with Rough Surfaces)", 《中国博士论文全文数据库(电子期刊)》 (CHINA DOCTORAL DISSERTATIONS FULL-TEXT DATABASE (ELECTRONIC JOURNALS)), no. 7, 15 July 2019 (2019-07-15), ISSN: 1674-022X *
高岩峰 等 (GAO, YANFENG ET AL.): "Kramers-Kronig变换在介电响应分析中的数值计算方法、意义及应用 (Numerical Computational Method, Application and Significance of the Kramers-Kronig Transform in the Analysis of Dielectric Response)", 中国电机工程学报 (PROCEEDINGS OF THE CSEE), vol. 40, no. 1, 5 January 2020 (2020-01-05), ISSN: 0258-8013 *

Also Published As

Publication number Publication date
CN111103511A (en) 2020-05-05
CN111103511B (en) 2021-06-15

Similar Documents

Publication Publication Date Title
CN106855605B (en) Frequency domain test analysis system and method for evaluating overall aging life of cable
CN116111727B (en) Comprehensive distribution box abnormity monitoring method based on dynamic temperature threshold
CN111079329B (en) Fatigue life assessment method based on similar structure test
CN108387609B (en) Nondestructive testing method, device and readable storage medium of equipment
CN108872820B (en) Method and system for evaluating insulation aging state of oil impregnated paper in high-voltage current transformer
WO2021088453A1 (en) Dielectric state analysis method and system, computer and storage medium
CN110095697B (en) Cable current-carrying capacity adjusting method, device, equipment and readable storage medium
CN111751199B (en) Fatigue life prediction method based on EIFS distribution
CN108828416A (en) A kind of crosslinking polyethylene-insulated cable insulation ageing state appraisal procedure based on to derivative analysis method
CN108872059B (en) Composite insulator aging state assessment method and terminal equipment
CN109142865B (en) Frequency domain spectrum identification method considering polarization equivalent circuit parameters of oiled paper insulation interface
CN110598254B (en) Translation factor-based temperature correction method for dielectric spectrum of solid insulation frequency domain of transformer
CN105445625A (en) Method for classifying alternating-current conductivity frequency domain spectrums of transformer oil clearance under difference temperatures into same reference temperature
CN110095696B (en) Cable current-carrying capacity adjusting method, device, equipment and readable storage medium
CN110321661B (en) Fractal theory-based double-pore structure soil pore dividing method
Fan et al. Polarization loss analysis and ageing characterisation of transformer oil‐immersed insulation by using decoupled frequency domain spectroscopy
CN109870635B (en) Oiled paper insulation aging state evaluation method and system based on activation energy iterative correction
Goldschmidt et al. Modelling and numerical investigations of the mechanical behavior of polyurethane under the influence of moisture
CN111638429A (en) Temperature correction method and device for insulating material state evaluation
CN108169013B (en) Fatigue life calibration method suitable for multi-nail connection test piece
CN116297229A (en) Method for detecting aging degree of crosslinked polyethylene insulating material based on yellow degree evolution
CN115840084A (en) Impedance testing method, device and equipment based on coaxial cable and storage medium
JP2015227804A (en) Remaining life calculation method, deterioration diagnosis method, deterioration diagnostic device and program
CN108375537A (en) Cable for ship insulation life detection method based on biggest quality minimizing speed
CN113792475A (en) Moisture content evaluation method considering transformer aging effect based on weighted KNN algorithm

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20885317

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 20885317

Country of ref document: EP

Kind code of ref document: A1