CN109577064A - Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine - Google Patents
Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine Download PDFInfo
- Publication number
- CN109577064A CN109577064A CN201811532221.8A CN201811532221A CN109577064A CN 109577064 A CN109577064 A CN 109577064A CN 201811532221 A CN201811532221 A CN 201811532221A CN 109577064 A CN109577064 A CN 109577064A
- Authority
- CN
- China
- Prior art keywords
- air
- gas hood
- page
- formula
- drying
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- D—TEXTILES; PAPER
- D21—PAPER-MAKING; PRODUCTION OF CELLULOSE
- D21F—PAPER-MAKING MACHINES; METHODS OF PRODUCING PAPER THEREON
- D21F5/00—Dryer section of machines for making continuous webs of paper
- D21F5/02—Drying on cylinders
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Paper (AREA)
Abstract
The invention discloses a mechanism modeling method for predicting energy consumption and evaporation capacity of a drying part of a toilet paper machine, which comprises the following steps: exporting input and output variable data of the drying part and preprocessing the input and output variable data; establishing a drying part mechanism model which comprises an air hood ventilation model, a paper sheet drying dynamic model and a drying energy consumption model; determining model parameters, including manual measurement and parameter adjustment by using the processed partial data; and (4) model verification, namely performing simulation operation by using the residual data, comparing a simulation result with actual output, and evaluating the prediction precision and the prediction performance of the optimized model. The invention obtains higher simulation precision on the numerical value and the variation trend of the drying energy consumption, and simultaneously can predict the evaporation dehydration capacity of the drying part, thereby providing convenience for energy conservation, consumption reduction and quality control of toilet paper factories.
Description
Technical field
The present invention relates to papermaking dry fields, and in particular to one kind is predicted for toilet paper machine drying section energy consumption and evaporation capacity
Modelling by mechanism method.
Background technique
Paper for daily use production capacity quicklys increase in recent years, and energy cost goes up, and how persistently to reduce production energy consumption, improves efficiency
Level is the common issue that service workers face.Drying section is the maximum workshop section of energy consumption in pulping and papermaking processes, accounts for about production
67% or more of process total energy consumption, and the operation of drying section has strong influence to its energy consumption.It is found by investigation, it is practical raw
In production, often when page just leaves drying cylinder, it is whether up to standard that web dryness is judged by feel, to adjust drying cylinder steaming
The temperature of pressure (saturation temperature) and the gas hood air-supply of vapour.Adjust which parameter, adjust to how many there is no specific schemes and
Research, therefore will appear excessively dry and energy the waste of page, and repeatedly quality fluctuation caused by feedback regulation.
Theoretical foundation is provided to the operation to drying section, modeling is important means.In paper: Drying Process modeling
In research, Lin Zhizuo establishes the drying dynamic model of toilet paper, but mode input and non-operational data, it is difficult to obtain, and mould
Type is not related to drying section energy consumption.In paper: in the paper for daily use machine drying section energy optimization research based on improved adaptive GA-IAGA, poplar
Profit coral establishes the energy consumption model of toilet paper machine drying section based on energy and material balance, but is not related to dry dynamic study, no
It can predict into the mass dryness fraction of paper and the dehydration by evaporation amount of drying section.In order to provide theoretical direction, drying section model to practical operation
It not only needs to carry out good prediction to energy consumption, while must realize to the dehydration by evaporation energy at paper mass dryness fraction index or drying section
The prediction of power could constrain energy optimization, and the energy consumption prediction model for lacking constraint condition is not no directive function.
Summary of the invention
In view of the deficiencies of the prior art, it is an object of the present invention to provide a kind of for toilet paper machine drying section energy consumption and steams
The modelling by mechanism method of hair amount prediction, the method had both considered the evaporation capacity of page compared to previous modelling by mechanism method
Prediction, it is also considered that dry energy consumption sufficiently meets the demand of Instructing manufacture operation.
The purpose of the present invention can be achieved through the following technical solutions:
A kind of modelling by mechanism method for toilet paper machine drying section energy consumption and evaporation capacity prediction, the method includes following
Step:
S1, export drying section input/output variable data are simultaneously pre-processed: choosing drying cylinder vapor (steam) temperature, the paper of drying section
Air mass flow, the supply air temperature of gas hood and the ambient temperature and humidity that locomotive speed, page quantification, blower pass through are mode input;Choosing
Drying cylinder steam flow, gas hood steam flow and drying process dehydration by evaporation amount is taken to export as model;It is led from historical data base
Data out;The rejecting of abnormal data is carried out according to the parameter area under nominal situation;Evidence of fetching at a certain time interval is put down
Mean value eliminates data fluctuations caused by valve normal regulating;
S2, it establishes drying section mechanism model: establishing the gas hood ventilation model of drying section based on energy and the conservation of mass, utilize
Mode input data carry out hard measurement to the humidity of air, enthalpy and mass flow in page surrounding air and ventilation shaft;Base
Drying Process is modeled in heat and mass rule, the gas hood blasting humidity acquired in gas hood ventilation model is substituted into paper
Page drying kinetic model solves temperature and moisture distribution of the page on drying cylinder;Drying cylinder is solved according to page Temperature Distribution to steam
Steam flow amount solves consumed by the heating of gas hood air according to air enthalpy and flow in the pipeline acquired in gas hood ventilation model
Steam flow;
S3, model parameter is determined: for the preset parameter of paper machine, including breadth, dry section length, the balanced ratio of gas hood,
Page initial temperature and mass dryness fraction, at the design mass dryness fraction of paper, directly progress manual measurement;It conducts heat for undetermined parameter, including drying cylinder
Coefficient, drying cylinder heat loss factor and gas hood heat loss factor carry out tune ginseng using 80% mode input output data, according to mould
Quasi- result and the comparison for exporting result, adjust the size of undetermined parameter;
S4, model verifying: simulation trial is carried out using remaining 20% mode input output data and by analog result and reality
Border output comparison, evaluates the precision of prediction and estimated performance of Optimized model, using average absolute percentage error as standard, carries out mould
Pattern intends the analysis and evaluation of effect.
Further, shown in such as formula (1) of the paper page drying position kinetic model formula in step S2:
Wherein, XpWater entrained by unit mass over dry page for each point on drying cylinder, kg water/kg bone dry fiber;TpFor
The temperature of each point page on drying cylinder, DEG C;L is the fore-and-aft distance that page passes through on drying cylinder, m;KmBetween page and air
Convective transfer coefficient, kg/ (m2·s);G is that the over dry of page is quantitative, kg/m2;V is paper machine speed, m/s;yV0For page surface
Steam quality score;y′V∞For the steam quality score in page local environment air;hs,pFor out of drying cylinder steam to paper
The overall heat-transfer coefficient of page, W/ (m2·℃);TsFor vapor (steam) temperature in drying cylinder, DEG C;hconvConvective heat transfer between page and air
Coefficient, W/ (m2·℃);CAFor the correction factor that conducts heat;TfFor the temperature of page local environment air, DEG C;ΔHv,ΔHsRespectively
The evaporative phase-change heat and heat of adsorption of moisture in page, kJ/kg;Cpf,CpwThe respectively specific heat of fiber and water, kJ/ (kg
℃)。
Further, in step S2, one kind has been used during the gas hood ventilation model for establishing toilet paper machine drying section
Iterative algorithm, to blasting humidity progress it is assumed that solving the temperature and humidity and mass flow of each point in circulating line on this basis, with
Resulting blasting humidity value is calculated instead of assumption value, is iterated calculating until humidity value is stablized.
Further, in step S2, establishing the gas hood ventilation model of toilet paper machine drying section, detailed process is as follows:
According to gas hood ventilation pipeline configuration carry out modelling by mechanism: hot wind is sent into gas hood by circulating fan, blow to page into
Row forces impingement drying, arranges after mixing with the moisture evaporated in the surrounding air and page in gas hood gap sucking workshop from gas hood
Out, two parts are divided into after discharge, a part is extracted out by exhaust blower, and another part humid air reuse is mixed with fresh wind, mixing
Wind is heated using steam heater, is entered gas hood by circulating fan, is formed circulation;It is right based on energy and the conservation of mass
Gas hood air-supply, the temperature of air draft and mixing wind in ventilating system, humidity, enthalpy, flow carry out hard measurement;
Since gas hood blasting humidity is unknown, cause its dependent variable can not direct solution, using a kind of building for iterative calculation
Mould method is shown in formula (2) using fresh wind moisture as the initial value of gas hood blasting humidity:
X1 (1)=X6 (2)
Wherein, X indicates air absolute humidity, kg water/kg dry air;X1 (1)The air-supply of gas hood is wet when for first time iteration
Degree, subscript (1) represent first time iteration, X6For the humidity of fresh wind, it is not involved in iteration;
According to the disengaging gas hood air energy conservation of mass, the temperature and humidity of gas hood air draft and flow blown by gas hood and leak out with
And evaporated water influences, and sees formula (3):
Wherein, X1、X2、X4To be respectively gas hood air-supply, the air humidity leaked out with gas hood air draft, kg water/kg dry air;The respectively gas hood mass flow of blowing, leak out, evaporating water and gas hood air draft, kg/s;H1、H2、H3、
H4The respectively gas hood enthalpy of blowing, leak out, evaporating water and gas hood air draft, kJ/kg;
Shown in the density of humid air needed for energy conservation of matter calculating process and the calculation formula such as formula (4) of enthalpy:
Wherein, pw, p be respectively humid air stagnation pressure and practical water vapor pressure, pa;ρ is density, kg/m3;T is air themperature,
℃;X is air absolute humidity, g/kg;H is air enthalpy, kJ/kg dry air;
The flow that gas hood leaks out can determine by gas hood balanced ratio, see formula (5):
Wherein, HB is gas hood balanced ratio, i.e. the ratio between the dry air mass flow from pipeline disengaging gas hood, passes through manual measurement
The mode to convert afterwards determines that conversion method is shown in formula (6):
Evaporated water is calculated according to page part water content, sees formula (7):
Wherein, W is width of the page on drying cylinder, m;Xa、XdRespectively page is aqueous when entering and leaving dry section
Amount, kg water/kg dry fibers;G is that the over dry of page is quantitative, kg/m2;V is paper machine speed, m/s;
Formula (8) are seen with the mass flow of the fresh wind of determination according to the entire ventilating system dry air conservation of mass:
Wherein, X5、X6The humid air that is respectively discharged from circulating line and the fresh air for entering circulating line it is wet
Degree, kg water/kg dry fibers;m5、m6The fresh air of the humid air and entrance circulating line that are respectively discharged from circulating line
Mass flow, kg/s;
According to air energy mass-conservation equation (9) in pipeline, the mass flow m of mixing wind is solved7, humidity X7With enthalpy
H7:
Heating front and back humidity is constant, sees formula (10):
The X that will be acquired1Substitution formula (3), (5), (8), (9), (10) are iterated calculating again, after iteration is multiple, ventilation
The numerical value of gas hood air-supply, the humidity of air draft and mixing wind, enthalpy, flow in system reaches stable, and divulges information mould as gas hood
The foundation of toilet paper machine drying section gas hood ventilation model is completed in the output of type.
Further, in step S2, after establishing drying section gas hood ventilation model, by the gas hood blasting humidity acquired and initially
Mode input data: air mass flow that drying cylinder vapor (steam) temperature, paper machine speed, page quantification, blower pass through, gas hood air-supply temperature
Degree substitutes into paper page drying position kinetic model together, solves temperature T of the page on drying cylinderpWith moisture XpDistribution, paper page drying position are dynamic
Shown in the main body of mechanical model such as formula (1);
Water vapor in air mass fraction in formula (1) converts according to blasting humidity, as shown in formula (11):
Convective heat-transfer coefficient in formula (1) is calculated as shown in formula (12):
Wherein Nu is nusselt number, λfFor the air thermal conductivity of page local environment, W/ (mk);LconvFor pair of gas hood
Spread thermal characteristics length, m;
The calculation method of nusselt number uses different formula in different situations, paper and air heat and mass outside gas hood
Situation is similar to the situation that turbulent flow flows through plate, is calculated using standard Nu Saier relational expression, as shown in formula (13):
Wherein, Pr is the Prandtl number of air;Re is the Reynolds number for describing air fluid properties;
Rule of thumb formula calculates the nusselt number of air in gas hood, as shown in formula (14):
Wherein f ' is gas hood nozzle bore rate, and h is distance of the nozzle to page, m;D ' is nozzle diameter, m;TaFor gas hood
Supply air temperature, DEG C;
Correction factor in formula (1) is in the evaporation mass transport process of page moisture, and evaporation water generates white heat
The scale factor that is quantified of influence, the calculating of correction factor is shown in formula (15):
Wherein, E is the dimensionless factor of heat transfer correction factor;mevFor the evaporated water in dry section, kg/s being capable of basis
Evaporated water in formula (7) calculates;A is dry section area, m2;CpVFor the specific heat at constant pressure of water vapour, J/ (kg DEG C);
Convective transfer coefficient in formula (1) is determined by the thermophysical property of air at heat transfer coefficient and mass transfer film, sees formula
(16):
Wherein, pTFor atmospheric pressure, pa;TfFor page local environment air themperature, DEG C;λfFor page local environment air
Thermal conductivity, W/ (mk);ρfFor atmospheric density, kg/m3;CpfFor the specific heat of fiber, kJ/ (kg DEG C);D is the diffusion of air
Coefficient, m2/s;
So far, longitudinally above the temperature and moisture distribution of every bit are solved page, and drying kinetic model, which derives, to be completed;
Finally, solving drying cylinder energy consumption according to page Temperature Distribution, the enthalpy variation before and after heating of being blown according to gas hood solves heater energy
Consumption;
Drying cylinder energy consumption is divided into the heat loss of the efficient drying energy consumption and drying cylinder shell that are transmitted to page, the steam stream of drying cylinder consumption
Amount, which solves, sees formula (17):
Wherein, TEFor environment temperature;W is page breadth, m;L is stroke of the page on drying cylinder, m;ACi、hCiRespectively dry
Cylinder does not cover the heat transfer coefficient and area in i-th piece of region in the partial shell of page, ∑ ACihCiFor each point heat transfer coefficient on drying cylinder
With the sum of products of area, as undetermined parameter;ΔHCFor the enthalpy change of drying cylinder steam to condensed water, kJ/kg;For drying cylinder consumption
Steam flow, kg/h;A, position when d is respectively page contact and leaves dry section;
Gas hood air heating energy consumption includes that air heating energy consumption and heat exchanger conduct heat to environment, sees formula (18):
Wherein, Δ HHFor latent heat of the steam in air heater;For the steam flow of gas hood air heating consumption, kg/
h;∑AHihHiFor the sum of products of each point heat transfer coefficient and area on gas hood, as undetermined parameter.
Further, in step S3, for can not undetermined parameter measured directly, carried out using mode input output data
Ginseng is adjusted, by comparing the analogue value and actual value of drying section evaporation capacity, adjusts heat transfer coefficient hs,pSize;It is gentle according to drying cylinder
The steam flow analogue value of cover and the comparison of actual value adjust separately heat loss factor's ∑ A of drying cylinderCihCiWith the heat waste of gas hood
Lose coefficient ∑ AHihHi。
Further, in step S4, shown in the calculation formula of the average absolute percentage error such as formula (19):
Wherein,For model output data, yiFor actual value, n is sample number, and MAPE is average absolute percentage error.
Compared with the prior art, the invention has the following advantages and beneficial effects:
Provided by the present invention for the modelling by mechanism method of toilet paper machine drying section energy consumption and evaporation capacity prediction, compared to previous
Method it is more comprehensive, both considered page evaporation capacity prediction, it is also considered that dry energy consumption sufficiently meets Instructing manufacture
The demand of operation.The present invention is while predicting result, to more intermediate variables such as gas hood air humidity, page
Temperature Distribution has carried out hard measurement, provides convenience for the efficiency research of drying section.
Detailed description of the invention
Fig. 1 is modelling by mechanism Method And Principle of the embodiment of the present invention for toilet paper machine drying section energy consumption and evaporation capacity prediction
Figure.
Fig. 2 is the structural schematic diagram of toilet paper machine drying section.
Fig. 3 is the ventilating system schematic diagram of toilet paper machine.
Fig. 4 is the gas hood steam flow comparison diagram of toilet paper machine of the embodiment of the present invention.
Fig. 5 is the gas hood steam flow Error Graph of toilet paper machine of the embodiment of the present invention.
Fig. 6 is the drying cylinder steam flow comparison diagram of toilet paper machine of the embodiment of the present invention.
Fig. 7 is the drying cylinder steam flow Error Graph of toilet paper machine of the embodiment of the present invention.
Fig. 8 is the evaporation capacity comparison diagram of toilet paper machine of embodiment of the present invention drying section.
Fig. 9 is the evaporation capacity Error Graph of toilet paper machine of embodiment of the present invention drying section.
Specific embodiment
Present invention will now be described in further detail with reference to the embodiments and the accompanying drawings, but embodiments of the present invention are unlimited
In this.
Embodiment:
As shown in Figure 1, present embodiments providing a kind of mechanism for toilet paper machine drying section energy consumption and evaporation capacity prediction
Modeling method establishes mechanism model including the heat and mass transfer process for drying section, uses the number of operation data and manual measurement
According to come the parameter that determines model, contrast model output and reality output amount, the prediction effect of evaluation model.It specifically include following step
It is rapid:
Step 1: export drying section input/output variable data are simultaneously pre-processed: being chosen the operating parameter of drying section, wrapped
Including paper machine speed, the air mass flow that quantitative, drying cylinder vapor (steam) temperature, the supply air temperature of gas hood, blower pass through and ambient temperature and humidity is
Mode input;It chooses steam flow consumed by gas hood heater and drying cylinder and is exported at the evaporation capacity of paper as model;From
Initial data is exported in historical data base;All data are matched in temporal sequence, reject nominal situation range outside, stop
Data when machine and in start-up course reject the incomplete data of variable match, with being averaged for the time interval access evidence of 30min
Value eliminates the interference that valve fluctuation generates.
Step 2: the realization of drying section mechanism model.
The drying section of target basis machine has structure shown in Fig. 2, forms by Flying Dutchman and using the gas hood of circulating air.
L Water Paper page after squeezing is transferred to drying cylinder by carrying roller, and most of stroke of the page on drying cylinder is covered by gas hood, be evaporated de-
Water and discharge moisture.After page air leaving cover, drying cylinder is left by scraper.The circulating ventilation pipeline of gas hood is as shown in figure 3, circulated air
Hot wind is sent into gas hood by machine, is blowed to page and is carried out pressure impingement drying, and gas hood gap sucks the cold air in workshop, evaporates in page
Moisture, gas hood discharge humid air be divided into two parts.A part is extracted out by exhaust blower, another part humid air reuse, with
Fresh wind mixing, mixing wind are heated using steam heater, are entered gas hood by circulating fan, are formed circulation.
Firstly, according to gas hood ventilation model, to the humidity of gas hood air-supply, air draft and mixing wind in ventilating system, enthalpy,
The numerical value of flow carries out hard measurement.
According to input data, it is known that the supply air temperature T of gas hood1And mass flowThe air mass flow of circulating line dischargeThe temperature T of environmentEWith humidity XE, temperature T that gas hood leaks out2With humidity X2, the temperature T of fresh wind6With humidity X6, gas hood leakage
The temperature and humidity of wind and fresh wind is identical as environment.
The step of using solution by iterative method gas hood ventilation model: the humidity X of fresh wind is taken6For wind pushing air humidity X1Just
Initial value, as shown in formula (2).The humidity X of gas hood air draft is solved according to the disengaging gas hood air energy conservation of mass4, flowWith enthalpy
Value H4, as shown in formula (3).Shown in atmospheric density used in calculating process and enthalpy formula such as formula (4).Utilize gas hood balanced ratio
The flow that gas hood leaks out is determined according to formula (5)Shown in the measurement of gas hood balanced ratio and reduction formula such as formula (6).According to formula
(7) evaporated water is calculatedThe fresh air flow for entering circulating line is calculated according to formula (8)According in ventilation shaft
The humidity X of air energy conservation of mass solution mixing wind7, flowWith enthalpy H7, as shown in formula (9).By the humidity of mixing wind
X7Gas hood blasting humidity X after being assigned to heating1, as shown in formula (10).By X1Formula (3), (5), (8), (9), (10) are substituted into again
It is iterated calculating.In this way after iteration 20 times, the temperature of gas hood air-supply, air draft and mixing wind in ventilating system, humidity,
The numerical value of enthalpy, flow, the output as gas hood ventilation model.
Then, gas hood blasting humidity previous step acquired and initial mode input data (drying cylinder vapor (steam) temperature, paper machine
Air mass flow that speed, page quantification, blower pass through, gas hood supply air temperature) together, paper page drying position kinetic model is substituted into, is asked
Solve temperature T of the page on drying cylinderpWith moisture XpDistribution.Shown in the main body of paper page drying position kinetic model such as formula (1).In formula (1)
Water vapor in air mass fraction conversion such as formula (11) shown in.Convective heat-transfer coefficient in formula (1) calculate such as formula (12),
(13), shown in (14).Correction factor in formula (1) is such as shown in (15).Shown in convective transfer coefficient such as formula (16) in formula (1).
The realization step of formula (1) paper page drying position kinetic model: temperature when drying cylinder has just been transmitted to by manual measurement page
Degree and mass dryness fraction, determine the initial water content of pagePage initial temperature TpL0=85 DEG C.Assuming that undetermined parameter
hs,p=650W/ (m2DEG C), it substitutes into formula (1), obtains the moist gradient in the 1st step-length 0.001mThen L1=
Page humidity at 0.001mIt willSubstitution formula (1), obtains
Temperature gradient in 1st step-lengthThe then page temperature Tp at L1=0.001mL1=85+0.001 ×
(- 29.69440)=84.97030 DEG C.Then on this basis, the page temperature that original equation calculates second step-length is substituted into again
TpL2And water contentThe rest may be inferred, calculates page in the temperature Tp of drying cylinder longitudinal direction each pointLiAnd water content
Finally, the Temperature Distribution according to page in drying cylinder longitudinal direction each point solves drying cylinder energy consumptionAs shown in formula (17).Root
It blows according to gas hood and heats the enthalpy and flow solution heater energy consumption of front and backAs shown in formula (18).
Step 3: determining the undetermined parameter of model.
Using 80% mode input output data, tune ginseng is carried out to undetermined parameter.By the drying section evaporation of model output
The analogue value of water is compared with actual value, then adjusts heat transfer coefficient hs,pSize, be allowed to close.It is gentle according to drying cylinder
The analogue value of the steam flow of cover is compared with actual value, adjusts separately heat loss factor's ∑ A of drying cylinderCihCiWith gas hood
Heat loss factor's ∑ AHihHi, keep model output and actual value close.
Step 4: model verifying.
The verifying of model is carried out using remaining 20% data, using average absolute percentage error as standard to steam stream
The simulation precision of amount and evaporation capacity is analyzed and evaluated.Shown in the calculating of average absolute percentage error such as formula (19).
The variation tendency and numerical value of the steam flow of gas hood have all obtained preferable simulation, as shown in Figure 4.Gas hood steam stream
Measuring average absolute percentage error is 3.16%, and error distribution is and relatively stable as shown in figure 5, worst error is only 0.5t/h.
The variation tendency and numerical value of the steam flow of drying cylinder have all obtained preferable simulation, as shown in Figure 6.Drying cylinder steam flow is average exhausted
It is 3.32% to percent error, error distribution is and relatively stable as shown in fig. 7, worst error is only 0.6t/h, can satisfy
Instruct the demand of practical operation.The variation tendency and numerical value of the evaporated water of drying section have all obtained preferable simulation, such as Fig. 8 institute
Show.The evaporated water average absolute percentage error of drying section is 3.40%, and error is distributed as shown in figure 9, worst error is only
0.08kg/s has substantially met the demand for predicting drying section evaporability under each operating condition, provides for the operation during production
Guidance.
The above, only the invention patent preferred embodiment, but the scope of protection of the patent of the present invention is not limited to
This, anyone skilled in the art is in the range disclosed in the invention patent, according to the present invention the skill of patent
Art scheme and its patent of invention design are subject to equivalent substitution or change, belong to the scope of protection of the patent of the present invention.
Claims (7)
1. a kind of modelling by mechanism method for toilet paper machine drying section energy consumption and evaporation capacity prediction, which is characterized in that the side
Method the following steps are included:
S1, export drying section input/output variable data are simultaneously pre-processed: choosing drying cylinder vapor (steam) temperature, the paper locomotive of drying section
The supply air temperature and ambient temperature and humidity of air mass flow, gas hood that speed, page quantification, blower pass through are mode input;It chooses and dries
Cylinder steam flow, gas hood steam flow and drying process dehydration by evaporation amount are exported as model;Number is exported from historical data base
According to;The rejecting of abnormal data is carried out according to the parameter area under nominal situation;The average value for evidence of fetching at a certain time interval,
Eliminate data fluctuations caused by valve normal regulating;
S2, it establishes drying section mechanism model: establishing the gas hood ventilation model of drying section based on energy and the conservation of mass, utilize model
Input data carries out hard measurement to the humidity of air, enthalpy and mass flow in page surrounding air and ventilation shaft;Based on biography
Hot mass transfer rule models Drying Process, and the gas hood blasting humidity acquired in gas hood ventilation model substitution page is done
Dry kinetic model solves temperature and moisture distribution of the page on drying cylinder;Drying cylinder steam stream is solved according to page Temperature Distribution
Amount solves gas hood air and heats consumed steam according to air enthalpy and flow in the pipeline acquired in gas hood ventilation model
Flow;
S3, model parameter is determined: for the preset parameter of paper machine, including breadth, dry section length, the balanced ratio of gas hood, page
Initial temperature and mass dryness fraction, at the design mass dryness fraction of paper, directly progress manual measurement;For undetermined parameter, including drying cylinder heat transfer coefficient,
Drying cylinder heat loss factor and gas hood heat loss factor carry out tune ginseng using 80% mode input output data, are tied according to simulation
The comparison of fruit and output result, adjusts the size of undetermined parameter;
S4, model verifying: using remaining 20% mode input output data progress simulation trial and analog result and reality is defeated
It compares out, evaluates the precision of prediction and estimated performance of Optimized model, using average absolute percentage error as standard, carry out pattern die
The analysis and evaluation of quasi- effect.
2. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 1
Method, which is characterized in that shown in the paper page drying position kinetic model formula such as formula (1) in step S2:
Wherein, XpWater entrained by unit mass over dry page for each point on drying cylinder, kg water/kg bone dry fiber;TpFor drying cylinder
The temperature of upper each point page, DEG C;L is the fore-and-aft distance that page passes through on drying cylinder, m;KmConvection current between page and air
Mass tranfer coefficient, kg/ (m2·s);G is that the over dry of page is quantitative, i.e. the product of page quantification and mass dryness fraction, kg/m2;V is paper locomotive
Speed, m/s;yV0For page surface steam quality score;y′V∞For the steam quality score in page local environment air;
hs,pFor overall heat-transfer coefficient of the steam to page, W/ (m out of drying cylinder2·℃);TsFor vapor (steam) temperature in drying cylinder, DEG C;hconvFor paper
Convective heat-transfer coefficient between page and air, W/ (m2·℃);CAFor the correction factor that conducts heat;TfFor page local environment air
Temperature, DEG C;ΔHv,ΔHsThe evaporative phase-change heat and heat of adsorption of moisture respectively in page, kJ/kg;Cpf,CpwRespectively fiber
With the specific heat of water, kJ/ (kg DEG C).
3. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 1
Method, it is characterised in that: in step S2, used a kind of iteration during the gas hood ventilation model for establishing toilet paper machine drying section
Algorithm carries out blasting humidity it is assumed that solving the temperature and humidity and mass flow of each point in circulating line, on this basis to calculate
Resulting blasting humidity value replaces assumption value, is iterated calculating until humidity value is stablized.
4. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 2
Method, which is characterized in that in step S2, establishing the gas hood ventilation model of toilet paper machine drying section, detailed process is as follows:
Carry out modelling by mechanism according to the pipeline configuration of gas hood ventilation: hot wind is sent into gas hood by circulating fan, is blowed to page and is carried out by force
Impingement drying processed is discharged after mixing with the moisture evaporated in the surrounding air and page in gas hood gap sucking workshop from gas hood, row
It is divided into two parts after out, a part is extracted out by exhaust blower, and another part humid air reuse mixes, mixing wind makes with fresh wind
It is heated with steam heater, gas hood is entered by circulating fan, forms circulation;Based on energy and the conservation of mass, to ventilation
Gas hood air-supply, the temperature of air draft and mixing wind in system, humidity, enthalpy, flow carry out hard measurement;
Since gas hood blasting humidity is unknown, cause its dependent variable can not direct solution, using a kind of modeling side of iterative calculation
Method is shown in formula (2) using fresh wind moisture as the initial value of gas hood blasting humidity:
X1 (1)=X6 (2)
Wherein, X indicates air absolute humidity, kg water/kg dry air;X1 (1)The blasting humidity of gas hood, subscript when for first time iteration
(1) first time iteration, X are represented6For fresh wind moisture, it is not involved in iteration;
According to the disengaging gas hood air energy conservation of mass, the temperature and humidity and flow of gas hood air draft are blown by gas hood and are leaked out and steam
Shampoo amount influences, and sees formula (3):
Wherein, X1、X2、X4To be respectively gas hood air-supply, the air humidity leaked out with gas hood air draft, kg water/kg dry air;The respectively gas hood mass flow of blowing, leak out, evaporating water and gas hood air draft, kg/s;H1、H2、H3、
H4The respectively gas hood enthalpy of blowing, leak out, evaporating water and gas hood air draft, kJ/kg;
Shown in the density of humid air needed for energy conservation of matter calculating process and the calculation formula such as formula (4) of enthalpy:
Wherein, pw, p be respectively humid air stagnation pressure and practical water vapor pressure, pa;ρ is density, kg/m3;T is air themperature, DEG C;X
For air absolute humidity, g/kg;H is air enthalpy, kJ/kg dry air;
The flow that gas hood leaks out can determine by gas hood balanced ratio, see formula (5):
Wherein, HB is gas hood balanced ratio, i.e. the ratio between the dry air mass flow from pipeline disengaging gas hood, by changing after manual measurement
The mode of calculation determines that conversion method is shown in formula (6):
Evaporated water is calculated according to page part water content, sees formula (7):
Wherein, W is width of the page on drying cylinder, m;Xa、XdRespectively water content of the page when entering and leaving dry section,
Kg water/kg dry fibers;G is that the over dry of page is quantitative, kg/m2;V is paper machine speed, m/s;
Formula (8) are seen with the mass flow of the fresh wind of determination according to the entire ventilating system dry air conservation of mass:
Wherein, X5、X6The humidity of the humid air and the fresh air into circulating line that are respectively discharged from circulating line, kg
Water/kg dry fibers;m5、m6The quality of the humid air and the fresh air into circulating line that are respectively discharged from circulating line
Flow, kg/s;
According to air energy mass-conservation equation (9) in pipeline, the mass flow m of mixing wind is solved7, humidity X7With enthalpy H7:
Heating front and back humidity is constant, sees formula (10):
The X that will be acquired1Substitution formula (3), (5), (8), (9), (10) are iterated calculating again, after iteration is multiple, in ventilating system
Gas hood air-supply, the humidity of air draft and mixing wind, enthalpy, flow numerical value reach stable, and as the defeated of gas hood ventilation model
Out, that is, the foundation of toilet paper machine drying section gas hood ventilation model is completed.
5. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 4
Method, which is characterized in that in step S2, after establishing drying section gas hood ventilation model, by the gas hood blasting humidity acquired and initially
Mode input data: air mass flow that drying cylinder vapor (steam) temperature, paper machine speed, page quantification, blower pass through, gas hood supply air temperature,
Paper page drying position kinetic model is substituted into together, solves temperature T of the page on drying cylinderpWith moisture XpDistribution, paper page drying position dynamics
Shown in the main body of model such as formula (1);
Water vapor in air mass fraction in formula (1) converts according to blasting humidity, as shown in formula (11):
Convective heat-transfer coefficient in formula (1) is calculated as shown in formula (12):
Wherein Nu is nusselt number, λfFor the air thermal conductivity of page local environment, W/ (mk);LconvIt is gas hood to spreading
Thermal characteristics length, m;
The calculation method of nusselt number uses different formula in different situations, paper and air heat and mass situation outside gas hood
The situation that plate is flowed through similar to turbulent flow is calculated using standard Nu Saier relational expression, as shown in formula (13):
Wherein, Pr is the Prandtl number of air;Re is the Reynolds number for describing air fluid properties;
Rule of thumb formula calculates the nusselt number of air in gas hood, as shown in formula (14):
Wherein f ' is gas hood nozzle bore rate, and h is distance of the nozzle to page, m;D ' is nozzle diameter, m;TaFor the air-supply of gas hood
Temperature, DEG C;
Correction factor in formula (1) is the shadow that evaporation water generates white heat in the evaporation mass transport process of page moisture
Quantified scale factor is rung, formula (15) are shown in the calculating of correction factor:
Wherein, E is the dimensionless factor of heat transfer correction factor;mevFor the evaporated water in dry section, kg/s can be according to formula
(7) evaporated water in calculates;A is dry section area, m2;CpVFor the specific heat at constant pressure of water vapour, J/ (kg DEG C);
Convective transfer coefficient in formula (1) is determined by the thermophysical property of air at heat transfer coefficient and mass transfer film, sees formula (16):
Wherein, pTFor atmospheric pressure, pa;TfFor page local environment air themperature, DEG C;λfFor the heat of page local environment air
Conductance, W/ (mk);ρfFor atmospheric density, kg/m3;CpfFor the specific heat of fiber, kJ/ (kg DEG C);D is the diffusion system of air
Number, m2/s;
So far, longitudinally above the temperature and moisture distribution of every bit are solved page, and drying kinetic model, which derives, to be completed;Most
Afterwards, drying cylinder energy consumption is solved according to page Temperature Distribution, the enthalpy and flow before and after heating of blowing according to gas hood solve heater energy
Consumption;
Drying cylinder energy consumption is divided into the heat loss of the efficient drying energy consumption and drying cylinder shell that are transmitted to page, and the steam flow of drying cylinder consumption is asked
Solution is shown in formula (17):
Wherein, TEFor environment temperature;W is page breadth, m;L is stroke of the page on drying cylinder, m;ACi、hCiRespectively drying cylinder is not
Cover the heat transfer coefficient and area in i-th piece of region in the partial shell of page, ∑ ACihCiFor each point heat transfer coefficient and face on drying cylinder
Long-pending sum of products, as undetermined parameter;ΔHCFor the enthalpy change of drying cylinder steam to condensed water, kJ/kg;For the steam of drying cylinder consumption
Flow, kg/h;A, position when d is respectively page contact and leaves dry section;
Gas hood air heating energy consumption includes that air heating energy consumption and heat exchanger conduct heat to environment, sees formula (18):
Wherein, Δ HHFor latent heat of the steam in air heater;For the steam flow of gas hood air heating consumption, kg/h;∑
AHihHiFor the sum of products of each point heat transfer coefficient and area on gas hood, as undetermined parameter.
6. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 5
Method, it is characterised in that: in step S3, for can not undetermined parameter measured directly, adjusted using mode input output data
Ginseng adjusts heat transfer coefficient h by comparing the analogue value and actual value of drying section evaporation capacitys,pSize;According to drying cylinder and gas hood
The steam flow analogue value and actual value comparison, adjust separately heat loss factor's ∑ A of drying cylinderCihCiWith the heat loss of gas hood
Coefficient ∑ AHihHi。
7. a kind of modelling by mechanism side for toilet paper machine drying section energy consumption and evaporation capacity prediction according to claim 1
Method, which is characterized in that in step S4, shown in the calculation formula of the average absolute percentage error such as formula (19):
Wherein,For model output data, yiFor actual value, n is sample number, and MAPE is average absolute percentage error.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811532221.8A CN109577064B (en) | 2018-12-14 | 2018-12-14 | Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811532221.8A CN109577064B (en) | 2018-12-14 | 2018-12-14 | Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109577064A true CN109577064A (en) | 2019-04-05 |
CN109577064B CN109577064B (en) | 2020-06-19 |
Family
ID=65928577
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811532221.8A Active CN109577064B (en) | 2018-12-14 | 2018-12-14 | Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109577064B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111241717A (en) * | 2020-03-11 | 2020-06-05 | 广州博依特智能信息科技有限公司 | Method for optimizing operation parameters of drying part of toilet paper machine based on mechanism model |
CN111241754A (en) * | 2020-01-17 | 2020-06-05 | 广州博依特智能信息科技有限公司 | Soft measurement method for key process parameters of paper sheet drying |
CN111893791A (en) * | 2020-07-17 | 2020-11-06 | 广州博依特智能信息科技有限公司 | Method for optimizing operation of drying part of domestic paper making machine based on intelligent algorithm |
CN111893793A (en) * | 2019-05-06 | 2020-11-06 | 维美德技术有限公司 | Method and apparatus for controlling energy consumption in a fiber web manufacturing process |
CN111926606A (en) * | 2020-07-30 | 2020-11-13 | 衢州学院 | Online monitoring method for energy efficiency of drying part of paper machine |
CN111967152A (en) * | 2020-08-10 | 2020-11-20 | 中国华能集团有限公司 | Method for correcting standard k-model by adding k-source term |
CN112949162A (en) * | 2021-01-25 | 2021-06-11 | 广州博依特智能信息科技有限公司 | Operation optimization method for drying part energy system of boxboard paper machine based on data driving |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020124982A1 (en) * | 2001-01-23 | 2002-09-12 | Takashi Sasaki | Method for estimating dryer vapor pressure in papermaking machine and apparatus therefor |
JP2002302890A (en) * | 2001-04-04 | 2002-10-18 | Yokogawa Electric Corp | Method for controlling paper machine and apparatus therefor |
CN101526804A (en) * | 2009-04-17 | 2009-09-09 | 华南理工大学 | Analysis method of heat exchange system of drying section of paper machine |
CN101760980A (en) * | 2010-01-22 | 2010-06-30 | 华南理工大学 | Monitoring control system of paper-making machine drying part energy system and working method thereof |
CN101792981A (en) * | 2010-03-23 | 2010-08-04 | 华南理工大学 | Optimization control system for energy system of drying part of paper machine |
CN201716569U (en) * | 2009-11-30 | 2011-01-19 | 华南理工大学 | Optimization system of heat and power cogeneration energy system in paper mill |
-
2018
- 2018-12-14 CN CN201811532221.8A patent/CN109577064B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020124982A1 (en) * | 2001-01-23 | 2002-09-12 | Takashi Sasaki | Method for estimating dryer vapor pressure in papermaking machine and apparatus therefor |
JP2002302890A (en) * | 2001-04-04 | 2002-10-18 | Yokogawa Electric Corp | Method for controlling paper machine and apparatus therefor |
CN101526804A (en) * | 2009-04-17 | 2009-09-09 | 华南理工大学 | Analysis method of heat exchange system of drying section of paper machine |
CN201716569U (en) * | 2009-11-30 | 2011-01-19 | 华南理工大学 | Optimization system of heat and power cogeneration energy system in paper mill |
CN101760980A (en) * | 2010-01-22 | 2010-06-30 | 华南理工大学 | Monitoring control system of paper-making machine drying part energy system and working method thereof |
CN101792981A (en) * | 2010-03-23 | 2010-08-04 | 华南理工大学 | Optimization control system for energy system of drying part of paper machine |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111893793A (en) * | 2019-05-06 | 2020-11-06 | 维美德技术有限公司 | Method and apparatus for controlling energy consumption in a fiber web manufacturing process |
CN111241754A (en) * | 2020-01-17 | 2020-06-05 | 广州博依特智能信息科技有限公司 | Soft measurement method for key process parameters of paper sheet drying |
CN111241754B (en) * | 2020-01-17 | 2023-08-11 | 广州博依特智能信息科技有限公司 | Soft measurement method for key process parameters of paper drying |
CN111241717A (en) * | 2020-03-11 | 2020-06-05 | 广州博依特智能信息科技有限公司 | Method for optimizing operation parameters of drying part of toilet paper machine based on mechanism model |
CN111893791A (en) * | 2020-07-17 | 2020-11-06 | 广州博依特智能信息科技有限公司 | Method for optimizing operation of drying part of domestic paper making machine based on intelligent algorithm |
CN111926606A (en) * | 2020-07-30 | 2020-11-13 | 衢州学院 | Online monitoring method for energy efficiency of drying part of paper machine |
CN112900136A (en) * | 2020-07-30 | 2021-06-04 | 衢州学院 | On-line monitoring method for key process parameters of paper machine drying cylinder |
CN111967152A (en) * | 2020-08-10 | 2020-11-20 | 中国华能集团有限公司 | Method for correcting standard k-model by adding k-source term |
CN112949162A (en) * | 2021-01-25 | 2021-06-11 | 广州博依特智能信息科技有限公司 | Operation optimization method for drying part energy system of boxboard paper machine based on data driving |
Also Published As
Publication number | Publication date |
---|---|
CN109577064B (en) | 2020-06-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109577064A (en) | Mechanism modeling method for predicting energy consumption and evaporation capacity of drying part of toilet paper machine | |
Heo et al. | Dynamic modeling of paper drying processes | |
CN101792981B (en) | Optimization control system for energy system of drying part of paper machine | |
CN111241717A (en) | Method for optimizing operation parameters of drying part of toilet paper machine based on mechanism model | |
CN106949598A (en) | Network center's machine room energy-saving optimization method when network traffic load changes | |
Di Marco et al. | Mathematical modelling and energy performance assessment of air impingement drying systems for the production of tissue paper | |
CN113609660B (en) | Numerical simulation method for predicting heat transfer rule of clothes in clothes dryer | |
WO1996005371A1 (en) | Method and apparatus for controlling moisture of product during paper changing time in paper machine | |
CN103939941A (en) | Method for boiler combustion optimization with combination of irreversible thermodynamics | |
CN110096845A (en) | A kind of mixed flow type closed cooling tower design and calculation method | |
CN104964527B (en) | Lumber kiln | |
Shao et al. | Investigation on weights setting rule for weights-based fuzzy logic control algorithm utilized in direct expansion air-conditioning systems | |
Ghodbanan et al. | Optimization for energy consumption in drying section of fluting paper machine | |
CN111926606B (en) | Online monitoring method for energy efficiency of drying part of paper machine | |
CN207527938U (en) | Combined drying system | |
Ahn et al. | Experimental study on drying time and energy consumption of a vented dryer | |
Jankes et al. | Waste heat potentials in the drying section of the paper machine in umka cardboard mill | |
Laskowski | A mathematical model of a steam condenser in off-design operation | |
Gelir et al. | The effect of the heat recovery on fuel consumption in the stenter machine | |
CN108669616A (en) | A kind of drum-type cut-tobacco drier moisture content of cut tobaccos control method based on water balance | |
Laskowski | A mathematical model of a steam condenser in the changed conditions | |
JP2018044719A (en) | Prediction method, device and program for inlet temperature at crusher of circulatory system pulverization plant, and computing method, device and program for raw material supply amount | |
Godin et al. | Model-based framework for the evaluation of energy-efficiency measures in multi-cylinder paper drying | |
CN205388414U (en) | Acetate drying machine | |
CN205512310U (en) | Tobacco leaf drying -machine and tobacco leaf baking house |
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 | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220516 Address after: 510663 room 1021, No.2, Tengfei 1st Street, Zhongxin Guangzhou Knowledge City, Huangpu District, Guangzhou City, Guangdong Province Patentee after: GUANGZHOU POI INTELLIGENT INFORMATION TECHNOLOGY CO.,LTD. Address before: 510640 No. five, 381 mountain road, Guangzhou, Guangdong, Tianhe District Patentee before: SOUTH CHINA University OF TECHNOLOGY |
|
TR01 | Transfer of patent right |