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 PDF

Info

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
Application number
CN201510336215.5A
Other languages
Chinese (zh)
Other versions
CN104950960A (en
Inventor
马海涛
徐生欣
尤元
李丙林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changchun University of Technology
Original Assignee
Changchun University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changchun University of Technology filed Critical Changchun University of Technology
Priority to CN201510336215.5A priority Critical patent/CN104950960B/en
Publication of CN104950960A publication Critical patent/CN104950960A/en
Application granted granted Critical
Publication of CN104950960B publication Critical patent/CN104950960B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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

A kind of control method of supercritical extract process temperature and pressure
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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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>&amp;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.
CN201510336215.5A 2015-06-17 2015-06-17 A kind of control method of supercritical extract process temperature and pressure Expired - Fee Related CN104950960B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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