Processing math: 0%

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

公路客运门到门运营模式下的末端线网优化

杨忠振 王文娣 TheoNOTTEBOOM

杨忠振, 王文娣, TheoNOTTEBOOM. 公路客运门到门运营模式下的末端线网优化[J]. 交通运输工程学报, 2017, 17(1): 119-128.
引用本文: 杨忠振, 王文娣, TheoNOTTEBOOM. 公路客运门到门运营模式下的末端线网优化[J]. 交通运输工程学报, 2017, 17(1): 119-128.
YANG Zhong-zhen, WANG Wen-di, Theo NOTTEBOOM. Optimization of end line network with door to door operation mode of highway passenger transport[J]. Journal of Traffic and Transportation Engineering, 2017, 17(1): 119-128.
Citation: YANG Zhong-zhen, WANG Wen-di, Theo NOTTEBOOM. Optimization of end line network with door to door operation mode of highway passenger transport[J]. Journal of Traffic and Transportation Engineering, 2017, 17(1): 119-128.

公路客运门到门运营模式下的末端线网优化

基金项目: 

国家自然科学基金项目 71431001

中央高校基本科研业务费专项资金项目 3132016303

详细信息
    作者简介:

    杨忠振(1964-), 男, 辽宁凌海人, 大连海事大学教授, 工学博士, 从事交通运输规划与管理研究

  • 中图分类号: U492.431

Optimization of end line network with door to door operation mode of highway passenger transport

