CN107609732B - Wind power consumption potential determination method and system - Google Patents

Wind power consumption potential determination method and system Download PDF

Info

Publication number
CN107609732B
CN107609732B CN201710638461.5A CN201710638461A CN107609732B CN 107609732 B CN107609732 B CN 107609732B CN 201710638461 A CN201710638461 A CN 201710638461A CN 107609732 B CN107609732 B CN 107609732B
Authority
CN
China
Prior art keywords
temperature
model
wind power
time
power
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.)
Active
Application number
CN201710638461.5A
Other languages
Chinese (zh)
Other versions
CN107609732A (en
Inventor
何桂雄
覃剑
唐艳梅
蒋利民
钟鸣
梁琛
王维洲
刘铠诚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
State Grid Corp of China SGCC
China Electric Power Research Institute Co Ltd CEPRI
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 STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE, State Grid Corp of China SGCC, China Electric Power Research Institute Co Ltd CEPRI filed Critical STATE GRID GASU ELECTRIC POWER RESEARCH INSTITUTE
Priority to CN201710638461.5A priority Critical patent/CN107609732B/en
Publication of CN107609732A publication Critical patent/CN107609732A/en
Application granted granted Critical
Publication of CN107609732B publication Critical patent/CN107609732B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P20/00Technologies relating to chemical industry
    • Y02P20/10Process efficiency
    • Y02P20/133Renewable energy sources, e.g. sunlight

Landscapes

  • Ceramic Products (AREA)
  • Carbon And Carbon Compounds (AREA)

Abstract

A method and a system for determining wind power consumption potential comprise the following steps: determining a temperature rise stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power; calculating the wind power consumption potential of the silicon carbide enterprise according to the wind power consumption potential model; the wind power consumption potential model comprises a resistance furnace power time model and a power load model of an auxiliary material system. According to the technical scheme provided by the invention, a wind power consumption potential model is established, the possibility of wind power consumption under normal production is provided, and the electric quantity consumed by the whole silicon carbide smelting enterprise is maximized through calculation, so that the enterprise can reasonably and efficiently use wind power, and the cost is reduced.

Description

Wind power consumption potential determination method and system
Technical Field
The invention relates to the fields of automation and power electronics, in particular to a method and a system for determining wind power consumption potential.
Background
Wind energy is regarded as renewable energy, and people pay more and more attention to and use the renewable energy in recent years, but wind power generation is abandoned because of the characteristics of large fluctuation, unbalanced output and the like, and therefore the problem of wind abandonment is increasingly serious in recent years. Silicon carbide smelting enterprises as high-energy-carrying enterprises have great significance in participating in wind power consumption at load ends.
The wind abandonment problem is essentially the problem of wind power consumption capability of a power system, the wind abandonment problem needs the joint efforts of all links including power generation, power transmission and distribution and power utilization and all interest relevant parties, the wind power large base is planned in a system, and the responsibility main body is implemented.
In the face of the problem of abandoned wind, a lot of researches are made on the problem of wind power consumption recently, and the problems can be classified into three types of source network loads in summary. The source, namely, the problem of wind abandon is solved from the power generation end; the grid, namely the wind power integration is realized in a large range, and the problem is solved from a transmission and distribution end; and loading, namely a wind power using end.
Therefore, a method and a system for constructing a wind power consumption potential model are found, and the method and the system are the problems that technical personnel in the field need to solve urgently in order to carry out wind power consumption under normal production and efficiently utilize wind power.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a method and a system for determining wind power digestion potential.
A wind power consumption potential determination method, the method comprising:
determining a temperature rise stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power;
calculating the wind power consumption potential of the silicon carbide enterprise according to the wind power consumption potential model;
the wind power consumption potential model comprises a resistance furnace power time model and a power load model of an auxiliary material system.
Preferably, the power time model of the resistance furnace comprises:
constructing a silicon carbide smelting two-dimensional heat transfer model;
carrying out differential processing on the two-dimensional heat transfer model, and combining boundary conditions to obtain a temperature gradient model on the premise of meeting stability conditions;
and obtaining a temperature gradient function from the temperature gradient model, and obtaining a power time model of the resistance furnace based on a time-temperature relation.
Preferably, the silicon carbide smelting two-dimensional heat transfer model is shown as the following formula:
Figure BDA0001364849270000021
wherein λ is the thermal conductivity of the material, ρ is the density of the material, C p Is the specific heat capacity of the material, q v To generate heat of reduction reaction of silicon carbide.
Preferably, the resistance furnace temperature gradient model is constructed according to a silicon carbide smelting two-dimensional heat transfer model, and comprises the following steps: and establishing a coordinate axis by taking the outer surface of the right side of the furnace core as a boundary, dividing the coordinate axis into rectangular grids in the X-axis direction and the Y-axis direction by the distances of delta X and delta Y respectively, performing differential processing on the two-dimensional heat transfer model, and combining boundary conditions to obtain a temperature gradient model on the premise of meeting stability conditions.
Preferably, the stability conditions are represented by the following formula:
Figure BDA0001364849270000022
preferably, the boundary condition is represented by the following formula:
T i,j 0 =30℃ (3)
T i,j 0 =T w n (4)
Figure BDA0001364849270000023
Figure BDA0001364849270000024
wherein P (t) is work at time t; l, h and b are respectively the length, the height and the width of the furnace core of the resistance furnace; t is w n The temperature of the outer surface of the furnace core at any moment can be measured; t is i,j 0 : when not electrified initially, the temperature of the whole resistance furnace; λ: the coefficient of thermal conductivity of the material; ρ: the density of the material; c p : material specific heat capacity; Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; q. q.s w : the density of the heat flux passing over the outer surface of the core; Δ t: the time period is divided in time at intervals Δ t.
Preferably, the temperature gradient model is represented by the following formula:
Figure BDA0001364849270000031
in the formula, λ: the coefficient of thermal conductivity of the material; Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; p (t): a function of power over time; t is 1 ......T r : and (3) temperature.
Preferably, the temperature raising stage of the resistance furnace divided by the temperature gradient model includes:
the first stage is as follows: heating the material, and when the highest temperature of the material is 1450 ℃, generating no silicon carbide;
and a second stage: heating the material, and only generating beta-SiC when the maximum temperature of the material is 1800 ℃;
and a third stage: heating the material, wherein 4H-alpha-SiC is generated when the maximum temperature of the material is 2200 ℃;
a fourth stage: heating the material, and generating 6H-alpha-SiC when the maximum temperature of the material is 2600 ℃;
the fifth stage: the material is heated and when the maximum temperature of the material exceeds 2600 ℃, the SiC decomposes into carbon and silicon.
Preferably, the obtaining a temperature gradient function from the temperature gradient model and obtaining a power time model of the resistance furnace based on a time-temperature relationship includes:
obtaining a temperature gradient function according to the temperature gradient model, wherein the temperature gradient function is as follows: t (x,0, T), wherein: x is a horizontal axis coordinate; t is the time;
arbitrarily giving a power value, and obtaining a time-temperature relationship through the temperature gradient function, wherein the time-temperature relationship expression is as follows:
Figure BDA0001364849270000032
in the formula: t (t) is: a function of temperature time; p is power, p k Is taken as a power value, x is a horizontal axis coordinate, x n Taking values for the abscissa;
and (4) simultaneously solving to obtain the power-time relation of the resistance furnace under the condition of any power value, and further obtaining a power-time model of the resistance furnace.
Preferably, any given one of the powers is set by a minimum regulation time of the transformer and a resistance furnace reaction;
the minimum regulating time of the transformer is the minimum unit of power change;
the resistance furnace reaction includes a lance pressure and a thermal reduction reaction rate.
Preferably, the constructing of the power load model of the auxiliary material system includes: grouping the auxiliary material systems according to whether the auxiliary material systems are transferable loads, and constructing a power load model of the auxiliary material systems for the transferable loads by using the following formula:
W auxiliary device =W P +W S +W H +W L
(8)
In the formula, W P : the amount of electricity consumed by the group with the crusher as the primary equipment; w is a group of S : the amount of electricity consumed by the group using the screening machine as the main equipment; w H : the electric quantity consumed by the group taking the mixer as main equipment; w L : there is no power consumed by the group of primary devices.
Preferably, the grouping of the auxiliary material systems includes: and (4) grouping the belt conveyor, the bucket elevator, the crusher, the screening machine and the mixer.
Preferably, the step of determining the temperature rise stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power includes: and determining the temperature rising stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power and the temperature rising stage divided by combining the temperature of the resistance furnace when the wind power arrives and a temperature gradient model.
A wind power consumption potential determination system, the system comprising: the system comprises a first construction module, a second construction module, a third construction module and a calculation module;
the first construction module is used for constructing a wind power consumption potential model;
the second construction module is used for constructing a power time model of the resistance furnace;
and the third construction module is used for constructing a power load model of the auxiliary material system.
And the calculation module is used for calculating the wind power consumption potential of the silicon carbide enterprise by combining a wind power consumption potential model according to the arrival time and duration of wind power and the temperature of the resistance furnace when the wind power arrives.
Compared with the prior art, the invention has the beneficial effects that:
according to the technical scheme provided by the invention, a resistance furnace power time model is constructed by combining a temperature gradient model constructed by a silicon carbide smelting two-dimensional heat transfer model with a power constraint condition, when wind power comes, different power combination modes are determined through the resistance furnace power time model, and the wind power absorption potential of a silicon carbide enterprise is calculated, so that the wind power utilization can be realized to the maximum extent, and the cost is saved;
according to the technical scheme provided by the invention, the wind power consumption potential model is established, and the wind power consumption potential of the silicon carbide enterprise is calculated within the wind power incoming time t informed by the dispatching station, so that the enterprise production can be guided, the electric quantity consumed by the whole silicon carbide smelting enterprise is maximized, the wind power is reasonably and efficiently used, and the cost is reduced;
according to the technical scheme provided by the invention, the power load model of the auxiliary material system is established according to the power consumption, the production time and the rated power of the batching system, so that the possibility of wind power consumption under normal production is provided.
Drawings
FIG. 1 is a flow chart of the construction of a wind power consumption potential model according to the present invention;
FIG. 2 is a process flow diagram of the silicon carbide smelting batching system of the present invention;
FIG. 3 is a schematic cross-sectional view of a silicon carbide of the present invention;
FIG. 4 is a schematic view of the present invention showing the establishment of the furnace coordinate axes;
FIG. 5 is a schematic diagram of a discrete node of the present invention;
wherein, 1-heat preservation material, 2-second material, 3-first material and 4-middle point of right boundary of the furnace core.
Detailed Description
For a better understanding of the present invention, reference is made to the following description taken in conjunction with the accompanying drawings and examples.
As shown in fig. 1, the method for determining wind power consumption potential provided by the present invention includes:
determining the working stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power;
calculating the wind power consumption potential of the silicon carbide enterprise according to the wind power consumption potential model;
the wind power consumption potential model comprises a resistance furnace power time model and a power load model of the auxiliary material system.
As shown in fig. 2, the main equipment categories are screen, crusher, blender and resistance furnace. According to the process flow, the method is divided into 10 groups.
Wherein the transferable load is defined as: not only material conveying equipment (belts, bucket elevators and the like) but also material storage equipment such as bins exist among the production equipment, so that some equipment can independently operate according to the capacity of the bins and the normal production process is not influenced. These consumers, which can operate independently and can be turned on and off at any time, are referred to as transferable loads.
The non-transferable load is defined as: the operation of the plant requires a certain time and material coupling relationship to be met, and therefore the plant without transferable additional conditions (e.g. a silo) is called a non-transferable load.
When the bin meets the size, 10 packets can be considered as transferable loads.
The power consumption of the whole silicon carbide enterprise can be approximately considered as the sum of the power consumption of the auxiliary material system and the power consumption of the resistance furnace.
∑W General assembly =∑W Auxiliary device +∑W Furnace with a heat exchanger (5)
The ten auxiliary material systems can be independently used as transferable loads to be started and stopped at any time.
The power consumption of the auxiliary material system is expressed as follows:
W auxiliary device =W P +W S +W H +W L (6)
In the formula, W P : the amount of electricity consumed by the group with the crusher as the primary equipment; w S : the amount of electricity consumed by the group using the screening machine as the main equipment; w H : the electric quantity consumed by the group taking the mixer as main equipment; w L : there is no power consumed by the group of primary devices.
When the system runs stably, the inlet and outlet flow of the materials reaches a balance,
therefore, we can get the corresponding power load relationship:
W auxiliary device =W P +W S +W H +W L =P 1 (M 1 )t+P 2 (M 2 )t+P 3 (M 3 )t (7)
Assuming that t is the given wind power coming time, when the wind power comes temporarily, the maximum electric quantity W which can be consumed within the time t can be determined only by determining the load quantities which need to be completed and correspond to the ten groups of equipment under the condition of meeting the process and normal production conditions.
The power time model of the resistance furnace comprises: constructing a silicon carbide smelting two-dimensional heat transfer model;
carrying out differential processing on the two-dimensional heat transfer model, and combining boundary conditions to obtain a temperature gradient model on the premise of meeting stability conditions;
and obtaining a temperature gradient function through the temperature gradient model, and obtaining a power time model of the resistance furnace based on the time-temperature relation.
As shown in FIG. 3, the silicon carbide smelting furnace can be approximately seen as a rectangular parallelepiped trough-shaped smelting device, the outermost of the cross-sectional view is the heat insulating material 1, the upper position inside the heat insulating material is two materials 2, and the lower position is one material 3, so that any cross section can be approximately seen as heat insulation along the longitudinal direction in the longitudinal length direction. Thus, a two-dimensional heat transfer model of a silicon carbide smelting furnace can be established:
Figure BDA0001364849270000061
wherein, lambda is the material heat conductivity coefficient, rho is the material density, C p Is the specific heat capacity of the material, q v Is an internal heat source (i.e. the heat of the reduction reaction for generating silicon carbide). Because the reaction heat can be approximately ignored compared with the internal heat source, the silicon carbide heat transfer model can be approximated to a second-order unsteady heat conduction model without the internal heat source.
The stability conditions are shown below:
Figure BDA0001364849270000071
in the formula, Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; Δ t: dividing the time period by the interval Δ t; c p : material specific heat capacity; λ: the coefficient of thermal conductivity of the material; ρ: the density of the material.
The boundary conditions are shown below:
T i,j 0 =30℃ (3)
T i,j 0 =T w n (4)
Figure BDA0001364849270000072
Figure BDA0001364849270000073
wherein P (t) is the power at time t; l, h and b are respectively the length, the height and the width of the furnace core of the resistance furnace; t is w n The temperature of the outer surface of the furnace core at any moment can be measured; t is i,j 0 : when not electrified initially, the temperature of the whole resistance furnace; λ: the coefficient of thermal conductivity of the material; q. q.s w : the density of the heat flux passing over the outer surface of the core.
Referring to FIG. 4, the cross-section is defined by the outer surface of the right side of the furnace core, and coordinate axes are established as shown in FIG. 4. And carrying out differential processing on the heat transfer model to obtain the temperature distribution condition of the whole furnace body at each moment, wherein the temperature distribution condition is represented as a temperature-time model by the following formula:
Figure BDA0001364849270000074
t (x,0, T) is a temperature gradient function, wherein: x is the abscissa: t is the time.
At any given power value, the time-temperature relationship of any point on the x-axis can be obtained through the temperature gradient function
Figure BDA0001364849270000081
In the formula: t (t) is: a function of temperature over time; p is power, p k Is taken as a power value, x is a horizontal axis coordinate, x n Taking values for the abscissa; simultaneous solution can obtain the power-time relation of the resistance furnace under the condition of power value, and finally obtain a data set of power,and obtaining the power time model of the resistance furnace.
Any given one of the powers is set by the minimum adjustment time of the transformer and the resistance furnace reaction;
the minimum regulating time of the transformer is the minimum unit of power change;
the resistance furnace reaction includes the lance pressure and the thermal reduction reaction rate.
Therefore, when wind power comes, the maximum electric quantity which can be consumed within the time t can be determined only by determining different power combination modes under the condition of meeting power constraint and normal production conditions.
And forming a wind power consumption potential model according to the resistance furnace power time model when the maximum power is combined and the power load model of the auxiliary material system when the wind power consumption is maximum.
According to a resistance furnace power time model in maximum power combination: determining the temperature rise stages of the resistance furnace according to the arrival time and duration of the wind power, wherein the resistance furnace power time model is in a power combination mode when each temperature rise stage reaches the maximum power;
the temperature rising stage of the resistance furnace comprises the following steps:
the first stage is as follows: when in use
Figure BDA0001364849270000082
The temperature of the materials from the furnace core to the outside is 1450-30 ℃, only the materials are heated, and no silicon carbide is generated;
and a second stage: when in use
Figure BDA0001364849270000083
The temperature of the materials from the furnace core to the outside is 1800-1450 ℃, 1450-30 ℃ in two temperature sections, only beta-SiC is generated, SiO2+3C → SiC +2 CO;
and a third stage: when in use
Figure BDA0001364849270000084
The temperature of materials from the furnace core to the outside is 2200-1800 ℃, 1800-1450 ℃ and 1450-30 ℃, and at the time, beta-SiC → 4H-alpha-SiC is obtained during the temperature of 2200-18000 ℃;
a fourth stage: when in use
Figure BDA0001364849270000085
The temperature of the material outside the furnace core is 2600-2200 ℃, 2200-1800 ℃, 1800-1450 ℃, 1450-30 ℃, and 4H-alpha-SiC → 6H-alpha-SiC is carried out during the temperature is 2600-2200 ℃;
the fifth stage: when in use
Figure BDA0001364849270000086
The temperature of the material from the furnace core to the outside is>2600 deg.C, 2600 deg.C-2200 deg.C, 2200 deg.C-1800 deg.C, 1800 deg.C-1450 deg.C, 1450 deg.C-30 deg.C, 2600 deg.C of silicon carbide and decomposition of SiC into carbon and silicon.
Correspondingly, the invention also provides a wind power consumption potential determining system, which comprises: the system comprises a first construction module, a second construction module, a third construction module and a calculation module;
the first construction module is used for constructing a wind power consumption potential model;
the second construction module is used for constructing a power time model of the resistance furnace;
and the third construction module is used for constructing a power load model of the auxiliary material system.
And the calculation module is used for calculating the wind power consumption potential of the silicon carbide enterprise by combining a wind power consumption potential model according to the arrival time and duration of wind power and the temperature of the resistance furnace when the wind power arrives.
The wind power consumption potential determining system provided by the invention substantially corresponds to the wind power consumption potential determining method, the same technical effect can be achieved, and specific implementation processes can refer to the contents of the embodiment of the wind power consumption potential determining method, and are not repeated herein.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The present invention is not limited to the above embodiments, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention are included in the scope of the claims of the present invention which are filed as the application.

Claims (4)

1. A wind power consumption potential determination method is characterized by comprising the following steps:
determining the temperature rise stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power;
calculating the wind power consumption potential of the silicon carbide enterprise according to the wind power consumption potential model;
the wind power consumption potential model comprises a resistance furnace power time model and a power load model of an auxiliary material system;
the power time model of the resistance furnace comprises:
constructing a silicon carbide smelting two-dimensional heat transfer model;
carrying out differential processing on the two-dimensional heat transfer model, and combining boundary conditions to obtain a temperature gradient model on the premise of meeting stability conditions;
obtaining a temperature gradient function from the temperature gradient model, and obtaining a power time model of the resistance furnace based on a time-temperature relation;
the silicon carbide smelting two-dimensional heat transfer model is shown as the following formula:
Figure FDA0003658127040000011
wherein λ is the thermal conductivity of the material, ρ is the density of the material, C p Is the specific heat capacity of the material, q v To generate the heat of reduction reaction of silicon carbide;
the stability conditions are shown below:
Figure FDA0003658127040000012
in the formula, Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; Δ t: dividing the time period by the interval Δ t; c p : material specific heat capacity; λ: the coefficient of thermal conductivity of the material; ρ: a material density;
the boundary condition is shown as follows:
T i,j 0 =30℃ (3)
T i,j 0 =T w n (4)
Figure FDA0003658127040000013
Figure FDA0003658127040000021
wherein P (t) is the power of time t; l, h and b are respectively the length, the height and the width of the furnace core of the resistance furnace; t is w n The temperature of the outer surface of the furnace core at any moment can be measured; t is i,j 0 : when not electrified initially, the temperature of the whole resistance furnace; λ: the thermal conductivity of the material; q. q.s w : the heat flux density passing through the outer surface of the furnace core;
the temperature gradient model is shown as follows:
Figure FDA0003658127040000022
in the formula, λ: the coefficient of thermal conductivity of the material; Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; p (t): a function of power over time; t is 1 ......T r : (ii) temperature;
the step of obtaining a temperature gradient function through the temperature gradient model and obtaining a power time model of the resistance furnace based on a time-temperature relation comprises the following steps:
obtaining a temperature gradient function according to the temperature gradient model, wherein the temperature gradient function is as follows: t (x,0, T), wherein: x is the abscissa and t is the time;
arbitrarily setting a power value, obtaining a time-temperature relation through the temperature gradient function, and obtaining a time-temperature relation tableThe expression is as follows:
Figure FDA0003658127040000023
in the formula: t (t) is: a function of temperature over time; p is power, p k Is taken as a power value, x is the abscissa, x n Taking values for the abscissa;
simultaneously solving to obtain the power-time relation of the resistance furnace under the condition of any power value so as to obtain a power-time model of the resistance furnace;
the arbitrary given power value is set by the minimum regulating time of the transformer and the reaction of the resistance furnace;
the minimum regulating time of the transformer is the minimum unit of power change;
the resistance furnace reaction comprises furnace spraying pressure and thermal reduction reaction rate;
the construction of the power load model of the auxiliary material system comprises the following steps: grouping the auxiliary material systems according to whether the auxiliary material systems are transferable loads, and constructing a power load model of the auxiliary material systems for the transferable loads by using the following formula:
W auxiliary device =W P +W S +W H +W L (8)
In the formula, W P : the amount of electricity consumed by the group with the crusher as the primary equipment; w S : the amount of electricity consumed by the group using the screening machine as the main equipment; w H : the electric quantity consumed by the group taking the mixer as main equipment; w L : no amount of power consumed by the group of primary devices;
the method for determining the temperature rise stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power comprises the following steps:
and determining the temperature rising stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power and the temperature rising stage of the resistance furnace divided by combining the temperature of the resistance furnace when the wind power arrives and a temperature gradient model.
2. The method of determining wind power absorption potential of claim 1 wherein the resistive furnace temperature ramp phase divided by the temperature gradient model comprises:
the first stage is as follows: heating the material, wherein the highest temperature of the material is 1450 ℃, and no silicon carbide is generated;
and a second stage: heating the material, wherein the highest temperature of the material is 1800 ℃, and only beta-SiC is generated;
and a third stage: heating the material, wherein the highest temperature of the material is 2200 ℃, and 4H-alpha-SiC is generated;
a fourth stage: heating the material, wherein the highest temperature of the material is 2600 ℃, and 6H-alpha-SiC is generated;
the fifth stage: the material is heated, the highest temperature of the material exceeds 2600 ℃, and SiC is decomposed into carbon and silicon.
3. The method for determining wind power consumption potential of claim 1, wherein the grouping of auxiliary material systems comprises: and (4) grouping the belt conveyor, the bucket elevator, the crusher, the screening machine and the mixer.
4. A wind power consumption potential determination system, the system comprising: the system comprises a first construction module, a second construction module, a third construction module and a calculation module;
the first construction module is used for constructing a wind power consumption potential model;
the second construction module is used for constructing a power time model of the resistance furnace;
the third construction module is used for constructing a power load model of the auxiliary material system;
the calculation module is used for calculating the wind power consumption potential of the silicon carbide enterprise according to the arrival time and the duration of wind power and the temperature of the resistance furnace when the wind power arrives by combining with a wind power consumption potential model;
the second building block is specifically configured to:
constructing a silicon carbide smelting two-dimensional heat transfer model;
carrying out differential processing on the two-dimensional heat transfer model, and combining boundary conditions to obtain a temperature gradient model on the premise of meeting stability conditions;
obtaining a temperature gradient function from the temperature gradient model, and obtaining a power time model of the resistance furnace based on a time-temperature relation;
the silicon carbide smelting two-dimensional heat transfer model is shown as the following formula:
Figure FDA0003658127040000041
wherein λ is the thermal conductivity of the material, ρ is the density of the material, C p Is the specific heat capacity of the material, q v To generate the heat of reduction reaction of silicon carbide;
the stability conditions are shown below:
Figure FDA0003658127040000042
in the formula, Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; Δ t: dividing the time period by the interval Δ t; c p : material specific heat capacity; λ: the coefficient of thermal conductivity of the material; ρ: the density of the material;
the boundary condition is shown as follows:
T i,j 0 =30℃ (3)
T i,j 0 =T w n (4)
Figure FDA0003658127040000043
Figure FDA0003658127040000044
wherein P (t) is the power at time t; l, h and b are respectively the length, the height and the width of the furnace core of the resistance furnace; t is w n The temperature of the outer surface of the furnace core at any momentA measurable value; t is i,j 0 : when not electrified initially, the temperature of the whole resistance furnace; λ: the coefficient of thermal conductivity of the material; q. q.s w : the density of the heat flux passing over the outer surface of the core;
the temperature gradient model is shown as follows:
Figure FDA0003658127040000051
in the formula, λ: the coefficient of thermal conductivity of the material; Δ x and Δ y: respectively representing the spacing distances of the rectangular grids divided in the X-axis direction and the Y-axis direction; p (t): a function of power over time; t is 1 ......T r : (ii) temperature;
obtaining a temperature gradient function according to the temperature gradient model, wherein the temperature gradient function is as follows: t (x,0, T), wherein: x is the abscissa and t is the time;
arbitrarily giving a power value, and obtaining a time-temperature relationship through the temperature gradient function, wherein the time-temperature relationship expression is as follows:
Figure FDA0003658127040000052
in the formula: t (t) is: a function of temperature over time; p is power, p k Is taken as a power value, x is the abscissa, x n Taking values for the abscissa;
simultaneously solving to obtain the power-time relation of the resistance furnace under the condition of any power value so as to obtain a power-time model of the resistance furnace;
the arbitrary given power value is set by the minimum regulating time of the transformer and the reaction of the resistance furnace;
the minimum regulating time of the transformer is the minimum unit of power change;
the resistance furnace reaction comprises furnace spraying pressure and thermal reduction reaction rate;
the third building block is specifically configured to:
grouping the auxiliary material systems according to whether the auxiliary material systems are transferable loads, and constructing a power load model of the auxiliary material systems for the transferable loads by using the following formula:
W auxiliary device =W P +W S +W H +W L (8)
In the formula, W P : the amount of electricity consumed by the group with the crusher as the primary equipment; w S : the amount of electricity consumed by the group using the screening machine as the main equipment; w H : the electric quantity consumed by the group taking the mixer as main equipment; w L : no amount of power consumed by the group of primary devices;
the calculation module is specifically configured to:
and determining the temperature rising stage of the resistance furnace when the wind power arrives according to the arrival time and duration of the wind power and the temperature rising stage of the resistance furnace divided by combining the temperature of the resistance furnace when the wind power arrives and a temperature gradient model.
CN201710638461.5A 2017-07-31 2017-07-31 Wind power consumption potential determination method and system Active CN107609732B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710638461.5A CN107609732B (en) 2017-07-31 2017-07-31 Wind power consumption potential determination method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710638461.5A CN107609732B (en) 2017-07-31 2017-07-31 Wind power consumption potential determination method and system

Publications (2)

Publication Number Publication Date
CN107609732A CN107609732A (en) 2018-01-19
CN107609732B true CN107609732B (en) 2022-09-02

Family

ID=61063833

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710638461.5A Active CN107609732B (en) 2017-07-31 2017-07-31 Wind power consumption potential determination method and system

Country Status (1)

Country Link
CN (1) CN107609732B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103219748A (en) * 2012-01-20 2013-07-24 华北电力科学研究院有限责任公司 Determination method and system of power grid wind electricity digestion capability
CN103810535A (en) * 2014-01-07 2014-05-21 广西大学 Wind power absorption capacity assessment method for electric power system
KR20150005185A (en) * 2013-07-05 2015-01-14 현대중공업 주식회사 Apparatus and method for reducing peak power consumption of wind power generation system
CN105281362A (en) * 2015-09-28 2016-01-27 国家电网公司 Assessment method for wind power consumption capability of electric power system
CN105932711A (en) * 2016-04-26 2016-09-07 国网电力科学研究院 Evaluation, measurement and calculation method for wind power generation-load real-time coordination of curtailed wind power accommodation
CN106099991A (en) * 2016-06-20 2016-11-09 广西电网有限责任公司电力科学研究院 A kind of power grid wind electricity digestion capability appraisal procedure

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103219748A (en) * 2012-01-20 2013-07-24 华北电力科学研究院有限责任公司 Determination method and system of power grid wind electricity digestion capability
KR20150005185A (en) * 2013-07-05 2015-01-14 현대중공업 주식회사 Apparatus and method for reducing peak power consumption of wind power generation system
CN103810535A (en) * 2014-01-07 2014-05-21 广西大学 Wind power absorption capacity assessment method for electric power system
CN105281362A (en) * 2015-09-28 2016-01-27 国家电网公司 Assessment method for wind power consumption capability of electric power system
CN105932711A (en) * 2016-04-26 2016-09-07 国网电力科学研究院 Evaluation, measurement and calculation method for wind power generation-load real-time coordination of curtailed wind power accommodation
CN106099991A (en) * 2016-06-20 2016-11-09 广西电网有限责任公司电力科学研究院 A kind of power grid wind electricity digestion capability appraisal procedure

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Optimal Control for Minimizing Power Consumption During Holding Patterns for Airborne Wind Energy Pumping System;G. Licitra,S. Sieberling,S. Engelen,P. Williams,;《2016 European Control Conference》;20160701;1574-1579 *
高载能负荷参与的电网消纳风/光电能力定量预评估方法;刘聪等;《电网技术》;20150105(第01期);223-229 *
高载能负荷参与调节消纳受阻风电的控制策略;吕良,刘文颖,吕思琦,朱丹丹,梁琛,周强;《电网与清洁能源》;20170625;103-115 *

Also Published As

Publication number Publication date
CN107609732A (en) 2018-01-19

Similar Documents

Publication Publication Date Title
Meng et al. Flexible system integration and advanced hierarchical control architectures in the microgrid research laboratory of Aalborg University
CN105846456B (en) A kind of AC-DC interconnecting power network wind, fire coordinate dynamic economic dispatch optimization method
JP6187463B2 (en) Grid integrated control device, grid control system, grid control device, program, and control method
CN103471176B (en) Energy-saving linkage control system and energy-saving linkage control method
Awad et al. Integrated optimal power flow for electric power and heat in a microgrid
CN103617552A (en) Power generation cost optimization method for iron and steel enterprise
CN107408816A (en) It is variable to feed out energy management
CN111667136A (en) Clearing method and device for regional power market and storage medium
CN107609732B (en) Wind power consumption potential determination method and system
CN107748963A (en) A kind of electric network coordinated operation method based on the optimization of OCD distributings
Kharitonov et al. Justifying the actuator type selection for the flow-line conveyor system with vibration-assisted transportation of bulk material
CN102790401B (en) Method for inhibiting fluctuation of random power of alternative-current junctor
CN107611977A (en) Industrial thermic load group participates in the control system and method for primary frequency regulation of power network
CN112330493A (en) Energy system modeling and comprehensive analysis method, device and storage medium
CN111064231A (en) New energy graded interactive consumption method and system
CN203718923U (en) Energy-saving coordinated control system
Huang et al. Numerical simulation of multiple polysilicon CVD reactors connected in series using CFD method
CN106326568A (en) Method for simulating and analyzing silicon carbide synthesis furnace based on CFD (computational fluid dynamics) technique
CN103033065B (en) Control method and device of afterheat power generation system
Sayed et al. A tight misocp formulation for the integrated electric-gas system scheduling problem
CN115096033A (en) Refrigeration strategy generation method, system, terminal and medium based on body temperature change
CN107632520B (en) Wind power consumption potential model optimization control method and system
CN107482663B (en) Polycrystalline silicon load-based frequency control method for wind power isolated power grid with high permeability
CN108876196B (en) Alternate iterative optimization scheduling method based on non-protrusion force characteristic of cogeneration unit
Xia et al. Numerical simulation of chemical vapor deposition reaction in polysilicon reduction furnace

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
GR01 Patent grant
GR01 Patent grant