大型离心式通风机性能预测方法.pdf
第45卷 第2期 2 0 1 9年2月 工矿自动化 In d u s t r y a n d M i n e A u t o m a t i o n Vo l . 4 5N o . 2 F e b. 2 0 1 9 文章编号1 6 7 1-25 1 X(2 0 1 9)0 2-00 7 0-05 D O I1 0. 1 3 2 7 2/ j . i s s n. 1 6 7 1-25 1 x. 2 0 1 8 1 0 0 0 1 4 大型离心式通风机性能预测方法 孙涛1, 代邦武2, 褚菲2, 马小平2 ( 1. 云南能投煤业有限公司,云南 昆明 6 5 0 2 2 8; 2. 中国矿业大学 信息与控制工程学院,江苏 徐州 2 2 1 1 1 6) 摘要 针对现有离心式通风机性能预测方法不能充分利用离心式通风机历史运行数据、 建模周期长等问 题, 提出了基于最小二乘支持向量机(L S S VM) 与拉丁超立方采样( LH S) 的大型离心式通风机性能预测方 法。选取出口压力作为衡量通风机性能的指标, 利用 L S S VM 建立离心式通风机性能预测模型; 通过 LH S 方法采集离心式通风机的入口温度、 入口压力、 入口流量和转速, 将采集的数据进行归一化处理后用于 L S S VM 模型的训练; 通过测试数据验证所建立模型的有效性。仿真结果表明, 基于 L S S VM 与 LH S的大 型离心式通风机性能预测方法能够充分利用已有通风机数据信息快速准确地预测通风机性能。 关键词 离心式通风机;性能预测;最小二乘支持向量机;拉丁超立方采样;L S S VM;LH S 中图分类号T D 6 3 5 文献标志码A 网络出版地址h t t p / /k n s . c n k i . n e t /k c m s/d e t a i l/32. 1 6 2 7. T P. 2 0 1 9 0 1 1 8. 1 6 3 9. 0 0 3. h t m l 收稿日期2 0 1 8-10-11; 修回日期2 0 1 9-01-08; 责任编辑 胡娴。 基金项目 江苏省重点研发计划资金支持项目(B E 2 0 1 6 0 4 6) ; 江苏省普通高校研究生创新计划项目(KY L X 1 6 0 5 3 3) 。 作者简介 孙涛(1 9 6 1-) , 男, 四川苍溪人, 高级工程师, 主要研究方向为矿山机电设备设计与安装、 矿井自动化监控等, E-ma i ls t c u m t@ 1 2 6. c o m。 引用格式 孙涛, 代邦武, 褚菲, 等.大型离心式通风机性能预测方法[J].工矿自动化,2 0 1 9,4 5(2) 7 0-74. S UN T a o,D A I B a n g w u,CHU F e i,e t a l . P e r f o r m a n c e p r e d i c t i o n m e t h o d f o r l a r g e-ce n t r i f u g a l v e n t i l a t o r[J]. I n d u s t r y a n d M i n e A u t o m a t i o n,2 0 1 9,4 5(2) 7 0-74. Pe r f o r m a n c e p r e d i c t i o n m e t h o d f o r l a r g e-ce n t r i f u g a l v e n t i l a t o r SUN T a o 1, DA I B a n g w u 2, C HU F e i 2, M A X i a o p i n g 2 ( 1 . Y u n n a n N e n g t o u C o a l I n d u s t r y C o . ,L t d .,K u n m i n g 6 5 0 2 2 8,C h i n a;2. S c h o o l o f I n f o r m a t i o n a n d C o n t r o l E n g i n e e r i n g,C h i n a U n i v e r s i t y o f M i n i n g a n d T e c h n o l o g y,X u z h o u 2 2 1 1 1 6,C h i n a) A b s t r a c tI n v i e w o f p r o b l e m s t h a t e x i s t i n g p e r f o r m a n c e p r e d i c t i o n m e t h o d s f o r c e n t r i f u g a l v e n t i l a t o r c a n n o t f u l l y u t i l i z e h i s t o r i c a l o p e r a t i o n d a t a o f c e n t r i f u g a l v e n t i l a t o r a n d h a v e l o n g m o d e l i n g p e r i o d, p e r f o r m a n c e p r e d i c t i o n m e t h o d f o r l a r g e-ce n t r i f u g a l v e n t i l a t o r b a s e d o n L S S VM a n d LH S w a s p r o p o s e d . O u t l e t p r e s s u r e i s s e l e c t e d a s i n d e x t o m e a s u r e p e r f o r m a n c e o f v e n t i l a t o r,a n d p e r f o r m a n c e p r e d i c t i o n m o d e l o f c e n t r i f u g a l v e n t i l a t o r i s e s t a b l i s h e d b y u s i n g L S S VM. I n l e t t e m p e r a t u r e,i n l e t p r e s s u r e ,i n l e t f l o w r a t e a n d r o t a t i o n a l s p e e d o f t h e c e n t r i f u g a l v e n t i l a t o r a r e c o l l e c t e d b y LH S m e t h o d,a n d t h e c o l l e c t e d d a t a a r e n o r m a l i z e d f o r t r a i n i n g o f L S S VM m o d e l .V a l i d i t y o f t h e e s t a b l i s h e d m o d e l i s v e r i f i e d b y t e s t i n g d a t a . T h e s i m u l a t i o n r e s u l t s s h o w t h a t t h e p e r f o r m a n c e p r e d i c t i o n m e t h o d f o r l a r g e-ce n t r i f u g a l v e n t i l a t o r b a s e d o n L S S VM a n d LH S c a n m a k e f u l l u s e o f e x i s t i n g v e n t i l a t o r d a t a i n f o r m a t i o n t o q u i c k l y a n d a c c u r a t e l y p r e d i c t p e r f o r m a n c e o f v e n t i l a t o r . K e y w o r d sc e n t r i f u g a l v e n t i l a t o r;p e r f o r m a n c e p r e d i c t i o n;l e a s t s q u a r e s s u p p o r t v e c t o r m a c h i n e;L a t i n h y p e r c u b e s a m p l i n g;L S S VM;LH S 0 引言 离心式通风机的主要作用是保证井下空气供 应, 稀释有害气体, 降低井下煤尘浓度, 在煤矿安全 生产中至关重要。稳定的通风机性能直接关系到井 下设备的可靠运转与人员的安全[ 1-2]。而对通风机 性能的预测控制与运行优化是以准确的性能预测模 型为基础的, 因而建立准确的通风机性能预测模型 十分重要。 建立离心式通风机性能预测模型的方法主要有 3类 ①以应用数学、 流体力学和流场等相关理论 为基础建立离心式通风机模型, 对通风机性能进行 预测[ 3]。② 实验法, 即借助先进的测量技术建立各 种工况下离心式通风机实验模型[ 4-5]。③ 以计算机 技术为依托, 利用各种 C F D(C o m p u t a t i o n a l F l u i d D y n a m i c s) 数值模拟技术建立离心式通风机性能预 测模型[ 6-9]。这些方法往往需要复杂的数学计算和 重复的实验设计, 建模周期长, 成本高, 并且存在不 能充分利用通风机历史运行数据等问题, 造成信息 资源的浪费。近年来, 随着人工智能算法的发展, 数 据驱动的建模方法逐渐被用于通风机性能预测。 张玮[ 10]利用离心式通风机历史运行数据提出了一 种基于模糊R B F神经网络的离心式通风机建模方 法, 该方法取得了一定效果, 但神经网络建模所需数 据量大, 建模周期长, 且未对建模数据分布进行优 化, 可能导致建模数据过度集中, 模型容易陷入局部 最优。 本文提出了基于最小二乘支持向量机(L e a s t S q u a r e s S u p p o r t V e c t o r M a c h i n e,L S S VM) 与拉丁 超立方采样(L a t i n H y p e r c u b e S a m p l i n g,LH S ) 的 大型离心式通风机性能预测方法, 利用L S S VM算 法和离心式通风机历史运行数据建立性能预测模 型, 并采用LH S方法保证建模数据均匀分布在建 模区间, 提高模型泛化性。 1 L S S VM模型 支持向量机(S u p p o r t V e c t o r M a c h i n e, S VM) 以统 计 学 习 理 论 为 基 础, 利 用 V C 维 (V a p n i k - Ch e r v o n e n k i s D i m e n s i o n)和 结 构 风 险 最 小 化 ( S t r u c t u r e R i s k M i n i m i z a t i o n,S RM) 原理的基本 思想, 在有限样本的情况下, 从学习能力和模型的复 杂度之间找到一个最优权衡, 从而获得一个强泛化 性的学习机器[ 11], 可以有效避免“ 过学习” 和“ 维数 灾难” 等问题, 在自然语言处理、 文本识别、 图像识别 等领域获得广泛的应用[ 12]。在此基础上, J . A. K. S u y k e n s等 [13]将损失函数直接定义为误差平方和, 提出了一种改进的 S VM 算法 L S S VM。其主 要优势是转化S VM 中的不等式约束为等式约束, 从而 降 低 算 法 计 算 的 复 杂 性, 加 快 算 法 求 解 的 速度[ 14-16]。 设有训练集(x i,yi) , i=1,2, ,l,l为样本数, xi∈R d , 是 第i个样本的输入, yi∈R , 是 第i个样本 的期望输出, 构造回归函数 f(xi)=ω T φ( xi)+b+ei ( 1) 式中 ω 为 L S S VM 的权系数; b为偏差;e为误差; φ( ) 为非线性情况下将输入从低维空间变换到高 维空间的核函数。 定义损失函数为 p(yi,f(xi) )=( yi-f(xi) ) 2 =e 2 ( 2) 根 据 结 构 风 险 最 小 原 理,L S S VM 风 险 函 数为[ 17] J= 1 2 ω 2 +r 2∑ l i=1 e 2 i ( 3) 式中r为惩罚系数。 据此可以得到最优化问题 m i n1 2 ω 2 +r 2∑ l i=1 e 2 i s. t .f(xi)=ω T φ( xi)+b+ei ( 4) 构造L a g r a n g e函数 L(ω,b,ei,ai)= 1 2 ω 2 +r 2∑ l i=1 e 2 i- ∑ l i=1 ai{ω T φ( xi)+b+ei-yi}(5) 式中a i 为L a g r a n g e乘子,ai≥0。 根据KKT条件[ 18] , 最 优条件为 L ω =ω-∑ l i=1 aiφ(xi)=0→ω= ∑ l i=1 aiφ(xi) L b =-∑ l i=1 ai=0→ ∑ l i=1 ai=0 L ei =r∑ l i=1 ei-∑ l i=1 ai=0→ai=rei L ai =0→ω T φ( xi)+b+ei-yi= 烅 烄 烆 0 ( 6) 消去ω和e , 求 解式( 6) 可得 01v 1TvΩ+D- 1 [] r b [] a = 0 [] y ( 7) 式中1v=[1,1,,1] ; 对角阵Dr=d i a g[r,r,, r];y=[y1,y2,,yl] T; a=[a1,a2,,al] T ; 矩 阵Ω 的元素Ω i j =φ T( xi) φ( xj)=K(xi,xj),K(xi,xj) 为满足 M e r c e r条件的对称函数, 也称核函数。 求解式( 7) 得到ai与b , 则L S S VM 回归预测模 型为 17 2 0 1 9年第2期 孙涛等 大型离心式通风机性能预测方法 f(x)=∑ l i=1 aiK(xi,xj)+b ( 8) L S S VM 的回归预测能力主要取决于核函数与 具体 参 数 的 选 取。 实 际 应 用 中 通 常 选 择 满 足 M e r c e r条件的核函数, 一般分为两大类 局部核和 全局核[ 19]。若将两者线性组合, 就可以将两者的优 点集于一身, 使得到的模型具有更强的泛化能力。 常用的核函数主要有高斯径向基核函数、 多项式核 函数、 多层感知器核函数等。 2 离心式通风机性能预测模型 2. 1 离心式通风机数据采集 数据是建立离心式通风机模型的基础, 因而需 要设计实验来采集必要的离心式通风机建模数据。 影响离心式通风机性能的输入变量有很多, 忽略次 要变量的影响, 选取入口压力、 入口温度、 入口流量 和转速作为影响离心式通风机性能的主要变量, 选 取出口压力作为衡量通风机性能的指标。为了提高 所建模型的泛化性, 避免局部建模, 采集的训练与测 试数据应该均匀分布在通风机的整个运行区间内。 LH S采用分层抽样, 将采样区间均匀划分为若干等 份, 在每个部分中随机采集数据, 保证了数据分布的 均匀 性, 避 免 了 数 据 过 度 集 中[ 20-21]。因此, 选 用 LH S方法采集离心式通风机的实验数据。G.E. P. B o x等 [22]指出, 在实验的初始阶段, 采集的数据 不宜超过实验数据总量的2 5%, 假设采集数据总量 n=1 0d(d为输入变量维度) , 那么初始实验中采集 的实验数据量n0应满足n0≤0. 25n=2. 5d, 因而, 本文采用n0=0. 2 5n作为初始试验阶段实验数据 量。数据采集硬件实现方案如图1所示, 首先通过 传感器对被测通风机的入口压力、 温度、 流量及转速 进行测量, 然后测量数据通过总线传输到 D AQ数 据采集系统, 最终 D AQ 数据采集系统将数据打包 通过I /O设备传入主机。 图1 离心式通风机数据采集硬件实现方案 Fi g . 1 H a r d w a r e i m p l e m e n t a t i o n o f d a t a a c q u i s i t i o n o f c e n t r i f u g a l v e n t i l a t o r 2. 2 数据处理 由于变量之间存在量纲的差异, 采集到的数据 并不直接应用于模型训练, 需要对数据做归一化处 理, 即把有量纲数据变为无量纲数据, 将采集的数据 映射到[ 0,1] 范围内, 从而提高模型的收敛速度, 提 升模型的精度。本文采用最大-最小值归一化方 法, 将采集的数据映射到[0,1] 的范围内, 具体表达 式为 S= s-smi n sma x-sma x ( 9) 式中S为归一化后的数据;s为采集数据;smi n为采 集数据最小值; sma x为采集数据最大值。 2. 3 模型训练与模型验证 离心式通风机性能预测模型训练结构如图2所 示。该结构可以分为2个部分, 即数据采集处理和 模型训练。前者主要完成实验数据的采集与处理, 后者实现性能预测模型的建立与验证。首先通过 LH S方法采集离心式通风机的实验数据( 入口温 度、 压力、 流量和通风机转速) , 将采集的数据进行归 一化处理后用于L S S VM模型的训练[ 23]。模型训 练完成后, 将测试数据应用到所建立的模型, 验证模 型的有效性。如果所建立的模型满足建模的停止条 件, 则应用此模型。若所建立的模型无法满足建模 的停止条件, 则需要采集更多的数据进行模型训练。 本文选取R B F核函数作为L S S VM的核函数, 核参 数通过网格搜索法获得[ 13]。 图2 离心式通风机性能预测模型训练结构 Fi g . 2 T r a i n i n g s t r u c t u r e o f p e r f o r m a n c e p r e d i c t i o n m o d e l f o r c e n t r i f u g a l v e n t i l a t o r 此 外, 选 取 均 方 根 误 差 (R o o t-Me a n-S q ua r e E r r o r,RMS E) 和 平 均 相 对 误 差 (M e a n R e l a t i v e E r r o r,MR E) 来评估模型的预测性能, 同 时选取 RMS E作为模型的训练停止条件, 二者具体表达 式为 27 工矿自动化 2 0 1 9年第4 5卷 RMS E= 1 l∑ l i=1 ( ti-t - i) 槡 2 ( 10) MR E= 1 l∑ l i=1 ti-t - i ti 10 0%(1 1) 式中 ti为真实值;t - i为模型预测值。 3 仿真分析 某煤矿主要通风机采用离心式通风机, 本文以 该离心式通风机为研究对象, 利用 L S S VM 算法建 立通风机性能预测模型, 验证所提方法的有效性。 模型训练与测试样本均从现场的集散控制系统中获 得。利用LH S方法从离心式通风机稳定运行区域 内选取1 0 0组数据用于模型训练, 选取5 0组测试数 据用于模型验证, 并确定模型训练的停止条件为 RMS E<0. 0 5。利用 M a t l a b实现上述模型。图3 给出了在不同数量训练样本的情况下, 预测模型的 RMS E情况。从图3可以看出, 随着训练样本的增 加, 预测模型的RMS E值不断下降, 最终趋于稳定, 在训练样本数为3 0时, 模型满足训练停止条件。 图3 LS S VM方法均方根误差 F i g . 3 R o o t-me a n-s q ua r e e r r o r o f L S S VM m e t h o d 图4展示了在满足模型停止条件, 即使用3 0个 训练样本时, 模型的预测值与真实值的对比。从 图4可以看出, 所建立的模型能够很好地预测离心 式通风机的输出压力, 预测值与真实数据较吻合。 图4 通风机性能预测模型预测值与真实值对比 Fi g . 4 C o m p a r i s o n o f p r e d i c t e d a n d a c t u a l v a l u e s o f p e r f o r m a n c e p r e d i c t i o n m o d e l f o r c e n t r i f u g a l v e n t i l a t o r 为了进一步体现所提方法的优越性, 表1给出 了在 满 足 停 止 条 件 后, 不 同 训 练 样 本 下 模 型 的 RMS E和 MR E。从表1中可以看出, 误差随着训练 数据的增加而呈下降趋势, 并最终逐渐收敛。 表1 不同训练样本下LS S VM方法的MR S E与MR E值 T a b l e 1 RM S E a n d MR E o f L S S VM m e t h o d u n d e r d i f f e r e n t t r a i n i n g s a m p l e s 训练样本数 RM S E/1 0- 2 MR E/% 3 0 1. 9 0 2. 3 0 4 0 1. 4 0 2. 0 0 5 0 0. 5 0 1. 8 0 6 0 0. 3 6 1. 5 0 7 0 0. 3 6 0. 5 3 8 0 0. 3 3 0. 5 5 9 0 0. 3 2 0. 5 2 10 0 0. 3 2 0. 5 3 4 结语 基于L S S VM算法建立矿用离心式通风机性能 预测模型, 采用LH S方法采集通风机的实验数据, 进一步降低建模成本, 提高建模精度, 通过仿真实例 验证所提方法的有效性。然而实际生产中存在大量 相似的离心式通风机, 尽管它们在大小、 结构、 转速 等方面存在差异, 但所遵循的机理相似, 因而如何利 用已有相似离心式通风机数据建立当前离心式通风 机模型成为下一步研究方向。 参考文献(R e f e r e n c e s) [1] 牛青.矿井通风机性能测试方法与系统的研究[D]. 西安 西安科技大学,2 0 1 4. [2] 万雷.单级离心风机气动性能模拟与实验验证[D]. 哈尔滨 哈尔滨工程大学,2 0 1 3. [3] 李超, 伍继浩.离心通风机整机损失数学模型的研究 及优化[J]. 风机技术, 2 0 0 0(4) 1 5-17. [4] 姜华, 宫武旗, 张炜, 等. 带分流叶片离心叶轮非定常 流场的 实 验 研 究 [J].西 安 交 通 大 学 学 报,2 0 0 9, 4 3(9) 1 4-18. J I ANG H u a,G ONG Wu q i,Z HANG W e i,e t a l . E x p e r i m e n t a l i n v e s t i g a t i o n o f t h e u n s t e a d y d i s c h a r g e f l o w i n a c e n t r i f u g a l i m p e l l e r w i t h s p l i t t e r b l a d e s[J]. J o u r n a l o f X i ' a n J i a o t o n g U n i v e r s i t y,2 0 0 9,4 3( 9) 1 4-18. [5] 聂波, 张俊林, 满建楠. 进口弯管对前向离心通风机性 能影响的研究[J].流体机械, 2 0 1 8,4 6(1) 6-9. N I E B o,Z HANG J u n l i n,MAN J i a n n a n . E f f e c t o n t h e 37 2 0 1 9年第2期 孙涛等 大型离心式通风机性能预测方法 p e r f o r m a n c e o f t h e f o r w a r d c e n t r i f u g a l f a n w i t h i n l e t b e n t p i p e[J]. F l u i d M a c h i n e r y,2 0 1 8,4 6(1) 6-9. [6] ONMA P,CHAN T R A S M I T. C o m p a r i s o n o f t w o m e t h o d s t o d e t e r m i n e f a n p e r f o r m a n c e c u r v e s u s i n g c o m p u t a t i o n a l f l u i d d y n a m i c s[C] / /I O P C o n f e r e n c e S e r i e sM a t e r i a l s S c i e n c e a n d E n g i n e e r i n g,2 0 1 8. [7] MA DHWE S H N,KA R AN TH K V,S HA RMA N Y. E f f e c t o f i n n o v a t i v e c i r c u l a r s h r o u d f e n c e s o n a c e n t r i f u g a l f a n f o r a u g m e n t e d p e r f o r m a n c e-A nu m e r i c a l a n a l y s i s[J]. J o u r n a l o f M e c h a n i c a l S c i e n c e a n d T e c h n o l o g y,2 0 1 8,3 2(1) 1 8 5-19 7. [8] 黄忠文, 王培, 韩海燕.离心通风机叶轮的有限元建模 与应力分析[J].流体机械, 2 0 1 5,4 3(1 0) 2 7-30. HUANG Z h o n g w e n,WANG P e i,HAN H a i y a n . F i n i t e e l e m e n t m o d e l i n g a n d s t r e s s a n a l y s i s o f c e n t r i f u g a l f a n i m p e l l e r[J]. F l u i d M a c h i n e r y,2 0 1 5, 4 3(1 0) 2 7-30. [9] 邵国建, 陈帮, 谢云川.进气箱对离心风机性能影响的 研究[J].风机技术, 2 0 1 8,6 0(1) 3 9-43. S HAO G u o j i a n,CHE N B a n g,X I E Y u n c h u a n . E f f e c t o f i n l e t b o x o n t h e p e r f o r m a n c e o f c e n t r i f u g a l f a n[J]. C h i n e s e J o u r n a l o f T u r b o m a c h i n e r y,2 0 1 8,6 0(1) 3 9-43. [1 0] 张玮.基于神经网络的数据统计建模[D].杭州 浙江 工业大学,2 0 0 9. [1 1] VA P N I K V N. A n o v e r v i e w o f s t a t i s t i c a l l e a r n i n g t h e o r y[J]. I E E E t r a n s a c t i o n s o n n e u r a l n e t w o r k s , 1 9 9 9,1 0(5) 9 8 8-99 9. [1 2] CHE R KA S S KY V,MA Y. P r a c t i c a l s e l e c t i o n o f S VM p a r a m e t e r s a n d n o i s e e s t i m a t i o n f o r S VM r e g r e s s i o n [J]. N e u r a l N e t w o r k s, 2 0 0 4,1 7(1) 1 1 3-12 6. [1 3] S UYK E N S J A K,VAN D EWA L L E J . L e a s t s q u a r e s s u p p o r t v e c t o r m a c h i n e c l a s s i f i e r s[J].N e u r a l P r o c e s s i n g L e t t e r s, 1 9 9 9,9(3) 2 9 3-30 0. [1 4] 周欣然.基于最小二乘支持向量机的在线建模与控制 方法研究[D]. 长沙 湖南大学, 2 0 1 2. [1 5] 陈爱军.最小二乘支持向量机及其在工业过程建模中 的应用[D].杭州 浙江大学, 2 0 0 6. [1 6] L U C h a o,CHE N J i e,HONG R o n g j i n g,e t a l . D e g r a d a t i o n t r e n d e s t i m a t i o n o f s l e w i n g b e a r i n g b a s e d o n L S S VM m o d e l[J]. M e c h a n i c a l S y s t e m s a n d S i g n a l P r o c e s s i n g,2 0 1 6,7 6-773 5 3-36 6. [1 7] WANG H,HU D. C o m p a r i s o n o f S VM a n d L S-SVM f o r r e g r e s s i o n[C] / /I n t e r n a t i o n a l C o n f e r e n c e o n N e u r a l N e t w o r k s a n d B r a i n,2 0 0 5. [1 8] 张檑, 李宏光.基于KK T条件选择被控变量的自优 化控制方法[J].北京化工大学学报( 自然科学版) , 2 0 1 3,4 0( 增刊1) 6 7-71. Z HANG L e i,L I H o n g g u a n g . S e l f-o p ti m i z i n g c o n t r o l a p p r o a c h w i t h K a r u s h K u h n T u c k e r(KK T)b a s e d c o n t r o l l e d v a r i a b l e a l t e r n a t i v e s[J]. J o u r n a l o f B e i j i n g U n i v e r s i t y o f C h e m i c a l T e c h n o l o g y( N a t u r a l S c i e n c e E d i t i o n) ,2 0 1 3,4 0(S 1) 6 7-71. [1 9] 郑丽媛, 孙朋, 张素君.煤矿瓦斯突出预测的P S O- LS S VM模 型 [J].仪 表 技 术 与 传 感 器,2 0 1 4(6) 1 3 8-14 0. Z HE NG L i y u a n,S UN P e n g,Z HE NG S u j u n . P r e d i c t i o n o n g a s e m i s s i o n v a l u e b a s e d o n P S O- LS S VM[J]. I n s t r u m e n t T e c h n i q u e a n d S e n s o r,2 0 1 4 (6) 1 3 8-14 0. [2 0] MC KAY M D,B E C KMAN R J,C ONOV E R W J . A c o m p a r i s o n o f t h r e e m e t h o d s f o r s e l e c t i n g v a l u e s o f i n p u t v a r i a b l e s i n t h e a n a l y s i s o f o u t p u t f r o m a c o m p u t e r c o d e[J]. T e c h n o m e t r i c s,2 0 0 0,2 1(2) 2 3 9-24 5. [2 1] 吴振君, 王水林, 葛修润. LH S方法在边坡可靠度分 析中的应用[J].岩土力学, 2 0 1 0,3 1(4) 1 0