发明名称 一种超临界萃取过程温度和压力的精准控制方法
摘要 本发明公开了一种超临界萃取过程温度和压力的精准控制方法,该方法针对二氧化碳超临界萃取过程温度和压力具有非线性、强耦合及干扰不可测,且对温度和压力控制精度要求高的特点,在深入分析萃取过程机理的基础上,建立了萃取过程的部分数学模型,基于模型,采用多变量推理控制方法控制萃取过程中萃取釜内的温度和压力,该方法不但实现了萃取过程温度和压力之间的解耦,而且实现了设定值扰动下的完全跟踪和不可测扰动下的完全补偿,提高了萃取过程温度和压力控制精度和物质的萃取率,提高了经济效益。
申请公布号 CN104950960A 申请公布日期 2015.09.30
申请号 CN201510336215.5 申请日期 2015.06.17
申请人 长春工业大学 发明人 马海涛;徐生欣;尤元;李丙林
分类号 G05D27/02(2006.01)I 主分类号 G05D27/02(2006.01)I
代理机构 长春市四环专利事务所(普通合伙) 22103 代理人 张建成
主权项 一种超临界萃取过程温度和压力的精准控制方法,该方法是:(一)、针对萃取釜的温度和压力进行控制,首先对其进行精确提取;根据工艺要求,测量萃取釜内温度的传感器安装在萃取釜循环热水夹套内;测量萃取釜内压力的传感器安装在超临界二氧化碳流体进入萃取釜入口处管道内;(二)、萃取斧内温度的升降所需的热量是通过萃取釜夹套里的循环热水提供的,而循环热水则由热水箱里电热管加热,再通过热力循环泵提供萃取釜夹套里的热水并进行循环的,加热管用固态继电器控制通断,根据温度传感器测量的实时数据,通过控制器对加热管的加热时间进行自动调节;萃取釜内压力是通过变频器控制柱塞泵转速对二氧化碳加压来进行控制,当压力传感器测量得到压力数据后,将这些数据输送到控制器中,与预先设定好的压力值比较后做分析运算处理,经过运算输出信号来控制变频器,进而控制柱塞泵电机转速,从而达到增压或减压的目的;(三)、整个系统由计算机自动控制,系统由多变量推理控制部分、过程部分、过程部分数学模型及不可测扰动部分组成,其中的多变量推理控制部分采用V规范型推理控制器G<sub>iv</sub>(s);过程部分是由加热管通断时间m<sub>1</sub>对萃取釜温度的影响关系通道、加热管通断时间m<sub>1</sub>对萃取釜压力的影响关系通道、泵转速m<sub>2</sub>对萃取釜内温度的影响关系通道以及泵转速m<sub>2</sub>对萃取釜内压力的影响关系通道构成;过程部分数学模型是由加热管通断时间m<sub>1</sub>对萃取釜温度的影响关系通道的数学模型、加热管通断时间m<sub>1</sub>对萃取釜压力的影响关系通道的数学模型、泵转速m<sub>2</sub>对萃取釜内温度的影响关系通道的数学模型以及泵转速m<sub>2</sub>对萃取釜内压力的影响关系通道的数学模型构成;不可测扰动部分由二氧化碳汽化吸收热量或液化放出的热量d<sub>1</sub>和电网电压波动对泵转速的影响d<sub>2</sub>构成。系统引入萃取过程中加热管通断时间m<sub>1</sub>对萃取釜温度和压力的影响关系通道及增压泵转速m<sub>2</sub>对萃取釜温度和压力的影响关系通道的数学模型传递函数矩阵为:<maths num="0001" id="cmaths0001"><math><![CDATA[<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>]]></math><img file="FDA0000739868220000021.GIF" wi="932" he="196" /></maths>设计多变量V规范型推理控制器为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>i</mi><mover><mi>&nu;</mi><mo>&OverBar;</mo></mover></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>]]></math><img file="FDA0000739868220000022.GIF" wi="788" he="182" /></maths>由图1得,<maths num="0003" id="cmaths0003"><math><![CDATA[<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></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><msub><mi>m</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo><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></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><msub><mi>m</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>)</mo><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>]]></math><img file="FDA0000739868220000023.GIF" wi="1846" he="226" /></maths>整理后,有<maths num="0004" id="cmaths0004"><math><![CDATA[<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><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><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>]]></math><img file="FDA0000739868220000024.GIF" wi="1926" he="226" /></maths>由图1及式(2)可得Y(s)=G(s)G<sup>‑1</sup><sub>iv</sub>(s)e(s)+d(s)   (3)其中<maths num="0005" id="cmaths0005"><math><![CDATA[<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>]]></math><img file="FDA0000739868220000025.GIF" wi="445" he="204" /></maths><maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>m</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>]]></math><img file="FDA0000739868220000031.GIF" wi="443" he="204" /></maths><maths num="0007" id="cmaths0007"><math><![CDATA[<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>]]></math><img file="FDA0000739868220000032.GIF" wi="461" he="210" /></maths>由图1和式(3)可得:当所建数学模型完全准确,即<img file="FDA0000739868220000033.GIF" wi="346" he="94" />并且<img file="FDA0000739868220000034.GIF" wi="508" he="90" />时,在设定值r<sub>1</sub>和r<sub>2</sub>作用下输出Y(s)=r(s),即,y<sub>1</sub>=r<sub>1</sub>、y<sub>2</sub>=r<sub>2</sub>,实现了设定值扰动作用下的完全跟踪;在过程不可测扰动d<sub>1</sub>和d<sub>2</sub>作用下,输出Y(s)=0,即,y<sub>1</sub>=0,y<sub>2</sub>=0,实现了不可测扰动作用下的完全补偿;此时,推理控制器<maths num="0008" id="cmaths0008"><math><![CDATA[<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><mover><mi>G</mi><mo>^</mo></mover><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000739868220000035.GIF" wi="411" he="86" /></maths>即<maths num="0009" id="cmaths0009"><math><![CDATA[<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>]]></math><img file="FDA0000739868220000036.GIF" wi="655" he="450" /></maths>因此,针对二氧化碳超临界萃取生产工艺过程,采用多变量推理控制方法控制萃取釜内温度和压力,不但实现了超临界萃取过程中温度和压力的解耦,而且可以实现设定值扰动作用下的完全跟踪和不可测扰动作用下的完全补偿;其中:<img file="FDA0000739868220000039.GIF" wi="129" he="64" />——多变量V规范型推理控制器矩阵(2*2矩阵);<img file="FDA0000739868220000037.GIF" wi="111" he="75" />——超临界萃取过程数学模型矩阵(2*2矩阵);G(s)——超临界萃取过程矩阵(2*2矩阵);<img file="FDA0000739868220000038.GIF" wi="188" he="86" />——加热管的通断时间m<sub>1</sub>对萃取釜内温度的影响关系通道数学模型;<img file="FDA0000739868220000041.GIF" wi="192" he="86" />——增压泵泵转速m<sub>2</sub>对萃取釜内温度的影响关系通道数学模型;<img file="FDA0000739868220000042.GIF" wi="188" he="83" />——加热管的通断时间m<sub>1</sub>对萃取釜内压力的影响关系通道数学模型;<img file="FDA0000739868220000043.GIF" wi="194" he="79" />——增压泵转速m<sub>2</sub>对萃取釜内压力的影响关系通道数学模型;<img file="FDA0000739868220000044.GIF" wi="195" he="80" />——加热管的通断时间m<sub>1</sub>对萃取釜内温度的影响关系通道;<img file="FDA0000739868220000045.GIF" wi="198" he="80" />——增压泵泵转速m<sub>2</sub>对萃取釜内温度的影响关系通道<img file="FDA0000739868220000046.GIF" wi="200" he="84" />——加热管的通断时间m<sub>1</sub>对萃取釜内压力的影响关系通道;<img file="FDA0000739868220000047.GIF" wi="199" he="82" />——增压泵转速m<sub>2</sub>对萃取釜内压力的影响关系通道r<sub>1</sub>——系统温度设定值;r<sub>2</sub>——系统压力设定值;e<sub>1</sub>——系统温度偏差值;e<sub>2</sub>——系统压力偏差值;y<sub>1</sub>——萃取釜内实际温度;y<sub>2</sub>——萃取釜内实际压力;d<sub>1</sub>——二氧化碳气化吸收或液化放出的热量;d<sub>2</sub>——电网电压波动对泵转速的影响;m<sub>1</sub>——加热管通断时间;m<sub>2</sub>——增压泵转速。
地址 130012 吉林省长春市延安大街2055号
您可能感兴趣的专利