More Information
Article Text (Baidu Translation)
  • 摘要: 为解决高铁时代公路客运难的问题, 提出门到门的公路客运运营模式, 并构建了双层规划模型, 上层模型以公路客运服务商的日运营总利润最大为目标, 优化门到门模式下公路客运的末端线网, 下层模型计算均衡状态下交通网络上的客流量, 最后基于大连至沈阳运输区段的实际数据进行了案例分析。计算结果表明: 公路客运的日运营总利润随末端线路数的增加先升后降, 运营4条末端线路、日发车25班最佳; 高铁开通后原有运营模式下公路客运服务商的日运营总利润为7 270元, 票价为99元, 日均客流量为130人次; 门到门运营模式下公路客运服务商的日运营总利润为37 755元, 票价区间为[88, 105]元, 日均客流量将达到778人次。可见, 公路客运门到门运营模式可以有效提高公路客运的市场份额和运营利润, 在保证公路客运服务商利润的前提下, 大幅度提高整体客运服务水平。

     

  • 空中交通流量的增大导致了空中交通网络的拥挤, 为解决飞行冲突日益严峻的问题, 缩小飞行间隔是重要手段之一。为确保飞行安全, 有必要对空中飞行碰撞问题进行研究。国外从20世纪60年代就开始着手空中飞行碰撞风险研究工作, 其中以对侧向碰撞风险的研究居多。1964年, Filkins等用概率论理论分析了侧向碰撞风险问题[1]; 1965年, 北大西洋规划小组运用Reich模型对缩小侧向间隔标准进行了安全评估分析; 1994年国际民航组织间隔总概念专家组提交了基于自动相关监视的侧向间隔的危险碰撞模型, 计算出整个侧向重叠概率。

    国内对于飞行碰撞风险的研究较少, 2001年, 王欣、徐肖豪用Filkins的概率论方法在不考虑纵向和垂直方向的情况下, 分析了侧向碰撞风险问题[2], 但只求得了碰撞概率, 而非碰撞次数; 2002年, 应爱玲、徐肖豪运用Reich模型研究了洋区飞行的侧向碰撞风险问题[3], 但只适应于洋区飞行, 没有深入探讨国内非洋区飞行的碰撞风险问题; 2005年, 朱代武采用速度矢量三角分析法分析了航空器在低空空域飞行中发生相撞的问题, 提出了避让飞行冲突的两种解决方法[4]

    关于碰撞问题研究, 在公路领域也有涉及。在2002年, 徐杰等通过对一系列外部因素影响下驾驶员行为与跟随车与引导车行驶一致性的研究, 寻求可避免发生追尾碰撞事故的安全距离恰当值[5]; 2004年, 侯志祥等应用自适应神经模糊推理系统, 建立了高速公路汽车追尾的自适应神经模糊推理系统概率模型, 计算出在不同车速差和行车间距时的高速公路汽车追尾概率[6]

    本文对非洋区平行航路基于VOR导航的空中飞行的侧向碰撞率进行了研究, 以Reich模型为基础, 对Reich模型进行了改进, 推导出了平行航路侧向碰撞率的计算公式, 然后对侧向碰撞率的决定因素分别进行了分析, 最后用实例进行了仿真计算。

    1.1.1   碰撞模板

    设航路中每架飞机的尺寸是相同的, 飞机的长度、翼展、高度分别为λxλyλz, 以飞机A为中心, 在纵向、横向、垂向分别以2λx、2λy、2λz虚拟出一个长方体区域(设飞机B为一个质点), 见图 1中内部盒体, 即为碰撞模板。当飞机B碰撞飞机A时, 飞机B的重心正好在飞机A碰撞模板的边缘上, 故可将飞机B视为一个微粒, 则碰撞过程可被认为是微粒B“轰炸”飞机A的模板。

    图  1  碰撞模板和临近层
    Figure  1.  Collision slab and proximity shell
    1.1.2   临近层

    SxSySz分别为纵向、侧向和垂向的间隔标准, 以飞机A为中心做一长、宽、高分别为2Sx、2Sy、2Sz的长方体, 见图 1中外部盒体附近, 即为临近层。当飞机B在临近层外时, 两机之间被视为无碰撞风险; 当飞机B十分靠近该层时, 碰撞风险就会上升; 当B越过临近层边缘, 位于临近层内部时, 表示发生了危险接近, 并且有碰撞风险。

    Hsu在研究交叉航线碰撞概率估算中曾提出将飞机用圆柱体代替[8]。以圆柱体碰撞模板代替传统的长方体碰撞模板, 更为合理且计算准确, 因为当飞机发生航线偏离时, 若按图 2(a)处理, 偏离航线时的两机距离和相撞判据将很难求得; 若按照图 2(b)处理, 不管如何偏离航线, 两机距离和相撞判据均不受影响[9]

    图  2  碰撞模板
    Figure  2.  Collision slabs

    由于中国航路内大部分运输机的长度和翼展相差不大, 故假设翼展和长度相等, 均用飞机长度D表示, 飞机高度用λz表示, 则Reich模型中的碰撞模板可改进如下: 以飞机A为中心, 分别以2D为直径, 2λz为高虚拟出一个圆柱体区域, 当飞机B正好在飞机A圆柱体的边缘上时, 两架飞机就发生了机体接触, 相当于飞机B与飞机A进行绝对碰撞, 见图 3

    图  3  改进的碰撞模板
    Figure  3.  Improved collision slab

    碰撞率R即单位时间内的碰撞次数, 等于单位时间所预计的微粒B进入飞机A模板的次数。设两机纵向距离小于D的频率, 即两机纵向发生重叠的频率为Nx, 两机纵向距离小于D的概率, 即两机纵向发生重叠的概率为Px, 飞机B纵向穿过A的碰撞模板的平均时间为tx, A、B两机在纵向上的相对速度为VxNyNzPyPztytzVyVz定义均同纵向各参数。

    因在临近层之外不存在飞行冲突, 故以上的频率NxPx概率均是指在临近层内飞行, 其值分别等于B在临近层内与A纵向距离小于D的次数与B纵向穿越临近层所需时间的比值和飞机B纵穿飞机A碰撞模板的时间与其纵向穿越临近层所需时间的比值。由定义知

    tx=2DVx,ty=2DVy

    当2架飞机在平行航路上巡航飞行时, 某方向发生重叠是指两机中心之间的距离小于该方向上的间隔标准值, 它并不对两机另外两方向上的距离产生影响, 因此, 飞机在侧向、纵向和垂向上的重叠是相互独立的。只有当两机在3个方向上同时发生重叠时, 2架飞机才会发生碰撞。故纵向碰撞概率即为3个方向同时发生重叠的概率PyPz, 则纵向碰撞率为NxPyPz。同理可得, 侧向和垂向的碰撞率分别为NyPxPzNzPxPy。故总碰撞率等于3个方向碰撞率之和, 即

    R=ΝxΡyΡz+ΝyΡxΡz+ΝzΡxΡy=ΡxΡyΡz(1tx+1ty)+ΝzΡxΡy(3)

    假设飞机以匀速的巡航速度飞行在高空区, 飞机处在盲飞状态, 即飞行员和管制员对其管辖空域内的飞机即将发生的危险情况失去正确判断的能力。考虑侧向碰撞率, 则垂向间隔记为0, 侧向间隔标准记为Sy。设沿航路飞行时侧向重叠概率为Py(Sy), 2架飞机垂向发生重叠的概率为Pz(0), 飞机长度为D, 则侧向碰撞率Cy

    Cy=ΡxΡy(Sy)Ρz(0)(1tx+1ty)+Νz(0)ΡxΡy(Sy)=ΡxΡy(Sy)Ρz(0)(Vx+Vy2D)+Νz(0)ΡxΡy(Sy)(4)

    根据参考文献[10], 取Pz(0)值为0.39, 得

    Cy=0.3912SxΡy(Sy)(Vx+Vy)+DSxΡy(Sy)Νz(0)(5)

    式中: Sx为已给定的间隔标准; VxVy可由飞机机型、巡航速度及偏航角度计算得出; D可由飞机的技术参数得知; Nz(0)可由历史数据分析求得; 故式(5)主要考虑Py(Sy)的计算。

    通过对雷达显示器进行拍摄, 测量飞机偏离航路的距离以及飞机到所使用的导航台的距离来研究飞机侧向偏航距离的分布规律。根据参考文献[1], 通过分析知航路上的每一个点, 在其上飞机的偏航距离都是服从正态分布的随机变量。设μ0σ0分别是当飞机通过VOR台上空时的偏航距离平均值和标准差, μρσρ分别是当飞机到VOR台的距离为ρ时飞机的偏航距离的平均值和标准差, bμbσ是各自的增长系数, 根据参考文献[1], 有

    μρ=μ0+bμρ,σρ=σ0+bσρ(6)

    设飞机A、B的侧向偏航分别为AB, 则

    {A~Ν[μ01+bμ1ρ1,(σ01+bσ1ρ1)2]B~Ν[μ02+bμ2ρ2,(σ02+bσ2ρ2)2](7)

    式中: μ01μ02分别为飞机A、B通过VOR台上空时的偏航均值的初始值; σ01σ02分别为飞机A、B通过VOR台上空时偏航均方差的初始值; bμ1bμ2分别为飞机A、B偏航均值随距离的增长率; bσ1bσ2分别为飞机A、B偏航均方差随距离的增长率; ρ1ρ2分别为飞机A、B距VOR台的距离。

    根据参考文献[2], 令μ01=μ02=μ0, bμ1=bμ2=bμ, σ01=σ02=σ0, bσ1=bσ2=bσ。设d为2条平行航路间的水平距离, L为2架飞机的距离, 两VOR台的水平位置相同, 以航路1为基准线(图 4), 则飞机A和B的侧向距离为

    L=(d+B)-A=d+(B-A)(8)

    图  4  基于VOR导航的平行航路
    Figure  4.  Parallel routes based on VOR navigation

    由以上分析可知

    L~N[d+(μ02-μ01)+bμ2ρ2-bμ1ρ1, (σ02+bσ2ρ2)2+(σ01+bσ1ρ1)2]~N[d+bμ(ρ2-ρ1), (σ0+bσρ2)2+(σ0+bσρ1)2] (9)

    L < D时, 两机发生侧向重叠, 侧向重叠概率Py(Sy)为L在(-D, D)内的积分。在已知航路距离dbμbσσ0及飞机A、B距导航台距离时, 可由式(9)得出侧向重叠概率。

    设飞机在航路飞行阶段, Sx为3 km; 以波音777为例, 机身长度D为63.73 m; 航路距离d为25 km, ρ1为50 km, ρ2为30 km; 参考文献[2], 取bμ为0.015, bσ为0.02, σ0为2.4 km, 则

    Py(Sy)=∫-0.063 730.063 73 1σ2πe-(x-μ)22σ2dx=4.0424×10-9

    即侧向重叠概率为4.042 4×10-9

    由飞机技术参数可得, 波音777的巡航速度为934.58 km·h-1, 当两机偏航角度为3°和5°时, 两机在纵向上的相对速度Vx和侧向上的相对速度Vy分别为2.17和124 km·h-1。参考文献[3], 根据北大西洋规划小组对北大西洋区域的雷达数据分析数据取Nz(0)为60架·h-1, 则

    Cy=0.3912SxΡy(Sy)(Vx+Vy)+DSxΡy(Sy)Νz(0)=4.0424×10-9×[0.39×2.17+1242×3+0.063733×60]=3.8304×10-8

    即在以上给定的条件下, 每108飞行小时内的可能碰撞次数为3.830 4次。中国民航总局的规定是飞机在空中发生碰撞的可能性不大于10-6级别[11], 而3.830 4×10-8小于10-6, 由此可知所给参数满足安全要求。

    针对平行航路上的侧向间隔安全性评估问题, 本文将Reich模型和概率论理论相结合, 给出了非洋区VOR导航下的两机侧向碰撞率的计算方法。该方法可以用于平行航路间隔的安全评估、基于VOR导航的平行航路纵向和垂向碰撞率的研究与非洋区基于CNS性能的平行航路碰撞率分析。

  • 图  1  城际客运网络

    Figure  1.  Intercity passenger transport network

    图  2  放射型线网结构

    Figure  2.  Radial line network structures

    图  3  公路客运末端线网的走行方向和站点设置

    Figure  3.  Travel directions and site setting of end line network for highway passenger transport

    图  4  研究区域

    Figure  4.  Study area

    图  5  公路客运末端线网的优化方案

    Figure  5.  Optimized schemes of end line network for highway passenger transport

    图  6  不同F值下的线网日运营总利润

    Figure  6.  Total daily operating profits of line networks with different F values

    表  1  各小区前往沈阳的日均客流量

    Table  1.   Daily passenger flow to Shenyang for each area

    下载: 导出CSV

    表  2  公路客运末端线网的优化结果

    Table  2.   Optimized result of end line network for highway passenger transport

    下载: 导出CSV

    表  3  不同F值下的日发车班数

    Table  3.   Daily departure numbers with different F values

    下载: 导出CSV

    表  4  大连某公路客运公司三种运营情况下各项指标的变化

    Table  4.   Index variations of one highway passenger transport company in Dalian under three operation situations

    下载: 导出CSV
  • [1] LAMPKIN W, SAALMANS P D. The design of routes, service frequencies, and schedules for a municipal bus undertaking: a case study[J]. Operational Research Quarterly, 1967, 18 (4): 375-397. doi: 10.1057/jors.1967.70
    [2] HIRSCH W M, DANTZIG G B. The fixed charge problem[J]. Naval Research Logistics Quarterly, 1968, 15 (3): 413-424. doi: 10.1002/nav.3800150306
    [3] SILMAN L A, BARZILY Z, PASSY U. Planning the route system for urban buses[J]. Computers and Operations Research, 1974, 1 (2): 201-211. doi: 10.1016/0305-0548(74)90046-X
    [4] BAAJ M H, MAHMASSANI H S. An AI-based approach for transit route system planning and design[J]. Journal of Advanced Transportation, 1991, 25 (2): 187-209. doi: 10.1002/atr.5670250205
    [5] XIONG Yi-hua, SCHNEIDER J B. Transportation network design using a cumulative genetic algorithm and neural network[J]. Transportation Research Record, 1992 (1364): 37-44.
    [6] PATTNAIK S B, MOHAN S, TOM V M. Urban bus transit route network design using genetic algorithm[J]. Journal of Transportation Engineering, 1998, 124 (4): 368-375. doi: 10.1061/(ASCE)0733-947X(1998)124:4(368)
    [7] BIELLI M, CARAMIA M, CAROTENUTO P. Genetic algorithms in bus network optimization[J]. Transportation Research Part C: Emerging Technologies, 2002, 10 (1): 19-34. doi: 10.1016/S0968-090X(00)00048-6
    [8] 张启人, 熊桂林. 公共交通大系统建模与优化[J]. 系统工程, 1986, 4 (6): 25-39. https://www.cnki.com.cn/Article/CJFDTOTAL-GCXT198606004.htm

    ZHANG Qi-ren, XIONG Gui-lin. Modeling and optimization of large public transport system[J]. Systems Engineering, 1986, 4 (6): 25-39. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-GCXT198606004.htm
    [9] 吴稼豪. 城市公共交通车线网络优化模型[J]. 上海机械学院学报, 1987, 9 (2): 7-13. https://www.cnki.com.cn/Article/CJFDTOTAL-HDGY198702003.htm

    WU Jia-hao. A model for the urban transit network optimization and evaluation[J]. Journal of Shanghai Institute of Mechanical Engineering, 1987, 9 (2): 7-13. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HDGY198702003.htm
    [10] 林柏梁, 杨富社, 李鹏. 基于出行费用最小化的公交网络优化模型[J]. 中国公路学报, 1999, 12 (1): 79-83. doi: 10.3321/j.issn:1001-7372.1999.01.012

    LIN Bo-liang, YANG Fu-she, LI Peng. Designing optimal bus network for minimizing trip times of passenger flows[J]. China Journal of Highway and Transport, 1999, 12 (1): 79-83. (in Chinese). doi: 10.3321/j.issn:1001-7372.1999.01.012
    [11] 韩印, 李维斌, 李晓峰. 城市公交线网调整优化PSO算法[J]. 中国公路学报, 1999, 12 (3): 100-104. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL199903014.htm

    HAN Yin, LI Wei-bin, LI Xiao-feng. The algorithm PSO of adjustment and optimization for public traffic network[J]. China Journal of Highway and Transport, 1999, 12 (3): 100-104. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL199903014.htm
    [12] 单连龙, 高自友. 城市公交系统连续平衡网络设计的双层规划模型及求解算法[J]. 系统工程理论与实践, 2000, 20 (7): 85-93. https://www.cnki.com.cn/Article/CJFDTOTAL-XTLL200007017.htm

    SHAN Lian-long, GAO Zi-you. A bilevel programming model for continuous equilibrium network design and its solution algorithm for urban transit system[J]. Systems Engineering Theory and Practice, 2000, 20 (7): 85-93. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-XTLL200007017.htm
    [13] 冯树民, 陈洪仁. 公共交通线网优化研究[J]. 哈尔滨工业大学学报, 2005, 37 (5): 691-693. doi: 10.3321/j.issn:0367-6234.2005.05.033

    FENG Shu-min, CHEN Hong-ren. Study of public transit network optimization method[J]. Journal of Harbin Institute of Technology, 2005, 37 (5): 691-693. (in Chinese). doi: 10.3321/j.issn:0367-6234.2005.05.033
    [14] 汤可夫, 吴大为. 基于改进遗传算法的公交线网整体优化方法[J]. 重庆交通学院学报, 2004, 23 (6): 97-101. https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT200406024.htm

    TANG Ke-fu, WU Da-wei. Study on the optimization methods of bus network based on adapted genetic algorithm[J]. Journal of Chongqing Jiaotong University, 2004, 23 (6): 97-101. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-CQJT200406024.htm
    [15] 常玉林, 胡启洲. 城市公交线网优化的线性模型[J]. 中国公路学报, 2005, 18 (1): 95-98. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200501021.htm

    CHANG Yu-lin, HU Qi-zhou. Optimal line model on urban public traffic line network[J]. China Journal of Highway and Transport, 2005, 18 (1): 95-98. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200501021.htm
    [16] 于滨, 杨永志, 杨忠振, 等. 基于直达客流密度最大的公交线网优化[J]. 哈尔滨工业大学学报, 2009, 41 (2): 205-207. https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX200902050.htm

    YU Bin, YANG Yong-zhi, YANG Zhong-zhen, et al. Transit network optimization based on direct passenger flow density maximization[J]. Journal of Harbin Institute of Technology, 2009, 41 (2): 205-207. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-HEBX200902050.htm
    [17] 于滨, 刘鸿婷, 闫博, 等. 公交线路网优化的双层模型及其解法[J]. 吉林大学学报: 工学版, 2010, 40 (2): 402-405. https://www.cnki.com.cn/Article/CJFDTOTAL-JLGY201002020.htm

    YU Bin, LIU Hong-ting, YAN Bo, et al. Bi-level model for bus route network optimization and its solution[J]. Journal of Jilin University: Engineering and Technology Edition, 2010, 40 (2): 402-405. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-JLGY201002020.htm
    [18] 刘伟铭, 姜山, 付凌峰. 多车型高速公路离散平衡网络设计的双层规划模型[J]. 中国公路学报, 2008, 21 (1): 94-99. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200801018.htm

    LIU Wei-ming, JIANG Shan, FU Ling-feng. Bi-level program model for multi-type freeway discrete equilibrium network design[J]. China Journal of Highway and Transport, 2008, 21 (1): 94-99. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL200801018.htm
    [19] 马嘉琪, 白雁, 韩宝明. 城市轨道交通线网基本单元与复杂网络性能分析[J]. 交通运输工程学报, 2010, 10 (4): 65-70, 102. http://transport.chd.edu.cn/article/id/201004011

    MA Jia-qi, BAI Yan, HAN Bao-ming. Characteristic analysis of basic unit and complex network for urban rail transit[J]. Journal of Traffic and Transportation Engineering, 2010, 10 (4): 65-70, 102. (in Chinese). http://transport.chd.edu.cn/article/id/201004011
    [20] 蒋阳升, 罗孝羚. 考虑首末站约束和站间客流强度的公交线网优化[J]. 长安大学学报: 自然科学版, 2017, 37 (1): 106-111. https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL201701014.htm

    JIANG Yang-sheng, LUO Xiao-ling. Optimization of public transit network considering initial and terminal stations location requirements and passenger flow intensity[J]. Journal of Chang'an University: Natural Science Edition, 2017, 37 (1): 106-111. (in Chinese). https://www.cnki.com.cn/Article/CJFDTOTAL-XAGL201701014.htm
    [21] YAO Jia, SHI Feng, AN Shi, et al. Evaluation of exclusive bus lanes in a bi-modal degradable road network[J]. Transportation Research Part C: Emerging Technologies, 2015, 60: 36-51.
    [22] PANASYUK M V, PUDOVIK E M, SABIROVA M E. Optimization of regional passenger bus traffic network[J]. Procedia Economics and Finance, 2013, 5: 589-596.
    [23] AMIRIPOUR S M M, CEDER A A, MOHAYMANY A S. Designing large-scale bus network with seasonal variations of demand[J]. Transportation Research Part C: Emerging Technologies, 2014, 48: 322-338.
    [24] YU Bin, KONG Lu, SUN Yao, et al. A bi-level programming for bus lane network design[J]. Transportation Research Part C: Emerging Technologies, 2015, 55: 310-327.
    [25] LI Qiang, WANG Yun-jing, ZHOU Yang. Readjustment effect analysis of public bus network in Beijing[J]. Journal of Transportation Systems Engineering and Information Technology, 2008, 8 (2): 27-33.
    [26] ROCA-RIU M, ESTRADA M, TRAPOTE C. The design of interurban bus network in city centers[J]. Transportation Research Part A: Policy and Practice, 2012, 46 (8): 1153-1165.
  • 加载中
图(6) / 表(4)
计量
  • 文章访问数:  702
  • HTML全文浏览量:  99
  • PDF下载量:  736
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-08-01
  • 刊出日期:  2017-02-25

目录

/

返回文章
返回