CN104933474B - 危险化学品运输的模糊双层优化方法 - Google Patents
危险化学品运输的模糊双层优化方法 Download PDFInfo
- Publication number
- CN104933474B CN104933474B CN201510268940.3A CN201510268940A CN104933474B CN 104933474 B CN104933474 B CN 104933474B CN 201510268940 A CN201510268940 A CN 201510268940A CN 104933474 B CN104933474 B CN 104933474B
- Authority
- CN
- China
- Prior art keywords
- mrow
- mtd
- msub
- mtr
- mtable
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000000126 substance Substances 0.000 title claims abstract description 15
- 238000005457 optimization Methods 0.000 title claims abstract description 8
- 239000002245 particle Substances 0.000 claims description 16
- 238000004088 simulation Methods 0.000 claims description 10
- 238000004422 calculation algorithm Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 230000003044 adaptive effect Effects 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 230000008929 regeneration Effects 0.000 claims description 2
- 238000011069 regeneration method Methods 0.000 claims description 2
- 239000002355 dual-layer Substances 0.000 abstract description 4
- 230000002123 temporal effect Effects 0.000 abstract description 2
- 239000000383 hazardous chemical Substances 0.000 description 10
- 238000005516 engineering process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 208000015181 infectious disease Diseases 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种危险化学品运输的模糊双层优化方法，该方法包括以下步骤：确定运输网络中的风险参数；搜集仓库供给量与需求量信息；上层决策者制定客户分配方案；下层决策者寻找最优路径；根据上下层决策，得出最优路径。该方法采用双层规划充分反映所需的决策层次，利用其分层决策的特点，通过规划者和运输者的双方相互影响的决策，以求得满意解。本方法提出的模型具有有效性以及在策略和时间上的最优性。
Description
技术领域
本发明提供了一种处理危险化学品运输问题的模糊优化方法，属于危险化学品运输技术领域；特别是对多仓库车辆路径问题使用了双层优化模型，并考虑了危险化学品运输风险的模糊性。
背景技术
伴随中国经济的快速发展，危险化学品在工业、农业、国防和人类日常生活已经变成必不可少的材料。危险化学品具有爆炸、毒害、感染、腐蚀、放射性等特性。全世界每年化学品的产量巨大，化学品的种类繁多。在中国，95％以上的危险化学品涉及异地运输。危险化学品运输占年货运总量的30％以上，并呈上升趋势。长距离、大吨位已成为危险化学品道路运输的普遍状况。
危险化学品运输事故后果严重危害人民群众的生命、财产安全和环境安全，往往造成重大的经济损失和恶劣的政治影响。危险化学品运输事故是典型的“低概率-高后果”事件。因此，为了提高危险化学品运输的安全性、可靠性，应加强危险化学品运输的风险度量与分析，优化危险化学品运输的路径。
车辆路径问题是物流优化中的关键问题。车辆路径优化问题的调度方法成为降低风险、优化路径的关键因素。因此，风险的度量和车辆路径优化调度方法成为提高运输质量的有效保证。
对于风险的度量，高估或低估都将产生较大的负面影响。一方面，过高的风险估计会导致社会恐慌、运输成本上升或过分严格的运输制度。另一方面，过低的风险估计，会给公众的人身、财产以及环境带来巨大隐患。由于历史数据缺乏，测量方法选择以及测量精度差异的影响，风险具有一定的不确定性。虽然对于多仓库车辆路径问题的研究十分广泛[1]，但均没有考虑运输风险的模糊性。
危险化学品多仓库运输是指多个仓库使用一些车辆寻求服务多个客户的组合优化问题，在满足供给和需求约束条件下最小化风险、成本等。多仓库车辆路径问题比单仓库车辆路径问题更符合现实也更复杂。多仓库车辆路径问题旨在处理的问题有：(1)客户的分配，(2)每个仓库服务客户的路径选择。对于多仓库车辆路径问题研究者建立了不同的数学模型[2-4],但对于多层决策问题，均没有考虑双层规划。
为了解决这些问题，本发明充分考虑了运输风险的模糊性以体现风险的不确定性。同时，本发明使用了双层优化方法以充分反映所需的决策层次。基于以上问题，设计了一个模糊双层规划模型最小化危险化学品运输风险。同时，本发明使用了基于模糊模拟的粒子群算法以求得满意结果。
参考文献
[1]Jairo,R.M.,Julián,L.F.,Santiago,N.I.,Heriberto,F.J.,Nilson,H.,2015.A literature review on the vehicle routing problem with multipledepots.Computers &Industrial Engineering.79,115-129.
[2]Gulczynski,D.,Golden,B.L.,Wasil,E.,2011.The multi-depot splitdelivery vehicle routing problem:An integer programming-based heuristic,newtest problems,and computational results.Computers&Industrial Engineering.61(3),794-804.
[3]Sitek,P.,Wikarek,J.,Grzybowska,K.,2014.A multi-agent approach tothe multiechelon capacitated vehicle routing problem.Highlights of PracticalApplications of Heterogeneous Multi-Agent Systems.The PAAMS Collection,Communications in ComputerandInformation Science.430,121-132.
[4]Venkatasubbaiah,K.,Acharyulu,S.G.,ChandraMouli,K.V.V.,2011.Fuzzygoal programmingmethodfor solvingmulti-objectivetransportationproblems.Global Journal of Researchin Engineering.11(3).
发明内容
本发明的目的在于提供一种基于危险化学品运输的多仓库车辆路径模糊双层优化方法。基于可信性理论，本发明将风险模糊化，更客观地体现了不确定性。本发明设计了一种基于双层规划的均衡策略规划。该方法最能反映所需的决策层次。该方法的基本思想是：在运输网络中存在上层规划者和下层运输者，上层规划者对客户进行分配，下层运输者针对服务客户路径制定最优策略。通过双层的策略，以达到最小化风险的期望值的目的。首先，为了满足客户的需求，规划者预先制定了一些客户分配策略。其次，运输者对于每种策略制定出风险最小的运输路线。最终，通过双方的决策制定出满意策略。
为实现上述目的，本发明提出了一种危险化学品运输的模糊双层优化方法，该方法包括以下步骤：1、确定运输网络中的风险参数；2、搜集仓库供给量与需求量信息；3、上层决策者制定客户分配方案；4、下层决策者寻找最优路径；5、根据上下层决策，得出最优路径。
各个详细步骤如下：
符号系统：
I 仓库个数；
J 客户个数；
i 仓库，i＝1,2,…,I；
j 客户，j＝1,2,…,J；
Qi 仓库i的容量，i＝1,2,…,I；
Dj 客户j的需求，j＝1,2,…,J；
ki分配给仓库i的客户个数；
Ai 仓库i和其客户的集合，例：Ai＝{i,j1,j2,…,jki}；
从m到n的旅行风险，m,n∈Ai,i＝1,2,…,I,m≠n，决策变量；
xij 仓库i到仓库j的供给量，i＝1,2,…,I,j＝1,2,…,J，决策变量；
从m到n的弧是可行的，则取值为1，否则为0。
m、n表示仓库和其所服务的顾客
S1、确定运输网络中的风险参数。
风险是对事故发生率和结果的度量，这种度量方式是危险化学品运输不同于一般运输问题的关键。风险表达式如下，
Pmn表示从m到n的事故发生率，κmn表示事故影响区域，τmn表示平均人口密度。由于人口密度一般是不确定的，将τmn作为模糊变量，则风险也为模糊变量。
S2、搜集仓库供给量Qi和客户需求量的信息Dj。
S3、上层决策者制定客户分配方案。
j＝1,2,…,J和i＝1,2,…,I表示客户和仓库，随机从j＝1,2,…,J生成C1,C2,…,CJ，其需求D1,D2,…,DJ。计算前j个客户的累积需求Cumj。Q1,Q2,…,QI表示仓库的供给。客户的分配过程如下，定义仓库1和仓库2的分配过程。定义Uj＝Cumj-Q1，j＝1,2,…,J。找到Um的第一个非负值。分配C1,C2,…,Cm给仓库1。如果Um＞0，设Um作为客户m的需求且k＝m。否则k＝m+1。对于仓库2，重新计算Ck,Ck+1,…,CJ的累计值为Cumk,Cumk+1,…,CumJ。定义Uj＝Cumj-Q2，j＝k,k+1,…,J。重复以上过程直到所有客户分配给仓库。
以上的分配方案的模型表示如下：
目标最小化总体运输风险，第一个约束确保分配不超过仓库的供给能力，第二个约束保证分配满足客户的需求，第三个约束表示仓库i到客户j的供给量是非负的。
S4、下层决策者寻找最优路径。下层规划构建了一个旅行商问题。模型如下：
为求解下层规划，使用了模糊模拟、二分法和数值积分法求得风险的期望值。
S4.1模糊模拟
为求解下层规划，引入模糊模拟方法以模拟以下映射
U:x→E[f(x,ξ)]
其中f是实值函数。模糊向量ξ＝(ξ1,ξ2,…,ξm)存在一个联合可信函数v。随机生成向量y1,y2,…,yN并计算其可信值vk＝v(yk),k＝1,2,…,N。对于任何实数r，Cr{f(x,ξ)≥r}的可信值如下：
Cr{f(x,ξ)≤r}的可信值为如下：
α悲观值是满足L(r)＝Cr{f(x,ξ)≤r}≥α的最小值。由于L(r)是增函数，通过使用二分法计算最小值。
S4.2二分法
步骤1初始化一个足够小的正数ε；
步骤2随机生成yi并计算vi，i＝1,2,…,I；
步骤3计算最小值a＝min{f(x,yi)|1≤i≤I}和最大值
b＝max{f(x,yi)|1≤i≤I}；
步骤4设r＝(a+b)/2；
步骤5如果Cr{f(x,ξ)≤r}≥α，则b＝r。否则，a＝r；
步骤6如果b-a＞ε，进行第4步骤；
步骤7返回(a+b)/2作为α悲观值。
S4.3数值积分算法
步骤1初始化积分点数N；
步骤2设e＝0和n＝1；
步骤3通过使用平分法模拟悲观值，计算βi＝(ξi)inf(n/N)，1≤i≤m；
步骤4设e＝e+f(β1,β2,…,βm)/N和n＝n+1；
步骤5如果n≤N，继续步骤3；
步骤6返回模拟值e；
S5、根据上下层决策，得出最优路径。
通过上下层决策相互作用，使用粒子群算法以求得最优解。
S5.1混合粒子群优化：
步骤1随机生成初始种群；
步骤2通过模糊模拟算法计算每个粒子的适应值；
步骤3通过选择、交叉、变异更新种群；
步骤4重复步骤2-步骤5最大代数次；
步骤5返回目标值最小的粒子作为最优解。
与现有技术相比，本发明具有如下有益效果。
该方法采用双层规划充分反映所需的决策层次，利用其分层决策的特点，通过规划者和运输者的双方相互影响的决策，以求得满意解。现有的基于危险化学品多仓库车辆路径问题。没有考虑运输过程中风险的不确定，本方法基于可信理论充分考虑了不确定性。利用模糊期望值模型，解决了不确定性条件下风险的度量。并设计了基于模糊模拟的粒子群算法求解最优策略满足了客户的需求，而且最小化风险。实验结果表明提出模型的有效性以及在策略和时间上的最优性。
附图说明
图1是三角模糊数的可信函数
图2是多仓库车辆路径问题的例图
图3是混合粒子群算法的流程图
具体实施方式
以下结合附图和具体实例对本发明内容作详细说明。
在以10个仓库和20个客户的车辆路径问题为例，对本发明进行详细分析。首先，搜集关于运输图中的风险信息如：事故发生率，影响区域和平均人口密度，根据以上信息设置不同路段的风险模糊值。假定运输风险是三角模糊变量。搜集历史信息设置仓库和客户的供给量以及需求量。然后，将以上搜集的信息通知给上层规划者和下层运输者，上层规划者根据客户分配方案对客户进行分配，满足了模型中的约束。
根据以上的分配结果作为初始化的粒子群，然后对每个粒子进行评估即得出每个粒子的风险值。模糊模拟的过程与评估相结合，具体过程：对三角模糊风险值进行模糊模拟，将其表示为一个矩阵，矩阵中的值都是在相应的三角模糊值中取得，采用了200个样本点即200个矩阵，根据模糊模拟的过程得出每个粒子的评估值。对粒子进行选择，交叉和变异的过程，最后得出满意解。置信水平为0.98。混合粒子群的最大迭代次数为200，粒子数设为20。得出的最优路径表示在表1。
表1
Claims (1)
1.一种危险化学品运输的模糊双层优化方法，其特征在于：
符号系统：
I仓库个数；
J客户个数；
i仓库，i＝1,2,…,I；
j客户，j＝1,2,…,J；
Qi仓库i的容量，i＝1,2,…,I；
Dj客户j的需求，j＝1,2,…,J；
ki分配给仓库i的客户个数；
Ai仓库i和其客户的集合，例：Ai＝{i,j1,j2,…,jki}；
从m到n的旅行风险，m,n∈Ai,i＝1,2,…,I,m≠n，决策变量；
xij仓库i到仓库j的供给量，i＝1,2,…,I,j＝1,2,…,J，决策变量；
从m到n的弧是可行的，则取值为1，否则为0；
m、n表示仓库和其所服务的顾客；
该方法的实现步骤如下：
S1、确定运输网络中的风险参数；
风险是对事故发生率和结果的度量，这种度量方式是危险化学品运输不同于一般运输问题的关键；风险表达式如下，
<mrow>
<msubsup>
<mi>R</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>&kappa;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>&times;</mo>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
</msub>
</mrow>
Pmn表示从m到n的事故发生率，κmn表示事故影响区域，τmn表示平均人口密度；由于人口密度一般是不确定的，将τmn作为模糊变量，则风险也为模糊变量；
S2、搜集仓库供给量Qi和客户需求量的信息Dj；
S3、上层决策者制定客户分配方案；
j＝1,2,…,J和i＝1,2,…,I表示客户和仓库，随机从j＝1,2,…,J生成C1,C2,…,CJ，其需求D1,D2,…,DJ；计算前j个客户的累积需求Cumj；Q1,Q2,…,QI表示仓库的供给；客户的分配过程如下，定义仓库1和仓库2的分配过程；定义Uj＝Cumj-Q1，j＝1,2,…,J；找到第一个非负值Um；分配C1,C2,…,Cm给仓库1；如果Um＞0，设Um作为客户m的需求且k＝m；否则k＝m+1；对于仓库2，重新计算Ck,Ck+1,…,CJ的累计值为Cumk,Cumk+1,…,CumJ；定义Uj＝Cumj-Q2，j＝k,k+1,…,J；重复以上过程直到所有客户分配给仓库；
以上的分配方案的模型表示如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>J</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&le;</mo>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>I</mi>
</munderover>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>D</mi>
<mi>j</mi>
</msub>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>J</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>I</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>J</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
</mtable>
</mfenced>
目标最小化总体运输风险，第一个约束确保分配不超过仓库的供给能力，第二个约束保证分配满足客户的需求，第三个约束表示仓库i到客户j的供给量是非负的；
S4、下层决策者寻找最优路径；下层规划构建了一个旅行商问题；模型如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>i</mi>
<mi>J</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>min</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>E</mi>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>&Element;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>&NotEqual;</mo>
<mi>n</mi>
</mrow>
</munder>
<msubsup>
<mi>R</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>&Element;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>&NotEqual;</mo>
<mi>n</mi>
</mrow>
</munder>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>n</mi>
<mo>&Element;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>m</mi>
<mo>&NotEqual;</mo>
<mi>n</mi>
</mrow>
</munder>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>&Element;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mtable>
<mtr>
<mtd>
<mrow>
<munder>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>,</mo>
<mi>n</mi>
<mo>&Element;</mo>
<mi>S</mi>
</mrow>
</munder>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>&le;</mo>
<mo>|</mo>
<mi>S</mi>
<mo>|</mo>
<mo>-</mo>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>S</mi>
<mo>&Subset;</mo>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mn>2</mn>
<mo>&le;</mo>
<mo>|</mo>
<mi>S</mi>
<mo>|</mo>
<mo>&le;</mo>
<msub>
<mi>k</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msubsup>
<mi>Z</mi>
<mrow>
<mi>m</mi>
<mi>n</mi>
</mrow>
<mi>i</mi>
</msubsup>
<mo>&Element;</mo>
<mo>{</mo>
<mn>0</mn>
<mo>,</mo>
<mn>1</mn>
<mo>}</mo>
<mo>.</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
为求解下层规划，使用了模糊模拟、二分法和数值积分法求得风险的期望值；
S4.1模糊模拟
为求解下层规划，引入模糊模拟方法以模拟以下映射
U:x→E[f(x,ξ)]
其中f是实值函数；模糊向量ξ＝(ξ1,ξ2,…,ξm)存在一个联合可信函数v；随机生成向量y1,y2,…,yN并计算其可信值vk＝v(yk),k＝1,2,…,N；对于任何实数r，Cr{f(x,ξ)≥r}的可信值如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo><</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo><</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0.5</mn>
<mo>,</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
Cr{f(x,ξ)≤r}的可信值为如下：
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo><</mo>
<mn>0.5</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>></mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>max</mi>
<mrow>
<mo>{</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>k</mi>
</msub>
<mo>|</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<mrow>
<mi>x</mi>
<mo>,</mo>
<msub>
<mi>y</mi>
<mi>k</mi>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<mi>r</mi>
</mrow>
<mo>}</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<mn>0.5.</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
α悲观值是满足L(r)＝Cr{f(x,ξ)≤r}≥α的最小值；由于L(r)是增函数，通过使用二分法计算最小值；
S4.2二分法
步骤4.2.1初始化一个足够小的正数ε；
步骤4.2.2随机生成yi并计算vi，i＝1,2,…,I；
步骤4.2.3计算最小值a＝min{f(x,yi)1≤i≤I}和最大值
b＝max{f(x,yi)|1≤i≤I}；
步骤4.2.4设r＝(a+b)/2；
步骤4.2.5如果Cr{f(x,ξ)≤r}≥α，则b＝r；否则，a＝r；
步骤4.2.6如果b-a＞ε，进行步骤4.2.4；
步骤4.2.7返回(a+b)/2作为α悲观值；
S5、根据上下层决策，得出最优路径；
通过上下层决策相互作用，使用粒子群算法以求得最优解；
S5.1混合粒子群优化：
步骤5.1.1随机生成初始种群；
步骤5.1.2通过模糊模拟算法计算每个粒子的适应值；
步骤5.1.3通过选择、交叉、变异更新种群；
步骤5.1.4重复步骤5.1.2-步骤5.1.5最大代数次；
步骤5.1.5返回目标值最小的粒子作为最优解；
最终得出最优路径。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510268940.3A CN104933474B (zh) | 2015-05-24 | 2015-05-24 | 危险化学品运输的模糊双层优化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510268940.3A CN104933474B (zh) | 2015-05-24 | 2015-05-24 | 危险化学品运输的模糊双层优化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104933474A CN104933474A (zh) | 2015-09-23 |
| CN104933474B true CN104933474B (zh) | 2018-04-27 |
Family
ID=54120634
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510268940.3A Active CN104933474B (zh) | 2015-05-24 | 2015-05-24 | 危险化学品运输的模糊双层优化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104933474B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107451693A (zh) * | 2017-08-02 | 2017-12-08 | 南京工业大学 | 多点多目标的危化品运输路径优化方法 |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110533238B (zh) * | 2019-08-22 | 2022-08-26 | 杭州电子科技大学 | 二型模糊环境下的危化品车辆路径规划方法 |
| CN111091329B (zh) * | 2019-12-18 | 2022-12-16 | 北京化工大学 | 一种危化品多车型运输的半开放式车辆路径优化方法 |
| CN112270517B (zh) * | 2020-11-17 | 2022-05-10 | 北京化工大学 | 一种基于变权风险的危险化学品多仓库车辆路径优化方法 |
| CN113139774B (zh) * | 2021-04-25 | 2023-07-11 | 广州大学 | 一种面向多仓库运输的车辆路径优化方法 |
| CN116342197B (zh) * | 2023-01-05 | 2023-11-14 | 上海朗晖慧科技术有限公司 | 一种结合地理位置信息的实时价格方案生成方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103761588A (zh) * | 2014-02-18 | 2014-04-30 | 张家港美核电子科技有限公司 | 基于多目标建模优化的危化品运输调度方法 |
| CN104463379A (zh) * | 2014-12-31 | 2015-03-25 | 广东易富网络科技有限公司 | 一种带时变需求的关联物流运输优化调度方法 |
| CN104537446A (zh) * | 2015-02-09 | 2015-04-22 | 四川大学 | 二层带模糊随机时间窗口车辆路径优化方法 |
-
2015
- 2015-05-24 CN CN201510268940.3A patent/CN104933474B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103761588A (zh) * | 2014-02-18 | 2014-04-30 | 张家港美核电子科技有限公司 | 基于多目标建模优化的危化品运输调度方法 |
| CN104463379A (zh) * | 2014-12-31 | 2015-03-25 | 广东易富网络科技有限公司 | 一种带时变需求的关联物流运输优化调度方法 |
| CN104537446A (zh) * | 2015-02-09 | 2015-04-22 | 四川大学 | 二层带模糊随机时间窗口车辆路径优化方法 |
Non-Patent Citations (1)
| Title |
|---|
| Exact and meta-heuristic approach for a general heterogeneous dial-a-ride problem with multiple depots;Kris Braekers，An Caris，Gerrit K.Janssens;《transpotation research part B:methodological》;20140930;全文 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107451693A (zh) * | 2017-08-02 | 2017-12-08 | 南京工业大学 | 多点多目标的危化品运输路径优化方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104933474A (zh) | 2015-09-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104933474B (zh) | 危险化学品运输的模糊双层优化方法 | |
| Tian et al. | Multi-objective scheduling of priority-based rescue vehicles to extinguish forest fires using a multi-objective discrete gravitational search algorithm | |
| Ahmed et al. | Knowledge graph based trajectory outlier detection in sustainable smart cities | |
| Duddu et al. | Principle of demographic gravitation to estimate annual average daily traffic: Comparison of statistical and neural network models | |
| Xu et al. | Hazmats transportation network design model with emergency response under complex fuzzy environment | |
| Jootoo et al. | Bridge type classification: Supervised learning on a modified NBI data set | |
| CN112270517B (zh) | 一种基于变权风险的危险化学品多仓库车辆路径优化方法 | |
| Ma et al. | A multiobjective route robust optimization model and algorithm for hazmat transportation | |
| Yedavalli et al. | Microsimulation analysis for network traffic assignment (MANTA) at metropolitan-scale for agile transportation planning | |
| Lee et al. | Dynamic OD estimation using three phase traffic flow theory | |
| Guzmáan et al. | Fuzzy models and resolution methods for covering location problems: an annotated bibliography | |
| Zhang et al. | A distributionally robust optimisation model for last mile relief network under mixed transport | |
| Wong et al. | Estimation of multiclass origin-destination matrices from traffic counts | |
| Sirbiladze et al. | Uncertainty modeling in multi-objective vehicle routing problem under extreme environment | |
| Xu et al. | Multi-objective bilevel construction material transportation scheduling in large-scale construction projects under a fuzzy random environment | |
| Zhang et al. | Bilevel optimization for the hazmat transportation problem with lane reservation | |
| Shaw | GIS for urban travel demand analysis: requirements and alternatives | |
| Jiang et al. | Nonlinear interval optimization for uncertain problems | |
| CN112926800B (zh) | 一种考虑复杂路网的物流配送区域划分方法 | |
| Guo et al. | Max-flow rate priority algorithm for evacuation route planning | |
| Khaleghi et al. | Robust sustainable multi-period hub location considering uncertain time-dependent demand | |
| Lu et al. | A geographically partitioned cellular automata model for the expansion of residential areas | |
| Huang et al. | Modelling uncapacitated facility location problem with uncertain customers’ positions | |
| Huo et al. | [Retracted] Research on Solving Postdisaster Material Distribution and Scheduling with Improved NSGA‐II Algorithm | |
| CN116703008A (zh) | 一种用于新建公路的交通出行量预测方法、设备及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |