首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 203 毫秒
1.
福州市PM2.5污染过程中大气边界层和区域传输研究   总被引:5,自引:0,他引:5  
以福建省会福州市2013年1月空气质量变化为对象,分析大气边界层变化和周边区域污染物传输对福州市大气颗粒物PM2.5的影响.利用福州市2013年1月逐日地面和探空观测资料以及NCEP提供的2013年1月FNL分析资料,通过大气边界层要素与PM2.5浓度之间的相关性,对PM2.5污染过程的大气边界层特征进行分析;同时采用HYSPLIT后向轨迹模拟及区域风场相关矢分析对影响福州雾霾的污染物区域传输路径进行探讨.结果表明:地面气温与PM2.5浓度呈正相关,地面风速与PM2.5浓度呈负相关,近地面边界层条件有利于霾颗粒物的形成和累积.但不同于我国东部主要污染源区霾污染过程中存在大气边界层逆温,福州PM2.5污染过程中并未出现大气边界层逆温结构,这一边界层结构的垂直混合可有利于区域传输的污染物从上层大气到达近地面从而加重福州霾污染,福州是华东地区一个PM2.5污染物的主要接受区,PM2.5污染物主要以外源输送为主.2013年1月份福州市清洁日近地面风向为海洋吹向大陆的东南风,霾污染日则为大陆吹向海洋的偏北风,PM2.5污染物主要从长三角地区、苏北以及安徽河南一带通过东北和西北方向的传输路径影响福州的空气质量.  相似文献   

2.
珠江三角洲秋季典型气溶胶污染的过程分析   总被引:2,自引:0,他引:2  
为了解大气中各物理和化学过程对气溶胶浓度的贡献情况,利用Models-3/CMAQ模式系统对珠江三角洲(以下简称珠三角)秋季典型气溶胶污染进行研究.模拟时间是2012年10月,期间珠三角主要受高压系统的控制,在17日冷锋过境前后高压天气形势发生转变,风向从东北风转为偏东风.结果表明,珠三角秋季PM2.5浓度呈现西高东低的水平分布特征,随着高度的上升浓度高值中心也向西南方向偏移;受大气边界层高度的影响,陆地上PM2.5输送高度呈现白天高夜晚低的变化特征;过程分析结果表明源排放,水平输送和垂直输送是影响近地面PM2.5浓度变化的主要过程;本地污染物排放是城市中心(广州站)PM2.5浓度升高的主要原因,而在下风向位置(江门站)外来污染物的水平输送过程是PM2.5的最主要来源.  相似文献   

3.
石家庄大气污染物输送通道及污染源区研究   总被引:4,自引:1,他引:3  
为探索石家庄的区域输送规律,确定主要污染源区,利用HYSPLIT(Hybrid Single Particle Lagrangian Integrater Trajectory)后向轨迹模式和NCEP的GDAS全球气象要素数据,对2013—2016年从不同高度上抵达石家庄地区的逐日72 h气流后向轨迹进行聚类分析,并结合石家庄逐小时颗粒物污染物浓度数据,分析石家庄PM2.5的潜在源贡献因子(WPSCF)和浓度权重轨迹(WCWT).结果表明,(1)石家庄PM2.5浓度具有明显单峰谷日变化,秋冬季与春夏季峰谷值出现时间不同;(2)近地层大气污染输送路径以近距离,移速慢的轨迹为主,轨迹较短的路径所占比例在40%以上.除夏季外,近距离输送路径均存在螺旋转向,在后向48~36 h内轨迹端点到达河北省内,转为东向和南向输送.(3)大气污染输送通道的垂直分布特征表明,输送轨迹中低于500 m高度的轨迹点占28.7%,高于1000 m低于3000 m高度的轨迹点占36.1%,高于3000 m高度的轨迹点占25.3%.低层多以近距离输送为主,高度越高,近距离输送轨迹的频率越低.500 m高度输送通道仍以近距离输送为主,并存在螺旋转向,1500 m高度以上多远距离输送.(4)石家庄PM2.5的主要污染源区范围较小.途径河北中南部、河南北部、山东西部和山西中北部地区的轨迹对石家庄PM2.5的污染贡献最大.  相似文献   

