CN109634105B - 一种面向加工方案切换的原油动态调合方法 - Google Patents
一种面向加工方案切换的原油动态调合方法 Download PDFInfo
- Publication number
- CN109634105B CN109634105B CN201811509467.3A CN201811509467A CN109634105B CN 109634105 B CN109634105 B CN 109634105B CN 201811509467 A CN201811509467 A CN 201811509467A CN 109634105 B CN109634105 B CN 109634105B
- Authority
- CN
- China
- Prior art keywords
- crude oil
- blending
- period
- buffer tank
- switching
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/0205—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
- G05B13/024—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Lubricants (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明提出一种面向加工方案切换的原油动态调合方法,首先利用原油近红外快速评价技术获取切换前加工原油和待切换原油的性质数据,然后基于原油数据以及工艺参数,以设定的调合周期对调合比进行优化,并以一定的跟踪周期进行滚动计算缓冲罐内混合原油的实时性质值,动态调整下一调合周期的目标性质值,进而完成原油切换。该方法可以解决一类具有缓冲罐的复杂原油加工工艺,在原油加工方案切换下的调合问题,可有效避免方案切换时,因人工操作对常减压装置生产稳定性的影响,实现装置自动平稳过渡,进而稳定产品质量。
Description
技术领域
本发明涉及石油化工领域的原油调合方面,具体是一种实现原油加工方案切换的动态调合技术。
背景技术
受市场和效益影响,原油加工方案切换是炼化企业的常见操作,即采用同一套常减压装置来加工不同产品,例如,炼化企业常采用同一套常减压装置进行焦化料及沥青料生产,这就要求常减压装置定期或不定期地在不同原油加工方案之间来回切换。
现有的方案切换,仍多为人工操作,受经验影响,不同员工操作的切换效果(过渡过程时间、波动剧烈程度等)有较大差异,特别是在常减压装置前还具有缓冲罐这一类复杂工艺时,原油切换难度更大。
现有的原油调合技术均面向正常加工情况,其本质上还属于一种静态或稳态调合。可见,企业还需要一种加工方案切换下的原油自动动态调合技术,保证常减压装置加工方案实现平稳过渡。
发明内容
本发明提出一种面向加工方案切换的原油动态调合方法:首先利用原油近红外快速评价技术获取切换前加工原油和待切换原油的实时性质数据,然后基于原油数据以及工艺参数,在每个调合周期优化原油配方,并以一定的跟踪周期进行滚动计算,动态调整下一调合周期的目标值,进而完成原油切换。
技术方案:
一种面向加工方案切换的原油动态调合方法,采用动态调合方法完成原油加工方案切换,具有如下步骤:
1)获取切换前加工原油和待切换原油的实时性质数据以及各组分油占比,所述切换前加工原油是加工方案切换前用于调合的组分油,所述待切换原油是指加工方案切换后用于调合的组分油;
2)初始化参数,包括原油切换操作开始时间、跟踪周期、调合周期;一次原油切换操作开始到结束的阶段分为多个调合周期,一个调合周期分多个跟踪周期进行;
3)基于设定的目标性质值以及工艺约束条件,优化计算当前调合周期内,进入缓冲罐中各组分原油的调合比;
4)根据组分油的调合比,控制调合头,输出当前调合周期的调合流量;
5)由组分油配比,按照跟踪周期滚动计算缓冲罐内当前的混合原油性质值,单次调合周期中最后一个跟踪周期结束时缓冲罐内混合原油的性质值pm即为跟踪计算得到的性质值pk;
7)判断混合原油性质是否到达待切换原油性质目标,是则完成原油切换过程,否则返回步骤3)。
具体的,步骤2)中设定有跟踪周期与调合周期,调合周期内假定缓冲罐输出的原油性质等于此调合周期前缓冲罐内的原油性质值,再以跟踪周期滚动计算当前罐内混合原油的实际性质值,跟踪周期t与调合周期T的关系如下:
t=T/u (1)
其中,u为一个调合周期内的跟踪周期个数,1≤u≤3,u为整数。
具体的,步骤3)中以缓冲罐中的原油性质数据偏差最小为目标,建立目标函数和约束条件,优化求解得到该调合周期的各组分油占比:
其中,Ψ(Xk)为第k个调合周期的目标函数,k=1,2...N,N为整个调合过程需要的总周期个数;i表示参与调合原油的组分个数,i=1,2,...,n;λj表示原油第j个性质的权重,j表示组分油性质个数,j=1,2...,m;Xk为第k个调合周期的组分油占比,Xk=[xk1,xk2…xkn]T,Xk≥0,xki表示第i种组分油在第k个调合周期中的占比;ai为第i种组分油的性质数据;dk为第k个调合周期开始前实际测得的缓冲罐中原油性质数据,dk≥0;A为缓冲罐的截面积,l为缓冲罐内原油液位高度;k=1时,为加工方案切换后混合原油的性质值;k≥2时,为动态调整后第k个调合周期的目标性质值;F为原油流量向量,组分原油流量fi≥0;M为原油管线的最大流速,M>0。
具体的,缓冲罐内当前的混合原油性质值的计算如式(6)所示:
其中,为第v个跟踪周期进入缓冲罐的混合原油质量占比,为第v个跟踪周期开始时混合原油在第v个跟踪周期中质量占比,pv为第v个跟踪周期结束后缓冲罐内混合原油的性质值,p为当前调合周期内进入缓冲罐内混合原油的性质值,pv-1为第v-1个跟踪周期结束后缓冲罐内混合原油的性质值,v=1,2,...,m;
缓冲罐内当前的混合原油的非线性调合性质密度采用式(7)计算:
其中,为第v个跟踪周期进入缓冲罐的混合原油质量占比,为第v个跟踪周期开始时混合原油在第v个跟踪周期中质量占比,ρv为第v个跟踪周期结束后缓冲罐内混合原油的非线性值,ρ为当前调合周期内进入缓冲罐内混合原油的密度值,ρv-1为第v-1个跟踪周期结束后缓冲罐内混合原油的密度值,v=1,2,...,m。
具体的,缓冲罐内当前的混合原油质量调合比通过下式计算:
式中,F为原油流量,A为缓冲罐的截面积,l为缓冲罐内原油液位高度,ρ为当前调合周期内进入缓冲罐内混合原油的密度值,ρv-1为第v-1个跟踪周期结束后缓冲罐内混合原油的密度值,t为跟踪周期。
具体的,该方法通过跟踪周期滚动计算出第k个调合周期结束时,缓冲罐内混合原油的实际性质值pk,计算得到混合原油实际性质值pk与第k个调合周期目标性质值的偏差,再对第k+1个调合周期目标性质值进行动态调整,具体计算式如下:
本方法中切换前加工原油和待切换原油的性质数据是通过近红外快速评价技术得到。
本方法原油调合环节中存在线性与非线性调合关系的性质,在动态调合的过程中,具有线性调合关系的原油性质主要包括酸值、硫含量和密度,具有非线性调合关系的原油性质主要包括原油密度。
有益效果:
本发明提出的面向加工方案切换的原油动态调合方法,在加工方案切换状态下,基于当前加工原油和待切换原油的性质数据,通过对调合周期对内调合比的优化,以及调合周期内罐内混合原油性质的滚动计算,动态调整下一调合周期的目标性质值,可有效避免方案切换时,因人工操作对常减压装置生产稳定性的影响,实现装置自动平稳过渡,进而稳定产品质量。
附图说明
图1原油加工方案切换下的动态调合流程图
具体实施案例
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述实施例。
以某炼化企业为例,在正常操作时,该企业首先依据加工方案以及各组分原油数据,计算出最佳的原油调合配方,然后将各组分原油输送至原油静态混合器进行调合,并将调合后的原油经缓冲罐送至常减压装置进行加工。该炼化企业在2018年4月3日进行方案切换,从焦化料生产方案更改为沥青料生产方案。方案切换前调合组分原油分别为伊轻和丹佛,在7时起开始进行原油切换操作,调合组分原油改为伊重和玛雅。
图1为原油加工方案切换下的动态调合流程图。下面以原油加工过程中最为重要的密度、硫含量和酸值3个性质为例,给出原油加工方案切换下动态调合的详细过程。具体的实施步骤如下:
1)获取切换前加工原油和待切换原油的性质数据
当前加工原油为伊轻和丹佛,这2种组分油占比均为50%。经过调合头调合后得到当前加工的混合原油,其中密度ρstart=816.960kg/m3,硫含量Sstart=0.741Wt%,酸值Cstart=0.079mgKOH/g。
待切换原油为伊重和玛雅,整个切换过程结束时,最终混合原油中这2种组分油所占比例分别为85%和15%,经过调合头调合后的混合原油性质为:密度ρfinal=858.972kg/m3,硫含量Sfinal=2.088Wt%,酸值Cfinal=0.202mgKOH/g。
2)初始化参数
原油开始切换时间为2018年4月3日7时00分,此时循环迭代标识符k=0,混合原油流量F=744.18m3/h,缓冲罐液位l=12.4m,其截面积A=400m2,调合周期T=30min,跟踪周期t=15min。
3)优化当前调合周期组分油的调合比
基于调合的目标值以及原油的性质数据,计算当前调合周期内,待切换原油中各组分油的调合比。组分原油的密度、硫含量以及酸值数据如表1所示。
表1组分原油的性质数据
伊重 | 伊轻 | 丹佛 | 玛雅 | |
密度(kg/m3) | 874.400 | 853.700 | 787.700 | 921.500 |
硫含量(Wt%) | 1.900 | 1.360 | 0.140 | 3.410 |
酸值(mgKOH/g) | 0.200 | 0.130 | 0.020 | 0.180 |
k=0时,表示切换操作开始前,此时缓冲罐内为50%伊轻和50%丹佛的混合原油,当前缓冲罐内的混合原油的密度ρ0=ρstart=816.960kg/m3,硫含量S0=Sstart=0.741Wt%,酸值C0=Cstart=0.079mgKOH/g。当k=1时,表示处于第一个调合周期,第一个调合周期的目标值取切换操作结束时的性质,即密度ρ′1=858.972kg/m3,硫含量S′1=2.088Wt%,酸值C′1=0.202mgKOH/g。根据式(2)计算当前调合周期内原油的调合比。
XF≤M
计算出当前周期内待切换组分原油伊重和玛雅的调合比为0%:100%。该比例对应的组分原油性质分别为:密度ρ1=921.500kg/m3,硫含量S1=3.410Wt%,酸值C1=0.180mgKOH/g。
4)根据组分油的调合比,控制调合头,输出当前调合周期的调合流量;
5)计算缓冲罐内当前的混合原油性质
由式(4)和式(5)得到第v个跟踪周期进入缓冲罐的混合原油质量比:
第v-1个跟踪周期结束后的混合原油,在第v个跟踪周期中,混合原油质量占比:
依据计算出的调合比以及式(6)以及式(7)可以计算得跟踪周期内的原油密度,硫含量以及酸值,计算结果如表2所示。
表2跟踪周期内调合比以及性质数据
6)动态调整下一个调合周期的目标值
基于当前加工原油性质值和待加工原油性质值的偏差,依据式(3)对下一个调合周期的目标值进行动态调整,其中修正系数取w=0.1,下一个调合周期的目标值分别为:
密度ρ′2=858.972+0.1*(824.655-858.972)=855.540kg/m3,
硫含量S′2=2.088+0.1*(0.960-2.088)=1.975Wt%,
酸值C′2=0.202+0.1*(0.087-0.202)=0.191mgKOH/g。
7)判断混合原油性质是否到达待切换原油性质目标
经过计算,在2018年4月3日14:00时得到的缓冲罐内混合原油性质已达到待切换原油性质目标,计算结束。缓冲罐内跟踪周期的混合原油的密度、硫含量以及酸值数据如表3所示。
表3缓冲罐内的混合原油的酸值、密度、倾点数据
由表3可见,该方法可以解决一类具有缓冲罐的复杂原油加工工艺,在原油加工方案切换下的调合问题,可有效避免方案切换时,因人工操作对常减压装置生产稳定性的影响,实现装置自动平稳过渡,进而稳定产品质量。
Claims (5)
1.一种面向加工方案切换的原油动态调合方法,其特征在于采用动态调合方法完成原油加工方案切换,具有如下步骤:
1)获取切换前加工原油和待切换原油的实时性质数据以及各组分油占比,所述切换前加工原油是加工方案切换前用于调合的组分油,所述待切换原油是指加工方案切换后用于调合的组分油;
2)初始化参数,包括原油切换操作开始时间、跟踪周期、调合周期;一次原油切换操作开始到结束的阶段分为多个调合周期,一个调合周期分多个跟踪周期进行;
3)基于设定的目标性质值以及工艺约束条件,优化计算当前调合周期内,进入缓冲罐中各组分原油的调合比;
4)根据组分油的调合比,控制调合头,输出当前调合周期的调合流量;
5)由组分油配比,按照跟踪周期滚动计算缓冲罐内当前的混合原油实际性质值pk;
通过跟踪周期滚动计算出第k个调合周期结束时缓冲罐内混合原油的实际性质值pk,计算得到混合原油实际性质值pk与第k个调合周期目标性质值的偏差,再对第k+1个调合周期目标性质值进行动态调整,具体计算式如下:
7)判断混合原油性质是否到达待切换原油性质目标,是则完成原油切换过程,否则返回步骤3)。
2.根据权利要求1所述的一种面向加工方案切换的原油动态调合方法,其特征在于步骤2)中设定有跟踪周期与调合周期,跟踪周期t与调合周期T的关系如下:
t=T/u (1)
其中,u为一个调合周期内的跟踪周期个数,1≤u≤3,u为整数。
3.根据权利要求1所述的一种面向加工方案切换的原油动态调合方法,其特征在于步骤3)中以缓冲罐中的原油性质数据偏差最小为目标,建立目标函数和约束条件,优化求解得到该调合周期的各组分油占比:
其中,Ψ(Xk)为第k个调合周期的目标函数,k=1,2...N,N为整个调合过程需要的总周期个数;i表示参与调合原油的组分个数,i=1,2,...,n;λj表示原油第j个性质的权重,j表示组分油性质个数,j=1,2...,m;Xk为第k个调合周期的组分油占比,Xk=[xk1,xk2…xkn]T,Xk≥0,xki表示第i种组分油在第k个调合周期中的占比;ai为第i种组分油的性质数据;dk为第k个调合周期开始前实际测得的缓冲罐中原油性质数据,dk≥0;A为缓冲罐的截面积,l为缓冲罐内原油液位高度;N=1时,为加工方案切换后混合原油的性质值;N≥2时,为动态调整后第k个调合周期的目标性质值;F为原油流量向量,组分原油流量fi≥0;M为原油管线的最大流速,M>0。
4.根据权利要求1所述的一种面向加工方案切换的原油动态调合方法,其特征在于:
缓冲罐内当前的混合原油性质值的计算如式(6)所示:
其中,为第v个跟踪周期进入缓冲罐的混合原油质量占比,为第v个跟踪周期开始时混合原油在第v个跟踪周期中质量占比,pv为第v个跟踪周期结束后缓冲罐内混合原油的性质值,p为当前调合周期内进入缓冲罐内混合原油的性质值,pv-1为第v-1个跟踪周期结束后缓冲罐内混合原油的性质值,v=1,2,...,m;
缓冲罐内当前的混合原油的非线性调合性质密度采用式(7)计算:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509467.3A CN109634105B (zh) | 2018-12-11 | 2018-12-11 | 一种面向加工方案切换的原油动态调合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811509467.3A CN109634105B (zh) | 2018-12-11 | 2018-12-11 | 一种面向加工方案切换的原油动态调合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109634105A CN109634105A (zh) | 2019-04-16 |
CN109634105B true CN109634105B (zh) | 2022-04-08 |
Family
ID=66072541
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811509467.3A Active CN109634105B (zh) | 2018-12-11 | 2018-12-11 | 一种面向加工方案切换的原油动态调合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109634105B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110726830B (zh) * | 2019-10-23 | 2022-03-04 | 汉谷云智(武汉)科技有限公司 | 一种多进多出混合油性质预测方法及装置 |
CN110703713B (zh) * | 2019-11-05 | 2022-05-06 | 青岛大学 | 一种提高单装置多产品加工方案切换效率的方法 |
CN110825045B (zh) * | 2019-11-05 | 2022-03-22 | 青岛大学 | 一种提高多装置多产品加工方案切换效率的方法 |
CN111120873B (zh) * | 2019-12-06 | 2021-11-30 | 浙江中控技术股份有限公司 | 用于油品在线调合的管道补偿方法 |
CN112947325B (zh) * | 2021-01-28 | 2022-11-08 | 南京富岛信息工程有限公司 | 一种面向原油调合的储罐动态调度方法 |
CN113110336B (zh) * | 2021-04-20 | 2022-07-15 | 南京富岛信息工程有限公司 | 一种考虑调度约束的原油动态调合方法 |
CN116343933B (zh) * | 2022-11-22 | 2023-09-01 | 中石油云南石化有限公司 | 对原油油段性质进行跟踪、控制并重定义的方法及用途 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102375411A (zh) * | 2010-08-21 | 2012-03-14 | 中国石油化工股份有限公司 | 一种油品在线调合方法及装置 |
CN102393636A (zh) * | 2011-09-22 | 2012-03-28 | 任季明 | 一种对石油炼制过程的控制系统及方法 |
CN102591369A (zh) * | 2012-02-15 | 2012-07-18 | 南京富岛信息工程有限公司 | 一种原油调合自动控制方法 |
CN102643662A (zh) * | 2012-04-25 | 2012-08-22 | 南京富岛信息工程有限公司 | 一种原油调合优化方法 |
CN102768513A (zh) * | 2012-07-02 | 2012-11-07 | 清华大学 | 基于智能决策的炼油生产过程调度优化方法 |
CN103324084A (zh) * | 2013-07-11 | 2013-09-25 | 南京富岛信息工程有限公司 | 一种面向沥青质控制的原油在线调合方法 |
CN103365208A (zh) * | 2013-07-12 | 2013-10-23 | 南京富岛信息工程有限公司 | 一种具有底油性质偏差修正的原油调合优化方法 |
CN104298198A (zh) * | 2014-09-26 | 2015-01-21 | 南京富岛信息工程有限公司 | 一种原油调合优化控制系统中的可靠通信方法 |
CN106124448A (zh) * | 2016-06-12 | 2016-11-16 | 南京富岛信息工程有限公司 | 一种原油切换状态下的常减压装置进料性质预测方法 |
-
2018
- 2018-12-11 CN CN201811509467.3A patent/CN109634105B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102375411A (zh) * | 2010-08-21 | 2012-03-14 | 中国石油化工股份有限公司 | 一种油品在线调合方法及装置 |
CN102393636A (zh) * | 2011-09-22 | 2012-03-28 | 任季明 | 一种对石油炼制过程的控制系统及方法 |
CN102591369A (zh) * | 2012-02-15 | 2012-07-18 | 南京富岛信息工程有限公司 | 一种原油调合自动控制方法 |
CN102643662A (zh) * | 2012-04-25 | 2012-08-22 | 南京富岛信息工程有限公司 | 一种原油调合优化方法 |
CN102768513A (zh) * | 2012-07-02 | 2012-11-07 | 清华大学 | 基于智能决策的炼油生产过程调度优化方法 |
CN103324084A (zh) * | 2013-07-11 | 2013-09-25 | 南京富岛信息工程有限公司 | 一种面向沥青质控制的原油在线调合方法 |
CN103365208A (zh) * | 2013-07-12 | 2013-10-23 | 南京富岛信息工程有限公司 | 一种具有底油性质偏差修正的原油调合优化方法 |
CN104298198A (zh) * | 2014-09-26 | 2015-01-21 | 南京富岛信息工程有限公司 | 一种原油调合优化控制系统中的可靠通信方法 |
CN106124448A (zh) * | 2016-06-12 | 2016-11-16 | 南京富岛信息工程有限公司 | 一种原油切换状态下的常减压装置进料性质预测方法 |
Non-Patent Citations (4)
Title |
---|
A compound rapid assay method for crude oil properties based on similarity discrimination and partial least squares;Xisong Chen 等;《2016 35th Chinese Control Conference (CCC)》;20161231;全文 * |
原油调合优化及控制系统设计与实现;张锐;《东南大学硕士学位论文》;20150520;全文 * |
原油调合优化控制系统设计与实现;周家庆 等;《炼油与化工》;20131231;第24卷(第3期);正文第1.1-1.2节,图1和图2 * |
基于粒子群算法实现原油加工配方优化;蒋立沫 等;《炼油与化工》;20180507;第29卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109634105A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109634105B (zh) | 一种面向加工方案切换的原油动态调合方法 | |
CN108333943B (zh) | 基于增量模式的原油调合在线优化方法 | |
WO2018218870A1 (zh) | 一种伏特加酒智能全自动勾兑系统 | |
CN102643662B (zh) | 一种原油调合优化方法 | |
CN109077342B (zh) | 一种多等级烟叶模块打叶复烤投料加工方法 | |
WO2018218869A1 (zh) | 一种智能全自动药酒生产方法 | |
WO2018218868A1 (zh) | 一种智能全自动兰姆酒勾兑方法 | |
CN106406099B (zh) | 基于模糊匹配和率值反馈的水泥配料系统及方法 | |
CN103365208A (zh) | 一种具有底油性质偏差修正的原油调合优化方法 | |
Yanita et al. | Analysis of factors affecting the competitiveness of indonesian crude palm oil (cpo) export in the global market | |
Hristov et al. | Improvement of wheat quality in cultivars released in Serbia during the 20 th century | |
CN102591369A (zh) | 一种原油调合自动控制方法 | |
CN106124448B (zh) | 一种原油切换状态下的常减压装置进料性质预测方法 | |
CN105139084A (zh) | 一种油品在线全流程优化方法 | |
JP2018514184A5 (zh) | ||
CN111915232B (zh) | 减少>60mm粒度焦炭质量百分比的配煤调整方法 | |
KR101859167B1 (ko) | 가축분뇨 액비의 부숙도 측정방법 | |
CN109439357A (zh) | 肥煤参与的配煤炼焦的方法 | |
CN110866635B (zh) | 一种提高装置加工方案切换预测精度的方法 | |
CN104007658B (zh) | 一种具有双重循环优化的汽油在线调合方法 | |
Ma et al. | Chinese Dairy Farm Performance and Policy Implications in the New Millennium | |
KR101450505B1 (ko) | 고기능성 중합용 첨가제 | |
Prasada et al. | The potential loss of rice production due to wetland conversion in East Java | |
Zhang et al. | Quality differences between NILs of wheat variety Long 97-586 possessing HMW-GS 7+ 8 and 7 | |
Taymaz et al. | Estimating plant-level marginal costs and mark-ups in the Turkish manufacturing industry |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |