期刊文献+
共找到14篇文章
< 1 >
每页显示 20 50 100
RTTOV散射模块在FY-3MWHS资料同化中的应用 被引量:5
1
作者 苏捷 朱克云 张杰 《气象》 CSCD 北大核心 2013年第11期1461-1472,共12页
一直以来,受同化系统中观测算子的限制,云区的卫星微波观测资料在使用时都常会做晴空假设,只考虑大气的吸收作用而忽略散射作用,从而会造成对大气中水汽含量的高估。为了克服这一缺点,文中设计了RTTOV微波散射模块在WRF-3DVAR系统中的... 一直以来,受同化系统中观测算子的限制,云区的卫星微波观测资料在使用时都常会做晴空假设,只考虑大气的吸收作用而忽略散射作用,从而会造成对大气中水汽含量的高估。为了克服这一缺点,文中设计了RTTOV微波散射模块在WRF-3DVAR系统中的应用方案,使得计算辐射能衰减的同时考虑大气的散射作用,并在此基础上通过4组试验进行了风云3-A星微波湿度计亮温以及多普勒雷达反射率因子的同化研究对比。结果表明:RTTOV微波散射模块的引入,能有效抑制晴空辐射模块对大气水汽含量的高估。若先同化雷达反射率因子后再继续同化卫星资料,会发现同化雷达后,同化区域内大气水凝物含量明显增加,受此影响同化卫星资料时大气的散射作用也会随之增大。并且雷达同化区域覆盖率越大,对降水预报的改善效果也愈明显。 展开更多
关键词 资料同化 风云三号 微波湿度计 RTTOV—SCATT 雷达反射率因子
在线阅读 下载PDF
基于Arase卫星数据的地球电子外辐射带数据同化建模
2
作者 罗陶荣 倪彬彬 +2 位作者 曹兴 王晓宇 雷源 《地球物理学报》 北大核心 2025年第4期1187-1200,共14页
Arase卫星是一颗专门探测地球辐射带的科学卫星,于2016年12月20日成功发射.在2019年10月范艾伦卫星终止运行后,Arase卫星仍然继续正常运行并提供高质量的观测数据.本文结合Arase卫星高能电子探测器的电子通量观测数据和基于Fokker-Planc... Arase卫星是一颗专门探测地球辐射带的科学卫星,于2016年12月20日成功发射.在2019年10月范艾伦卫星终止运行后,Arase卫星仍然继续正常运行并提供高质量的观测数据.本文结合Arase卫星高能电子探测器的电子通量观测数据和基于Fokker-Planck扩散方程的三维辐射带数值模型,利用卡尔曼滤波方法实现对2017年4月至2022年3月期间地球外辐射带(L^(*)=3~7)电子通量时空演化过程的数据同化建模.结果显示,同化模型有效弥补了Arase卫星数据对于高赤道投掷角电子的观测缺失,重构了不同能量、投掷角辐射带电子在不同L^(*)处的长期动态演化过程及其对地磁活动的响应特征.通过计算电子通量的数据同化结果与卫星观测结果之间的皮尔逊相关系数、平均误差、平均绝对误差和均方误差,发现电子通量同化结果和观测结果整体一致.通过计算同化过程中的新息矢量,分析卫星观测结果对辐射带数值模型模拟结果的修正程度及其对电子能量和L^(*)的依赖性,这为分析辐射带数值模型所缺失的物理机制以及后续模型的改进与完善提供了重要参考. 展开更多
关键词 Arase卫星 地球辐射带 电子通量 三维数据同化 模型性能
在线阅读 下载PDF
能见度观测同化对天津地区雾霾低能见度数值预报影响研究
3
作者 杨旭 樊文雁 +3 位作者 蔡子颖 朱玉强 唐颖潇 董琪如 《中国环境科学》 北大核心 2025年第11期5959-5968,共10页
以天津地区两次典型雾霾低能见度过程为例,通过发展适用于雾霾能见度观测的资料同化方法,实现PM_(2.5)浓度、相对湿度及雾水含量等影响要素的模拟优化,提升模式能见度数值预报效果.研究表明,典型霾过程能见度变化主要受PM_(2.5)浓度影响... 以天津地区两次典型雾霾低能见度过程为例,通过发展适用于雾霾能见度观测的资料同化方法,实现PM_(2.5)浓度、相对湿度及雾水含量等影响要素的模拟优化,提升模式能见度数值预报效果.研究表明,典型霾过程能见度变化主要受PM_(2.5)浓度影响,最低值1km左右,随相对湿度日变化幅度为2~3km,雾霾混合过程相对湿度高,夜间雾显著影响导致局地能见度低于50m,能见度日变幅大、过程极值低是其典型特征;通过能见度观测资料同化应用,雾霾低能见度模拟效果有显著提升,两次过程的能见度均方根误差分别由7km和4km减小至2km左右,间隔6h同化可以及时修正模拟偏差,部分时段提升效果更优;通过同化实现PM_(2.5)浓度和相对湿度模拟优化,霾能见度模拟偏差下降幅度最大达到6~8km;通过增加雾水含量修正,结合相对湿度同化改进,能够进一步改善雾霾混合过程夜间雾水含量考虑不足的问题,能见度模拟偏差由1~2km减小至300~500m,低能见度影响程度、影响范围均有明显提升;一个月模拟试验显示,低于10km的能见度模拟均方根误差由4.4km减小至3.2km. 展开更多
关键词 雾霾 能见度 资料同化 三维变分 数值预报
在线阅读 下载PDF
多普勒天气雷达资料对中尺度模式短时预报的影响 被引量:68
4
作者 盛春岩 浦一芬 高守亭 《大气科学》 CSCD 北大核心 2006年第1期93-107,共15页
利用中尺度模式ARPS(The Advanced Regional Prediction System)及其资料分析系统ADAS(ARPS Data Analysis System),将国内新一代多普勒雷达(CINRAD)反射率及径向风资料直接用于中尺度数值模拟,通过一次华北地区暴雨过程的模拟... 利用中尺度模式ARPS(The Advanced Regional Prediction System)及其资料分析系统ADAS(ARPS Data Analysis System),将国内新一代多普勒雷达(CINRAD)反射率及径向风资料直接用于中尺度数值模拟,通过一次华北地区暴雨过程的模拟对比试验,分析了雷达资料对初始场的改进效果及其对模拟结果的影响,结果表明:(1)利用雷达径向风资料对初始风场进行调整后,自近地面到对流层顶的u,υ,ω都发生了变化,调整后的初始风场在对流层中层变化最大。(2)利用雷达反射率进行微物理调整和云分析能调整初始场中的云水信息,使得雷达回波附近3km以下的水汽混合比(qv)增加,4km以下的雨水混合比(qr)增加,对流层(约10km以下)的云水混合比(qc)增加,4~9km的对流层上部云冰混合比(qi)和雪混合比(qs)增加。ADAS通过非绝热初始化调整温度场,从而得到了一个动力和热力上平衡的初始场。(3)模拟的1h雨量与实况的对比表明,同时利用雷达反射率和径向风改进过的初始场能明显增强3h内的降水强度和落区预报,改善中尺度数值模式短时定量降水预报。模拟的1h流场对比分析表明,经雷达径向风调整后,能够在初始场中增加气旋性涡旋等中小尺度风信息,明显减少模式的spin-up时间。(4)通过对雷达径向风和反射率对模式初始场和模拟结果影响的对比分析发现,雷达径向风主要是改进初始风场,而雷达反射率主要是改进初始场中的湿度参数,增加初始场中云水等的含量,调整温度场。通过模拟的6h降水对比发现,利用雷达径向风调整初始场后,对降水模拟有一定的改进,但效果不甚明显,而雷达反射率资料对定量降水预报改进效果明显,同时使用雷达径向风和反射率资料改进初始场后对降水的模拟效果最明显。 展开更多
关键词 ARPS模式 雷达资料同化 初始场 3h短时预报
在线阅读 下载PDF
最优插值方法在西北太平洋海温同化中的应用研究 被引量:10
5
作者 李云 刘钦政 +2 位作者 张建华 王旭 李燕 《海洋预报》 2008年第2期25-32,共8页
本文开发出一套与三维温盐流数值模式POM相匹配的最优插值数值同化模块,具有将船舶报资料和Argo海温数据加入数值模拟的能力。利用该同化系统,本文对2001~2002年进行了数值模拟试验,通过模拟结果和未加入同化的模拟结果的对比表明,该... 本文开发出一套与三维温盐流数值模式POM相匹配的最优插值数值同化模块,具有将船舶报资料和Argo海温数据加入数值模拟的能力。利用该同化系统,本文对2001~2002年进行了数值模拟试验,通过模拟结果和未加入同化的模拟结果的对比表明,该同化系统能有效地改进三维海温模拟结果,使之更为接近于观测值。 展开更多
关键词 资料同化 最优插值 船舶报资料 ARGO资料 三维海洋数值模式
在线阅读 下载PDF
淮河流域东北部一次异常特大暴雨的数值模拟研究Ⅰ:结果检验和β中尺度对流系统的特征分析 被引量:15
6
作者 王亦平 陆维松 +1 位作者 潘益农 王元 《气象学报》 CSCD 北大核心 2008年第2期167-176,共10页
在2000年12号台风(Prapiroon)影响期间,其外围对淮河流域东北部造成了一次罕见的特大暴雨,暴雨中心响水24h降水量达到800mm。文中所用的模式是俄克拉荷马大学风暴分析和预测中心研制的一个三维非静力可压缩数值区域预报模式ARPS(V5.2)... 在2000年12号台风(Prapiroon)影响期间,其外围对淮河流域东北部造成了一次罕见的特大暴雨,暴雨中心响水24h降水量达到800mm。文中所用的模式是俄克拉荷马大学风暴分析和预测中心研制的一个三维非静力可压缩数值区域预报模式ARPS(V5.2)。采用了3层单向嵌套网格,Domain1中心取为(27.5°N,117.5°E),格距45km,格点数为75×75;Do-main2中心取为(31.5°N,119.5°E),格距15km,格点数为140×140;Domain3中心取为(33.5°N,119.5°E),格距5km,格点数为180×180;垂直方向分为35层,垂直格距为625m。所利用的资料为:2000年8月29—31日每日4个时次(00、06、12、18时)1°×1°的NCEP/NCAR再分析资料及713雷达资料、GMS-5红外云图、探空报、地面加密资料。物理过程选用简单冰相方案,Kain和Fritsch积云参数化方案。对逐时的卫星云图、713雷达图像进行数值反演,结合探空资料反演出暴雨中深对流系统的水汽三维分布情况,通过三维同化系统ADAS处理,同化初始场和侧边界条件,再用ARPS模式进行数值积分,从8月29日08时开始到31日08时结束,积分48h。并结合雷达资料、红外云图、探空报、地面加密资料等对数值模拟结果进行了对比分析。结果表明:ARPS较好地模拟了在台风移动的左前方、在高空槽与副高之间出现的中尺度强暴雨区,模拟降水区及暴雨中心位置与实况较为一致。利用模拟大气中的水物质模拟了雷达回波,与实际雷达回波进行了对比分析,揭示了该过程中β中尺度对流系统的演变特征,4条对流带的交汇点在响水附近。并将模式模拟的T-lgp图与实际T-lgp图进行了对比分析,揭示了该次暴雨过程的不稳定性。持续的高空风垂直切变为对流系统的发展提供动能,造成对流系统斜压发展,有利于降水集中在某一固定的地点。由于模拟结果与实况较为接近,因此可以利用模拟结果作为对该暴雨过程作进一步研究的基础。 展开更多
关键词 特大暴雨 数值模拟 三维同化 雷达回波 T-logp图
在线阅读 下载PDF
国家气候中心全球海洋资料四维同化系统在热带太平洋的结果初步分析 被引量:9
7
作者 刘益民 李维京 张培群 《海洋学报》 CAS CSCD 北大核心 2005年第1期27-35,共9页
介绍了国家气候中心(NCC)“全球海洋资料四维同化系统”(简称NCC-GODAS).该系统包含观测资料预处理系统、插值分析系统和所应用的动力模式.插值分析系统采用四维同化技术方案,在时间上设置一个4周的窗口,将此窗口之内的观测资料以一定... 介绍了国家气候中心(NCC)“全球海洋资料四维同化系统”(简称NCC-GODAS).该系统包含观测资料预处理系统、插值分析系统和所应用的动力模式.插值分析系统采用四维同化技术方案,在时间上设置一个4周的窗口,将此窗口之内的观测资料以一定的权重插入插值分析系统, 在空间上采用三维变分方案.海洋动力模式为“九五”期间LASG研制的L30T63OGCM 1 0版本.文中考察了该系统从1982年到2003年3月在热带太平洋的部分同化分析结果,并与NCEP的再分析资料和EMC/NCEP的太平洋区域海洋资料同化系统的结果进行了对比分析.结果显示,该系统的同化结果(如SST,SSTA,Nino指数,次表层海温变化等)与NCEP的同期同化结果具有很好的一致性.同时,该系统得出的海洋同化资料气候场与海洋模式气候场相比,有显著的改进.表明该系统具有较好的同化能力,其同化结果可为海-气耦合模式进行季节和跨季节的气候预测业务提供可信的海洋初始场资料,同时还可以为相关研究提供海洋分析场. 展开更多
关键词 海洋资料同化 资料预处理 海洋动力模式 三维变分 时间窗口 SST Nino指数
在线阅读 下载PDF
一个用于气候模式的三维变分资料同化系统——方案设计及检验 被引量:9
8
作者 官元红 周广庆 陆维松 《大气科学学报》 CSCD 北大核心 2009年第3期392-404,共13页
介绍了一种适合格点模式的三维变分同化方案,并应用于中国科学院大气物理研究所的9层大气环流格点模式(IAP9L2°×2.5°-AGCM)中。同化方案考虑质量场和风场的动力约束关系,并通过变量变换将模式变量转化为互不相关的控制变... 介绍了一种适合格点模式的三维变分同化方案,并应用于中国科学院大气物理研究所的9层大气环流格点模式(IAP9L2°×2.5°-AGCM)中。同化方案考虑质量场和风场的动力约束关系,并通过变量变换将模式变量转化为互不相关的控制变量:流函数、非平衡速度势、非平衡位势和比湿;同时采用共轭梯度法求解代价函数避免了直接计算背景误差协方差矩阵逆的困难。通过理想试验、模式连续同化积分试验以及对东亚地区特殊寒潮天气的同化试验,结果表明:该系统能够正确反应多变量之间的相互关系,并且长时间运行稳定,更重要的是同化不仅可以明显改善直接参与同化的物理量,而且经过模式的动力调整也改善了未参与同化的物理量(如海平面气压、近地面气温、降水等),这为大气的持续性信息存储于陆面物理量提供了可能,从而可为IAP9L2°×2.5°-AGCM提供陆—气相协调的初始场。 展开更多
关键词 资料同化 气候模式 三维变分 短期气候预测
在线阅读 下载PDF
COSMIC资料在GRAPES全球三维变分同化系统的初步研究 被引量:17
9
作者 唐细坝 薛纪善 《热带气象学报》 CSCD 北大核心 2009年第5期521-531,共11页
利用COSMIC获得的GPS无线电掩星数据所反演得到的大气温湿度廓线资料,在GRAPES全球三维变分同化系统中做连续循环同化试验,研究了加入COSMIC资料后得到的GRAPES分析场及其5天中期预报是否有所改善。研究结果表明:通过加入COSMIC资料连... 利用COSMIC获得的GPS无线电掩星数据所反演得到的大气温湿度廓线资料,在GRAPES全球三维变分同化系统中做连续循环同化试验,研究了加入COSMIC资料后得到的GRAPES分析场及其5天中期预报是否有所改善。研究结果表明:通过加入COSMIC资料连续循环同化得到的分析场,相对原来没有加COSMIC资料连续同化出来的GRAPES分析场有明显改善,由此GRAPES全球模式全球5天中期预报水平有明显提高,其中加入没有稀疏化后COSMIC资料对南半球GRAPES分析场的改善以及GRAPES全球模式5天中期预报水平的提高都比较显著;加入稀疏化后的COSMIC资料对北半球GRAPES分析场的改善以及GRAPES全球模式5天中期预报水平的提高都有比较显著的效果。 展开更多
关键词 数值预报 COSMIC资料同化 GPS无线电掩星 GRAPES全球三位变分同化系统
在线阅读 下载PDF
海洋模态在海洋水声学中的潜在应用前景 被引量:3
10
作者 崔茂常 乔方利 +1 位作者 莫军 郭炳火 《海洋学报》 CAS CSCD 北大核心 2002年第6期127-134,共8页
关键词 海洋水声学 海洋模态 Carton海洋四维同化数据库 声层析技术 反演 海洋调查
在线阅读 下载PDF
基于Fortran语言的地球外辐射带电子三维数据同化建模 被引量:1
11
作者 罗陶荣 曹兴 +4 位作者 马新 付松 倪彬彬 王晓宇 雷源 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2024年第4期1285-1298,共14页
开发地球电子辐射带的数据同化模型,对于理解辐射带电子的动态演化过程和辐射带空间天气预报具有重要意义.结合范阿伦卫星的辐射带电子观测数据和外辐射带三维扩散模型,采用卡尔曼滤波算法,本文开发了基于Fortran语言的外辐射带电子三... 开发地球电子辐射带的数据同化模型,对于理解辐射带电子的动态演化过程和辐射带空间天气预报具有重要意义.结合范阿伦卫星的辐射带电子观测数据和外辐射带三维扩散模型,采用卡尔曼滤波算法,本文开发了基于Fortran语言的外辐射带电子三维数据同化模型(Three-dimensional Data Assimilative Model of Outer Radiation belt Electrons,简称TDAMORE),实现对L^(*)=3~7、能量范围为0.1~5 MeV、投掷角范围为5°~90°的外辐射带电子时空变化过程的三维重构.通过对2018年8月期间外辐射带电子通量演化过程的重构,证实TDAMORE模型可以较好地重现不同能量和不同投掷角电子通量在磁暴前后的演化特征.通过分析电子通量的观测和同化结果之间的相关系数、平均误差、平均绝对误差和均方误差,发现对于能量低于4 MeV的电子,观测与同化结果之间的相关系数基本大于0.8且误差相对较低.而对于更高能量的电子,观测与同化结果之间的误差相对较高,这可能是同化模型忽略了电磁离子回旋波对电子的散射损失导致的. 展开更多
关键词 地球外辐射带 电子通量 三维数据同化 范阿伦卫星 FORTRAN语言
在线阅读 下载PDF
2018年中国中东部地区臭氧浓度分布的时空特征 被引量:5
12
作者 毛毛 彭珍 《南京大学学报(自然科学版)》 CAS CSCD 北大核心 2021年第3期460-465,共6页
近几年,受人为排放的挥发性有机物逐年增加和大气中PM2.5浓度逐年降低的影响,近地面的O_(3)已成为影响我国夏季空气质量的首要污染物,全国很多地区每年5月至9月会出现O_(3)污染超标的现象和O_(3)大范围蔓延的趋势.利用2018年逐时的O_(3... 近几年,受人为排放的挥发性有机物逐年增加和大气中PM2.5浓度逐年降低的影响,近地面的O_(3)已成为影响我国夏季空气质量的首要污染物,全国很多地区每年5月至9月会出现O_(3)污染超标的现象和O_(3)大范围蔓延的趋势.利用2018年逐时的O_(3)浓度分析场数据对我国中东部尤其是京津冀、长三角和珠三角地区O_(3)浓度的时空变化趋势进行研究.结果表明,随着时间变化,月平均O_(3)浓度在京津冀地区呈倒“V”型,在长三角和珠三角地区呈“M”型,这是因为O_(3)生成与地表接受的短波辐射通量有较好的线性关系.在京津冀地区6月太阳辐射强度最大且云量较小,因此O_(3)污染严重,O_(3)超标天数接近一半;而长三角和珠三角地区在7月受雨季的影响,云量较大,因此O_(3)月平均浓度相对较小,月平均浓度变化呈现出“M”型.而在一天内O_(3)浓度在上午逐渐上升,午后达到最大且稍滞后于地表所受到的短波辐射通量的日变化. 展开更多
关键词 O_(3)污染 资料同化 “M”型双峰特征 地表短波辐射通量
在线阅读 下载PDF
基于气溶胶三维变分同化天津PM_(2.5)数值预报研究 被引量:9
13
作者 杨旭 唐颖潇 +5 位作者 蔡子颖 韩素芹 董琪如 杨健博 朱玉强 樊文雁 《中国环境科学》 EI CAS CSCD 北大核心 2021年第12期5476-5484,共9页
基于气溶胶三维变分同化技术,建立天津空气质量数值模式气溶胶同化模块,通过天津地区两次重污染过程同化模拟敏感性试验,分析了观测资料范围对同化结果的影响,并结合一个月的滚动预报试验,分析了气溶胶同化对天津地区PM_(2.5)数值预报... 基于气溶胶三维变分同化技术,建立天津空气质量数值模式气溶胶同化模块,通过天津地区两次重污染过程同化模拟敏感性试验,分析了观测资料范围对同化结果的影响,并结合一个月的滚动预报试验,分析了气溶胶同化对天津地区PM_(2.5)数值预报效果的影响,以期为提升天津空气质量预报能力提供支撑.结果表明:气溶胶同化各控制变量背景误差水平相关系数的衰减尺度约50km,垂直方向上400m高度与模式底层的相关系数衰减至0.6左右;观测资料范围对同化结果影响显著,仅采用天津地区观测数据进行同化,对重污染过程期间天津地区PM_(2.5)浓度模拟的影响时效约12h,采用模拟区域内所有观测数据进行同化影响时效可持续24h以上,且模拟效果更优;采用三维变分同化技术,实现地面PM_(2.5)观测资料同化,天津地区PM_(2.5)数值预报效果显著提升,预报值和实况值之间的相关系数由0.74增加到0.87,均方根误差由32.3μg/m^(3)减小为22.4μg/m^(3),平均相对误差由39.9%减小为27.1%;同化对模式初始时刻的改进效果最明显,随时间同化效果衰减,14h内改进效果最佳,对24h PM_(2.5)浓度预报也有明显改进. 展开更多
关键词 气溶胶 资料同化 三维变分 数值预报
在线阅读 下载PDF
地面PM_(2.5)和气溶胶激光雷达联合变分同化对天津冬季一次重污染过程数值模拟改进研究 被引量:1
14
作者 杨旭 樊文雁 +3 位作者 蔡子颖 朱玉强 唐颖潇 董琪如 《中国环境科学》 EI CAS CSCD 北大核心 2023年第10期5188-5197,共10页
针对输送型重污染过程,大气化学模式垂直模拟偏差较大,地面观测同化对污染垂直初始场改进不明显,高空输送模拟偏低的问题,以2023年1月5~6日天津地区一次重污染过程为例,通过PM_(2.5)地面观测和气溶胶激光雷达垂直观测联合应用,发展相应... 针对输送型重污染过程,大气化学模式垂直模拟偏差较大,地面观测同化对污染垂直初始场改进不明显,高空输送模拟偏低的问题,以2023年1月5~6日天津地区一次重污染过程为例,通过PM_(2.5)地面观测和气溶胶激光雷达垂直观测联合应用,发展相应的气溶胶变分同化技术,改进大气化学模式污染垂直模拟效果,服务重污染天气预报预警和成因解析.研究表明:此次重污染过程由于大气化学模式垂直模拟能力不足,导致S1发展阶段地面PM_(2.5)增长速率和S2维持阶段的峰值浓度明显低估,垂直分布上污染层高度模拟偏低200~300m,白天输送阶段高空PM_(2.5)浓度高于地面、夜间稳定边界层以上污染高值特征均未能有效反映;通过同化地面PM_(2.5)观测有效实现PM_(2.5)空间分布模拟优化,过程期间上游城市石家庄和衡水的模拟误差分别下降35.3μg/m^(3)和77.3μg/m^(3),天津地面PM_(2.5)质量浓度模拟误差下降39.7μg/m^(3),但地面同化对垂直分布模拟改进有限,尤其是对于输送型污染过程中PM_(2.5)垂直分布的复杂变化,基本无法准确反映;通过增加气溶胶激光雷达PM_(2.5)垂直观测同化,模式PM_(2.5)垂直模拟效果进一步提升,450m左右均方根误差减小幅度最大,约19.8μg/m^(3).在同化时刻及临近2~3h正效果明显,后续时效中S2阶段表现更优,夜间稳定边界层以上污染高值得到有效模拟,污染层高度模拟差异减小至50~100m,低于未同化CTRL试验(200~300m)和地面同化DASFC试验(100~200m);S1阶段随着时效延长垂直模拟偏差明显增加,但同化的正作用依旧存在,污染层高度模拟误差比DASFC试验降低幅度不低于50m. 展开更多
关键词 气溶胶同化 三维变分 激光雷达 数值模拟
在线阅读 下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部