4.
北京地区秋季雾霾天PM2.5污染与气溶胶光学特征分析   总被引:15,自引:9,他引:6  
利用北京城区和郊区2011年9月1日~12月7日PM2.5质量浓度、气溶胶散射系数(σsca)和黑碳浓度观测资料,研究了雾霾天气条件下北京地区PM2.5污染与气溶胶光学参数的变化特征,并讨论了气象条件的作用.结果表明,北京地区PM2.5污染和气溶胶光学特性受雾霾天气的影响非常明显.PM2.5浓度、σsca和气溶胶吸收系数(σabs)在雾霾期均明显高于非雾霾期,雾霾期日均PM2.5浓度在城区和郊区分别达到97.6μg·m-3和64.4μg·m-3,为非雾霾期日均浓度的3.3和4.8倍.城区高PM2.5浓度造成雾霾类天气出现频率明显高于郊区.轻雾天城区PM2.5浓度、σsca和σabs明显高于郊区,区域输送的影响相对较弱,轻雾和霾天城郊差异较小,区域性特征明显,而雾天σsca城郊非常接近且在各雾霾类天气中相对最高,气溶胶散射能力最强,区域性特征较为明显.气象条件的不同造成各雾霾过程PM2.5浓度、σsca和σabs的空间分布、PM2.5污染及气溶胶消光强度上呈现不同的特点.边界层以上偏南风将南部地区污染物向北京输送,在整层下沉气流作用下使得边界层内污染物浓度增加,加之边界层高度持续稳定在600 m左右,边界层内风速很低,污染物水平、垂直扩散均很弱,造成局地污染物的累积,形成了PM2.5污染和气溶胶消光强度最强的一次雾霾过程.  相似文献   

5.
长江三角洲地区大气污染物水平输送场特征分析   总被引:24,自引:3,他引:21  
以长江三角洲地区为研究对象,探讨城市群环境污染规律.利用区域中尺度大气数值预报模式MM5模拟2004年1,4,7和10月长江三角洲地区的气象场.结果表明:冬季,长江三角洲地区近地面及500 m高度层主要为西北风控制,输送气流主要来自西部内陆地区;春季,盛行东南风和偏东风,存在明显的东南向西北方向的输送气流;夏季,则以偏南输送气流为主,杭州湾地区海面向内陆方向以及太湖湖面风速较大,输送扩散能力较强;秋季则转为东北风,近地面杭州湾以北盛行北风,以南主要受海面东北风的影响.结合HYSPLIT-4三维轨迹模式计算分析该地区典型污染过程时污染物的输送气流轨迹,证实了污染过程伴随500 m高度处东北主频气流的“外源”输入现象.   相似文献   

6.
准确的空气质量数值预报模式依赖于精确的气象条件模拟,尤其依赖于大气边界层的准确模拟.为理解边界层过程如何影响空气污染物的传输与混合,利用WRF-Chem模式不同边界层方案(YSU和MYJ)进行敏感性试验,针对山西冬季典型静稳天气,对地面温度场、地面风场、PM_(2.5)浓度及边界层内部的动力和热力层结进行模拟分析,并与观测资料进行对比,分析不同PBL方案对于气象要素和PM_(2.5)浓度分布的模拟能力,探讨边界层内部热力层结和湍流输送差异对PM_(2.5)浓度模拟的影响.结果表明:2种边界层方案均能较好模拟出冬季静稳天气背景下地面温度、风速及PM_(2.5)浓度的空间分布和日变化特征,气温模拟的较大误差主要出现在夜间,而地面风速和PM_(2.5)浓度的模拟结果在午后误差较大;相对于YSU方案,局地MYJ方案模拟的温度、风场和PM_(2.5)浓度的误差更小,模拟结果更接近于实况观测.地面PM_(2.5)浓度的模拟误差可能与近地面逆温层、混合层及地面风速等的模拟误差有关;不同边界层参数化方案导致的边界层内热力层结和湍流输送的模拟差异,可能是影响近地面PM_(2.5)浓度模拟差异的主要原因;夜间MYJ方案逆温层厚度较厚,地面PM_(2.5)模拟浓度较低;午后MYJ方案混合层高度较低,加之地面风速较弱,导致地面PM_(2.5)模拟浓度较高.  相似文献   

