CN115659659A - Method and system for calculating winding temperature field under natural convection heat dissipation of transformer - Google Patents
Method and system for calculating winding temperature field under natural convection heat dissipation of transformer Download PDFInfo
- Publication number
- CN115659659A CN115659659A CN202211345515.6A CN202211345515A CN115659659A CN 115659659 A CN115659659 A CN 115659659A CN 202211345515 A CN202211345515 A CN 202211345515A CN 115659659 A CN115659659 A CN 115659659A
- Authority
- CN
- China
- Prior art keywords
- oil
- transformer
- winding
- temperature
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004804 winding Methods 0.000 title claims abstract description 105
- 230000017525 heat dissipation Effects 0.000 title claims abstract description 51
- 238000000034 method Methods 0.000 title claims abstract description 51
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims abstract description 45
- 238000012932 thermodynamic analysis Methods 0.000 claims abstract description 18
- 238000004364 calculation method Methods 0.000 claims description 18
- 239000012530 fluid Substances 0.000 claims description 10
- 230000005855 radiation Effects 0.000 claims description 6
- 230000004907 flux Effects 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 7
- 238000009529 body temperature measurement Methods 0.000 abstract description 3
- 239000000463 material Substances 0.000 description 9
- 230000008878 coupling Effects 0.000 description 7
- 238000010168 coupling process Methods 0.000 description 7
- 238000005859 coupling reaction Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004146 energy storage Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Housings And Mounting Of Transformers (AREA)
Abstract
The embodiment of the invention discloses a method and a system for calculating a winding temperature field under natural convection heat dissipation of a transformer, wherein the method comprises the following steps: performing thermodynamic analysis on oil, windings and an iron core part of a target transformer to obtain oil parameters, winding parameters and iron core parameters corresponding to the target transformer, and establishing a constraint equation according to the thermodynamic analysis result and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition; acquiring a temperature value of any point in the space where the target transformer is located; and solving to obtain a temperature field corresponding to the target transformer by using the temperature model of the target transformer based on the obtained temperature value. The temperature model provided by the invention can realize the acquisition of the winding temperature spatial distribution under the natural convection heat dissipation of the transformer under the condition of only acquiring a point of temperature measurement data, thereby greatly improving the efficiency of calculating and monitoring the temperature of the transformer winding.
Description
Technical Field
The invention relates to the technical field of temperature detection of transformers of power systems, in particular to a method and a system for calculating a winding temperature field of a transformer under natural convection heat dissipation.
Background
With the development of power systems, such as on-load capacity-regulating and voltage-regulating distribution transformers are increasingly applied to the power systems, but due to the complex application environment of the transformers, when a fault or an emergency occurs, the transformer fault reflected by the temperature change of the transformers cannot be timely sensed, the transformers are easily damaged, power failure occurs, and economic loss is caused, so that in order to ensure the safety of devices, the temperature of the transformers needs to be monitored; however, the existing transformer temperature monitoring method needs to monitor and fuse temperature data of multiple points in space and construct a corresponding finite element temperature field, and has the problems of need to measure multidimensional parameters, complex calculation and low monitoring efficiency.
The digital twin is one of core technical means for solving the problems of distributed energy with high-proportion access randomness, intermittence and fluctuation characteristics of a novel power system and interactive energy coordination control of an energy storage device, a V2G and the like in the future. The digital power grid is initially built in 2021 year no matter in the national power grid or in the south, namely advanced digital map and digital twin technology are utilized to realize digital twin of a main grid of 110kV or more, which marks that at the end of 2021 year, the national power grid and the south power grid basically realize digital twin of unit level and system level defined in digital twin power grid white paper, and the power grid develops towards equipment level digital twin in the future, so that the construction of the digital twin power grid defined by the national power grid white paper is really realized. According to the 'digital twin application white paper' of the Ministry of industry and belief, the digital twin senses, diagnoses and predicts the state of a physical entity by means including simulation, actual measurement and data analysis, namely compared with the traditional on-line monitoring, one of the core characteristics of the digital twin technology of the electrical equipment is how to realize the holographic acquisition of the state of the equipment based on the minimum acquisition of sensor data and simulation calculation.
Disclosure of Invention
In view of the above, the invention provides a method and a system for calculating a winding temperature field under natural convection heat dissipation of a transformer based on a digital twinning technology, which are used for solving the problems that the existing transformer temperature detection method needs to measure multidimensional parameters, is complex in calculation and is low in monitoring efficiency. To achieve one or a part of or all of the above or other objects, the present invention provides a method for calculating a winding temperature field under natural convection heat dissipation of a transformer, comprising:
performing thermodynamic analysis on oil, winding and iron core parts of the target transformer to obtain oil parameters, winding parameters and iron core parameters corresponding to the target transformer, and establishing a constraint equation according to the thermodynamic analysis result and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition;
acquiring a temperature value of any point in the space where the target transformer is located; and solving to obtain a temperature field corresponding to the target transformer by using the temperature model of the target transformer based on the obtained temperature value.
According to a specific embodiment, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the oil parameters include: oil density, oil constant pressure heat capacity parameters, oil flow rate and oil pressure; the winding parameters include: the coil, the number of turns of the coil, the length of each cake of the coil, and the constant-voltage heat capacity and density of the coil; the iron core parameters include: the geometric structure of the iron core and the geometric structure of the box body.
According to a specific embodiment, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the constraint equation includes: the fluid control equation of oil, and the thermodynamic equation of oil, winding and iron core.
According to a specific embodiment, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the fluid control equation of the oil is as follows:
wherein ρ oil The oil density around the wire,
The flow rate of oil around the wire, p is the pressure of oil around the wire,
Is a unit vector, C p,oil Is the specific heat capacity of the oil around the wire, and T is the temperature of the oil or the winding.
According to a specific embodiment, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, thermodynamic equations of the oil, the winding and the iron core are as follows:
wherein the content of the first and second substances,
Q＝Q I -Q con -Q rad ；
where ρ is oil The oil density around the wire,
The flow rate of oil around the wire, p is the pressure of oil around the wire,
Is a unit vector, C p,oil Is the specific heat capacity of oil around the lead, T is the temperature of the oil or the winding, C p winding Is the specific heat capacity, Q of the wire I For passing heat generated by the current itself through the winding, Q con Heat for convection heat dissipation of the winding, Q rad Represents the heat dissipation of the winding radiation and satisfies the Stefin-Boltz radiation law.
According to a specific embodiment, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the boundary condition includes:
a slip boundary condition and a wall boundary condition corresponding to a fluid control equation of the oil;
and the sum of the heating power of the winding and the core, and the dissipation of the boundary heat flux under thermal equilibrium conditions.
According to a specific implementation manner, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the slip boundary condition is a geometric structure boundary of the iron core; the wall boundary condition is the boundary where the oil intersects with the geometric structures of the winding and the iron core.
According to a specific implementation manner, in the method for calculating the winding temperature field of the transformer under natural convection heat dissipation, solving and obtaining the temperature field corresponding to the target transformer by using the temperature model of the target transformer based on the obtained temperature value includes:
substituting the obtained oil parameters, winding parameters and iron core parameters into the transformer temperature model to obtain a model to be solved;
and solving the model to be solved by taking the temperature value of any point in the space where the target transformer is located as a constraint point to obtain a temperature field corresponding to the target transformer.
According to a specific implementation manner, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the solving the model to be solved by using the temperature value of any point in the space where the target transformer is located as a constraint point includes:
solving the model to be solved by taking the uncertainty of the temperature calculation value of the model to be solved and the actually obtained temperature value smaller than or equal to a threshold value as a target;
calculating the uncertainty by:
wherein, T (r) 0 ,z 0 ) Calculating the value, T, for the temperature 0 (r 0 ,z 0 ) Obtaining the temperature value of any point in the space; the threshold is 1%.
In another aspect of the present invention, a system for calculating a winding temperature field under natural convection heat dissipation of a transformer is provided, including:
the acquisition unit is used for acquiring the temperature value of any point in the space where the target transformer is located and outputting the temperature value to the calculation unit;
the storage unit is used for storing the constructed temperature model of the target transformer; the method comprises the steps of obtaining oil parameters, winding parameters and iron core parameters corresponding to a target transformer by performing thermodynamic analysis on oil, winding and iron core parts of the target transformer, and establishing a constraint equation according to thermodynamic analysis results and the obtained parameters; determining a boundary condition according to the obtained parameters, and obtaining a temperature model of the target transformer based on the constraint equation and the boundary condition;
and the calculating unit is used for solving to obtain a temperature field corresponding to the target transformer by utilizing the temperature model of the target transformer based on the obtained temperature value.
The embodiment of the invention has the following beneficial effects:
according to the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, provided by the embodiment of the invention, the oil parameter, the winding parameter and the iron core parameter corresponding to the target transformer are obtained by carrying out thermodynamic analysis on the oil, the winding and the iron core part of the target transformer, and a constraint equation is established according to the thermodynamic analysis result and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition; therefore, when the temperature value of any point in the space where the target transformer is located is obtained, the temperature field corresponding to the target transformer can be obtained by solving by using the temperature model of the target transformer; the temperature model provided by the embodiment of the invention is constructed based on the thermodynamic principle, the corresponding accuracy is higher, meanwhile, the space distribution of the winding temperature under the natural convection heat dissipation of the transformer can be obtained under the condition that only one point of temperature measurement data can be obtained based on the temperature model provided by the embodiment of the invention, and the efficiency of calculating and monitoring the temperature of the transformer winding is greatly improved.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
Wherein:
FIG. 1 is a schematic flow chart illustrating a method for calculating a winding temperature field under natural convection heat dissipation of a transformer according to an embodiment;
FIG. 2 is a schematic flow chart illustrating a method for calculating a winding temperature field under natural convection heat dissipation of a transformer according to an embodiment;
FIG. 3 is a finite element temperature field calculation model constructed according to certain transformer parameters in one embodiment;
FIG. 4 is a schematic diagram illustrating a spatial distribution of transformer temperature calculated according to a transformer parameter in one embodiment;
FIG. 5 is a block diagram of a system for calculating a winding temperature field under natural convection cooling of a transformer according to an embodiment.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Example 1
Fig. 1 illustrates a method for calculating a winding temperature field under natural convection heat dissipation of a transformer according to an exemplary embodiment of the present invention, including:
performing thermodynamic analysis on oil, windings and an iron core part of a target transformer to obtain oil parameters, winding parameters and iron core parameters corresponding to the target transformer, and establishing a constraint equation according to the thermodynamic analysis result and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition;
acquiring a temperature value of any point in the space where the target transformer is located; and solving to obtain a temperature field corresponding to the target transformer by using the temperature model of the target transformer based on the obtained temperature value.
In the method provided by this embodiment, thermodynamic analysis is performed on oil, winding and iron core parts of a target transformer to obtain oil parameters, winding parameters and iron core parameters corresponding to the target transformer, and a constraint equation is established according to thermodynamic analysis results and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition; the temperature model is constructed based on the thermodynamic principle, the corresponding accuracy is high, meanwhile, the space distribution of the winding temperature under the natural convection heat dissipation of the transformer can be obtained under the condition that only one point of temperature measurement data can be obtained based on the temperature model provided by the embodiment of the invention, and the temperature calculation and monitoring efficiency of the transformer winding is greatly improved.
Example 2
In a possible implementation manner, as shown in fig. 2, the method for calculating a winding temperature field under natural convection heat dissipation of a transformer specifically includes:
s1, obtaining oil material parameters (constant-voltage heat capacity (J/(kg & K)), density (kg/m ^ 3), heat conductivity coefficient (W/(m & K)), dynamic viscosity (Pa & S)), winding parameters and iron core parameters of a target transformer; according to the geometric structure of a transformer iron core, a coil and a box body, the number of coil turns, the length of each coil cake, coil material parameters (constant-voltage heat capacity (J/(kg & K)) of a lead and insulating paper, the density (kg/m & ltlambda > 3) of the lead and the insulating paper, the heat conductivity (W/(m & K)) of the lead and the insulating paper, the thickness (mm) of the insulating paper), oil material parameters (constant-voltage heat capacity (J/(kg & K)), the density (kg/m & ltlambda > 3), the heat conductivity (W/(m & K)), dynamic viscosity (Pa & s)), and rated current, establishing a constraint equation and a constraint condition, and obtaining a finite element temperature field calculation model (namely the temperature model of the target transformer) of the two-dimensional flow-heat strong coupling of the transformer;
s2, measuring the temperature of any point in the space near the transformer winding in practice, and taking the temperature value of the point as a constraint point in a two-dimensional flow heat strong coupling temperature field calculation model;
s3, solving the model to be solved according to the condition that the uncertainty of the temperature calculation value of the model to be solved and the actually acquired temperature value is less than or equal to a threshold value;
calculating the uncertainty by:
wherein, T (r) 0 ,z 0 ) Calculating the value, T, for the temperature 0 (r 0 ,z 0 ) Obtaining the temperature value of any point in the space; the threshold is 1%.
In a possible implementation manner, in the method for calculating a winding temperature field under natural convection heat dissipation of a transformer, the oil parameter in S1 includes: oil material parameters (constant pressure heat capacity (J/(kg. K)), density (kg/m ^ 3), thermal conductivity (W/(m. K)), dynamic viscosity (Pa. S)), oil flow rate and oil pressure; the winding parameters include: the coil, the number of turns of the coil, the length of each cake of the coil, the coil material parameters (constant voltage heat capacity (J/(kg.K)) of the wire and the insulating paper, the density (kg/m ^ 3) of the wire and the insulating paper, the thermal conductivity (W/(m.K)) of the wire and the insulating paper and the thickness (mm)) of the insulating paper); the iron core parameters include: the geometric structure of the iron core and the geometric structure of the box body.
In a possible implementation manner, in the method for calculating the winding temperature field under natural convection heat dissipation of the transformer, the constraint equation established in S1 includes: the fluid control equation of oil, and the thermodynamic equation of oil, winding and iron core.
In a possible implementation manner, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the fluid control equation of the oil is as follows:
where ρ is oil The oil density around the wire,
The flow rate of oil around the wire, p is the pressure of oil around the wire,
Is a unit vector, C p,oil Is the specific heat capacity of the oil around the wire, and T is the temperature of the oil or the winding.
In a possible implementation manner, in the method for calculating the winding temperature field under natural convection heat dissipation of the transformer, thermodynamic equations of the oil, the winding, and the iron core portion are as follows:
wherein
Q＝Q I -Q con -Q rad (3)
Where ρ is oil The oil density around the wire,
The flow rate of oil around the wire, the pressure of oil around the wire p,
Unit vector, C p,oil Specific heat capacity of oil around the lead, temperature of T oil or winding,
Specific heat capacity, Q of the wire I For passing heat generated by the current itself through the winding, Q con Heat, Q, for convective heat dissipation of the winding rad Represents the heat dissipation of the winding radiation and satisfies the Stefin-Boltz radiation law.
In a possible implementation manner, in the method for calculating the winding temperature field under natural convection heat dissipation of the transformer, the boundary condition includes: a slip boundary condition and a wall boundary condition corresponding to a fluid control equation of the oil; and the sum of the heating power of the winding and the core, and the dissipation of the boundary heat flux under thermal equilibrium conditions. The slip boundary condition is a geometric structure boundary of the iron core; the wall boundary condition is the boundary where the oil intersects with the geometric structures of the winding and the iron core.
Specifically, by adopting the method, a finite element temperature field calculation model of the transformer two-dimensional flow thermal strong coupling established according to the geometric structure, the number of turns of the coil, the coil material parameter, the oil material parameter and the rated current of a certain transformer iron core, the coil and the box body is shown in fig. 3, wherein ABCF is an iron core part:
the boundary conditions corresponding to the constraint equation are specifically as follows:
the fluid control equations CD, DE, EF for the oil section are slip boundary conditions, and the other high, medium, low and core-to-oil boundaries are no-slip wall boundary conditions.
Thermodynamic equation of oil, winding, core part:
the sum of the heating power of the high voltage, the medium voltage, the low voltage and the iron core is P 0 The dissipation of the heat flux through the boundaries AB, BC, CD, DE, EF and FA at thermal equilibrium is-P 0 。
In a possible implementation manner, in the method for calculating a winding temperature field under natural convection heat dissipation of a transformer, the step S2 specifically includes:
measuring the temperature T of any point in the space near the transformer winding 0 (r 0 ,z 0 ) The point is positioned at any position of oil in a calculation model area, even the condition on the boundaries of CD, DE and EF or the boundaries of high voltage, medium voltage, low voltage and iron core and oil can be met, and the temperature value T of the point is calculated by a two-dimensional flow heat strong coupling temperature field calculation model 0 (r 0 ,z 0 ) As a constraint point.
In a possible implementation manner, in the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer, the step S3 specifically includes:
when the two-dimensional temperature field calculation model with strong heat coupling of the flow is used for solving the equations (1), (2) and (3), the value T (r) is calculated 0 ,z 0 ) Approximately equal to the actual measured value T 0 (r 0 ,z 0 ) (the uncertainty of the two is less than or equal to the threshold, wherein the threshold can be set to 1% -0.01% according to the requirement, and the threshold is equal to 1% in this embodiment):
and the calculated T (r, z) is the spatial distribution of the winding temperature under the natural convection heat dissipation of the transformer.
Specifically, by using the method for calculating the winding temperature field under the natural convection heat dissipation of the transformer provided by the embodiment of the invention, a finite element temperature field calculation model of the two-dimensional heat strong coupling of the transformer, which is established according to the geometric structures of an iron core, a coil and a box body of the transformer, the number of turns of the coil, the coil material parameters, the oil material parameters and the rated current, is shown in fig. 3, and the spatial distribution of the winding temperature under the natural convection heat dissipation of the corresponding transformer is shown in fig. 4.
Example 3
In another aspect of the present invention, a system for calculating a winding temperature field under natural convection heat dissipation of a transformer is provided, including:
the acquisition unit is used for acquiring the temperature value of any point in the space where the target transformer is located and outputting the temperature value to the calculation unit;
the storage unit is used for storing the constructed temperature model of the target transformer; the method comprises the steps of obtaining oil parameters, winding parameters and iron core parameters corresponding to a target transformer by performing thermodynamic analysis on oil, winding and iron core parts of the target transformer, and establishing a constraint equation according to thermodynamic analysis results and the obtained parameters; determining a boundary condition according to the obtained parameters, and obtaining a temperature model of the target transformer based on the constraint equation and the boundary condition;
and the calculating unit is used for solving and obtaining the temperature field corresponding to the target transformer by utilizing the temperature model of the target transformer based on the obtained temperature value.
Specifically, as shown in fig. 5, the storage unit and the computing unit may be provided by an electronic device, where the electronic device may be one of a terminal, a tablet computer, a computer, and the like, and the embodiment is not limited thereto. The acquiring unit may be a temperature sensor disposed at any point in the space.
Specifically, the electronic device includes a processor, a network interface, and a memory, where the processor, the network interface, and the memory are connected to each other, where the memory is used to store a computer program, and the computer program includes program instructions, and the processor is configured to call the program instructions to perform the method for calculating the winding temperature field under natural convection heat dissipation of the transformer.
In another aspect of the present invention, a computer storage medium is further provided, in which program instructions are stored, and the program instructions, when executed by at least one processor, are used to implement the above method for calculating the winding temperature field under natural convection heat dissipation of a transformer.
In an embodiment of the invention, the processor may be an integrated circuit chip having signal processing capabilities. The Processor may be a general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component.
The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The processor reads the information in the storage medium and completes the steps of the method in combination with the hardware.
The storage medium may be a memory, for example, which may be volatile memory or nonvolatile memory, or which may include both volatile and nonvolatile memory.
The nonvolatile Memory may be a Read-Only Memory (ROM), a Programmable ROM (PROM), an Erasable PROM (EPROM), an Electrically Erasable PROM (EEPROM), or a flash Memory.
The volatile Memory may be a Random Access Memory (RAM) which serves as an external cache. By way of example, and not limitation, many forms of RAM are available, such as Static Random Access Memory (SRAM), dynamic RAM (DRAM), synchronous DRAM (SDRAM), double Data Rate SDRAM (DDRSDRAM), enhanced SDRAM (ESDRAM), SLDRAM (SLDRAM), and Direct Rambus RAM (DRRAM).
The storage media described in connection with the embodiments of the invention are intended to comprise, without being limited to, these and any other suitable types of memory.
It should be understood that the disclosed system may be implemented in other ways. For example, the division of the modules into only one logical function may be implemented in another way, for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not implemented. In addition, the communication connection between the modules may be an indirect coupling or communication connection of the server or the unit through some interfaces, and may be an electrical or other form.
In addition, functional modules in the embodiments of the present invention may be integrated into one processing unit, or each module may exist alone physically, or two or more modules are integrated into one processing unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a removable hard disk, a magnetic or optical disk, and other various media capable of storing program codes.
The above disclosure is only for the purpose of illustrating the preferred embodiments of the present invention, and it is therefore to be understood that the invention is not limited by the scope of the appended claims.
Claims (10)
1. A method for calculating a winding temperature field under natural convection heat dissipation of a transformer is characterized by comprising the following steps: the method comprises the following steps:
performing thermodynamic analysis on oil, windings and an iron core part of a target transformer to obtain oil parameters, winding parameters and iron core parameters corresponding to the target transformer, and establishing a constraint equation according to the thermodynamic analysis result and the obtained parameters; determining a boundary condition according to the obtained parameters, and constructing a temperature model of the target transformer based on the constraint equation and the boundary condition;
acquiring a temperature value of any point in the space where the target transformer is located; and solving to obtain a temperature field corresponding to the target transformer by using the temperature model of the target transformer based on the obtained temperature value.
2. The method for calculating the winding temperature field under the natural convection heat dissipation of the transformer according to claim 1, wherein: the oil parameters include: oil density, oil constant pressure heat capacity parameters, oil flow rate and oil pressure; the winding parameters include: the coil, the number of turns of the coil, the length of each cake of the coil, the constant-voltage heat capacity and the density of the coil; the iron core parameters include: the geometric structure of the iron core and the geometric structure of the box body.
3. The method for calculating the winding temperature field under the natural convection heat dissipation of the transformer according to claim 1, wherein: the constraint equations include: the fluid control equation of oil, and the thermodynamic equation of oil, winding and iron core.
4. The method of calculating the winding temperature field under natural convection heat dissipation of a transformer of claim 2, wherein: the fluid handling equation for the oil is:
where ρ is oil Is the density of the oil around the wire,
is the flow rate of oil around the wire, p is the pressure of oil around the wire,
Is a vector of the unit,
t is the specific heat capacity of the oil surrounding the wire, and T is the temperature of the oil or the winding.
5. The method of calculating the winding temperature field under natural convection heat dissipation of a transformer of claim 2, wherein: the thermodynamic equation of the oil, the winding and the iron core is as follows:
wherein the content of the first and second substances,
Q＝Q I -Q con -Q rad ；
where ρ is oil The oil density around the wire,
The flow rate of oil around the wire, p is the pressure of oil around the wire,
Is a unit vector,
The specific heat capacity of oil around the lead, T the temperature of the oil or the winding,
Is the specific heat capacity, Q of the wire I For passing heat generated by the current itself through the winding, Q con For heat dissipated by convection of the winding, Q rad To represent the heat dissipation of the winding radiation and to satisfy the stefin-boltz radiation law.
6. The method of calculating the winding temperature field under natural convection heat dissipation of a transformer of claim 2, wherein: the boundary conditions include:
a slip boundary condition and a wall boundary condition corresponding to a fluid control equation of the oil;
and the sum of the heating power of the winding and the core, and the dissipation of the boundary heat flux under thermal equilibrium conditions.
7. The method of calculating the winding temperature field under natural convection heat dissipation of a transformer of claim 6, wherein: the slippage boundary condition is a geometric structure boundary of the iron core; the wall boundary condition is the boundary where the oil intersects with the geometric structures of the winding and the iron core.
8. The method for calculating the winding temperature field under the natural convection heat dissipation of the transformer according to any one of claims 1 to 7, wherein: the obtaining of the temperature field corresponding to the target transformer by solving the temperature model of the target transformer based on the obtained temperature value includes:
substituting the obtained oil parameters, winding parameters and iron core parameters into the transformer temperature model to obtain a model to be solved;
and solving the model to be solved by taking the temperature value of any point in the space where the target transformer is located as a constraint point to obtain a temperature field corresponding to the target transformer.
9. The method of calculating the winding temperature field under natural convection heat dissipation of a transformer of claim 8, wherein: solving the model to be solved by taking the temperature value of any point in the space where the target transformer is located as a constraint point, wherein the solving comprises the following steps:
solving the model to be solved by taking the uncertainty of the temperature calculation value of the model to be solved and the actually obtained temperature value smaller than or equal to a threshold value as a target;
calculating the uncertainty by:
wherein, T (r) 0 ,z 0 ) Calculating the value, T, for the temperature 0 (r 0 ,z 0 ) The temperature value of any point in the acquired space is obtained.
10. The utility model provides a system for winding temperature field under the heat dissipation of calculating transformer natural convection which characterized in that: the method comprises the following steps:
the acquisition unit is used for acquiring the temperature value of any point in the space where the target transformer is located and outputting the temperature value to the calculation unit;
the storage unit is used for storing the constructed temperature model of the target transformer; the method comprises the steps of obtaining oil parameters, winding parameters and iron core parameters corresponding to a target transformer by performing thermodynamic analysis on oil, winding and iron core parts of the target transformer, and establishing a constraint equation according to thermodynamic analysis results and the obtained parameters; determining a boundary condition according to the obtained parameters, and obtaining a temperature model of the target transformer based on the constraint equation and the boundary condition;
and the calculating unit is used for solving and obtaining the temperature field corresponding to the target transformer by utilizing the temperature model of the target transformer based on the obtained temperature value.
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211345515.6A CN115659659A (en) | 2022-10-31 | 2022-10-31 | Method and system for calculating winding temperature field under natural convection heat dissipation of transformer |
| PCT/CN2022/134842 WO2024092924A1 (en) | 2022-10-31 | 2022-11-29 | Method and system for calculating winding temperature field under natural convection heat dissipation of transformer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211345515.6A CN115659659A (en) | 2022-10-31 | 2022-10-31 | Method and system for calculating winding temperature field under natural convection heat dissipation of transformer |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN115659659A true CN115659659A (en) | 2023-01-31 |
Family
ID=84993324
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211345515.6A Pending CN115659659A (en) | 2022-10-31 | 2022-10-31 | Method and system for calculating winding temperature field under natural convection heat dissipation of transformer |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN115659659A (en) |
| WO (1) | WO2024092924A1 (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117910393A (en) * | 2024-03-20 | 2024-04-19 | 南方电网科学研究院有限责任公司 | Method and device for calculating internal pressure of transformer bushing |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| TWI256000B (en) * | 2004-12-31 | 2006-06-01 | Tatung Co | A method for analyzing temperature and fluid field on power transformers |
| CN104834781B (en) * | 2015-05-08 | 2018-04-27 | 国网河南省电力公司电力科学研究院 | A kind of Analysis of Transient Thermal Field emulation mode during commutation failure multiple based on smoothing reactor |
| CN113901686A (en) * | 2021-10-08 | 2022-01-07 | 广东电网有限责任公司广州供电局 | Transformer thermal circuit model construction method and transformer temperature field calculation method |
| CN114970259B (en) * | 2022-05-20 | 2024-05-31 | 中国南方电网有限责任公司超高压输电公司广州局 | Converter transformer modeling method, apparatus, device, storage medium and program product |
-
2022
- 2022-10-31 CN CN202211345515.6A patent/CN115659659A/en active Pending
- 2022-11-29 WO PCT/CN2022/134842 patent/WO2024092924A1/en unknown
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117910393A (en) * | 2024-03-20 | 2024-04-19 | 南方电网科学研究院有限责任公司 | Method and device for calculating internal pressure of transformer bushing |
| CN117910393B (en) * | 2024-03-20 | 2024-05-31 | 南方电网科学研究院有限责任公司 | Method and device for calculating internal pressure of transformer bushing |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2024092924A1 (en) | 2024-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115659659A (en) | Method and system for calculating winding temperature field under natural convection heat dissipation of transformer | |
| CN107192738A (en) | A kind of detection method and device of distribution transformer winding material | |
| JP7048284B2 (en) | Transformer diagnostic system, transformer diagnostic method, and transformer | |
| CN104697664B (en) | The method of adjustment coiling hot point of transformer on-line temperature monitoring | |
| CN111623884B (en) | Transformer hot spot temperature identification method and system based on improved heat network model | |
| US9419430B1 (en) | System for monitoring and modeling operation of a transformer | |
| Diaz‐Aguiló et al. | Adaptive soil model for real‐time thermal rating of underground power cables | |
| Wang et al. | Top‐oil temperature modelling by calibrating oil time constant for an oil natural air natural distribution transformer | |
| Li et al. | Test and analysis on extended temperature rise of 110 kV transformer based on distributed temperature sensing | |
| Tao et al. | Three‐dimensional thermal modelling of transformers in transformer room for spatial and temporal failure analysis | |
| Soltanbayev et al. | Automated dry-type transformer aging evaluation: A simulation study | |
| CN111177907A (en) | Method and device for automatically evaluating service life of reactor | |
| CN112710413B (en) | Temperature sensor checking method and device, storage medium and electronic equipment | |
| Correa et al. | Power transformer temperature–moisture dynamics modeling: an experimental validation | |
| Askari et al. | Analysis of thermal models to determine the loss of life of mineral oil immersed transformers | |
| CN114595528A (en) | Visualization system and method of digital twin-based transformer outgoing line device | |
| CN111859802A (en) | Method, system and equipment for calculating steady-state temperature of overhead conductor | |
| Melka et al. | Effective cooling of a distribution transformer using biodegradable oils at different climate conditions | |
| CN113237506B (en) | Early warning method and system for high-voltage bushing tap device and storage medium | |
| Zhang et al. | Improving IEC thermal model for oil natural air natural transformers using optimised parameters based on dynamic simulation | |
| CN117350093A (en) | Temperature-based fault early warning method, system, device and medium | |
| He et al. | Investigating temperature rise dynamics at hot-spots within dry-Type transformer windings: A comparative analysis across varied loading rates and an extrapolative computational model | |
| CN116184104B (en) | Fault detection method, device and equipment of oil immersed transformer and storage medium | |
| Simonson | Transformer ratings and transformer life | |
| CN118211435A (en) | Equivalent modeling method and system for sleeve coiled pipe and cable penetrating air interlayer |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |