基于Plackett-Burman设计的气井采出水储罐安全研究

本文选项
	PDF全文阅读

	本文摘要

	本文图片

	参考文献


扩展功能




	电子期刊订阅

	RSS

本文作者相关文章
	曾维国

	许晓伟

	李超

	王轩

	范峥

曾维国¹ , 许晓伟² , 李超³ , 王轩² , 范峥⁴

1. 中国特种设备检测研究院;
2. 中国石油长庆油田分公司第二采气厂;
3. 中国石油长庆油田分公司气田开发事业部;
4. 西安石油大学化学化工学院

收稿日期：2019-11-04

基金项目：中国国家留学基金“高能效化工过程中的气液污染控制”(201908610135)

作者简介：曾维国(1984-)，硕士，工程师，主要从事油气田压力管道检验及风险评估等工作。E-mail：godblesszwg@126.com.

通信作者：范峥(1982-)，博士，副教授，主要从事油气加工过程腐蚀控制等工作。E-mail：fanzheng@xsyu.edu.cn.

摘要：针对天然气气井采出水储罐可能存在易燃易爆安全风险的实际生产问题，首先利用X射线衍射仪、离子色谱仪、气质联用仪对气井采出水中的固体杂质、无机杂质和有机杂质进行了检测分析，并利用Plackett-Burman设计通过室内实验研究了不同杂质含量下气井采出水体系的气-液相平衡变化规律。结果表明：气井采出水中的杂质主要包括Fe₃O₄等固体杂质、K⁺等无机杂质以及甲酸等有机杂质几大类；与纯甲醇-水体系相比，杂质会对气井采出水的气-液相平衡产生不同程度的影响且总体上表现出一定的正偏差；在这些杂质中，FeS、MgCO₃、K⁺、Ca²⁺、F^-、甲酸、正癸烷、正庚烷、1, 2, 4-三甲基苯和乙基环己烷等对甲醇气相摩尔分数均存在显著影响，其中乙基环己烷影响最大，Ca²⁺次之，FeS最小。研究结果可为气井采出水储罐安全运行提供科学、可靠的理论依据和数据来源，同时对同类型设备也具有一定的借鉴作用。

关键词：气井采出水储罐杂质检测分析气-液相平衡 Plackett-Burman设计

Safety research on gas well produced water storage tank based on Plackett-Burman design

Zeng Weiguo¹ , Xu Xiaowei² , Li Chao³ , Wang Xuan² , Fan Zheng⁴

1. China Special Equipment Inspection and Research Institute, Beijing, China;
2. The Second Gas Production Plant, PetroChina Changqing Oilfiled Company, Xi'an, Shaanxi, China;
3. Gas Field Development Office, PetroChina Changqing Oilfiled Company, Xi'an, Shaanxi, China;
4. College of Chemistry & Chemical Engineering, Xi'an Shiyou University, Yulin, Shaanxi, China

Abstract: Aiming at the severe problems of possibilities for flammable and explosive risks during the actual production procedure, the X-ray diffractometer, ion chromatography and GC/MS were used to detect and analyze solid impurities, inorganic impurities and organic impurities in the gas well produced water initially. Then the changing law of vapor-liquid phase equilibrium of the gas well produced water system with diverse impurities under different contents was explored by laboratory experiments using the Plackett-Burman design. The results showed that the impurities in the gas well produced water mainly include solid impurities such as Fe₃O₄, inorganic impurities such as K⁺, and organic impurities such as formic acid. The impurities could affect the gas-liquid phase equilibrium of gas well produced water to different degrees. Compared with the pure methanol-water system, the impact shows a certain positive deviation overall. Moreover, FeS, MgCO₃, K⁺, Ca²⁺, F^-, formic acid, n-decane, n-heptane, 1, 2, 4-trimethylbenzene and ethylcyclohexane among these impurities have significant influences on the molar fraction of methanol in vapor phase, among which ethylcyclohexane has the greatest effect followed by Ca²⁺, and FeS is the smallest. The research results could provide a scientific and reliable theoretical basis and data source for the safe operation of gas well produced water storage tanks, and a reference to the same types of equipment.

Key words: gas well produced water storage tank impurities determination and analysis vapor-liquid phase equilibrium Plackett-Burman design