7.
冷空气过程对江苏持续性霾的影响研究   总被引:1,自引:0,他引:1  
利用地面气象观测资料、PM_(2.5)浓度监测资料和数值模式产品对2016年12月14—24日江苏遭遇的一次长时间霾天气过程进行分析.研究结果表明:过程期间有两次冷空气南下影响江苏省,两次冷空气均带来大风和降水,有效地清除了前期污染物,但随后全省PM_(2.5)浓度开始升高.第一次冷空气强度强于第二次,造成的污染也较重.WRF-Chem模式对本次过程的气象场和PM_(2.5)浓度模拟均较好,模拟观测相关系数分别达到0.52~0.99和0.40.模式能够较好地模拟出污染物的输送过程和时空分布.与第二次冷空气过程相比,第一次冷空气过程存在明显的污染物自北向南输送过程,100~500 m高空持续偏北气流(第2次过程为西北-偏西气流),期间全省平均边界层高度(PBLH)只有260 m(低于第2次过程的500 m),不利于污染物垂直扩散,造成地面浓度较高.利用HYSPLIT-4模式追踪了两次过程中淮安、泰州、无锡三站上空100 m处大气48 h后向轨迹,发现第一次过程中污染物来自山东中西部,第二次来自西部内陆地区.  相似文献   

8.
吴蒙  罗云  吴兑  范绍佳 《中国环境科学》2016,36(11):3263-3272
通过2004年10月在珠江三角洲(以下简称“珠三角”)开展大气边界层观测试验得到的垂直风温资料和逐时PM2.5浓度资料,利用局地环流指数(RF)等方法研究了珠三角海陆风特征及其对空气质量的影响.结果表明:局地环流指数是表征局地大气输送能力的有效指标;冷暖气团对峙导致珠三角污染日背景风场较弱,沿海海陆风活动活跃,空气质量指数与RF系数相关性颇高,珠三角沿海100~400m处RF系数值主要分布在0.5~0.8之间;在海陆风影响下低层风场的有效输送能力较弱,不利于污染物的输送扩散.而随着冷空气全面控制珠三角,垂直风场RF系数值高达0.9以上,海陆风难以发展,风场输送能力强,能够持续的将污染物输送出去.观测发现沿海观测点试验期间海陆风发生频率约为47.8%,其中72.7%的海陆风日出现了污染天气,海陆风日地面风向呈现出明显的随时间顺时针偏转特征,海风约从16:00时开始出现,并在20:00达到最大影响高度约为600~800m.夜间海风将污染物输送回内陆观测点,导致内陆PM2.5浓度在19:00~21:00时出现浓度峰值,呈现出明显的双峰结构.  相似文献   

9.
2014年春季山东省PM2.5跨界输送研究   总被引:2,自引:2,他引:2  
利用CAMx模式模拟了2014年春季(3~5月)山东省ρ(PM2.5)浓度空间分布、输送、沉降特征并建立了各市PM2.5跨界输送矩阵.结果表明,山东省ρ(PM2.5)分布呈现明显的空间差异特征,鲁西地区浓度较高,半岛地区浓度较低;垂直方向上ρ(PM2.5)主要集中在2 000 m以下.全省ρ(PM2.5)水平输送通量平均达110μg·(m2·s)-1,全省总沉降量为23.05×104t.全省近地面ρ(PM2.5)外来输送贡献率为21.08%±3.83%,各市外来输送贡献率平均为40.45%±5.96%;随着高度的升高,京津冀、背景及边界等条件对山东省ρ(PM2.5)的贡献率明显上升,分别上升了7.56%、6.18%.  相似文献   

10.
为了深入认识宁波市冬季细颗粒物(PM2.5)的污染特征和主要影响因素的作用规律,利用Models-3/CMAQ模式系统对2013年1月宁波市的PM2.5污染形成过程进行了模拟分析.结果表明,宁波市PM2.5的重点污染区主要分布在市区、北部地区及东部沿海,除了受到局地污染源排放的影响外,对比非污染的情况,大气输入和气溶胶生成作用的增强是引起PM2.5污染的主导因素,其中水平传输过程对PM2.5浓度升高的贡献最为突出.气溶胶过程的贡献在近地面(0~80 m)最显著,随着高度升高而逐渐减弱.硝酸盐在局地二次生产的细颗粒物中占主要份额(~70%).对于硫酸盐,局地二次生成所占的比例很低,主要来自宁波局地排放和宁波以外地区的大气传输(贡献比例分别为44%和40%).宁波市的PM2.5污染主要受到来自北向沿岸气团(占比54%)、西北向大陆气团(占比21%)和西向局地气团(占比25%)的传输影响.在西北方向短距离区域传输的作用下PM2.5浓度最高;在我国中东部大范围灰霾天气的影响下,西北向和北向的长距离传输作用也会导致宁波地区的PM2.5污染.  相似文献   

