
Citation: | LIN Zhibin, GAO Hongli, WU Yudong, TAN Yongwen. Denoising of Acoustic Emission of Diamond-Coated Mechanical Seals Wear Based on Empirical Wavelet Transform and Kullback-Leibler Divergence[J]. Journal of Southwest Jiaotong University, 2024, 59(1): 177-184. doi: 10.3969/j.issn.0258-2724.20210599 |
In order to obtain the pure wear acoustic emission of diamond-coated mechanical seal, the denoising method based on empirical wavelet transform (EWT) and Kullback-Leibler divergence (KLD) was proposed. Firstly, filter bank was calculated with empirical wavelet transform on acquired acoustic emission signal. Then the filter bank was applied to both the acquired acoustic emission signal and background noise acoustic emission signal. The Kullback-Leibler divergences were calculated between the corresponding bands of two signals. The cumulative sum algorithm was employed to find a threshold for determining whether the corresponding band is used for signal reconstruction. The results show that the proposed method can effectively suppress the noise of acoustic emission signals under different working conditions and wear states, and effectively improve the signal-to-noise ratio of wear acoustic emission signals, especially weak wear signals. Compared with the traditional denoising methods, the proposed EWT-KLD method has stronger adaptability and stability for denoising of wear acoustic emission signal under different working conditions, which is of great significance for the monitoring early seal wear and the cumulative wear process of seal.
随着城市范围的急剧扩张以及城市经济的快速发展,城市出行需求快速增长,随之而来的城市交通拥堵问题也开始加剧. 为缓解急速增长的出行需求与缓慢扩张的路网通行能力之间的矛盾,城市管理者在增加公交线网覆盖率的同时,也在努力提高公共交通服务水平,促使使用低乘坐率小汽车出行的居民改用高乘客率的公共交通出行.
在公交服务的优化中,许多研究已将一些实质性的现实世界细节明确地纳入其模型,例如不同的车队构成(如电动车和内燃车)、多种运营方案(如全停和跃站停靠)以及站点位置的地理限制[1-4]. 然而,现有的大多数模型都建立在固定的公交线路上. 现实中有些公交系统(例如地铁)由于其高昂的建设成本,在规划阶段就需要考虑到远期需求的变化[5-6]. 因此,传统的优化方法(公交线一次性建设完全线并投入使用)无法适用于这种需要兼顾远期需求变化的公交线设计.
多阶段设计在轨道交通建设中比较常见,例如成都地铁、上海轨道交通和深圳地铁的线路,都至少分两期完成全线建设,时间跨度在十年左右. 多阶段设计方法包括向两端延伸线路和线路内增设站点 2种方式,本文主要研究第1种方式. 针对公交系统多阶段设计,已有研究者作了相关研究. Matisziw等[7]提出在已有轨道路线不满足需求时如何延伸线路的方法;Cheng等[8]提出一个将预定轨道交通线最佳细分为多个部分以进行分期建设的模型;Sun等[9]通过双层规划方法提出一个轨道交通线多阶段设计的模型. 虽然上述模型都被应用到预定的轨道交通线路中,即给出一些预设车站的位置,并着重于解决何时延伸以及被延长公交线路分期实施问题,然而,忽略此时的公交站点选址优化问题. 因此,Xia等[10]提出解决上述问题的一种方法,但其仍假设出行需求是固定变化的,忽略了出行需求的随机特性.
随机优化和鲁棒性优化是2种研究随机需求对公交系统影响的方法. 前者旨在获得公交系统或乘客的最低期望成本[11-14],后者是最小化与不满足系统服务情况相关的成本[15-16]. 大多数研究人员倾向于使用随机优化方法[12,14,17],而忽略随机因素影响下方案的鲁棒性问题. Amiripour等[18]根据不同季节需求的波动,为鲁棒性网络设计提出一种基于期望系统成本的模型;An等[15]研究了随机需求下固定线路和“拨号叫车”的多方式服务优化问题,并构建一个基于糟糕情况下的成本最小化非线性模型;Hassannayebi等[11]通过引入过载情况下的惩罚成本,提出随机需求下列车时刻表优化模型;魏长钦等[19]以最大服务率和最小动态行程为目标,提出定制公交的站点及路径优化方法. 这些研究通常着重于追求目标函数(如系统成本)期望值的最优化,而忽略目标值的波动问题.
系统成本的方差是衡量系统鲁棒性能的重要属性,因此,同时考虑系统成本的均值和方差的优化方法更有利于决策者得到随机需求影响下的理想方案. 均值-方差(MV)理论在公交系统设计的应用方面,Yan等[20]提出一个鲁棒性的优化模型,用于固定公交路线的时间表设计,旨在最小化随机时刻表偏差的期望值及其波动的加权总和;Huang等[14]通过将乘客的出行时间、运营成本的期望值以及运营成本的方差线性加权成目标函数,并将其降至最低来优化发车时刻表. 上述相关研究虽然使用MV理论保障了方案在随机需求下的稳定性,但仍存在一些局限性:其求解过程由于利用离散模型而复杂;大多数研究都集中在发车频率或时刻表的优化上,而忽略了站点位置优化问题.
本文基于文献[10],利用MV理论构建随机增长需求下的多阶段公交走廊设计决策方法. 主要贡献如下:1) 构建多种风险态度下的决策模型,以表达不同投资者在随机增长需求下的不同决策态度;2) 提出适用于随机需求下的多阶段公交系统优化模型,以优化不同时期的站点布局和发车间隔;3) 使用连续近似方法将离散问题连续化,可将求解结果表达为解析解,简化求解过程. 最后,本文使用不同的设计策略(公交线路全覆盖、单阶段、多阶段设计)和分析决策态度分析模型效果.
本文假设公交走廊布局如图1所示. 为提供日常通勤服务,公交车辆往返于CBD与公交线路的边界,lt⩽Lt(lt、Lt分别为t时段的公交线路边界、居民分布边界),t为年度. 公交车辆在每个车站停下来接送乘客. 在t=0,1,⋯,T, t∈T 期间,住户人口沿着CBD和边界Lt的走廊连续分布,其中T是计划周期时间段的集合. 假设出行需求是由住户的通勤行为决定的,住户往返于其住所和预定的就业地点(即本文中的CBD). 在规划期间,假设走廊的总住户数Nt以每年gt的速度增长,gt为随机变量,服从特定概率分布. 假设公交线路的长度可随住户向外扩张而延伸,考虑到需求的增加和走廊的延伸,公交运营商会寻求优化公交系统的设计,以更好地满足出行需求.
为融入方差对决策结果的影响,使用MV函数(式(1))作为目标函数,随机需求下多阶段公交走廊设计模型可表示为式(1)~(5).
minlt,ρ(x),htU(E(Z),D(Z)), | (1) |
lt−1⩽lt,∀t∈T, | (2) |
(τe−∫ltlt−1ρ(x)dx)(lt−lt−1)⩽0,∀t∈T, | (3) |
¯λo,t(0)κht⩽CV,∀t∈T, | (4) |
{ρ(x)⩾0,∀t∈T,x∈[0,lt],ht⩾0,∀t∈T, | (5) |
式中:ρ(x)为站点分布密度(公交站点位置可由站点密度函数积分得到[4,21,22]);ht为发车间隔;E(Z)、D(Z)分别为Z的期望、方差,Z为所有标准化系统成本Z的集合;τe为线路延伸范围内的最小公交站点数(本文中取1);ˉλo,t(0)是t时段处于中心的总乘客数的期望值;CV为车辆的载客量(乘客/车);κ为高峰小时期间每小时出行需求占全天的比例.
式(2)为公交线路长度约束,防止公交线缩短;式(3)为线路延伸长度最短约束,避免公交线路延伸范围内无公交站点;式(4)为车容量约束;式(5)为非负约束.
根据风险态度定义3种决策者:中立风险、规避风险和寻求风险. 众所周知,具有不同风险态度的决策者会对同一问题做出不同的决策. 针对不同风险态度的公交系统设计问题,可以采用效用理论方法或MV方法. 本文中的“风险态度”是指在MV框架下对风险的态度. 因此,中立风险的决策者漠视收益的波动,规避风险的决策者将尝试使收益的方差最小化,而寻求风险的决策者将尝试使收益的方差最大化. 如果通过成本期望值量化收益,并通过成本方差量化收益的不确定性水平,则可以通过MV函数来表达决策者的风险态度,具体描述如下:
1) 中立风险决策
该类决策者对收益不确定性“中立”且漠不关心,其满意度只取决于预期的系统成本. 因此,其MV函数为
U(E(Z),D(Z))=E(Z). | (6) |
2) 规避风险决策
该类决策者保守且不享有任何系统成本不确定性的决策者,其满意度随着系统成本和其不确定性的增加而降低,其MV函数表示为
U(E(Z),D(Z))=E(Z)+ξ√D(Z), | (7) |
式中:ξ为方差的影响因子.
3) 寻求风险决策
这类决策者是赌徒,享受系统成本不确定性带来的兴奋. 随着系统成本的增加,其满意度会降低,而随着系统成本的不确定性的增加,其满意度则会增加. 这类风险决策态度的MV函数为
U(E(Z),D(Z))=E(Z)−ξ√D(Z). | (8) |
上述对风险决策的描述方法源自文献[23]. 综合3种风险决策态度,MV函数可表示为
U(E(Z),D(Z))=E(Z)+ζξ√D(Z), | (9) |
式中:ζ为决策分析态度;ζ=0,表示中立风险决策;ζ=1,表示规避风险决策;ζ=−1,表示寻求风险决策.
系统成本Z是乘客在规划范围内各阶段t的通勤成本Γt和公交运营商的成本Λt的总和,并将其标准化为每户的日均时间成本,即
Z=∑Tt=0(Γt+Λt)μTw∑Tt=0Nt, | (10) |
式中:μ为时间价值;Tw为每年工作时间.
1) 用户成本
在时间段t中的任何一天,通勤者都从其居住地x乘公交出行,其旅行时间用φt(x)表示. φt(x)由三部分组成:乘客前往最近公交站点的走行时间At(x)、在公交站的等待时间Wt(x)和从站点到目的地的在车时间It(x),如式(11)所示.
φt(x)=At(x)+Wt(x)+It(x). | (11) |
为建立简约模型,应用连续近似法求解.
步骤1 假设到站点的走行时间At(x)是走行距离d(x)除以步行速度vw,即d(x)/vw. 如果乘客住在公交线所覆盖的范围内(即x⩽lt),则其走行距离近似定义为2个连续站点之间距离的1/4[4,10,24],即1/(4ρ(x));否则,需要增加从其居住位置x到公交线路边界的距离,即x−lt+1/(4ρ(lt)).
步骤2 乘客在车站的等候时间Wt(x)近似为发车间隔时间的一半[20-21] .
步骤3 乘客的在车时间组成:1) 车辆的巡航时间,即∫0x1vcdu,vc为公交车辆的巡航速度;2) 公交车辆经过车站时由于加速、减速和滞留而导致的车站延误,可通过∫0xTdρ(u)du进行估算,其中,Td为每个站点的延迟,并且假定为常数[22,25].
因此,时间段t的年度通勤总成本Γt与年度通勤次数和单次通勤成本有关,表示为
Γt=Twμ∫lt0λt(x)φt(x)dx, | (12) |
式中:λt(x)为日出行需求密度,是住户密度nt(x)(具体见2.4节)和每户日均出行量η(假设为给定的)的乘积,即λt(x)=ηnt(x).
2) 运营商成本
公交运营商的成本取决于4个指标[26]:公交线的长度Λ1(t)、公交站点数量Λ2(t)、每辆车行驶的距离Λ3(t)、公交车辆的行驶时间Λ4(t).
Λ1(t)=2lt, | (13) |
Λ2(t)=∫lt0ρt(x)dx, | (14) |
Λ3(t)=TrΛ1(t)ht, | (15) |
Λ4(t)=Λ3(t)vc+TrtdΛ2(t)ht, | (16) |
式中:Tr为公交每天运营时长.
式(13)中的“2”表示2个运行方向;式(14)表示公交2个运行方向共享一个站点;式(15)中,Λ1(t)/ht表示1 h内出发的所有车辆的行驶里程;式(16)车辆行驶时间由巡航速度行驶的时间和由于停车而延迟时间所构成.
因此,公交商运营成本为
Λt=Tw∑4i=1πiΛi(t), | (17) |
式中:πi是与Λi相关的单位成本.
当住户总数已知时,居民以特定分布函数在走廊中,且其需满足总住户数平衡条件,即
∫Lt0nt(x)dx=Nt, | (18) |
每年走廊总住户数会有多种可能性,Nt是随机变量. 针对每年可能的总住户数Nt,存在一种住户分布情况nt(x|Nt)以及对应的概率密度 P(Nt). 因此,其均值和方差分别为ˉnt(x)=∑nt(x|Nt)P(Nt),σ2(nt(x))=1N−1∑(nt(x|Nt)−ˉnt(x))2P(Nt). 其中,N为Nt的可能情况数量. 需求密度函数λt(x)、累计需求密度函数λo,t(x)的均值ˉλt(x)、ˉλo,t(x)以及标准差σ(λt(x))、σ(λo,t(x))都可由此方法计算得到,并用于公交系统优化.
该问题假定确定了满足住宅分布模型的解,即住户密度和走廊的大小已知,通过式(1)求解公交系统优化模型. 利用局部分解法和最优化理论[10,21]可得
ρ∗(x,ζ)=12√f1(x,ζ)f2(x,ζ),∀x∈[lt−1,lt],t∈T, | (19) |
h∗t(ζ)=min(htmax,4√g1g2(ζ)),∀t∈T, | (20) |
式中:ρ∗(x,ζ)和h∗t(ζ)分别为风险态度决策ζ下ρ(x)和ht的最优解;htmax为最大发车间隔;f1(x,ζ)、f2(x,ζ)、g1、g2(ζ)的取值与风险决策态度有关,分别如式(21)~(24)所示.
f1(x,ζ)=μ∑Tt′=tˉλt'(x)+ζξμ∑Tt′=tσ(λt′(x)), | (21) |
f2(x,ζ)=μtdvw(∑Tt′=tˉλo,t'(x)+ζξ∑Tt′=tσ(λo,t′(x))+π2(T−t+1)+2π4TrTd∑Tt′=t1ht′), | (22) |
g1=π3Trlt+π4Tr(ltvc+Td2∫lt0ρ(x)dx), | (23) |
g2(ζ)=μˉλo,t(0)−ζξμσ(λo,t(x)). | (24) |
最优的公交线路长度l∗t可以通过比较3种情况下的MV函数值决定:1) 公交线路长度无变化;2) 不受式(2)、(3)(公交线路长度)约束的式(1)的最优解;3) 约束式(3)的临界值.
l∗t=argminlt∈{l∗t−1,ltmin,ˆlt}U(lt,ρ∗(x),h∗t),t∈T, | (25) |
式中:ltmin为由公交线路延长约束式(3)确定的最短公交线路长度;ˆlt为没有公交线路延长约束式(3)的目标函数式(1)的最优解,可以通过求解其斜率为0的位置来获得,且满足式(26)所示关系.
μ(ˉλo,t(ˆlt)+ζξσλt(ˆlt))vw=2π1+π2ρ(ˆlt)+2Trπ3ht+Trπ4ht(2vc+Tdρ(ˆlt))+μ(ˉλo,t(ˆlt)+ζξσ(λo,t(ˆlt)))(1vc+Tdρ(ˆlt)). | (26) |
首先,假设户数Nt=Nt−1+gt, 初始住户数为N0,用户数增长率gt∈{ga,gb},其中,ga=ˉg−σg,gb=ˉg+σg,对应的概率分别为pg、1−pg,ˉg、σg分别为gt的期望、偏差. 然后,以轨道交通为例,分析住户均匀分布下(住户分布密度为k),3种设计策略和3种风险决策态度对轨道公交系统的优化参数和各成本情况的影响(所提模型适用于任何公交方式). 3种设计策略分别为:1) 全覆盖设计,也就是公交线路全覆盖居民分布走廊,令lt=Lt;2) 单阶段设计,即在开始时修建的公交线路在之后所有阶段都使用,并不发生变化(可由本文所提模型简化得到);3) 多阶段设计,即本文所提模型. 最后,表1列出了模型的参数值[25-27].
参数 | 数值 | 参数 | 数值 | |
vw/(km·h−1) | 2 | ζ | 1 | |
vc/(km·h−1) | 60 | k/(户·km−1) | 4600 | |
Td/(h·站−1) | 0.0125 | ˉg/(户·年−1) | 0.03 | |
CV/(人·车−1) | 1800 | σg/(户·年−1) | 0.01 | |
π1/(h·(km·d)−1) | 14256 + 475.2μ | pg | 0.7 | |
π2/ (h·(站·d)−1) | 7056 + 235.2μ | κ | 0.1 | |
π3/(h·km−1) | 2.2 | Tw/(d·年−1) | 250 | |
π4/(h·h−1) | 101 + 5 μ | N0/ 户 | 30000 | |
Tr/ (h·d−1) | 16 | μ/ (h·h−1) | 20 | |
η/ (次·(户·d) −1) | 4 |
公交线路全覆盖设计是最常见的公交线路设计策略. 表2为公交线路全覆盖设计结果,表中,括号中两数值分别为初始和最终的公交线路/走廊长度,余表同. 图2为3种风险决策态度公交线路全覆盖设计中各阶段的发车间隔.
风险决策 态度 | 公交线路/ 走廊长度/km | 平均站点 密度/ (站·km−1) | 平均发 车间隔/ min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | (6.81, 14.93) | 1.0937 | 2.1989 | 1.1280 | 0.8123 | 1.9403 | 0.0794 |
规避风险 | (6.81, 14.93) | 1.1773 | 2.1721 | 1.1261 | 0.8154 | 1.9414 | 0.0761 |
寻求风险 | (6.81, 14.93) | 0.9814 | 2.2252 | 1.152 | 0.8076 | 1.9595 | 0.1057 |
与中立风险决策态度的结果相比,规避风险态度决策者通过提供更优质的公交服务(更小的发车间隔和更密集的站点分布),降低用户通勤成本(约0.17%),从而得到更稳定的收益方案(即系统成本的波动降低了约4.16%). 在这种情况下,尽管更优质的公交服务增加了运营商成本(约0.38%),但系统成本仅增加0.06%,这相对于降低的系统成本标准差而言,是可接受的. 对于寻求风险态度,其结果相反. 该结论说明所提出的模型在3种风险决策态度下对公交线路全覆盖设计都具有良好的效果. 此外,在3种风险决策态度中,随着住户数增加,发车间隔都随之增加(见图2). 这是因为在公交线路随住户分布延伸的同时,为避免运营商成本的过度增加,延长的公交线路限制了公交系统其他参数性能的提升,即发车间隔.
公交线路单阶段设计内容包括公交线路长度、站点分布和发车间隔. 这意味公交线路并不一定全覆盖范围居民分布范围. 表3和图3为3种风险决策态度公交线路单阶段设计的结果.
风险决策 态度 | 线路 长度/km | 走廊 长度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | 10.97 | (6.81, 14.93) | 1.1011 | 2.2979 | 1.208 | 0.8504 | 2.0584 | 0.1334 |
规避风险 | 11.40 | (6.81, 14.93) | 1.0976 | 2.3131 | 1.1813 | 0.8827 | 2.064 | 0.1169 |
寻求风险 | 10.48 | (6.81, 14.93) | 1.1099 | 2.2766 | 1.2499 | 0.8144 | 2.0644 | 0.1546 |
从表3可得,公交线路单阶段设计中,与中立风险风险决策态度对比,规避风险态度提高12.37%的系统成本稳定性,寻求风险方式降低了15.89%. 虽然规避风险方式仍然能获得更稳定的系统成本(表3最后一列),但其主要原因不再是更低发车间隔和更密集站点分布,而是更长的公交线路(3.92%). 由于规避风险决策提供更长的公交线路,缩短了公交线路覆盖范围外住户的通勤时间,从而降低了波动. 这暗示最佳的公交线路长度也受风险决策态度的影响. 此外,为降低因公交线路长度增加而上涨的运营成本,公交系统发车间隔有所增加,平均站点密度有所降低,与全覆盖设计方案相反. 该结论验证了所提模型对公交线路单阶段设计的有效性,并表明其作用机理与全覆盖设计有所差异. 与全覆盖设计类似,单阶段设计中,相较于中立风险决策,规避风险决策也会降低居民通勤成本,增加运营商成本,并导致系统成本及其稳定性增加.
与全覆盖设计相比,公交线路单阶段设计的系统成本虽然更低,但其波动更大. 这种现象是合理的,因为其固定的公交线路布局难以适应因住户随机增长所造成的住户分布范围随机波动.
由图3可知:在公交线路单阶段设计方案中,随着住户数的增加(即公交出行需求的增加),公交系统的发车间隔也随之降低;与中立风险设计相比,发车间隔在规避风险设计中更高,而在寻求风险设计中更低.
多阶段设计中,公交线路随需求的增加而延长,这能有效地避免过早投资修建所造成的浪费(主要是运营和维护成本). 表4列出了3种风险决策态度下公交线路多阶段设计的结果(阶段数为5个),图4呈现了均匀需求下3种风险决策态度公交线路多阶段设计中各阶段的发车间隔和公交线路长度.
风险决策 态度 | 线路长度/km | 走廊长 度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 阶段 数/个 | 期待用户 通勤成本/ (h· (户·d)−1) | 期待运营 商成本/(h· (户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/(h· (户·d)−1) | |
初始 | 最终 | |||||||||
中立风险 | 6.82 | 12.31 | (6.81, 14.93) | 1.2080 | 2.0977 | 5 | 1.1539 | 0.7398 | 1.8936 | 0.1095 |
规避风险 | 6.82 | 12.76 | (6.81, 14.93) | 1.2186 | 2.0930 | 5 | 1.1420 | 0.7580 | 1.9000 | 0.0947 |
寻求风险 | 6.82 | 11.88 | (6.81, 14.93) | 1.2209 | 2.0916 | 5 | 1.1856 | 0.7147 | 1.9003 | 0.1352 |
与公交线路单阶段设计中类似,规避风险决策态度促使公交线路更长,而寻求风险的更短(表4和图4(b)). 不同的是,多阶段设计在公交线路单阶段风险决策态度下不仅系统成本更低,其系统成本的稳定性也更高,这与公交线路的可延伸性有关. 多阶段的公交线路设计能适时地响应增长的出行需求,降低因过早修建公交线路和站点所带来的空载运营成本. 与公交线路全覆盖设计相比,多阶段设计的优势在于其系统成本更低. 与公交线路全覆盖和单阶段设计都不同的是,寻求风险方式与规避风险方式对公交系统设计的影响并非完全相反. 规避风险决策和寻求风险决策都会增加平均站点密度和降低平均发车间隔,该现象与寻求风险决策中第一次公交线路延伸实施时间被推迟有关. 综上,多阶段设计中风险决策态度对公交系统设计的影响比较复杂.
图4表明,发车间隔在未延伸公交线路时随需求的增加而降低,而在整体上随出行需求的增加而增加. 前者因为公交线路长度未变化,类似公交线路单阶段设计中的发车间隔变化特征,而后者因为公交线路长度在变化,这与全覆盖设计中的结果类似. 图4(a)还说明在寻求风险方式中,更高的站点密度推迟了公交线路各次延伸的实施时间,但未减少设计方案的阶段数. 该结论暗示当目标函数中方差影响系数ξ足够大时,其设计方案的阶段数可能会减少.
本文对不同设计策略(全覆盖设计、单阶段设计和多阶段非全覆盖设计)和不同的风险决策态度(中立风险决策、规避风险决策和寻求风险决策)下的公交走廊设计进行了对比分析. 实验结果表明:
1) 本文提出的基于均值-方差理论的风险态度决策方法在公交走廊多阶段设计中效果良好,具体表现为与中立风险相比, 规避风险方式能提供系统成本更稳定的方案,而追求风险方式则得到系统成本波动更大的方案.
2) 多阶段设计具有比单阶段设计和全覆盖设计更低的系统成本,比单阶段设计更稳定的系统成本.
3) 不同风险决策态度都是通过调整用户的通勤成本来控制系统成本的波动(表现为系统成本标准差),即在规避风险模型中降低通勤成本,寻求规避模型中增加通勤成本.
4) 在不同设计策略中,不同风险决策态度对公交系统设计参数(即公交线路长度、站点分布和发车间隔)的影响有所不同,具体表现为与中立风险决策方案相比:在全覆盖设计中,规避风险决策会增加平均站点密度,并降低平均发车间隔,而寻求风险决策相反;在单阶段设计中,规避风险决策会增加公交线路长度和发车间隔,降低平均站点密集度,而寻求风险决策相反;在多阶段设计中,风险决策态度对公交系统平均发车间隔和站点密度的影响较为复杂,这与线路延伸的实施时间有关.
致谢:太原科技大学博士科研启动基金(20202048).
[1] |
矫健. 高温油泵机械密封失效分析与对策[J]. 内燃机与配件,2020,321(21): 75-76. doi: 10.3969/j.issn.1674-957X.2020.21.034
JIAO Jian. Failure analysis and countermeasure of mechanical seal of high temperature oil pump[J]. Internal Combustion Engine and Parts, 2020, 321(21): 75-76. doi: 10.3969/j.issn.1674-957X.2020.21.034
|
[2] |
郝丽霞. 化工离心泵机械密封失效原因分析[J]. 中国设备工程,2020(3): 143-145.
HAO Lixia. Failure analysis of chemical centrifugal pump mechanical seal[J]. China Equipment Engineering, 2020(3): 143-145.
|
[3] |
王超. 高温油泵机械密封失效分析及改进[J]. 石油化工设备,2019,48(6): 62-67.
WANG Chao. Failure analysis and improvement of mechanical seal of high temperature oil pump[J]. Petro-demical Equipment, 2019, 48(6): 62-67.
|
[4] |
黄文涛. 机械密封在化工上的应用及泄漏的应对[J]. 化学工程与装备,2018(3): 185-186.
HUANG Wentao. Application of mechanical seal in chemical industry and response to leakage[J]. Chemical Engineering and Equipment, 2018(3): 185-186.
|
[5] |
袁明蛟. 石油化工生产常用机械密封形式及失效原因探讨[J]. 中国化工贸易,2017,9(10): 204-206. doi: 10.3969/j.issn.1674-5167.2017.10.182
YUAN Mingjiao. Discussion on commonly used mechanical seal forms and failure reasons in petrochemical production[J]. China Chemical Trade, 2017, 9(10): 204-206. doi: 10.3969/j.issn.1674-5167.2017.10.182
|
[6] |
陈志高. 某重水堆核电站主泵机械密封失效原因分析[J]. 设备管理与维修,2019(13): 161-164. doi: 10.16621/j.cnki.issn1001-0599.2019.07.79
CHEN Zhigao. Failure analysis of mechanical seal of main pump in a heavy water reaction nuclear power plant[J]. Equipment Management and Maintenance, 2019(13): 161-164. doi: 10.16621/j.cnki.issn1001-0599.2019.07.79
|
[7] |
郭洪涛. 反应器冷却水泵机械密封失效原因及对策[J]. 齐鲁石油化工,2019,47(1): 51-52,56. doi: 10.3969/j.issn.1009-9859.2019.01.014
GUO Hongtao. Failure causes and countermeasures of mechanical seal in cooling water pump for reactor[J]. Qilu Petrochemical Technology, 2019, 47(1): 51-52,56. doi: 10.3969/j.issn.1009-9859.2019.01.014
|
[8] |
赵远. 单端面机械密封故障的分析与处理分析[J]. 装备维修技术,2020,176(2): 119.
ZHAO Yuan. Fault analysis and treatment analysis of single end mechanical seal[J]. Equipment Maintenance Technology, 2020, 176(2): 119.
|
[9] |
尹源,黄伟峰,刘向锋,等. 机械密封智能化的技术基础和发展趋势[J]. 机械工程学报,2021,57(3): 116-128. doi: 10.3901/JME.2021.03.116
YIN Yuan, HUANG Weifeng, LIU Xiangfeng, et al. Technical Basis and Development Trend of Mechanical Seal Intelligence[J]. Journal of Mechanical Engineering, 2021, 57(3): 116-128. doi: 10.3901/JME.2021.03.116
|
[10] |
TOWSYFYAN H, GU F S, BALL A D, et al. Modelling acoustic emissions generated by tribological behaviour of mechanical seals for condition monitoring and fault detection[J]. Tribology International, 2018, 125: 46-58. doi: 10.1016/j.triboint.2018.04.021
|
[11] |
李晓晖. 流体动压型机械密封端面状态监测及健康评估[D]. 成都: 西南交通大学, 2016.
|
[12] |
张尔卿. 机械密封端面状态监测及寿命预测关键技术研究[D]. 成都: 西南交通大学, 2015.
|
[13] |
FAN Y E, GU F S, BALL A. A review of the condition monitoring of mechanical seals[C]//Proceedings of ASME 7th Biennial Conference on Engineering Systems Design and Analysis. Manchester: [s.n.], 2008: 179-184.
|
[14] |
LI R Y, HE D. Rotational machine health monitoring and fault detection using EMD-based acoustic emission feature quantification[J]. IEEE Transactions on Instrumentation and Measurement, 2012, 61(4): 990-1001. doi: 10.1109/TIM.2011.2179819
|
[15] |
HUANG W F, LIN Y B, LIU Y, et al. Face rub-impact monitoring of a dry gas seal using acoustic emission[J]. Tribology Letters, 2013, 52(2): 253-259. doi: 10.1007/s11249-013-0210-2
|
[16] |
刘忠,宋嘉城,邹淑云,等. 基于EMD的水轮机空化声发射信号阈值降噪方法[J]. 动力工程学报,2018,38(6): 501-507. doi: 10.3969/j.issn.1674-7607.2018.06.012
LIU Zhong, SONG Jiacheng, ZOU Shuyun, et al. Threshold denoising method for acoustic emission signals from hydraulic turbine cavitation based on EMD[J]. Journal of Chinese Society of Power Engineering, 2018, 38(6): 501-507. doi: 10.3969/j.issn.1674-7607.2018.06.012
|
[17] |
BACCAR D, SÖFFKER D. Wear detection by means of wavelet-based acoustic emission analysis[J]. Mechanical Systems and Signal Processing, 2015, 60/61: 198-207. doi: 10.1016/j.ymssp.2015.02.012
|
[18] |
BIANCHI D, MAYRHOFER E, GRÖSCHL M, et al. Wavelet packet transform for detection of single events in acoustic emission signals[J]. Mechanical Systems and Signal Processing, 2015, 64/65: 441-451. doi: 10.1016/j.ymssp.2015.04.014
|
[19] |
关山,彭昶. 刀具磨损声发射信号的混沌特性分析[J]. 农业工程学报,2015,31(11): 60-65. doi: 10.11975/j.issn.1002-6819.2015.11.009
GUAN Shan, PENG Chang. Chaotic characteristic analysis of tool wear acoustic emission signal[J]. Transactions of the Chinese Society of Agricultural Engineering, 2015, 31(11): 60-65. doi: 10.11975/j.issn.1002-6819.2015.11.009
|
[20] |
律方成,张瑜,董蒙,等. 基于复合评价指标的金属微粒声发射信号最优小波去噪及其特征提取[J]. 高压电器,2017,53(11): 1-8. doi: 10.13296/j.1001-1609.hva.2017.11.001
LV Fangcheng, ZHANG Yu, DONG Meng, et al. Optimal wavelet denoising and feature extraction of metal particulate acoustic emission signal based on composite evaluation index[J]. High Voltage Appliances, 2017, 53(11): 1-8. doi: 10.13296/j.1001-1609.hva.2017.11.001
|
[21] |
WEI N S, CHEN Z, GU F S, et al. Wavelet package denoising of acoustic emission signals for lubrication oil monitoring in engine systems[C]//Proceedings of IncoME-V & CEPE Net-2020. [S.l.]: Springer, 2021: 69-80.
|
[22] |
GILLES J. Empirical wavelet transform[J]. IEEE Transactions on Signal Processing, 2013, 61(16): 3999-4010. doi: 10.1109/TSP.2013.2265222
|
[23] |
JIANG W L, LI Z B, ZHANG S, et al. Hydraulic pump fault diagnosis method based on EWT decomposition denoising and deep learning on cloud platform[J]. Shock and Vibration, 2021, 2021(8): 6674351.1-6674351.18.
|
[24] |
张超,张晨,郭晓苏,等. 基于MCKD-EWT的滚动轴承故障诊断研究[J]. 轴承,2020(5): 43-48.
ZHANG Chao, ZHANG Chen, GUO Xiaosu, et al. Study on fault diagnosis for rolling bearings based on MCKD-EWT[J]. Bearing, 2020(5): 43-48.
|
[25] |
CHEGINI S N, BAGHERI A, NAJAFI F. Application of a new EWT-based denoising technique in bearing fault diagnosis[J]. Measurement, 2019, 144: 275-297. doi: 10.1016/j.measurement.2019.05.049
|
[26] |
王友仁,陈伟,孙灿飞,等. 基于能量聚集度经验小波变换的齿轮箱早期微弱故障诊断[J]. 中国机械工程,2017,28(12): 1484-1490. doi: 10.3969/j.issn.1004-132X.2017.12.015
WANG Youren, CHEN Wei, SUN Canfei, et al. Early weak fault diagnosis of gearboxes based on energy aggregation and EWT[J]. China Mechanical Engineering, 2017, 28(12): 1484-1490. doi: 10.3969/j.issn.1004-132X.2017.12.015
|
[27] |
张瑞,邓艾东,司晓东,等. 一种新的声发射信号消噪及故障诊断方法[J]. 振动与冲击,2018,37(4): 75-81.
ZHANG Rui, DENG Aidong, SI Xiaodong, et al. A new method for acoustic emission signal de-noised and fault diagnosis[J]. Journal of Vibration and Shock, 2018, 37(4): 75-81.
|
[28] |
DUBITZKY W, WOLKENHAUER O, CHO K, et al. Encyclopedia of systems biology[M]. New York: Springer, 2013.
|
[1] | WANG Min, CHEN Peng, LIU Yingfeng, FENG Gang, JIANG Yan. Numerical Simulation of In-plane Tensile Properties of High-Strength Steel Wire Mesh[J]. Journal of Southwest Jiaotong University, 2023, 58(2): 446-452. doi: 10.3969/j.issn.0258-2724.20210379 |
[2] | KANG Azhen, GU Yuhang, ZHANG Dongming, JIN Ke, ZHU Bing, YANG Bing. Wave Flume Test and Simplified Algorithm for Freak Wave Forces on a Dumbbell-Shaped Bridge Structure[J]. Journal of Southwest Jiaotong University, 2023, 58(5): 1017-1025. doi: 10.3969/j.issn.0258-2724.20210380 |
[3] | YANG Zhiying, HUANG Bo, DUAN Lunliang, ZHANG Jiawei, ZHU Bing. Study of Extreme Wave Forces on T-and Box-Types Girder[J]. Journal of Southwest Jiaotong University, 2021, 56(1): 12-19. doi: 10.3969/j.issn.0258-2724.20181069 |
[4] | QI Xin, YU Zhixiang, ZHANG Lijun, XU Hu, ZHAO Lei. Propagation Effect of Passive Flexible Protection System on Rockfall Impact[J]. Journal of Southwest Jiaotong University, 2020, 55(5): 1085-1093. doi: 10.3969/j.issn.0258-2724.20180442 |
[5] | DUAN Lunliang, ZHENG Dongsheng, WANG Shaohua, ZHANG Qibo. Numerical Study on Wave-Induced Oscillatory Liquefaction in Anisotropic Seabed[J]. Journal of Southwest Jiaotong University, 2019, 54(4): 741-747. doi: 10.3969/j.issn.0258-2724.20170810 |
[6] | DENG Shasha, SHEN Huoming, LIU Lang, TANG Huaiping. Research of Calculation Method for Wave Forces Acting on Large-Scale Bridge Piers Based on Diffraction Theory[J]. Journal of Southwest Jiaotong University, 2018, 53(2): 266-271. doi: 10.3969/j.issn.0258-2724.2018.02.007 |
[7] | WANG Yusuo, LI Junjie, LI Zhenghui, FENG Gaofei, WU Hao, HE Junnan. Assessment of Rockfall Impact Force by Particle Flow Code Numerical Simulation Based on Discrete Element Model[J]. Journal of Southwest Jiaotong University, 2016, 29(1): 22-29. doi: 10.3969/j.issn.0258-2724.2016.01.004 |
[8] | LIU Junqing, LI Qian, WANG Baoshi. Numerical Model of Asphalt Pavement Considering Transversely Isotropic Fatigue Damage[J]. Journal of Southwest Jiaotong University, 2016, 29(4): 670-676. doi: 10.3969/j.issn.0258-2724.2016.04.010 |
[9] | YE Kun, LI Renxian. Optimization Analysis of Height and Distance of Shelter Wind Wall for High-Speed Railway[J]. Journal of Southwest Jiaotong University, 2014, 27(2): 240-246. doi: 10.3969/j.issn.0258-2724.2014.02.009 |
[10] | KANG Azhen, ZHU Bing, XING Fan. Numerical Simulation of Wave Forces on Large Circular Cylinder Based on -Coordinate and Immersed Boundary Method[J]. Journal of Southwest Jiaotong University, 2014, 27(2): 233-239. doi: 10.3969/j.issn.0258-2724.2014.02.008 |
[11] | LUO Zhongxian, QIU Yanjun, YU Huaqian. Computational Modeling of Flows in Random Porous Media using Lattice Boltzmann Method[J]. Journal of Southwest Jiaotong University, 2014, 27(1): 93-96,118. doi: 10.3969/j.issn.0258-2724.2014.01.015 |
[12] | LIU Jiaxin, QIN Sicheng, XU Zhenyuan, ZHANG Ao, XI Yu, ZHANG Xuelin. Numerical Simulation of Heat Exchange Performance of Radiator Module in Construction Vehicles[J]. Journal of Southwest Jiaotong University, 2012, 25(4): 623-628. doi: 10.3969/j.issn.0258-2724.2012.04.014 |
[13] | , ZHU Bing, SONG Suidi, TAN Changjian, . Numerical Simulation for Diffraction around Large-Diameter Circular Cylinder Subjected to Three-Dimension Wave[J]. Journal of Southwest Jiaotong University, 2012, 25(2): 224-229. doi: 10.3969/j.issn.0258-2724.2012.02.009 |
[14] | CHEN Yun, CHEN Yun-Min, ZHOU Qun-Jian. Measurement of Water Content of Multiple Geomaterials by Time-Domain Reflectometry Technique[J]. Journal of Southwest Jiaotong University, 2011, 24(1): 42-48. doi: 10.3969/j.issn.0258-2724.2011.01.007 |
[15] | XU Rui, HUANG Dingfa, LI Chenggang, ZHOU Letao. Evaluation of GPS Baseline Quality Based on Probabilistic Characteristic of Integer Ambiguity[J]. Journal of Southwest Jiaotong University, 2007, 20(6): 739-742,747. |
[16] | TIAN Yun-de, QIN Shi-lun. ImprovedM ixed-M odeM ethod ofValid ElasticM odulus ofComposites[J]. Journal of Southwest Jiaotong University, 2005, 18(6): 783-787. |
[17] | XIAO Shi-guo, ZHOUDe-pei. Determination and Numerical Analysis Method of Relaxation Region for Cutting Slope[J]. Journal of Southwest Jiaotong University, 2003, 16(3): 318-322. |
[18] | LUOXiong, LI Yong-le. A Numerical Method for Simulating Natural Wind Based on Nonstationary Random Process[J]. Journal of Southwest Jiaotong University, 2002, 15(4): 367-370. |
[19] | TANWEIgang, TANGHUAIping. Adhesion Analysis of Wheel and Rail with Water Medium[J]. Journal of Southwest Jiaotong University, 2000, 13(5): 509-512. |
参数 | 数值 | 参数 | 数值 | |
vw/(km·h−1) | 2 | ζ | 1 | |
vc/(km·h−1) | 60 | k/(户·km−1) | 4600 | |
Td/(h·站−1) | 0.0125 | ˉg/(户·年−1) | 0.03 | |
CV/(人·车−1) | 1800 | σg/(户·年−1) | 0.01 | |
π1/(h·(km·d)−1) | 14256 + 475.2μ | pg | 0.7 | |
π2/ (h·(站·d)−1) | 7056 + 235.2μ | κ | 0.1 | |
π3/(h·km−1) | 2.2 | Tw/(d·年−1) | 250 | |
π4/(h·h−1) | 101 + 5 μ | N0/ 户 | 30000 | |
Tr/ (h·d−1) | 16 | μ/ (h·h−1) | 20 | |
η/ (次·(户·d) −1) | 4 |
风险决策 态度 | 公交线路/ 走廊长度/km | 平均站点 密度/ (站·km−1) | 平均发 车间隔/ min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | (6.81, 14.93) | 1.0937 | 2.1989 | 1.1280 | 0.8123 | 1.9403 | 0.0794 |
规避风险 | (6.81, 14.93) | 1.1773 | 2.1721 | 1.1261 | 0.8154 | 1.9414 | 0.0761 |
寻求风险 | (6.81, 14.93) | 0.9814 | 2.2252 | 1.152 | 0.8076 | 1.9595 | 0.1057 |
风险决策 态度 | 线路 长度/km | 走廊 长度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | 10.97 | (6.81, 14.93) | 1.1011 | 2.2979 | 1.208 | 0.8504 | 2.0584 | 0.1334 |
规避风险 | 11.40 | (6.81, 14.93) | 1.0976 | 2.3131 | 1.1813 | 0.8827 | 2.064 | 0.1169 |
寻求风险 | 10.48 | (6.81, 14.93) | 1.1099 | 2.2766 | 1.2499 | 0.8144 | 2.0644 | 0.1546 |
风险决策 态度 | 线路长度/km | 走廊长 度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 阶段 数/个 | 期待用户 通勤成本/ (h· (户·d)−1) | 期待运营 商成本/(h· (户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/(h· (户·d)−1) | |
初始 | 最终 | |||||||||
中立风险 | 6.82 | 12.31 | (6.81, 14.93) | 1.2080 | 2.0977 | 5 | 1.1539 | 0.7398 | 1.8936 | 0.1095 |
规避风险 | 6.82 | 12.76 | (6.81, 14.93) | 1.2186 | 2.0930 | 5 | 1.1420 | 0.7580 | 1.9000 | 0.0947 |
寻求风险 | 6.82 | 11.88 | (6.81, 14.93) | 1.2209 | 2.0916 | 5 | 1.1856 | 0.7147 | 1.9003 | 0.1352 |
参数 | 数值 | 参数 | 数值 | |
vw/(km·h−1) | 2 | ζ | 1 | |
vc/(km·h−1) | 60 | k/(户·km−1) | 4600 | |
Td/(h·站−1) | 0.0125 | ˉg/(户·年−1) | 0.03 | |
CV/(人·车−1) | 1800 | σg/(户·年−1) | 0.01 | |
π1/(h·(km·d)−1) | 14256 + 475.2μ | pg | 0.7 | |
π2/ (h·(站·d)−1) | 7056 + 235.2μ | κ | 0.1 | |
π3/(h·km−1) | 2.2 | Tw/(d·年−1) | 250 | |
π4/(h·h−1) | 101 + 5 μ | N0/ 户 | 30000 | |
Tr/ (h·d−1) | 16 | μ/ (h·h−1) | 20 | |
η/ (次·(户·d) −1) | 4 |
风险决策 态度 | 公交线路/ 走廊长度/km | 平均站点 密度/ (站·km−1) | 平均发 车间隔/ min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | (6.81, 14.93) | 1.0937 | 2.1989 | 1.1280 | 0.8123 | 1.9403 | 0.0794 |
规避风险 | (6.81, 14.93) | 1.1773 | 2.1721 | 1.1261 | 0.8154 | 1.9414 | 0.0761 |
寻求风险 | (6.81, 14.93) | 0.9814 | 2.2252 | 1.152 | 0.8076 | 1.9595 | 0.1057 |
风险决策 态度 | 线路 长度/km | 走廊 长度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 期待用户 通勤成本/ (h·(户·d)−1) | 期待运营 商成本/ (h·(户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/ (h·(户·d)−1) |
中立风险 | 10.97 | (6.81, 14.93) | 1.1011 | 2.2979 | 1.208 | 0.8504 | 2.0584 | 0.1334 |
规避风险 | 11.40 | (6.81, 14.93) | 1.0976 | 2.3131 | 1.1813 | 0.8827 | 2.064 | 0.1169 |
寻求风险 | 10.48 | (6.81, 14.93) | 1.1099 | 2.2766 | 1.2499 | 0.8144 | 2.0644 | 0.1546 |
风险决策 态度 | 线路长度/km | 走廊长 度/km | 平均站 点密度/ (站·km−1) | 平均发车 间隔/min | 阶段 数/个 | 期待用户 通勤成本/ (h· (户·d)−1) | 期待运营 商成本/(h· (户·d)−1) | 期待系统 成本/ (h·(户·d)−1) | 系统成本 标准差/(h· (户·d)−1) | |
初始 | 最终 | |||||||||
中立风险 | 6.82 | 12.31 | (6.81, 14.93) | 1.2080 | 2.0977 | 5 | 1.1539 | 0.7398 | 1.8936 | 0.1095 |
规避风险 | 6.82 | 12.76 | (6.81, 14.93) | 1.2186 | 2.0930 | 5 | 1.1420 | 0.7580 | 1.9000 | 0.0947 |
寻求风险 | 6.82 | 11.88 | (6.81, 14.93) | 1.2209 | 2.0916 | 5 | 1.1856 | 0.7147 | 1.9003 | 0.1352 |