版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
交通工程学(第4章交通流理论)
交通工程学4.1概述(了解)4.2交通流的统计分布特性(熟练掌握)4.3排队论模型(熟练掌握)4.4跟驰模型(熟练掌握)4.5流体模型(熟练掌握)第4章交通流理论4.1概述(了解)第4章交通流理论交通流理论是交通工程学的基本理论,是借助于物理、数学的定律与方法来阐明交通流基本特性的一种理论。4.1概述交通流理论是交通工程学的基本理论,历史沿革:诞生于20世纪30年代的概率论的方法.50年代以后成为交通工程的专题研究研究内容宏观稳态的交通流基本参数模型交通流统计分布特性交通流理论的模拟与仿真交通流模型的理论与方法:排队论、跟驰理论、流体力学理论、细胞源动机理论历史沿革:诞生于20世纪30年代的概率论的方法.50年代以后交通模型:描述交通流状态变量随时间与空间而变化的分布规律及其与交通控制变量之间关系的方程。参数模型:交通流参数之间的关系宏观模型:描述车队的运动规律微观模型:描述单个车辆的运动规律静态模型:不随时间改变的稳恒交通流随空间分布的规律动态模型:时间改变的稳恒交通流随空间分布的规律交通模型:描述交通流状态变量随时间与空间而变化的分布规律及其4.2.1交通流统计分布的含义4.2.2离散型分布4.2.3连续性分布4.2交通流的统计分布特性4.2.1交通流统计分布的含义4.2交通流的统计分布特性车辆的到达在某种程度上具有随机性,描述这种随机性的统计规律的方法称为交通流的统计分布。离散型分布:考察在一段固定长度的时间内到达某场所的交通数量或一定距离内分布的交通数量的波动性。
信号周期内到达的车辆数。连续型分布:描述事件之间时间间隔的连续型分布为工具,研究事件发生的间隔时间或距离的统计分布特性。车头时距分布、速度分布和可穿越空档分布。
4.2.1交通流统计分布的含义车辆的到达在某种程度上具有随机性,描述这种随机4.2.2离散型分布4.2.2.1泊松分布4.2.2.2二项分布4.2.2离散型分布4.2.2.1泊松分布4.2.2.24.2.2.1泊松分布(1)基本公式,k=0,1,2,…
Pk—在计数间隔t内到达k辆车或k个人的概率;λ—单位时间间隔的平均到达率(辆/s或人/s);t—每个计数间隔持续的时间(s)。若令m=λt为计数间隔t内平均到达的车辆(人)数,则,当m为已知时,可求出在计数间隔t内恰好有k辆车(人)到达的概率。4.2.2.1泊松分布(1)基本公式第4章交通流理论课件第4章交通流理论课件4.2.2.1泊松分布(续)(2)递推公式:,(3)适用条件:车流密度不大,车辆间相互影响较弱,其他外界干扰因素基本上不存在,即车流是随机的。(4)泊松分布的均值M和方差D都等于λt,而观测数据的均值m和方差S2均为无偏估计,因此,当观测数据表明S2/m显著地不等于1.0时,就是泊松分布不合适的表示。m—在某一给定时间间隔周期内到达车辆的平均数;S2—各车辆到达数与均值之差的平方和的平均数。4.2.2.1泊松分布(续)(2)递推公式:(3)适用条件4.2.2.1泊松分布(续)例4-1某路段每小时有120辆车通过,假设车辆到达服从泊松分布,问在指定的某一分钟内有3辆车通过的概率是多大,而一分钟内不超过3辆车的概率又是多大。(5)应用举例例4-2某信号灯交叉口的周期C=97s,有效绿灯时间g=44s,在有效绿灯时间内排队的车流以S=900(辆/h)的交通量通过交叉口,在有效绿灯时间外到达的车辆要停车排队。设信号灯交叉口上游车辆的到达率q=369(辆/h),服从泊松分布公式中,求到达车辆不致二次排队的周期数占周期总数的最大百分率。4.2.2.1泊松分布(续)例4-1某路段每小时有1204.2.2.1泊松分布(续)例4-2解:一个周期内能通过的最大车辆数A=gS=900×44/3600=11辆,当某周期到达的车辆数N≻11辆时,则最后到达的(N-11)辆车就不能在本周期内通过而发生二次排队。在泊松分布中,一个周期内平均到达的车辆数m=λt=369×97/3600=9.9辆。则可能到达车辆数大于11辆的周期出现的概率为即到达车辆不致两次排队的周期数最多占71%。4.2.2.1泊松分布(续)例4-2即到达车辆不致两次排队4.2.2.2二项分布(1)基本公式:,k=0,1,2,…Pk—在计数间隔t内到达k辆车或k个人的概率;λ—单位时间间隔的平均到达率(辆/s或人/s);t—每个计数间隔持续的时间(s)或距离(m);n—观测次数,正整数。通常记,则二项分布为:4.2.2.2二项分布(1)基本公式:通常记4.1.2.2二项分布(续)(2)递推公式:(3)适用条件:车辆比较拥挤、自由行驶机会不多的车流。(4)分布的均值M和方差D分别为M=np,D=np(1-p),显然有M>D。用观测数据计算出来的样本均值m和方差S2代替M和D,因此,S2/m应当小于1。当S2/m显著大于1.0时,就是二项分布不适的表示。4.1.2.2二项分布(续)(2)递推公式:(3)适用条件4.2.2.3负二项分布基本公式:4.2.2.3负二项分布基本公式:适用条件:车流受到干扰。车辆到达起伏幅度比较大的车流,而计数周期比较短的高方差分布分布的均值M和方差D分别为M=kp/p,D=kp/p2,显然有M<D。用观测数据计算出来的样本均值m和方差S2代替M和D,所以负二项分布的S2/m应当大于1,当S2/m显著小于1.0时,就是负二项分布不适的表示。适用条件:车流受到干扰。车辆到达起伏幅度比较大的车流,而计数4.2.3连续型分布4.2.3.1负指数分布4.2.3.2移位负指数分布4.2.3连续型分布4.2.3.1负指数分布4.2.3.1负指数分布(1)
基本公式:P(h>t)——到达的车头时距h大于t秒的概率;λ——车流的平均到达率(辆/s)。推导:由可知,在计数间隔t内没有车辆(k=0)到达的概率,这表明,在具体的时间间隔t内,无车辆到达,则上次车到达和下次车到达之间,车头时距至少有t,即。4.2.3.1负指数分布(1)基本公式:推导:由4.2.3.1负指数分布(续)(2)负指数分布的均值M和方差D分别为M=1/λ,D=1/λ2,用样本均值m代替M、样本的方差S2代替D,既可算出负指数分布的参数λ
。(3)适用条件:用于描述有充分超车机会的单列车流和密度不大的多列车流的车头时距分布,它常与计数的泊松分布相对应。(4)负指数分布的概率密度函数是单降的,车头时距越短,其出现的概率越大,但车头时距至少有一个车长,所以车头时距必有一个大于零的最小值τ。4.2.3.1负指数分布(续)(2)负指数分布的均值M和方负指数分布负指数分布4.2.3.1负指数分布(续)(5)应用举例例4-3某交通流属泊松分布,已知交通量为1200辆/h,求:a)车头时距t>
5s的概率;b)在1小时内,车头时距t>5s所出现的次数;
4.2.3.1负指数分布(续)(5)应用举例在次要车流通行能力研究中的应用在次要车流通行能力研究中的应用4.2.3.2移位负指数分布(1)基本公式为克服负指数分布的车头时距趋近于零其频率出现愈大这一缺点,可将负指数分布曲线从原点O沿t向右移一个最小间隔长度τ,得到移位负指数分布曲线:τ—大于零的一个最小车头时距,一般在1.0~1.5s之间。(2)移位负指数分布的均值M和方差D分别为M=1/λ+τ
,D=1/λ2,用样本均值m代替M、样本的方差S2代替D,则可算出移位负指数分布的参数λ和τ。4.2.3.2移位负指数分布(1)基本公式τ—大于零的一个4.2.3.2移位负指数分布(续)(3)适用条件用于描述不能超车的单列车流的车头时距分布和车流量低的车流的车头时距分布。(4)移位负指数分布的局限移位负指数分布的概率密度函数曲线是随t-τ单调递降的,车头时距愈接近τ,其出现的可能性愈大。这在一般情况下是不符合驾驶员的心理习惯和行车特点的。从统计角度看,车头时距分布的概率密度曲线一般总是先升后降的。4.2.3.2移位负指数分布(续)(3)适用条件其他常用分布形式爱尔兰分布:T:观测时间间隔的平均值T:车头时距(s)H:车头时距的观测值当k=1时,为负指数分布当k>1时,为爱尔兰分布K:确定分布曲线形状的参数K值4舍5入,取整数其他常用分布形式爱尔兰分布:K:确定分布曲线形状的参数K值4对数正态分布:对数正态分布:分布检验在实际观测中,不可能对观测值的全部取值的概率进行反复观测,往往是以局部观测数列的分析和观测值的算术平均值或方差为依据,推断其符合某种分布规律为了运用局部观测资料,即用样本推算总体的分布,需要先对总体的分布概率进行假设,然后用局部观测的数据来验证其符合程度拟合度检验:实际样本与理论样本之间总存在差异,且随机取样也存在样本之波动,其差异与变化程度究竟如何,即拟合度如何,只能通过拟合度检验来鉴别。常用的检验为x2检验(Chiquaretest)分布检验在实际观测中,不可能对观测值的全部取值的概率进行反复检验的原理:首先假设观测数列具有某种分布特性,于是建立实际频率与理论频率之间的差异,此差异用计算值X2表示。在确定的显著水平的条件下确定临界值x2。当计算值x2小于临界值x2时,假设分布被接受,否则,重新假设分布,重新进行计算检验的原理:首先假设观测数列具有某种分布特性,于是建立实际频检验计算过程:1、建立原假设H0计算p1、p2、p3。。。。。Pn计算F1、F2、F3………….Fn、2、选取统计量:检验计算过程:1、建立原假设H03、确定临界值:由显著水平与自由度DF确定DF=c-a-1由表4-2查出临界值x24、求统计检验结论:X2计算≤X2临界,假设成立,分布被接受否则,重新假设其分布,重新进行检验3、确定临界值:常用统计分布中的a值与DF值分布aDF泊松分布1C-2二项分布2C-3负二项分布2C-3正态分布2C-3均匀分布0C-1常用统计分布中的a值与DF值分布aDF泊松分布1C-2二项分X2检验中需要注意的事项:样本量较大,N≥50分组数应该连续,以7-9组为宜,一般不小于5组各组的理论频数不得小于5,如Ej=5,则应该将相邻的组项合并,直至Ej≥5为止。这时应以合并后的组数作为计算自由度的c值X2检验中需要注意的事项:样本量较大,N≥50例:下表为某观测现场的车流量数据,时间间隔为1min,试检验其分布规律是否服从泊松分布?显著水平为5%组序号每分钟到达的车辆数xi频数fi累计车辆数100021993261243927541144659457653087179≥800=50=174例:下表为某观测现场的车流量数据,时间间隔为1min,试检验解:1、根据泊松分布计算Pi2、计算理论频数Ej解:1、根据泊松分布计算Pi2、计算理论频数Ej组号fiPinFj100.0310501.550.57290.1079505.43260.1877509.391.234390.21775010.890.3354110.1894509.470.246590.1318506.590.87750.0765503.830.07810.0380501.90900.0196500.981503.31960.13890.13457.0716组号fiPinFj100.0310501.550.572903、确定DF=c-a-1=6-1-1=4查表的:4、判断分布是否成立:所以原假设成立,分布服从泊松分布3、确定4.3.1基本概念4.3.2基本原理4.3.3排队系统的表示4.3排队论模型4.3.1基本概念4.3排队论模型第4章交通流理论课件(1)排队论:是研究“服务”系统因“需求”拥挤而产生等待行列(即排队)的现象,以及合理协调“需求”与“服务”关系的一种数学理论。(2)排队:单指等待服务的车辆,不包括正在被服务的车辆。(3)排队系统:既包括了等待服务的,又包括了正在被服务的车辆。(4)排队论的应用:电话自动交换机;车辆延误、通行能力、信号灯配时以及停车场、加油站等交通设施的设计与管理;收费亭的延误估计。4.3.1基本概念(1)排队论:是研究“服务”系统因“需求”拥挤而产生等待行列(1)排队系统的3个组成部分输入过程:各种类型的“顾客(车辆或行人)”按怎样的规律到达。(到达时距符合什么样的分布)如定长输入D;泊松输入M;爱尔郎输入EK
排队规则:指到达的顾客按怎样的次序接受服务。如损失制;等待制;混合制。等待制:先到先服务:FIFO后到先服务:LIFO随机服务:SIRO4.3.2基本原理(1)排队系统的3个组成部分4.3.2基本原理服务方式:指同一时刻多少服务台可接纳顾客,每一顾客服务了多少时间。服务台:个数、排列方式(平行排列、成串排列)服务时间服从何种分布形式:如定长分布D;负指数分布M;爱尔朗分布Ek。服务方式:指同一时刻多少服务台可接纳顾客,每一顾客服务了多少(2)排队系统的主要数量指标队长和排队长:若排队系统中的顾客数为n,排队顾客数为q,正在被服务的顾客数位s,则n=q+s。队长是排队系统提供的服务水平的一种衡量。逗留时间和等待时间:逗留时间是指一个顾客逗留在排队系统中的总时间。等待时间是指从顾客到达时起到他开始接受服务时止这段时间。忙期和闲期:忙期是指服务台连续繁忙的时期,相对应的是闲期,这关系到服务台的工作强度。4.3.2基本原理(续)(2)排队系统的主要数量指标4.3.2基本原理(续)4.3.3排队系统的表示类别输入分布服务方式服务台数量符号含义M—泊松或负指数分布M—负指数分布1D—定长D—定长NEk—爱尔朗分布Ek—爱尔朗分布M/M/N——泊松输入、负指数分布服务、N个服务台M/D/1——泊松输入、定长服务、单个服务台a/b/c(L/Disc):a:车辆到达的分布b:服务时间分布c:服务台个数L:允许排队长度Disc:排队规则4.3.3排队系统的表示类别输入分布服务方式服务台数量符号M/M/1系统参数判别指标:顾客的平均到达率:系统的服务率:系统不稳定,队伍越来越长(或一直不消散)排队消散,系统稳定M/M/1系统参数判别指标:顾客的平均到达率:系统的服务率:M/M/1系统参数计算M/M/1系统参数计算M/M/1系统参数计算M/M/1系统参数计算M/M/N系统的计算公式判别指标:多路排队多通道系统:相当于N个M/M/1系统M/M/N系统的计算公式判别指标:多路排队多通道系统:相当于M/M/N系统的计算公式单路排队多通道服务系统:系统中没有车辆的概率:系统中有n辆车的概率:系统的平均车辆数:系统中排队的平均长度:系统中的平均消耗时间:系统中排队的平均等待时间:M/M/N系统的计算公式单路排队多通道服务系统:系统中没有车4.4.1车辆跟驰特性分析4.4.2线形跟驰模型4.4跟驰模型4.4.1车辆跟驰特性分析4.4跟驰模型4.3跟驰理论用动力学方法描述车队后车跟随前车行驶状态,通过描述车辆之间的行驶关联性描述交通流能够准确交通服务水平广泛应用于交通模拟与仿真4.3跟驰理论用动力学方法描述车队后车跟随前车行驶状态,通(1)跟驰理论的定义:运用动力学的方法,研究在无法超车的单一车道上车辆列队行驶时,后车跟随前车的行驶状态的一种理论。(2)车辆跟驰特性分析(非自由行驶状态的车队)
制约性:后车紧随前车前进。延迟性(滞后性):后车运行状态的改变在前车之后。传递性:前车的运行状态制约着后车的运行状态。4.4.1车辆跟驰特性分析车速条件距离条件(1)跟驰理论的定义:运用动力学的方法,研究在无法超车的单一4.4.2线形跟驰模型在t时刻,由于前车n的减速造成后车n+1的减速,由于车辆跟驰的延迟性,后车的减速滞后了T(驾驶员的反应时间)。在t时刻,前车和后车的位置分别为xn(t)和xn+1(t),两车此时的距离为S(t)=xn(t)-xn+1(t)。后车在反应时间T内行驶的距离。表示第i辆车在时刻t的速度。4.4.2线形跟驰模型在t时刻,由于前车n4.3.2线性跟驰模型微分方程:4.3.2线性跟驰模型微分方程:4.4.2线形跟驰模型(续)其中为后车在时刻(t+T)的加速度,称为后车的反应;为敏感度;为时刻t的刺激。所以:反应=敏感度×刺激。假定d2=d3,要使在时刻t两车的间距能保证在突然刹车事件中不发生碰撞,则有:即对t微分得:或4.4.2线形跟驰模型(续)其中非线性跟车理论称为灵敏度m和L为常数当m=0,L=0,即为线形模型非线性跟车理论称为灵敏度m和L为常数当m=0,L=0,即为线4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.3.3跟驰模型的讨论4.5.1理论概述4.5.2车流连续性方程4.5.3波动理论4.5.4交通波理论的应用举例4.5流体模型4.5.1理论概述4.5流体模型1955年,英国学者莱脱希尔和惠特汉提出。车流波动理论的定义:通过分析车流波的传播速度,以寻求车流流量和密度、速度之间的关系,并描述车流的拥挤——消散过程。适用条件:流体力学模拟理论假定在车流中各个单个车辆的行驶状态与它前面的车辆完全一样,这与实际不符,因此该模型运用于车辆拥挤路段较为合适。4.5.1理论概述1955年,英国学者莱脱希尔和惠特汉提出。4.5.1理论概交通流与流体流特性对比物理特性流体动力学系统交通流系统连续体形态单向不可压缩流体单车道不可压缩车流离散元素分子车辆运动方向一向性单方向变量质量m密度K速度V车速V压力P流量Q动量mVkm状态方程Q=km连续性方程运动方程交通流与流体流特性对比物理特性流体动力学系统交通流系统连续体假设车辆顺次通过断面Ⅰ和Ⅱ的时间间隔为dt,两断面的间距为dx。同时,车流在断面Ⅰ的流入量为q,密度为K。车流在断面Ⅱ的流出量为(q+△q),密度为(K-△K)。4.5.2车流连续性方程根据物质守恒定律,流入量-流出量=△x内车辆数的变化:
或取极限得:当车流量随距离而降低时,车流密度随时间而增大。
又因为q=Kv,交通流的运动方程为假设车辆顺次通过断面Ⅰ和Ⅱ的时间间隔为dt,两断面的间距为d车流的波动:车流中两种不同密度的分界面经过一辆辆车向后部传播的现象。波速:车流波动沿道路移动的速度。前进波:沿道路前进的波,波速为正。后退波;沿道路后退的波,波速为负。集结波:波阵面过后,车流密度变大。疏散波:波阵面过后,车流密度变小。集散波:包括集结波和疏散波。4.5.3车流波动理论(1)基本概念车流的波动:车流中两种不同密度的分界面经过一辆辆车向后部传播车队从速度Vl、密度K1(对应于车间距离l1)转变到速度V2,密度K2(对应于车间距离l2)。O为第一辆车的变速点,A为第二辆车的变速点、虚线OA的斜率就是集散波的波速。4.5.3车流波动理论(续)一个车队中前三辆车运行的时间-空间轨迹(2)交通波的基本方程车队从速度Vl、密度K1(对应于车间设变速点A的时刻为t,位置为x,则在时刻0到时刻t之间,两车车间距的变化为l2-l1,第一辆车行驶的距离为tV2,第二辆车行驶的距离为tV1,则l2-l1=tV2-tV1,t=(l2-l1)/(V2-V1)又因x=-l1+tV1,则可得波速公式:4.5.3车流波动理论(续)如果车流前后两行驶状态的流量和密度非常接近,则上式叫可演化为,这个公式是微弱波的波速公式,即车流中传播小紊流的速度公式。(2)交通波的基本方程设变速点A的时刻为t,位置为x,则在时刻0到4.5.3车流波动理论(续)设有一个交通波以速度w沿车道稳定地向右传播,波阵面s前车流密度为k1,速度为u1,波传过后车流密度变为k2,速度为u2。以波阵面s为界面,将看到的原车流以w—u1的速度向左流过波阵面,而以w—u2的速度从波阵面流出。假设为单车道,根据质量守恒定律,在波稳定传播的条件下,时间t内从波阵面右侧流入的车辆数应等于从左侧流出的车辆数,得到:(3)交通波的基本方程(简单证明)4.5.3车流波动理论(续)设有一个交通波以速度w沿车道稳4.5.3车流波动理论(续)集散波总是从前车向后车传播的,把单位时间内集散波所掠过的车辆数称为波流量。通常意义下的流量总是相对于道路的一个固定断面而言,而波流量则是相对于移动的波界面来计算的。可以证明,波流量的公式为Qw——车流波W的波流量;V2、V1——前后两种车流状态的车速;K2、K1——前后两种密度。(4)波流量4.5.3车流波动理论(续)集散波总是从前4.5.3车流波动理论(续)①k2>k1,且q2>q1——集结波、前进波相当于以较大的间距行驶的车队,后车催促前车依次不断加速逐步缩小间距的情况。②k2>k1,且q2<q1——集结波、后退波相当于车队中的头车减速或刹车,跟随车辆依次采取同样行为的情况,如车队驶进信号灯控制的交叉口而红灯启亮的情况。③k2<k1,且q2>q1——发散波、后退波相当于停在停车线后的车队,绿灯启亮后逐渐启动的情形。④k2<k1,且q2<q1——发散波、前进波相当于以较小间距行驶的车队,从队尾起,各车辆依次减速,逐步拉大车距的情况(5)波速公式的分析4.5.3车流波动理论(续)①k2>k1,且q2>q1—4.5.4交通波理论的应用举例例4-4车流在一条6车道的公路上畅通行驶,其速度为V=80km/h。路上有座4车道的桥,每车道的通行能力为1940辆/h,高峰时流量为4200辆/h(单向)。在过渡段的车速降至22km/h,这样持续了1.69h,然后车流量减到1956辆/h(单向)。试估计桥前车辆的排队长度和阻塞时间。(1)分析道路上瓶颈地段的车流状况(2)低速车插入高速车流产生的影响(3)车队在信号等交叉口处的排队长度4.5.4交通波理论的应用举例例4-4车流在一条6车
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 过渠涵施工方案(3篇)
- 金属厂房防腐施工方案(3篇)
- 韩国炸鸡店营销方案(3篇)
- 鱼塘边桥墩施工方案(3篇)
- 26年腭癌靶向禁忌症速记
- 眩晕与老年健康
- 胎儿窘迫的胎儿神经保护措施
- 柔性版印刷员保密意识知识考核试卷含答案
- 重轨加工工操作知识模拟考核试卷含答案
- 手术器械装配调试工安全生产能力水平考核试卷含答案
- 2025高考志愿第五轮学科评估(部分)+第四轮学科评估结果Excel表格
- 2025榆林能源集团有限公司招聘工作人员(473人)笔试参考题库附带答案详解析
- 云南省土地征收农用地转用审批管理细则 (2023年修订)
- 2024年中央司法警官学院招聘笔试真题
- 小红书运营:小红书账号运营培训课件
- 《运筹学(第3版)》 课件 第5章 整数规划;第6章 动态规划
- 全回转钻机在拔桩、清障中的应用
- 全国职业院校技能大赛高职组(商务数据分析赛项)备赛试题库(含答案)
- (正式版)QBT 2174-2024 不锈钢厨具
- 生态环境保护论文生态环境建设与水环境保护
- 建筑消防设施年度检测报告
评论
0/150
提交评论