11.
河南省冬季3次重污染过程的数值模拟及输送特征分析   总被引:1,自引:0,他引:1  
利用WRF-Chem模式模拟2015年11月27日—12月1日、12月5—14日、12月19—25日河南3次重污染过程,结合空气污染资料和ERA-Interim再分析资料,对比分析了这3次重污染过程的开始、持续和结束及污染物的输送特征.结果表明,静稳天气有利于污染的发展持续,3次重污染过程的结束均是由西路冷空气入侵造成的.第1次重污染过程平均风场上的风速均为小风或静风,从湖北到河南南部风向为偏南风;而第2和第3次重污染过程平均风场分别以偏东和偏北风为主.第2和第3次重污染过程中均存在明显的由北向南的污染物输送过程.3次重污染过程中,河南省本地排放对本省PM_(2.5)浓度的平均贡献率最大,而河南省周边区域对河南PM_(2.5)浓度的平均贡献率在这3次过程中不一样,第1次重污染过程,河南南部主要受偏南风影响,湖北对河南PM_(2.5)浓度的平均贡献率最大,为20.7%;第2和第3次重污染过程主要受偏东风影响,安徽和江苏对河南PM_(2.5)浓度的平均贡献率最大,分别为17.7%和18.5%.3次重污染过程中,安阳的主要污染输送源均不相同,分别来自河北、江苏和安徽、本省.  相似文献   

12.
利用2015~2019年山东省日照市PM2.5质量浓度和气象要素的小时数据,对日照市PM2.5季节污染特征和日照市海陆风特征进行了分析,并基于HYSPLIT模式计算了5年逐日02:00、08:00、14:00和20:00(BTC)的48h后向轨迹,不仅通过轨迹聚类分析和潜在源区分析探讨了日照市不同季节PM2.5主要传输路径和其轨迹污染特征及其潜在源区分布和贡献,也分析了海陆风对日照市污染物的影响.结果表明:日照市PM2.5呈现冬季最高、夏季最低的分布特征,监测站点颗粒物浓度在偏西北风影响下较高.日照市不同季节主要输送路径存在差异:春季主要受到偏东和偏北方向气流影响;夏季在副热带高压影响下主要受到来自海上的较为清洁的偏东气流影响;秋季主要受到西北和偏东气流影响;冬季主要受西北和偏北气流影响.整体而言,不同季节受偏西至偏南气流影响时,日照市对应的PM2.5浓度较高.日照市海陆风春秋季多,夏冬季少;在海陆风影响下,日照市PM2.5染和臭氧污染呈现不同的分布特征,且在不同PM2.5污染等级下,PM2.5浓度日变化特征也与其在非海陆风日的日变化有所差异.污染潜在源区分析结果表明,日照市最主要的潜在源区位于山东省临沂市、潍坊市、青岛市和江苏省连云港市.  相似文献   

13.
为探索后奥运时期京津冀区域大气本底污染状况及变化趋势,在中国科学院华北兴隆大气本底观测站,对夏季大气主要污染物(NOx、SO2、O3和PM2.5)进行了连续3年的在线观测,结合气流轨迹模式对大气污染物的传输路径及贡献率进行了分析. 结果表明:兴隆站大气本底夏季ρ(NOx)、 ρ(SO2)、 ρ(O3)和ρ(PM2.5)的平均值分别为(11.5±5.9)、(8.3±7.0)、(137.6±38.4)和(50.9±33.0)μg/m3;首要污染物为O3和PM2.5,ρ(O3)和ρ(PM2.5)日均值超过GB 3095─2012《环境空气质量标准》二级标准〔ρ(O3)为200 μg/m3,ρ(PM2.5)为75 μg/m3〕的天数分别为102和60 d,占观测期间有效天数的39%和23%,表现为大气氧化性增强、二次污染逐年上升;受偏南气流影响,太行山沿线区域和山东半岛—渤海湾地区是兴隆夏季大气本底污染的主要贡献区域,特别是京津冀城市区域夏季高浓度O3和PM2.5,对华北区域大气本底污染物浓度的整体上升具有重要影响.   相似文献   

14.
2011~2012年冬春期间(11月到翌年4月),通过设置在海伦农田生态系统国家野外科学观测研究站内的在线监测仪器获取了PM2.5和气态污染物(NOx、O3和SO2)质量浓度的时间变化,同时结合地面气象资料和HYSPLIT后向气团轨迹模型分析了该地大气污染物的污染水平、可能来源及传输过程.结果表明:观测期间PM2.5、NOx和SO2的24h均值(范围)分别为(54.7±45.7)(8.0~217.8),(23.0±11.5)(4.5~59.6), (10.0±10.3)(0.3~56.0)μg/m3, O3的日最大8h平均值(范围)为(62.4±18.7)(24.1~173.5)μg/m3,其中除O3在4月份超过国家一级标准8d外,其它气态污染物均未超过国家一级标准;PM2.5超过国家二级标准的天数为40d,占整个观测期间的22.5%.PM2.5和SO2各月质量浓度变化较大,最高值出现在12月份,是冬季采暖的高峰期.NOx、PM2.5和SO2日变化呈双峰型,峰值出现在07:00和17:00左右;O3为单峰型,峰值出现13:00~15:00.通过对海伦地区72h内HYSPLIT后向气团轨迹模拟结果和该站点的气象数据进行分析,表明该农业区大气污染受本地源和区域输送共同影响,偏南气流易造成污染物积累,而偏北气流有利于污染物扩散.  相似文献   

15.
利用2017~2019年晋城市和长治市冬季PM2.5逐时浓度资料、地面风场数据等,结合HYSPLIT轨迹模型和中尺度数值模式WRFV4.2分析了晋东南地区冬季PM2.5污染的特征和传输特点.结果表明,晋城市冬季PM2.5污染程度高于长治市.受地形影响,晋城市地面盛行偏南风、偏北风和西北风,污染方向主要为偏南风和偏北风;长治市近地面盛行偏南风,该风向污染频率最高.影响晋城市和长治市污染的潜在源区主要分布在偏西、东北和东南方向,偏西气流来自陕西省中部,东北气流来自河北省西南部,东南气流来自河南省中东部.污染经过晋东南地区主要影响山西省中南部和北京南部.通过数值模拟流场,结合潜在源区和影响区域的分析结果,在均压场或高压后部的天气形势下,晋东南地区污染输送路径包括来自东北方向(河北省西南部一带)的气流,沿长治市东北部的滏口陉向晋东南地区输送污染物及沿太行山东麓向南在晋豫交界处的太行陉发生转折向晋东南地区输送污染物;来自东南方向(河南北部及东部)的气流输送和来自偏西方向(陕西中南部)的气流输送.污染物经过晋东南地区向北输送至山西省中南部,部分经过山西省中东部的井陉输送至北京南部.  相似文献   

16.
为了探究边界层气象要素时空分布及其变化对银川市冬季持续污染天气过程污染物质量浓度变化的影响机制,利用2016年12月1日-2017年1月31日逐时空气质量以及地面和逐日定时探空气象观测数据,根据大气污染级别和过程持续时间,选取2016年12月9-21日(简称"1211过程")和2016年12月29日-2017年1月9日(简称"1231过程")为研究对象,采用统计和天气诊断相结合的方法,在分析比较银川市冬季两次典型持续污染过程演变特征及其与地面气象要素关系的基础上,探讨了大气环流、边界层要素变化对银川市冬季典型污染过程的可能影响机制.结果表明:①银川市冬季两次大气污染过程持续阶段,地面均以偏东或偏南风为主,风速较小,相对湿度较大,能见度较低;在污染清除阶段,地面风向转为西北或偏北风,风速较大,相对湿度较小,能见度较高.②当冬季欧亚大陆中纬度区域500 hPa高空盛行纬向气流,850 hPa高度上银川市受反气旋环流和暖温度脊控制,并且有弱暖平流从西南部向北输送时,银川市易出现静稳型持续污染天气.③冬季银川市持续大气污染过程中,ρ(PM2.5)与风速呈负相关(R平均值为-0.326),与相对湿度呈正相关(R平均值为0.688),与能见度呈显著负相关(R平均值为-0.905),与边界层高度呈较显著负相关(R平均值为-0.575).④银川市冬季静稳型持续污染天气主要分为弱西北和平直西风气流型两种,弱西北气流型具有近地面层逆温弱,污染物积累慢,清除快的特征;平直西风气流型具有近地面层逆温强,污染物积累快,清除慢的特征.研究显示,冬季银川市上空500 hPa高度盛行纬向气流,地面主导风向为偏东或偏南风时,随着地面相对湿度增大、近地层风速减小、大气垂直上升运动减弱、边界层高度降低,大气中ρ(PM2.5)将迅速升高,银川市易出现以PM2.5为首要污染物的静稳型持续污染天气.   相似文献   

