期刊文献+
共找到33篇文章
< 1 2 >
每页显示 20 50 100
厦门九龙江河口区无机氮和活性磷酸盐含量分布的计算 被引量:1
1
作者 潘伟然 骆智斌 +2 位作者 张国荣 王君 马腾 《台湾海峡》 CAS CSCD 北大核心 2009年第3期331-335,共5页
本文利用三维水动力和盐扩散方程,计算厦门湾九龙江河口区的盐度时空分布,并利用监测调查所得到的盐度与无机氮、活性磷酸盐的经验关系,对厦门九龙江河口区的无机氮、活性磷酸盐的三维时空分布进行计算,对监测站位的计算值和实测值进行... 本文利用三维水动力和盐扩散方程,计算厦门湾九龙江河口区的盐度时空分布,并利用监测调查所得到的盐度与无机氮、活性磷酸盐的经验关系,对厦门九龙江河口区的无机氮、活性磷酸盐的三维时空分布进行计算,对监测站位的计算值和实测值进行的比较和分析,结果说明在潮周变化尺度,物理混合及扩散对无机氮、活性磷酸盐的分布起主要作用. 展开更多
关键词 海洋环境科学 九龙江河口区 盐度 无机氮 活性磷酸盐 数值计算
在线阅读 下载PDF
湄洲湾海水交换率和半更换期的计算 被引量:15
2
作者 潘伟然 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 1992年第1期65-68,共4页
利用湄洲湾海域现场观测数据,采用单箱模型,二线数值模型分别计算了湄洲湾海水的平均交换率,平均半更换期和各区段海水的半更换期,计算结果表明:湄洲湾海水的平均交换率为0.127,平均半更换期为7.1d,湾顶附近水域的半更换期超过19d。
关键词 海水交换 半更换期 单箱模型
在线阅读 下载PDF
台风期间厦门湾潮流边界层特征研究
3
作者 潘伟然 骆智斌 +1 位作者 张国荣 马腾 《台湾海峡》 CAS CSCD 北大核心 2009年第4期540-545,共6页
本文基于2007年8月圣帕台风影响台湾海峡期间的厦门湾垂向海流流速剖面观测资料,运用Karman-Prandtl模型统计和分析了流速对数剖面的出现频率和边界层参数(摩阻流速与粗糙长度).结果表明流速剖面满足对数分布出现频率都超过50%,但进行... 本文基于2007年8月圣帕台风影响台湾海峡期间的厦门湾垂向海流流速剖面观测资料,运用Karman-Prandtl模型统计和分析了流速对数剖面的出现频率和边界层参数(摩阻流速与粗糙长度).结果表明流速剖面满足对数分布出现频率都超过50%,但进行内部一致性原则分析后,流速剖面满足对数分布情况均不足30%.造成流速剖面偏离对数分布的因素主要有:台风、径流、流速加减速、水体分层、底床结构、海底悬沙浓度等等.台风期间的摩阻流速和粗糙长度都大于非台风期间的数值. 展开更多
关键词 物理海洋学 潮流边界层 摩阻流速 粗糙长度 台风 厦门湾
在线阅读 下载PDF
厦门港二维盐度场的数值计算
4
作者 潘伟然 陈金泉 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 1991年第1期69-72,共4页
由二维流体动力学和盐扩散方程所组成的偏微分方程组出发,以有限差分法和欧拉—拉格朗日方法分别求解潮波运动方程和盐扩散方程,计算结果较理想地反映了厦门湾潮流场和盐度场的分布规律,此计算方法还可用于研究该海区的污染物迁移扩散... 由二维流体动力学和盐扩散方程所组成的偏微分方程组出发,以有限差分法和欧拉—拉格朗日方法分别求解潮波运动方程和盐扩散方程,计算结果较理想地反映了厦门湾潮流场和盐度场的分布规律,此计算方法还可用于研究该海区的污染物迁移扩散问题。 展开更多
关键词 盐度场 数值计算 污染物 厦门港
在线阅读 下载PDF
福建省台风灾害直接经济损失预评估模型 被引量:13
5
作者 于小龙 潘伟然 +1 位作者 张国荣 骆智斌 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2011年第6期1047-1052,共6页
基于2000—2009年对福建省造成灾害损失的22个台风数据,将台风分为气象型灾害和混合型灾害,利用数理统计的方法建立台风灾害直接经济损失预评估模型,在得到致灾因子和损失因子的多项式拟合关系基础上,进行了台风损失的悲观估计和乐观估... 基于2000—2009年对福建省造成灾害损失的22个台风数据,将台风分为气象型灾害和混合型灾害,利用数理统计的方法建立台风灾害直接经济损失预评估模型,在得到致灾因子和损失因子的多项式拟合关系基础上,进行了台风损失的悲观估计和乐观估计,并对台风灾害的影响等级进行了评估.研究结果显示该模型具有较高的拟合率,由误差分布和实际评估结果计算的模型的有效率均达到80%以上,能够较好地预评估台风可能造成的直接经济损失并能有效地分析台风灾害的趋势;灾害等级预报的准确性稍低于灾害经济损失预报的准确性,但总体属于同一水平.模型的建立可为防灾减灾提供科学依据,具有实际意义. 展开更多
关键词 台风 直接经济损失 灾害等级 数理统计 预评估模型
在线阅读 下载PDF
九龙江口-厦门湾三维潮流数值模拟 被引量:10
6
作者 骆智斌 潘伟然 +2 位作者 张国荣 马腾 王君 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2008年第6期864-868,共5页
采用浅水半隐式三维数值模型,考虑九龙江径流影响,运用垂向水平分层法,计算了包括九龙江河口、厦门西港、东海域、同安湾,以及大小金门的厦门湾全海域各固定垂直分层的水平和垂直流速分量.由于在浅滩处采用了干湿动边界处理方法,较好地... 采用浅水半隐式三维数值模型,考虑九龙江径流影响,运用垂向水平分层法,计算了包括九龙江河口、厦门西港、东海域、同安湾,以及大小金门的厦门湾全海域各固定垂直分层的水平和垂直流速分量.由于在浅滩处采用了干湿动边界处理方法,较好地重现了涨落潮造成的潮间带间歇性出没,提高了模型的计算精度和适用范围,计算结果与实测资料吻合良好.通过讨论各层的潮流场分布特征,显示厦门湾海域潮流以正规的半日潮为主,大部分海区潮流运动形式都为往复流,其中围头湾、大小金门水道、厦门东侧水道以及青屿水道为强流区. 展开更多
关键词 潮流 三维数值模拟 九龙江 厦门湾
在线阅读 下载PDF
福建东山湾海浪现场观测的统计特征 被引量:13
7
作者 叶雨颖 潘伟然 +2 位作者 张国荣 董剑 林道平 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2007年第3期386-389,共4页
为分析福建省南部港湾──东山湾的波浪特征,采用S4型压力式方向波浪仪观测获取现场一个月的风浪过程资料,对其进行了统计分析和频谱分析,给出了波高和周期分布,特征波要素与频谱关系,结果表明:实测波高分布与瑞利分布基本相符;波浪谱... 为分析福建省南部港湾──东山湾的波浪特征,采用S4型压力式方向波浪仪观测获取现场一个月的风浪过程资料,对其进行了统计分析和频谱分析,给出了波高和周期分布,特征波要素与频谱关系,结果表明:实测波高分布与瑞利分布基本相符;波浪谱的显著部分在狭窄的频率内,对应3、5、6、7级风谱峰密度值(m2s)分别为0.084 3、0.372 8、0.941 8、1.837 7,随着风速的增加,谱能也显著地增加,谱曲线下面的总面积增大,谱的显著部分波及的频率范围也扩大,且谱的显著部分沿低频率的方向推移,对应风浪的波高及周期范围增大. 展开更多
关键词 海浪 现场观测 海浪统计特征 频谱
在线阅读 下载PDF
浔江湾海水交换时间的计算 被引量:8
8
作者 许苏清 潘伟然 +1 位作者 张国荣 陈华 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2003年第5期629-632,共4页
采用箱式模型、二维ADI流体方程和对流 扩散模型计算浔江湾海水交换时间,分析浔江湾的水动力特点和湾内海水交换时空分布的特点.计算表明,浔江湾的海水平均交换时间为1.04d;湾内海水涨潮速度比落潮速度快;湾口的海水交换比湾内的其他地... 采用箱式模型、二维ADI流体方程和对流 扩散模型计算浔江湾海水交换时间,分析浔江湾的水动力特点和湾内海水交换时空分布的特点.计算表明,浔江湾的海水平均交换时间为1.04d;湾内海水涨潮速度比落潮速度快;湾口的海水交换比湾内的其他地方快. 展开更多
关键词 浔江湾 海水交换时间 水动力学 箱式模型 二维ADI流体方程 对流-扩散模型 水环境容量
在线阅读 下载PDF
厦门湾常风浪场数值模拟 被引量:11
9
作者 林毅辉 潘伟然 +1 位作者 张国荣 马腾 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2009年第2期298-301,共4页
通过采用目前国际上较为先进的第三代近岸海浪数值模型SWAN(Simulation waves nearshore),实现厦门湾常风浪场的数值模拟.先分别计算平潭、崇武、东山3个海洋观测站16个风向的多年平均风场,实现台湾海峡常风浪场的数值模拟,再将模拟结... 通过采用目前国际上较为先进的第三代近岸海浪数值模型SWAN(Simulation waves nearshore),实现厦门湾常风浪场的数值模拟.先分别计算平潭、崇武、东山3个海洋观测站16个风向的多年平均风场,实现台湾海峡常风浪场的数值模拟,再将模拟结果嵌套入厦门湾,根据实测风场,进行厦门湾常风浪场的数值模拟.模拟结果显示厦门湾常风浪场主要受两个因素的影响,一是风的作用,二是从台湾海峡传进港湾的波浪.厦门湾湾口的有效波高和有效周期均大于湾顶.在主风向风作用下,湾口区有效波高可达1.8 m,有效周期可达6 s,湾顶则约为0.2 m和2.1 s. 展开更多
关键词 SWAN 厦门湾 常风浪场 数值模拟
在线阅读 下载PDF
湄洲湾潮流特性的数值研究 被引量:9
10
作者 刘浩 潘伟然 +1 位作者 骆智斌 石晓枫 《台湾海峡》 CAS CSCD 北大核心 2009年第1期118-122,共5页
本文运用一个含动边界的二维河口海岸动力模型模拟了湄洲湾潮汐潮流的基本特征,模拟结果与实测数据吻合较好.在此基础上估算了湄洲湾大、小潮过程的纳潮量,并根据湾内保守示踪物的质量(浓度)随涨、落潮流周期性的变化,进一步估算了湄洲... 本文运用一个含动边界的二维河口海岸动力模型模拟了湄洲湾潮汐潮流的基本特征,模拟结果与实测数据吻合较好.在此基础上估算了湄洲湾大、小潮过程的纳潮量,并根据湾内保守示踪物的质量(浓度)随涨、落潮流周期性的变化,进一步估算了湄洲湾的水交换周期,其半交换和80%的交换周期分别为5d和15d.同时可以看出,主航道深水区的水交换特性明显强于湾顶浅水区. 展开更多
关键词 湄洲湾 模型 潮汐潮流 水交换
在线阅读 下载PDF
福清湾潮流场及污染物输运特性的模型研究 被引量:5
11
作者 刘浩 吴青松 +1 位作者 骆智斌 潘伟然 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2009年第1期139-143,共5页
运用一个含动边界的二维河口海岸动力模型模拟了福清湾潮汐潮流的基本特征,模拟结果与实测数据吻合较好.在此基础上估算了福清湾涨落潮流历时,分别为6 h 19 min 2 s和6 h 3 min 48 s;大小潮周期的纳潮量,分别为10.94×107和4.47... 运用一个含动边界的二维河口海岸动力模型模拟了福清湾潮汐潮流的基本特征,模拟结果与实测数据吻合较好.在此基础上估算了福清湾涨落潮流历时,分别为6 h 19 min 2 s和6 h 3 min 48 s;大小潮周期的纳潮量,分别为10.94×107和4.47×107m3;福清湾的半交换和80%的交换周期分别为3和6 d.最后,以COD为例,给出了在排放稳定后保守污染物的分布特征. 展开更多
关键词 福清湾 数值模型 潮汐潮流 污染物
在线阅读 下载PDF
台湾海峡1999年8月表层温盐的走航式观测 被引量:5
12
作者 陈华 胡建宇 +5 位作者 潘伟然 曾淦宁 陈照章 贺志刚 张彩云 李虹 《海洋通报》 CAS CSCD 北大核心 2001年第6期25-31,共7页
利用1999年8月台湾海峡的调查资料中走航式温盐仪的资料,对该海区的表层温盐分布进行分析?分析结果表明:调查期间在该海区存在多处上升流和冲淡水?台湾海峡调查海区的上升流区大致可分为两类:(1)海峡西部沿岸上升流区:平潭上升流区?东... 利用1999年8月台湾海峡的调查资料中走航式温盐仪的资料,对该海区的表层温盐分布进行分析?分析结果表明:调查期间在该海区存在多处上升流和冲淡水?台湾海峡调查海区的上升流区大致可分为两类:(1)海峡西部沿岸上升流区:平潭上升流区?东山上升流区和南澳上升流区;(2)台湾浅滩附近的上升流区?冲淡水主要有九龙江冲淡水?珠江冲淡水和闽江冲淡水,此外还有粤东沿岸水? 展开更多
关键词 台湾海峡 上升流 冲淡水 表层温盐分布 1998年
在线阅读 下载PDF
台湾海峡南部1999年8月温、盐、密度的分析 被引量:3
13
作者 曾淦宁 胡建宇 +5 位作者 潘伟然 陈照章 陈华 张彩云 贺志刚 李虹 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2001年第3期675-678,共4页
利用 1999年 8月份的 CTD资料 ,对台湾海峡南部温、盐、密度进行了分析 .分析结果表明 :1)粤东及闽南沿岸的上升流在调查期间较为明显 ,上升流的低温中心位于南澳岛以南区域 ;2 )台湾浅滩东侧在调查期间也存在上升流 ;3)在南澳与台湾浅... 利用 1999年 8月份的 CTD资料 ,对台湾海峡南部温、盐、密度进行了分析 .分析结果表明 :1)粤东及闽南沿岸的上升流在调查期间较为明显 ,上升流的低温中心位于南澳岛以南区域 ;2 )台湾浅滩东侧在调查期间也存在上升流 ;3)在南澳与台湾浅滩之间海域 ,存在一个来自西南方向的高温、低盐、低密水舌 ,显示此区域受来自粤东沿岸水的影响 .水平方向 ,粤东沿岸水与表层离岸的低温高盐水逐渐混合 ;垂直方向 ,粤东沿岸水影响范围可达 2 0~ 30 展开更多
关键词 台湾海峡南部 温度 盐度 密度 1999年 8月
在线阅读 下载PDF
泉州湾水体结构的潮周期变化 被引量:8
14
作者 刘浩 骆智斌 潘伟然 《台湾海峡》 CAS CSCD 北大核心 2009年第3期316-320,共5页
泉州湾6个站点的观测数据显示:内湾涨潮流历时由底层向表层逐渐变短,而落潮流历时则逐渐变长;同时,各个站点实测的盐度水深结构也具有明显的潮周期变化特征,这可能与局地水体的层化和混合机制的交替变化密切相关.进一步定量分析s2站位... 泉州湾6个站点的观测数据显示:内湾涨潮流历时由底层向表层逐渐变短,而落潮流历时则逐渐变长;同时,各个站点实测的盐度水深结构也具有明显的潮周期变化特征,这可能与局地水体的层化和混合机制的交替变化密切相关.进一步定量分析s2站位上影响水体结构变化的四种机制发现:外湾的河口环流作用与潮汐张力、风混合以及潮汐混合作用相比要小一个量级.在一般天气条件下的大潮周期,潮汐混合与潮汐张力相互竞争是导致水体结构交替变化的根本原因;小潮周期由于潮流混合作用减弱,水体层化现象得以持续较长时间,风的搅动在特定时刻可以起到削弱层化结构的作用. 展开更多
关键词 物理海洋学 层化和混合 潮汐张力 河口环流 泉州湾
在线阅读 下载PDF
基于三维调和分离的同安湾口门断面冬季潮流和余流特征的分析 被引量:5
15
作者 沈俊强 潘伟然 +2 位作者 张国荣 骆智斌 林毅辉 《台湾海峡》 CAS CSCD 2012年第3期338-344,共7页
利用同安湾口门断面走航ADCP观测数据,采用基于高斯基函数的Candela空间插值方法,对海流进行空间三维的调和分离,并对分离的各分潮流与余流进行流场的回归检验分析,F(α=0.01)检验显示流场的回归效果显著.观测期间同安湾口门断面平均潮... 利用同安湾口门断面走航ADCP观测数据,采用基于高斯基函数的Candela空间插值方法,对海流进行空间三维的调和分离,并对分离的各分潮流与余流进行流场的回归检验分析,F(α=0.01)检验显示流场的回归效果显著.观测期间同安湾口门断面平均潮差为5.15 m,最小潮差4.38m,最大潮差5.98 m,断面上涨落潮最大流速分别为92.3、80.3 cm/s,潮流特征分析表明,同安湾口门断面的潮流运动形式为往复流,以半日潮流为主,约占70%潮流信息,潮流流速从海表面向海底递减.余流最大值为12.5 cm/s,空间分布上将同安湾口门断面一分为二形成北进南出的余流进出通道,这与以往余流格局的认识相符. 展开更多
关键词 物理海洋学 同安湾 走航ADCP 空间插值 潮流分离 余流
在线阅读 下载PDF
营养盐负荷对浮游植物水华影响的模型研究 被引量:17
16
作者 刘浩 潘伟然 《水科学进展》 EI CAS CSCD 北大核心 2008年第3期345-351,共7页
运用一个耦合的生物物理模型模拟了渤海初级生产力,浮游植物生物量和氮磷含量的季节变化特征,模拟结果与实测值吻合较好。在模型验证的基础上,进一步考察了两种营养盐负荷——河流和沉积物对渤海生态系统动力过程的影响,发现河载营养盐... 运用一个耦合的生物物理模型模拟了渤海初级生产力,浮游植物生物量和氮磷含量的季节变化特征,模拟结果与实测值吻合较好。在模型验证的基础上,进一步考察了两种营养盐负荷——河流和沉积物对渤海生态系统动力过程的影响,发现河载营养盐对渤海生态体系的影响主要集中在河口水域,而限制沉积物中的营养盐进入水体则能显著抑制渤海范围内藻类水华的爆发。 展开更多
关键词 生态模型 浮游植物 营养盐负荷 水华 渤海
在线阅读 下载PDF
GA-BP神经网络在罗源湾口波浪模拟研究中的应用 被引量:4
17
作者 于小龙 潘伟然 +3 位作者 张国荣 崔文惠 林毅辉 骆智斌 《台湾海峡》 CAS CSCD 2012年第2期166-172,共7页
由于BP神经网络存在收敛速度慢和易于陷入极小值等缺点,引入遗传算法(GA)对网络的权值和阈值加以优化,并采用不同组合的输入因子和样本数,对福建省罗源湾口的波浪进行模拟研究.对输入因子的分析结果表明,研究区域的波浪主要受台湾海峡... 由于BP神经网络存在收敛速度慢和易于陷入极小值等缺点,引入遗传算法(GA)对网络的权值和阈值加以优化,并采用不同组合的输入因子和样本数,对福建省罗源湾口的波浪进行模拟研究.对输入因子的分析结果表明,研究区域的波浪主要受台湾海峡波浪传播影响,与局地气象因子(风速、风向、海气温差)的相关性较弱.训练样本试验表明,30 d以上的波浪历史数据可使GA-BP神经网络充分学习研究区域的波浪特征,从而实现对波浪要素的高精度模拟.模拟结果显示,对春、夏季实测波浪数据的模拟效果均很好,其中相关性分别为0.967和0.938,均方根误差分别为0.112 m和0.107 m,表明GA-BP神经网络在近岸波浪模拟预报中有较广阔的应用前景. 展开更多
关键词 物理海洋学 BP神经网络 遗传算法 波浪模拟 罗源湾
在线阅读 下载PDF
基于走航ADCP资料的厦门内湾东西口门海域潮流与潮致余流特征分析 被引量:3
18
作者 沈俊强 潘伟然 +2 位作者 张国荣 骆智斌 林毅辉 《应用海洋学学报》 CAS CSCD 2014年第3期395-403,共9页
应用厦门内湾东西3个口门断面(五通—下店、嵩屿—鼓浪屿、鼓浪屿—轮渡)一个潮周日走航ADCP资料的空间插值潮、流调和分离结果,分析占主导的M2分潮流与余流沿断面的空间分布特征.分析结果表明受地形约束的3个断面M2分潮流均以往复流为... 应用厦门内湾东西3个口门断面(五通—下店、嵩屿—鼓浪屿、鼓浪屿—轮渡)一个潮周日走航ADCP资料的空间插值潮、流调和分离结果,分析占主导的M2分潮流与余流沿断面的空间分布特征.分析结果表明受地形约束的3个断面M2分潮流均以往复流为主,断面深槽为主流通道,M2分潮最大可能流速均沿主流深槽向左右岸递减,三断面中以鼓浪屿—轮渡断面最为强劲,平均值为46 cm/s,五通—下店断面次之为43 cm/s,嵩屿—鼓浪屿断面最弱,仅为33 cm/s.该半日潮波进入厦门内湾首先在鼓浪屿—轮渡断面达到最大值,出现时间在月中天后2.1 h,其次为五通—下店断面2.7 h,最后为嵩屿—鼓浪屿断面3.3 h.再者地形变化剧烈的嵩屿—鼓浪屿断面M2分潮非线性运动较强,浅水分潮的幅度增大.潮波的非线性效应又衍生了非周期性的潮致余流,该余流在厦门内湾东西口门有较为稳定的存在,且表现为近岸区流幅大,潮流通道内小,嵩屿—鼓浪屿与鼓浪屿—轮渡断面均值为9 cm/s,五通—下店断面均值为5 cm/s,同时各个断面余流均存在使内湾海水正常更替,沿断面法向一分为二,一进一出的态势.基于本次观测结果计算出一潮周日厦门内湾东西口门余流通量为76.156×106m3,净潮通量为69.591×106m3,相对偏差小于8.6%. 展开更多
关键词 物理海洋学 厦门内湾 走航ADCP 空间插值 潮流分离 余流
在线阅读 下载PDF
厦门湾海上溢油风险的评估与分析 被引量:3
19
作者 王君 潘伟然 +1 位作者 张国荣 马腾 《台湾海峡》 CAS CSCD 北大核心 2009年第4期534-539,共6页
本文使用数理统计与概率的相关知识,并且根据厦门湾近几年进出船舶的统计资料,估算了未来几年中厦门湾海域船舶突发性溢油事故发生的概率.结果表明,厦门湾未来几年的溢油发生概率为0.27次/a.然后运用数值模拟方法对假设的突发溢油事故... 本文使用数理统计与概率的相关知识,并且根据厦门湾近几年进出船舶的统计资料,估算了未来几年中厦门湾海域船舶突发性溢油事故发生的概率.结果表明,厦门湾未来几年的溢油发生概率为0.27次/a.然后运用数值模拟方法对假设的突发溢油事故中溢油的扩展进行预测分析并结合建立的欧拉潮流场对油膜的漂移路径进行研究,分别从高潮、低潮、涨急、落急四个时刻开始计算溢油路径,结果显示油膜影响范围可从海沧南部延伸至浯村屿东部,根据模拟结果可对突发溢油的风险进行评估. 展开更多
关键词 物理海洋学 突发性溢油 概率估算 扩散 漂移
在线阅读 下载PDF
筼筜湖纳潮量与海水交换时间的计算 被引量:4
20
作者 姬厚德 潘伟然 +1 位作者 张国荣 骆智斌 《厦门大学学报(自然科学版)》 CAS CSCD 北大核心 2006年第5期660-663,共4页
通过在筼筜湖引水渠和出水渠断面设置流向、流速、水深等要素的测量断面,以获得筼筜湖纳潮情况下的水动力特征并由此计算筼筜湖的纳潮量和海水交换周期,在中潮条件下,筼筜湖的纳潮量为3.5×105m3,海水平均半交换周期为8.5 d.
关键词 海水交换时间 纳潮量 筼筜湖
在线阅读 下载PDF
上一页 1 2 下一页 到第
使用帮助 返回顶部