天然气气井采出水来自不同的地质层位，成分复杂且差异较大，水质普遍呈偏酸性胶体状态，具有高矿化度、高浊度、高铁含量、高腐蚀性和低pH值这“四高一低”的显著特点，极易出现腐蚀结垢现象^[1-3]。同时，为了防止地面管线设备节流后出现天然气水合物冻堵，生产现场还经常向气井内加注甲醇来抑制水合物的生成^[4-8]，而天然气中的凝析油和压裂液、泡排剂、洗井液等油田助剂则也随之混入气井采出水中^[9-11]。除此之外，虽然液体甲醇的自燃点高达463 ℃，一般情况下很难发生自燃和爆炸，但在气井采出水储罐这一相对封闭的体系内，由于甲醇属于挥发性组分，极易在罐内上部空间形成甲醇蒸气和空气混合物，而一旦甲醇蒸气的含量处于可燃可爆范围内，只需极小的点火便会发生爆炸并产生巨大能量，进而引燃下部液体^[12-15]。鉴于气井采出水储罐甲醇蒸气含量主要是由罐内气-液相平衡所决定，重点加强气井采出水中关键杂质的日常监测分析，并进一步明确不同杂质含量下气井采出水体系的气-液相平衡变化规律是十分必要的。

尽管目前对于纯甲醇-水体系的气-液相平衡关系已有较为准确、深入的研究，但对气井采出水杂质的组成与含量缺乏必要的检测数据，对上述杂质影响气-液相平衡的作用规律缺乏明确的科学认识。因此，本研究以陕北某气田的气井采出水为对象，分别利用X射线衍射仪、离子色谱仪、气质联用仪对气井采出水中的固体杂质、无机杂质和有机杂质进行检测，并借助Plackett-Burman设计通过室内实验研究了不同杂质含量下气井采出水体系的气-液相平衡变化规律，从而为气井采出水储罐的正常、平稳、长周期运行提供准确、可靠的理论支撑和科学依据，以有效杜绝气井采出水储罐安全事故的发生。

1 实验部分

1.1 气井采出水中不同杂质的检测与分析

采用D/MAX-2400型粉末X射线衍射仪对气井采出水中的固体杂质进行检测分析。钼靶片，锆滤波片，管压46 kV，管流150 mA，氮气氛，扫描速度8°/min，扫描范围3°~90°，步进宽度0.02°/步。

采用ICS-3000型离子色谱仪对气井采出水中的无机杂质进行检测分析。阴离子分析柱为AS11-HC，淋洗液为氢氧化钾溶液，阳离子分析柱为CS12A，淋洗液为甲基磺酸溶液，淋洗液流量1.5 mL/min，柱温30 ℃，进样量25 μL。

采用5973-6890型气质联用仪对气井采出水中的有机杂质进行检测分析。具体色谱条件如下：一阶保留时间为1 min，一阶升温速率为0.5 ℃/min，一阶温度为65 ℃；二阶保留时间为1 min，二阶升温速率为30 ℃/min，二阶温度为150 ℃；三阶保留时间为1 min，三阶升温速率为30 ℃/min，终温为250 ℃；汽化室温度为250 ℃；检测器温度为250 ℃。

1.2 气井采出水气-液相平衡数据的测定

采用循环法测定气井采出水的气-液相平衡数据。在沸腾器中盛入一定量的气井采出水，恒压下加热。待液体沸腾后，溢出的气相经完全冷凝后流入收集器，溢流后经回流管流回到沸腾器中。由于气相冷凝液组成与沸腾器中液相组成不同，所以随着沸腾过程的进行，沸腾器和收集器中的液相组成不断改变，直至达到两容器中的液相组成不再变化而保持恒定。此时即达到了气-液相平衡，分别从气、液相取样口中取样进行气相色谱分析，即可得到在某一平衡温度下的一组气-液相平衡数据^[16-18]。利用注射器向沸腾器中不断补加甲醇以改变其液相组成，如此反复操作，最终得到气井采出水的气-液相平衡关系。

当采用气相色谱法时，选择HP-50⁺毛细管色谱柱，分流比1:20，柱流量1.2 mL/min，程序升温，初始柱温65 ℃，保持1 min，以30 ℃/min升至250 ℃，继续保持1 min，通过面积归一化法对其甲醇含量进行分析。

1.3 Plackett-Burman设计

Plackett-Burman设计作为一种从多个因素中选取对实验指标有显著影响因素的方法^[19-20]，属于典型的二水平部分实验设计，它对每个因子取高和低两水平进行分析，其中上限用“1”表示，下限用“-1”表示，通过比较各个因子两水平之间的差异确定因子的显著性。一般来说，对于n次实验至多可研究n-1个因子，但实际因子应该不多于n-4个，保留3个以上虚拟变量，用以估计实验误差。同时，对实验结果还需要进行方差分析，一般选择可信度大于95%或者显著性水平达到0.05的因子作为重要因子^[21-23]。