17.
2014年京津冀地区PM2.5浓度时空分布及来源模拟   总被引:3,自引:0,他引:3  
采用模式(CAMx)模拟与污染物、气象观测资料相结合的方式,分析了2014年京津冀地区PM2.5时空分布及来源特征.结果表明:PM2.5具有较为明显的时间变化规律,呈秋冬高、春夏低的规律和双峰型分布的日变化特征;重污染日PM2.5高浓度(PM2.5>150μg/m3)主要分布在太行山前的华北平原区,特别是北京、保定、石家庄一线,而太行山、燕山等西部及北部山区PM2.5浓度明显低于平原区;重污染日京津冀地区PM2.5平均浓度在150μg/m3以上的面积约占总面积的73%;重污染日北京、天津、石家庄市的PM2.5外来输送率分别为58%、54%、39%;2014年10月6~12日京津冀地区发生的一次重污染过程中污染物由南向北输送,区域输送对于各地区PM2.5浓度有着十分重要的影响.  相似文献   

18.
珠江三角洲典型灰霾过程的边界层特征   总被引:3,自引:0,他引:3  
2014年1月初珠江三角洲(以下简称珠三角)地区出现严重的灰霾天气.利用广东省城市空气质量资料、中国气象局番禺大气成分站逐时PM10、PM2.5、PM1以及能见度数据、广东省946个自动站风向风速数据、广州基本站风向风速数据、2013~2014年东莞佛山大气边界层观测试验资料等,研究了珠三角2014年1月1~8日典型灰霾过程的边界层特征.结果表明:大陆冷高压变性出海是这次灰霾过程发生的主要天气系统.其使得珠三角地区受静稳天气控制,造成微风或静风状况以及低的边界层和逆温结构;6日珠三角主要城市空气质量呈现不同的变化,其原因主要是区域内风向的变化,5日上午珠三角为弱的偏北风,整个珠三角都为灰霾天气,而5日下午风向由偏北风转为偏南风,珠三角东南部地区空气质量改善,北部(内陆)地区的空气质量恶化;后向轨迹分析发现所有气流呈明显下沉作用,且5日20:00的后向轨迹簇来自陆面,6日20:00的后向轨迹簇来自海面,两者对珠三角的空气质量具有不同的影响.  相似文献   

19.
北京PM2.5背景值定值方法及其变化特征研究   总被引:1,自引:0,他引:1  
通过气团识别方式,确定了较强气团判别指标,定量评估了较强气团对北京上甸子区域大气本底站PM2.5浓度影响.结果显示,偏北风过程中PM2.5的背景值介于10.3~13.5μg/m3,平均值为(11.7±1.3)μg/m3;偏南风过程中PM2.5的背景值介于60.2~92.6μg/m3,平均值为(76.2±12.8)μg/m3.偏北风过程中,PM2.5背景值变化趋势不明显;偏南风过程中,PM2.5背景浓度呈线性下降趋势,下降速率为5.5μg/(m3×a).2009~2012年期间,区域本底站观测到的PM2.5年均值变化不大,这与北京地区的局地产生的二次PM2.5逐步增大有关,在一定程度上抵消了区域背景值下降的影响.当出现系统性偏南风时,北京以南区域输送对北京当地PM2.5浓度增加的平均贡献接近60%,但是自2009年开始逐步下降,截止2012年该贡献降低到44%.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号