CN104950960B - A kind of control method of supercritical extract process temperature and pressure - Google Patents
A kind of control method of supercritical extract process temperature and pressure Download PDFInfo
- Publication number
- CN104950960B CN104950960B CN201510336215.5A CN201510336215A CN104950960B CN 104950960 B CN104950960 B CN 104950960B CN 201510336215 A CN201510336215 A CN 201510336215A CN 104950960 B CN104950960 B CN 104950960B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- pressure
- 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.)
- Expired - Fee Related
Links
Abstract
The invention discloses a kind of supercritical extract process temperature and the accurate control method of pressure, this method has non-linear for carbon dioxide supercritical extraction process temperature and pressure, close coupling and interference can not survey, and to temperature and pressure control accuracy requirement it is high the characteristics of, on the basis of extraction process mechanism is analysed in depth, establish the part mathematical modeling of extraction process, based on model, temperature and pressure in extraction process in extraction kettle is controlled using Multivariable Inferential Control method, this method not only realizes the decoupling between extraction process temperature and pressure, and realize the perfect tracking under setting value disturbance and the not full remuneration under measurable disturbance, improve extraction process temperature and pressure control accuracy and the extraction yield of material, improve economic benefit.
Description
Technical field
The present invention relates to a kind of carbon dioxide supercritical extraction production process, more particularly to a kind of supercritical extract mistake
The control method of journey temperature and pressure.
Background technology
Supercritical carbon dioxide is widely used in shooting flow as a kind of fluid media (medium) of green wide material sources
In the research and application process of body abstraction technique.The critical pressure of pure carbon dioxide is 7.38Mpa, and critical-temperature is 31.06 DEG C,
Carbon dioxide in critical pressure and critical-temperature above state is referred to as supercritical carbon dioxide.Supercritical carbon dioxide extracts
The solvability to solute dependent on supercritical carbon dioxide fluid is taken, solubility of the solute in CO 2 fluid is faced with super
The density of boundary's carbon dioxide is relevant, and the density of supercritical carbon dioxide depends on its temperature and pressure.
Temperature and pressure is most important control parameter in carbon dioxide supercritical extraction production process, and temperature is to overcritical
Fluid solvent power influence is more complicated, under a certain pressure, the volatility increase of rise temperature extract, which adds
Concentration of the extract in overcritical gas phase, so that extraction quantity increases;But then, temperature raises, supercritical fluid
Density reduces, so that chemical constituent solubility reduces.Under a constant, it is close at a temperature of the height of extraction yield depends on this
The influence that the influence of degree or the volatility of solute increase to solubility is dominant.Influence of the pressure to extraction is also critically important, if
The pressure of carbon dioxide is too big, can cause carbon dioxide flow rate increase, carbon dioxide shorter residence time, with being extracted in extractor
Take thing time of contact to reduce, be unfavorable for the raising of extraction yield;But then, the pressure increase of carbon dioxide, it is possible to increase extraction
The mass transfer force of process, correspondingly increases mass tranfer coefficient, accelerates mass transfer rate, so as to improve supercritical extract ability.By
Carbon dioxide solubility ability can be just influenceed in the minor variations of the temperature and pressure near super critical point, and significant changes occur, this
Certain order of magnitude must be reached by requiring the control of temperature and pressure, otherwise can reduce extraction yield, influence to extract quality.
At present, the control for temperature and pressure during supercritical extract is most or uses simple PID control side
It method, can control very well for the PID controller of specific system design in some cases, but they still suffer from some and asked
Topic.A test signal, but this method meeting are inserted when control system is in closed-loop working state, it is necessary in control process
Cause disturbance, the influence that PID controller is interfered can produce overshoot.And during the supercritical extract of carbon dioxide in itself
A variety of survey is there is certain coupling between the two with immesurable interference, temperature and pressure to be present and interfere;
Influence of the vaporization or liquefaction absorption or liberated heat of carbon dioxide to temperature;Partial CO 2 fails to be in supercritical state
State and reduce quantity of solvent, reduce extraction yield, problems above is waited to need to solve.
Non-linear, close coupling and existing each is had based on temperature and pressure during above-mentioned carbon dioxide supercritical extraction
Kind disturbs immesurable problem, and the present invention is using the temperature and pressure of Multivariable Inferential Control method control extraction process, the party
Method not only realizes the decoupling between extraction process temperature and pressure, and realizes the perfect tracking and not under setting value disturbance
Full remuneration under measurable disturbance, therefore carry out supercritical extract process temperature and pressure essence on the basis of advanced control method
The research of quasi- control is significant.
The content of the invention
It is an object of the invention to provide a kind of carbon dioxide supercritical extraction process temperature and the control method of pressure, for
There is non-linear, close coupling and interference can not survey for carbon dioxide supercritical extraction process temperature and pressure, and to temperature and pressure
The characteristics of control accuracy requirement is high, on the basis of extraction process mechanism is analysed in depth, establish the part mathematics of extraction process
Model, based on model, the temperature and pressure in extraction process in extraction kettle, this method are controlled using Multivariable Inferential Control method
Not only realize the decoupling between extraction process temperature and pressure, and realize the lower perfect tracking of setting value disturbance with can not
The full remuneration surveyed under disturbance, improves extraction process temperature and pressure control accuracy and the extraction yield of material, improves economy
Benefit.
The method of the present invention is:
1st, it is controlled for the temperature and pressure of extraction kettle, it is accurately extracted first;According to technological requirement, survey
The temperature sensor of amount extraction temperature in the kettle is arranged in extraction kettle circulating hot water chuck;The pressure for measuring pressure in extraction kettle passes
Sensor is arranged on supercritical carbon dioxide fluid and entered in the pipeline of extraction kettle porch;
2nd, the heat extracted needed for the lifting of temperature in the kettle is by the cycling hot water extraction in extraction kettle circulating hot water chuck
Supply, and circulating hot water is then heated by heating tube in boiler, wherein the heating tube is electrothermal tube, then passes through thermal circulating pump
Circulating hot water in extraction kettle circulating hot water chuck is provided and circulated, heating tube controls break-make, root with solid-state relay
According to the real time data of temperature sensor measurement, the heat time of heating tube is automatically adjusted by controller;
Pressure is that carbon dioxide pressurization is controlled by Frequency Converter Control booster motor rotating speed in extraction kettle, when
After pressure sensor measurement obtains pressure data, these data are transported in controller, with pre-set pressure value ratio
After do analytic operation processing, frequency converter is controlled by computing output signal, and then control booster motor rotating speed, so as to reach
To the purpose of supercharging or decompression.
3rd, Multivariable Inferential Control System is automatically controlled by a computer, and Multivariable Inferential Control System is by Multivariable Inferential control
Make part, procedure division, procedure division mathematical modeling and disturbance unit packet can not be surveyed into Multivariable Inferential Control part therein
Using V canonical form reasoning and decision device matrixes Giv(s);Procedure division is by heating tube make-and-break time m1Influence to extraction kettle temperature
Relation passage, heating tube make-and-break time m1Influence relation passage, booster motor rotating speed m to extraction kettle pressure2To in extraction kettle
The influence relation passage and booster motor rotating speed m of temperature2The influence relation passage of pressure in extraction kettle is formed;Process portion
Point mathematical modeling is by heating tube make-and-break time m1Mathematical modeling, the heating tube break-make of influence relation passage to extraction kettle temperature
Time m1The mathematical modeling of influence relation passage to extraction kettle pressure, booster motor rotating speed m2Shadow to extracting temperature in the kettle
The mathematical modeling and booster motor rotating speed m of the relation of sound passage2To the mathematical modulo of the influence relation passage of pressure in extraction kettle
Type is formed;Disturbance part can not be surveyed by carbon dioxide vaporization absorption heat or liquefaction liberated heat d1With voltage ripple of power network pair
It is pressurized the influence d of revolution speed2Form.
Multivariable Inferential Control System is introduced into heating tube make-and-break time m in extraction process1To extraction kettle temperature and pressure
Influence relation passage and booster motor rotating speed m2The mathematical modeling transmission of influence relation passage to extraction kettle temperature and pressure
Jacobian matrix is:
Designing multivariable V canonical form reasoning and decision device matrixes is:
Obtained by Multivariable Inferential Control System,
After arrangement, have
Wherein,It can be obtained by Multivariable Inferential Control System and formula (2)
Y (s)=G (s) G-1 iv(s)e(s)+d(s) (3)
Wherein
It can be obtained by Multivariable Inferential Control System and formula (3):When built mathematical modeling entirely accurate, i.e.,AndWhen, in setting value r1And r2Lower output Y (the s)=r (s) of effect,
That is, y1=r1、y2=r2, realize the perfect tracking under setting value perturbation action;In process not measurable disturbance d1And d2Under effect,
Export Y (s)=0, i.e. y1=0, y2=0, realize the full remuneration that can not be surveyed under perturbation action;
Now, multivariable V canonical forms reasoning and decision device matrixI.e.
Therefore, for carbon dioxide supercritical extraction production process, controlled and extracted using Multivariable Inferential Control method
Temperature in the kettle and pressure are taken, not only realizes the decoupling of temperature and pressure during supercritical extract, and setting can be realized
Perfect tracking under value perturbation action and the full remuneration under perturbation action can not be surveyed.
The beneficial effects of the invention are as follows:
1st, the present invention has non-linear close coupling and interference not for carbon dioxide supercritical extraction process temperature and pressure
Measurable problem, using Multivariable Inferential Control method, realize the perfect tracking under setting value perturbation action and can not survey and disturb
The full remuneration under is acted, the control accuracy of extraction process temperature and pressure is improved, improves extraction yield.
2nd, using the computer automatic control system that process control and line solver are carried out based on mathematical modeling, life is improved
Efficiency is produced, extract quality is improved, reduces raw material and energy expenditure, improve economic benefit.
Brief description of the drawings
Fig. 1 is the composition frame chart of the Multivariable Inferential Control System of the present invention.
Embodiment
The method of the present invention is:
(1), it is controlled for extracting the temperature and pressure of axe, first has to accurately extract it.Will according to technique
Ask, the temperature sensor of measurement extraction temperature in the kettle is arranged in extraction axe circulating hot water chuck;Measure pressure in extraction kettle
Pressure sensor is arranged on supercritical carbon dioxide fluid and entered in the pipeline of extraction kettle porch.
(2) heat, extracted in axe needed for the lifting of temperature is by the circulating hot water in extraction kettle circulating hot water chuck
There is provided, and circulating hot water is then heated by heating tube in boiler, wherein the heating tube is electrothermal tube, then passes through thermodynamic cycle
What pump provided the circulating hot water in extraction kettle circulating hot water chuck and circulated, heating tube controls break-make with solid-state relay,
According to the real time data of temperature sensor measurement, the heat time of heating tube is automatically adjusted by controller.
Pressure is that carbon dioxide pressurization is controlled by Frequency Converter Control booster motor rotating speed in extraction kettle, when
After pressure sensor measurement obtains pressure data, these data are transported in controller, with pre-set pressure value ratio
After do analytic operation processing, frequency converter is controlled by computing output signal, and then control booster motor rotating speed, so as to reach
To the purpose of supercharging or decompression.
(3), Multivariable Inferential Control System is automatically controlled by a computer, and Multivariable Inferential Control System is pushed away by multivariable
Manage control section, procedure division, procedure division mathematical modeling and disturbance unit packet can not be surveyed into Multivariable Inferential Control therein
Part uses V canonical form reasoning and decision device matrixes Giv(s);Procedure division is by heating tube make-and-break time m1To extraction kettle temperature
Influence relation passage, heating tube make-and-break time m1Influence relation passage, booster motor rotating speed m to extraction kettle pressure2To extraction
The influence relation passage and booster motor rotating speed m of temperature in the kettle2The influence relation passage of pressure in extraction kettle is formed;Cross
Journey part mathematical modeling is by heating tube make-and-break time m1Mathematical modeling, the heating tube of influence relation passage to extraction kettle temperature
Make-and-break time m1The mathematical modeling of influence relation passage to extraction kettle pressure, booster motor rotating speed m2To extracting temperature in the kettle
Influence relation passage mathematical modeling and booster motor rotating speed m2To the number of the influence relation passage of pressure in extraction kettle
Model is learned to form;Disturbance part can not be surveyed by carbon dioxide vaporization absorption heat or liquefaction liberated heat d1With line voltage ripple
The dynamic influence d to revolution speed2Form.
System is introduced into heating tube make-and-break time m in extraction process1Influence relation passage to extraction kettle temperature and pressure and
Booster motor rotating speed m2The mathematical modeling transfer function matrix of influence relation passage to extraction kettle temperature and pressure is:
Designing multivariable V canonical form reasoning and decision device matrixes is:
Obtained by Multivariable Inferential Control System,
After arrangement, have
Wherein,It can be obtained by Multivariable Inferential Control System and formula (2)
Y (s)=G (s) G-1 iv(s)e(s)+d(s) (3)
Wherein
It can be obtained by Multivariable Inferential Control System and formula (3):When built mathematical modeling entirely accurate, i.e.,AndWhen, in setting value r1And r2Lower output Y (the s)=r (s) of effect,
That is, y1=r1、y2=r2, realize the perfect tracking under setting value perturbation action;In process not measurable disturbance d1And d2Under effect,
Export Y (s)=0, i.e. y1=0, y2=0, realize the full remuneration that can not be surveyed under perturbation action;
Now, multivariable V canonical forms reasoning and decision device matrixI.e.
Therefore, for carbon dioxide supercritical extraction production process, controlled and extracted using Multivariable Inferential Control method
Temperature in the kettle and pressure are taken, not only realizes the decoupling of temperature and pressure during supercritical extract, and setting can be realized
Perfect tracking under value perturbation action and the full remuneration under perturbation action can not be surveyed.
In Fig. 1:
Giv(s) --- multivariable V canonical form reasoning and decision device matrixes;
--- mathematical modeling transfer function matrix;
G (s) --- supercritical extract process matrix;
--- the make-and-break time m of heating tube1To the mathematical modulo for the influence relation passage for extracting temperature in the kettle
Type;
--- booster motor rotating speed m2To the mathematical modeling for the influence relation passage for extracting temperature in the kettle;
--- the make-and-break time m of heating tube1To the mathematical modulo of the influence relation passage of pressure in extraction kettle
Type;
--- booster motor rotating speed m2To the mathematical modulo of the influence relation passage of pressure in extraction kettle
Type;
--- the make-and-break time m of heating tube1Influence relation passage to extracting temperature in the kettle;
--- booster motor rotating speed m2Influence relation passage to extracting temperature in the kettle;
--- the make-and-break time m of heating tube1To the influence relation passage of pressure in extraction kettle;
--- booster motor rotating speed m2To the influence relation passage of pressure in extraction kettle;
r1--- Multivariable Inferential Control System desired temperature;
r2--- Multivariable Inferential Control System pressure set points;
e1--- Multivariable Inferential Control System temperature deviation value;
e2--- Multivariable Inferential Control System pressure divergence value;
y1--- actual temperature in extraction kettle;
y2--- actual pressure in extraction kettle;
d1--- carbon dioxide gasification absorbs or liquefaction liberated heat;
d2--- influence of the voltage ripple of power network to revolution speed;
m1--- heating tube make-and-break time;
m2--- booster motor rotating speed.
Claims (1)
1. the control method of a kind of supercritical extract process temperature and pressure, this method are:
(1), it is controlled for the temperature and pressure of extraction kettle, it is extracted first;According to technological requirement, measurement extraction
The temperature sensor of temperature in the kettle is taken to be arranged in extraction kettle circulating hot water chuck;Measure the pressure sensor of pressure in extraction kettle
Enter installed in supercritical carbon dioxide fluid in the pipeline of extraction kettle porch;
(2), extracting the heat needed for the lifting of temperature in the kettle is provided by the circulating hot water in extraction kettle circulating hot water chuck
, and circulating hot water is then heated by heating tube in boiler, wherein the heating tube is electrothermal tube, then is carried by thermal circulating pump
Circulating hot water for extraction kettle circulating hot water chuck is simultaneously circulated, and heating tube controls break-make with solid-state relay, according to temperature
The real time data of sensor measurement is spent, the heat time of heating tube is automatically adjusted by controller;
Pressure is that carbon dioxide pressurization is controlled by Frequency Converter Control booster motor rotating speed in extraction kettle, works as pressure
After sensor measurement obtains pressure data, these data are transported in controller, with pre-set pressure value relatively after
Analytic operation processing is done, frequency converter is controlled by computing output signal, and then controls booster motor rotating speed, is increased so as to reach
Pressure or the purpose of decompression;
(3), Multivariable Inferential Control System is automatically controlled by a computer, and Multivariable Inferential Control System is by Multivariable Inferential control
Make part, procedure division, procedure division mathematical modeling and disturbance unit packet can not be surveyed into Multivariable Inferential Control part therein
Using V canonical form reasoning and decision device matrixes Giv(s);Procedure division is by heating tube make-and-break time m1Influence to extraction kettle temperature
Relation passage, heating tube make-and-break time m1Influence relation passage, booster motor rotating speed m to extraction kettle pressure2To in extraction kettle
The influence relation passage and booster motor rotating speed m of temperature2The influence relation passage of pressure in extraction kettle is formed;Process portion
Point mathematical modeling is by heating tube make-and-break time m1Mathematical modeling, the heating tube break-make of influence relation passage to extraction kettle temperature
Time m1The mathematical modeling of influence relation passage to extraction kettle pressure, booster motor rotating speed m2Shadow to extracting temperature in the kettle
The mathematical modeling and booster motor rotating speed m of the relation of sound passage2To the mathematical modulo of the influence relation passage of pressure in extraction kettle
Type is formed;Disturbance part can not be surveyed by carbon dioxide vaporization absorption heat or liquefaction liberated heat d1With voltage ripple of power network pair
The influence d of booster motor rotating speed2Form;
Multivariable Inferential Control System is introduced into heating tube make-and-break time m in extraction process1Influence to extraction kettle temperature and pressure is closed
It is passage and booster motor rotating speed m2The mathematical modeling transmission function square of influence relation passage to extraction kettle temperature and pressure
Battle array be:
<mrow>
<mover>
<mi>G</mi>
<mo>^</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mn>11</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>^</mo>
</mover>
<mn>11</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>^</mo>
</mover>
<mn>12</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>^</mo>
</mover>
<mn>21</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mover>
<mi>g</mi>
<mo>^</mo>
</mover>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mover>
<mi>&tau;</mi>
<mo>^</mo>
</mover>
<mn>22</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Designing multivariable V canonical form reasoning and decision device matrixes is:
<mrow>
<msub>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>11</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>12</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>21</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>22</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
Obtained by Multivariable Inferential Control System,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>21</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>11</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>m</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>12</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<msub>
<mi>m</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>22</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
After arrangement, have
<mrow>
<mi>m</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>11</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>12</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>21</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<msup>
<mi>g</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>22</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<msubsup>
<mi>G</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
Wherein,It can be obtained by Multivariable Inferential Control System and formula (2)
Y (s)=G (s) G-1 iv(s)e(s)+d(s) (3)
Wherein
<mrow>
<mi>Y</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>y</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>d</mi>
<mn>2</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
It can be obtained by Multivariable Inferential Control System and formula (3):When built mathematical modeling entirely accurate, i.e.,And
AndWhen, in setting value r1And r2Lower output Y (the s)=r (s) of effect, i.e. y1=r1、y2=r2, it is real
The perfect tracking under setting value perturbation action is showed;In process not measurable disturbance d1And d2Under effect, Y (s)=0 is exported, i.e. y1
=0, y2=0, realize the full remuneration that can not be surveyed under perturbation action;
Now, multivariable V canonical forms reasoning and decision device matrixI.e.
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>11</mn>
</mrow>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<mo>=</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>g</mi>
<mn>11</mn>
</msub>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mn>11</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>12</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>g</mi>
<mn>12</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mn>12</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>21</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>g</mi>
<mn>21</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mn>21</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>g</mi>
<mrow>
<mi>i</mi>
<mi>v</mi>
<mn>22</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>/</mo>
<msub>
<mi>g</mi>
<mn>22</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>s</mi>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mn>22</mn>
</msub>
<mi>s</mi>
</mrow>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Therefore, for carbon dioxide supercritical extraction production process, extraction kettle is controlled using Multivariable Inferential Control method
Interior temperature and pressure, the decoupling of temperature and pressure during supercritical extract is not only realized, and can realize that setting value is disturbed
Act the perfect tracking under and the full remuneration under perturbation action can not be surveyed;
Wherein:
Giv(s) --- multivariable V canonical form reasoning and decision device matrixes;
--- mathematical modeling transfer function matrix;
G (s) --- supercritical extract process matrix;
--- the make-and-break time m of heating tube1To the mathematical modeling for the influence relation passage for extracting temperature in the kettle;
--- booster motor rotating speed m2To the mathematical modeling for the influence relation passage for extracting temperature in the kettle;
--- the make-and-break time m of heating tube1To the mathematical modeling of the influence relation passage of pressure in extraction kettle;
--- booster motor rotating speed m2To the mathematical modeling of the influence relation passage of pressure in extraction kettle;
--- the make-and-break time m of heating tube1Influence relation passage to extracting temperature in the kettle;
--- booster motor rotating speed m2Influence relation passage to extracting temperature in the kettle;
--- the make-and-break time m of heating tube1To the influence relation passage of pressure in extraction kettle;
--- booster motor rotating speed m2To the influence relation passage of pressure in extraction kettle;
r1--- Multivariable Inferential Control System desired temperature;
r2--- Multivariable Inferential Control System pressure set points;
e1--- Multivariable Inferential Control System temperature deviation value;
e2--- Multivariable Inferential Control System pressure divergence value;
y1--- actual temperature in extraction kettle;
y2--- actual pressure in extraction kettle;
d1--- carbon dioxide gasification absorbs or liquefaction liberated heat;
d2--- influence of the voltage ripple of power network to revolution speed;
m1--- heating tube make-and-break time;
m2--- booster motor rotating speed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510336215.5A CN104950960B (en) | 2015-06-17 | 2015-06-17 | A kind of control method of supercritical extract process temperature and pressure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510336215.5A CN104950960B (en) | 2015-06-17 | 2015-06-17 | A kind of control method of supercritical extract process temperature and pressure |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104950960A CN104950960A (en) | 2015-09-30 |
CN104950960B true CN104950960B (en) | 2017-12-29 |
Family
ID=54165677
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510336215.5A Expired - Fee Related CN104950960B (en) | 2015-06-17 | 2015-06-17 | A kind of control method of supercritical extract process temperature and pressure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104950960B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108803329B (en) * | 2018-06-19 | 2021-06-18 | 长春工业大学 | Solubility optimization method in supercritical extraction process |
CN112057891A (en) * | 2020-09-22 | 2020-12-11 | 南通市华安超临界萃取有限公司 | Self-adaptive control device for supercritical extraction separation device |
CN117225004B (en) * | 2023-11-14 | 2024-01-26 | 湖南禧又多油脂有限公司 | Automatic control system for vegetable oil low-temperature extraction treatment based on Internet of things |
CN117438123A (en) * | 2023-12-22 | 2024-01-23 | 中国核电工程有限公司 | Method and device for determining uranium extraction zone in extraction column and method for controlling position of uranium extraction zone |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4904514B2 (en) * | 2006-05-19 | 2012-03-28 | 独立行政法人産業技術総合研究所 | Highly efficient separation method and apparatus for solid-liquid mixture with supercritical carbon dioxide and / or liquid carbon dioxide |
CN100524140C (en) * | 2007-08-28 | 2009-08-05 | 中昊晨光化工研究院 | Polymer polymerization reaction temperature control method and device thereof |
CN101391156A (en) * | 2008-10-22 | 2009-03-25 | 南方医科大学珠江医院 | Fine particle preparation device using supercritical fluid quick expansion |
CN101642632B (en) * | 2009-09-03 | 2012-05-30 | 江苏大学 | Subcritical liquid extraction device and subcritical liquid extraction method of active ingredients of natural products |
RU2458053C1 (en) * | 2011-06-06 | 2012-08-10 | Открытое акционерное общество "КуйбышевАзот" | Method of controlling extraction of caprolactam |
CN103869696A (en) * | 2014-03-19 | 2014-06-18 | 长春工业大学 | Distillation parameter decoupling control method |
-
2015
- 2015-06-17 CN CN201510336215.5A patent/CN104950960B/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN104950960A (en) | 2015-09-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104950960B (en) | A kind of control method of supercritical extract process temperature and pressure | |
CN104297072B (en) | A kind of supercritical carbon dioxide closed cycle compression property experimental provision and method | |
CN105952424B (en) | A kind of supercritical water displacement of reservoir oil simulator and method | |
CN104934082B (en) | A kind of controlled ocean thermal environment analog systemss of temperature and pressure | |
CN105807632A (en) | Heating furnace temperature controller design method based on improved Dahlin algorithm | |
CN113656986A (en) | Method for rapidly calculating long-term operation heat exchange performance of intermediate-deep geothermal buried pipe | |
CN112307626B (en) | Operation and influence analysis method for multiple groups of passive heat export systems in containment | |
CN103776502B (en) | Fired power generating unit mesolow cylinder entrance reheat heat steam mass flow real time measure method | |
CN104888396A (en) | Method for detecting transient faults of fire pump in quick starting process | |
CN108086960A (en) | A kind of water erosion method exploitation of gas hydrates experimental simulation method and apparatus | |
CN104013089B (en) | A kind of pipe tobacco hot blast alternating temperature drying device and alternating temperature drying test method thereof | |
Vaganova et al. | Simulation of optimal exploitation of an open geotermal loop | |
CN110414174B (en) | Modeling method and system for steady-state characteristics of natural circulation of supercritical fluid | |
CN103697958B (en) | The real time measure method of coal unit drum outlet saturation steam mass rate | |
CN206111158U (en) | Linear physical simulation measurement and control device of viscous crude oil reservoir thermal recovery | |
CN106934134A (en) | A kind of method for characterizing the centrifugation unsteady intensity of pump cavitation | |
CN106769528A (en) | A kind of temperature control system and its temperature control method and tensile strength test device | |
CN103838961A (en) | Method for monitoring three-dimensional temperature and thermal stress of ultra-supercritical steam turbine rotor in real time | |
CN107451304B (en) | Mechanism modeling calculation method of reheater transfer function model | |
CN104806533B (en) | The hot endurance quality pilot system of a kind of nuclear power station main feed pump | |
CN109115532B (en) | Performance testing device and method for middle-deep concentric sleeve type heat exchanger | |
Wei et al. | Boiler drum water level control based on nonlinear valve characteristic compensation | |
CN105956329A (en) | Calculation method for mechanism modeling of each channel gain of heat exchanger | |
CN105402712B (en) | One kind is instantaneously without pressure steam generator | |
Wolszczak et al. | The control system of the yeast drying process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171229 Termination date: 20180617 |