本研究采用MINITAB 17.2.0.0软件对气井采出水气-液相平衡的影响因素进行分析。

2 结果与讨论

2.1 气井采出水中杂质的组成和含量

为了研究不同杂质对气井采出水气-液相平衡的作用规律，首先对气井采出水中杂质的组成和含量进行了检测，结果见表 1。

表 1 气井采出水水质检测结果

由表 1可知，气井采出水中的固体杂质主要由Fe₃O₄、FeCO₃、FeS、CaCO₃和MgCO₃等组成，尤以氧铁化合物、硫铁化合物居多，故气井井筒的Fe-H₂S-CO₂-H₂O电化学腐蚀产物是固体杂质的主要来源。同时，由于该气井采出水属于CaCl₂水型，Ca²⁺、Mg²⁺含量很高，极易造成钙镁沉积并出现结垢型腐蚀，且导致金属出现点蚀现象的Cl^-质量浓度亦高达34 800.01 mg/L，因此该气井采出水储罐存在一定的泄漏隐患。除此之外，气井采出水中的有机杂质包括甲酸、乙酸、C₇⁺长链饱和烷烃以及少量芳香烃等。其中，甲酸和乙酸属于甲醇经氧化、缩合后的降解产物，而C₇⁺等重烃则来自天然气凝析油和气井压裂液、泡排剂、洗井液等油田助剂。

2.2 气井采出水的气-液相平衡

天然气气井采出水的气-液相平衡关系图如图 1所示。

图 1 天然气气井采出水的气-液相平衡关系图

由图 1可知，与纯甲醇-水体系相比，由于气井采出水中的杂质会对其气-液相平衡产生不同程度的影响且总体表现出一定的正偏差，故它们较纯甲醇-水体系存在明显差异。这是因为在溶液依数性的作用下，气井采出水会出现蒸气压下降、沸点升高和凝固点降低等现象。对于蒸发过程而言，气井采出水表面能量较大的甲醇分子将会克服分子间作用力从表面逸出，同时蒸发出来的甲醇分子也有一部分将会重新返回液体表面，经过一段时间后，甲醇的蒸发和凝聚速率相等，两者处于动态平衡，此时气井采出水表面将形成稳定的甲醇饱和蒸气压。然而，需要注意的是，若气井采出水中存在难挥发物质，则有部分液体表面被这种难挥发分子占据，相当于减少了单位面积上甲醇分子的数量，从而导致单位面积上甲醇分子在单位时间内进入气相的数量减少，使气井采出水中的甲醇蒸气压较纯甲醇-水体系中的甲醇蒸气压低。

2.3 气井采出水中杂质对气-液相平衡的影响

根据气井采出水各杂质组成和含量的检测结果，本研究在常压下分别对甲醇液相摩尔分数为0.457 9的气井采出水中的固体杂质Fe₃O₄、FeCO₃、FeS、CaCO₃、MgCO₃，无机杂质K⁺、Na⁺、Ca²⁺、Mg²⁺、总铁、F^-、Cl^-、SO₄^2-和有机杂质甲酸、乙酸、甲基环己烷、正壬烷、1, 3-二甲基环己烷、正癸烷、正辛烷、正庚烷、正十一烷、1, 2, 4-三甲基苯、间二甲苯、3-甲基庚烷、乙基环己烷共26个因素进行了全面考察，设计了实验次数为32、虚拟变量数为5的Plackett-Burman设计(见表 2、表 3)。

表 2 Plackett-Burman设计的因素与水平

表 3 Plackett-Burman设计方案

序号	A	B	C	D	E	F	G	H	I	J	K	L	M	N	O	P	Q	R	S	T	U	V	W	X	Y	Z	y(甲醇气相)
1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	0.634 2
2	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	0.686 8
3	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	0.732 4
4	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	0.680 5
5	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	0.652 1
6	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	0.622 4
7	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	0.665 8
8	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	0.626 4
9	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	0.669 6
10	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	0.671 4
11	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	0.702 8
12	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	0.736 5
13	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	0.680 1
14	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	0.659 8
15	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	-1	0.723 6
16	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	1	0.689 8
17	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	1	0.727 8
18	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	1	0.695 2
19	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	-1	0.747 8
20	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	1	0.654 5
21	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	-1	0.688 9
22	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	1	0.720 2
23	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	-1	0.655 7
24	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	-1	0.723 5
25	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	-1	0.685 7
26	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	-1	0.695 2
27	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	1	0.714 2
28	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	-1	0.707 5
29	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	-1	0.685 2
30	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	1	0.723 6
31	-1	-1	-1	1	-1	1	-1	1	1	1	-1	1	1	-1	-1	-1	1	1	1	1	1	-1	-1	1	1	-1	0.696 2
32	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	-1	0.670 5
注：1.“-1”代表最小值。 2.“1”代表最大值。

表 3 Plackett-Burman设计方案

利用方差分析对Plackett-Burman设计进行研究，以期进一步确定气井采出水中各杂质对气-液相平衡的影响大小(见表 4)。

表 4 Plackett-Burman设计方差分析

由表 4可知，在显著性水平α=0.05(置信水平95%)下，由于所采用的实验模型的F值为6.63，大于临界值F_0.05(26, 5)=4.515 1且概率P小于临界值0.05，故它是十分显著的，其噪音出现的概率为2.20%，这表明此模型的拟合度和可靠性较好，适合于探讨气井采出水中各杂质因素对响应值的影响程度。同时，FeS、MgCO₃、K⁺、Ca²⁺、F^-、甲酸、正癸烷、正庚烷、1, 2, 4-三甲基苯和乙基环己烷等对甲醇气相摩尔分数均存在显著影响，而Fe₃O₄、FeCO₃、CaCO₃、Na⁺、Mg²⁺、总铁、Cl^-、SO₄^2-、乙酸、甲基环己烷、1, 3-二甲基环己烷、正辛烷、正十一烷、间二甲苯和3-甲基庚烷等对其影响不显著。需要说明的是，显著性水平作为一种有量度的或然性评价参数，其大小用于表征当原假设为真时，根据样本观察结果计算的检验统计量落入拒绝域的概率。因此，显著性水平越小，置信水平越高，发生这种情况的概率就越小。

在气井采出水中，影响甲醇气相摩尔分数的主要因素及其排序如下：

影响因素：FeS、MgCO₃、K⁺、Ca²⁺、F^-、甲酸、正癸烷、正庚烷、1, 2, 4-三甲基苯、乙基环己烷；

影响大小：乙基环己烷>Ca²⁺>正癸烷>MgCO₃>K⁺>甲酸>F^->1, 2, 4-三甲基苯>正庚烷>FeS。

3 结论

(1) 对于气井采出水储罐这一相对封闭的体系而言，由于气井采出水不仅表现出“四高一低”的显著特点，同时还含有大量的水合物抑制剂甲醇以及天然气中的凝析油和压裂液、泡排剂、洗井液等油田助剂，故一旦甲醇蒸气和空气在罐内上部空间形成爆炸性混合物，则极易发生爆炸并造成严重安全事故。

(2) 气井采出水中的固体杂质由Fe₃O₄、FeCO₃、FeS、CaCO₃和MgCO₃等组成，以氧铁化合物、硫铁化合物居多，其主要来自气井井筒的Fe-H₂S-CO₂-H₂O电化学腐蚀产物；无机杂质包括K⁺、Na⁺、Ca²⁺、Mg²⁺、总铁、F^-、Cl^-和SO₄^2-等，其中Ca²⁺、Mg²⁺、Cl^-含量很高，存在结垢型腐蚀和点蚀的风险；有机杂质中的甲酸和乙酸属于甲醇经氧化、缩合后的降解产物，而C₇⁺等重烃则来自天然气凝析油和气井压裂液、泡排剂、洗井液等油田助剂。

(3) 与纯甲醇-水体系相比，由于气井采出水中的杂质会对其气-液相平衡产生不同程度的影响且总体表现出一定的正偏差，故它们较纯甲醇-水体系存在明显差异。这是因为在溶液依数性的作用下，气井采出水会出现蒸气压下降、沸点升高和凝固点降低等现象。

(4) 在以上研究的基础上，根据Plackett-Burman设计的实验结果，通过方差分析研究了不同杂质含量下气井采出水体系的气-液相平衡变化规律。研究结果表明，FeS、MgCO₃、K⁺、Ca²⁺、F^-、甲酸、正癸烷、正庚烷、1, 2, 4-三甲基苯和乙基环己烷等对甲醇气相摩尔分数均存在显著影响，其中乙基环己烷影响最大，Ca²⁺次之，FeS最小。而Fe₃O₄、FeCO₃、CaCO₃、Na⁺、Mg²⁺、总铁、Cl^-、SO₄^2-、乙酸、甲基环己烷、1, 3-二甲基环己烷、正辛烷、正十一烷、间二甲苯和3-甲基庚烷等则对响应值影响不显著。

参考文献

[1]	刘宏彬. 某气田采出水水质特性及处理现状分析[J]. 油气田地面工程, 2019, 38(8): 55-59.
[2]	黄斌, 王捷, 傅程, 等. 油田采出水处理技术研究新进展[J]. 现代化工, 2018, 38(8): 52-57.
[3]	张永红, 金艳, 何化, 等. 页岩气采出水处理工艺试验研究[J]. 天然气与石油, 2019, 37(3): 88-93.
[4]	宋光春, 李玉星, 王武昌, 等. 油气混输管道中天然气水合物的形成和堵塞过程研究[J]. 石油与天然气化工, 2017, 46(2): 38-43.
[5]	李栋梁, 唐翠萍, 梁德青, 等. 天然气水合物抑制过程中甲醇用量的影响[J]. 石油化工, 2009, 38(12): 1292-1296.
[6]	李聪, 赵会军, 王桌芳, 等. 天然气水合物抑制剂甲醇回收中结垢机理研究[J]. 天然气化工-C1化学与化工, 2019, 44(1): 67-72.
[7]	张琦. 天然气水合物的生产储运技术及现状[J]. 石油与天然气化工, 2013, 42(3): 261-264.
[8]	邵聪, 诸林, 何阳东, 等. 公式法预测天然气水合物形成条件及其评价[J]. 石油与天然气化工, 2018, 47(4): 50-56.
[9]	李建国, 董杰, 李浩. 浅谈天然气开采过程中凝析油的环境风险[J]. 环境与发展, 2014, 26(8): 93-96.
[10]	叶喆. 天然气输送过程中凝析油形成机理研究[J]. 煤气与热力, 2016, 36(10): 32-36.
[11]	陈宏, 邓征宇, 宫兆波, 等. 化学驱油田采出水处理技术的研究进展[J]. 环境科学与技术, 2017, 40(5): 95-101.
[12]	MOLINO A, LAROCCA V, VALERIO V, et al. Supercritical water gasification of lignin solution produced by steam explosion process on Arundo Donax after alkaline extraction[J]. Fuel, 2018, 221: 513-517.
[13]	王志荣, 蒋军成, 姜慧. 室外池火灾火焰环境研究进展[J]. 石油与天然气化工, 2005, 34(4): 321-324.
[14]	宋毅. 甲醇储罐区火灾后果危险性评价[J]. 广州化工, 2015, 43(5): 233-235.
[15]	姜海燕, 潘国军. 池火灾模型在甲醇储罐定量评价中的应用[J]. 山东化工, 2010, 39(3): 41-43.
[16]	于雪敏, 崔鹏, 刘海, 等. 水-环己酮-乙酸丁酯液液相平衡数据的测定关联[J]. 化学工程, 2018, 46(10): 35-39.
[17]	曹玲, 李群生. 乙酸乙酯-乙腈-1-辛基-3-甲基咪唑六氟磷酸盐物系等压汽液相平衡数据的测定与关联[J]. 石油化工, 2016, 45(2): 194-199.
[18]	VALDERRAMA J O, ALVAREZ V H. A versatile thermodynamic consistency test for incomplete phase equilibrium data of high-pressure gas-liquid mixtures[J]. Fluid Phase Equilibria, 2004, 226: 149-159.
[19]	MOSS A F, CHRYSTAL P V, DERSJANT-LI Y, et al. The ranked importance of dietary factors influencing the performance of broiler chickens offered phytase-supplemented diets by the Plackett-Burman screening design[J]. British Poultry Science, 2019, 60(4): 439-448.
[20]	MAULVI F A, PARMAR R J, SHUKLA M R, et al. Plackett-Burman design for screening of critical variables and their effects on the optical transparency and swelling of gatifloxacin-Pluronic-loaded contact lens[J]. International Journal of Pharmaceutics, 2019, 566: 513-519.
[21]	马秉智, 李丽, 赫军, 等. Plackett-Burman设计结合响应面法优化苦豆子效应成分的提取工艺[J]. 中国药学杂志, 2019, 54(12): 953-959.
[22]	范峥, 姬盼盼, 李超, 等. 模糊神经网络-遗传算法优化丙烯酸苄酯合成工艺[J]. 化工学报, 2019, 70(11): 4315-4324.
[23]	BENNETT K, CHEN Y T. A two-level Plackett-Burman non-geometric experimental design for main and two factor interaction sensitivity analysis of zigzag-channel PCHEs[J]. Thermal Science and Engineering Progress, 2019, 11: 167-194.