期刊文献+
共找到190篇文章
< 1 2 10 >
每页显示 20 50 100
非均匀介质中地震波走时与射线路径快速计算技术 被引量:27
1
作者 赵爱华 张中杰 +1 位作者 王光杰 王辉 《地震学报》 CSCD 北大核心 2000年第2期151-157,共7页
介绍了一种基于惠更斯原理和费马原理求取地震波走时及其反射波射线路径的新方法 .该方法具有原理简单、易于实现、能适应较为复杂地质模型以及易于将其推广到各向异性介质等优点 .为了克服基本算法速度较慢的缺陷 ,提出了一种地震波走... 介绍了一种基于惠更斯原理和费马原理求取地震波走时及其反射波射线路径的新方法 .该方法具有原理简单、易于实现、能适应较为复杂地质模型以及易于将其推广到各向异性介质等优点 .为了克服基本算法速度较慢的缺陷 ,提出了一种地震波走时和反射波射线路径计算的改进方法 .在保证精度的条件下 ,该改进算法的计算速度显著提高 . 展开更多
关键词 非均匀介质 射线路径 计算技术 地震 走时
在线阅读 下载PDF
三维非均匀介质中真振幅地震偏移算子研究 被引量:42
2
作者 崔兴福 张关泉 吴雅丽 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2004年第3期509-513,共5页
利用三维非均匀介质中的波动方程 ,进行振幅保真波场偏移算子分解 ,得到用于真振幅偏移的单程波方程 .经过数学推理 ,导出裂步Fourier法真振幅偏移和Fourier有限差分法真振幅偏移的算子方程 ,并给出其具体的实现过程 .
关键词 三维非均匀介质 单程波方程 真振幅偏移 裂步Fourier偏移算子 Fourier有限差分 地震反射信号
在线阅读 下载PDF
错格实数傅里叶变换微分算子及其在非均匀介质波动传播研究中的应用 被引量:26
3
作者 赵志新 徐纪人 堀内茂木 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2003年第2期234-240,共7页
提出一种新的数值微分运算方法 ,即错格实数傅里叶变换微分法 .该方法的运算速度比错格复数傅里叶变换数值微分解法快 0 .33倍 ;因为该微分算法在整个微分运算过程中保留了奈奎斯特分量 ,使得它比普通分格的实数傅里叶变换数值微分算法... 提出一种新的数值微分运算方法 ,即错格实数傅里叶变换微分法 .该方法的运算速度比错格复数傅里叶变换数值微分解法快 0 .33倍 ;因为该微分算法在整个微分运算过程中保留了奈奎斯特分量 ,使得它比普通分格的实数傅里叶变换数值微分算法的精度高 ,稳定性好 .将该方法和Cagniard DeHoop解析法在求解半无限空间地震波动的问题中进行比较 ,结果表明 ,新微分法的精度和解析方法的精度相同 .在非均匀介质中的地震波传播数值模拟的结果表明 ,该方法是一种研究非均匀介质中地震波传播问题的有效的数值微分方法 . 展开更多
关键词 错格分割法 错格实数傅里叶变换微分 数值模拟 非均匀介质 地震波传播
在线阅读 下载PDF
非均匀介质条件下矿震震波三维传播模型构建及其应用 被引量:8
4
作者 贾宝新 赵培 +2 位作者 姜明 李国臻 王帅 《煤炭学报》 EI CAS CSCD 北大核心 2014年第2期364-370,共7页
为了研究速度模型对矿震定位误差的影响,基于惠更斯原理建立了矿震波在层状非均匀介质中的三维传播模型,通过现场爆破振动试验,采用矿山微震定位监测系统进行监测。对矿震震波的正演建立了矿震震波传播的三维波动方程,并给出了矿震波在... 为了研究速度模型对矿震定位误差的影响,基于惠更斯原理建立了矿震波在层状非均匀介质中的三维传播模型,通过现场爆破振动试验,采用矿山微震定位监测系统进行监测。对矿震震波的正演建立了矿震震波传播的三维波动方程,并给出了矿震波在介质分界面上传播速度计算公式,由此对经典线性定位方法进行合理改进,使得基于均匀介质的线性定位方法适合于非均匀分层介质。将理论推导结果应用于矿山微震监测定位系统中,进行实际矿震定位计算。试验表明,速度模型对定位的影响较大,运用所建立的三维传播模型应用于矿震定位,有效提高了定位精度。 展开更多
关键词 矿震震波 非均匀介质 三维传播模型 定位
在线阅读 下载PDF
用格子玻尔兹曼方法模拟非均匀介质中的电场响应 被引量:6
5
作者 李学民 曹俊兴 +1 位作者 何晓燕 王兴建 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2004年第2期349-353,共5页
介绍了用格子玻尔兹曼方法模拟非均匀介质中的电场响应的数值模拟方法 .格子玻尔兹曼方法是从微观领域出发进行数值计算的一种全新的正演模拟方法 ;从玻尔兹曼碰撞模型出发 ,利用泰勒展开和Chapman Enskog展开 ,在基本力学守恒条件和约... 介绍了用格子玻尔兹曼方法模拟非均匀介质中的电场响应的数值模拟方法 .格子玻尔兹曼方法是从微观领域出发进行数值计算的一种全新的正演模拟方法 ;从玻尔兹曼碰撞模型出发 ,利用泰勒展开和Chapman Enskog展开 ,在基本力学守恒条件和约束条件的限制下 ,导出了电场响应的扩散方程 ,得到了局部平衡分布函数的表达式 ,给出了若干正演模拟的结果 ;其结果表明 ,利用这种方法进行非均匀介质中的电场响应正演模拟具有灵活、方便和简单等优点 . 展开更多
关键词 格子玻尔兹曼方法 非均匀介质 电场响应 数值模拟 数值计算
在线阅读 下载PDF
核素在非均匀介质中的迁移预测 被引量:6
6
作者 郭敏丽 王金生 +1 位作者 滕彦国 田浩 《中国环境科学》 EI CAS CSCD 北大核心 2009年第3期325-329,共5页
针对非均匀性明显的某放射性废物处置场候选场址,应用分形理论方法,在介质水力性质参数与介质结构分维之间建立函数关系式来描述场址介质的非均匀性,并选择放射性废物中的代表性核素3H和90Sr,应用不同方法对核素迁移进行预测比较.结果表... 针对非均匀性明显的某放射性废物处置场候选场址,应用分形理论方法,在介质水力性质参数与介质结构分维之间建立函数关系式来描述场址介质的非均匀性,并选择放射性废物中的代表性核素3H和90Sr,应用不同方法对核素迁移进行预测比较.结果表明,采用分形理论方法预测的核素迁移浓度比参数直接回归法更接近实际情况;在应用不同参数确定方法预测核素迁移浓度时,分配系数的取值对预测结果影响较大. 展开更多
关键词 放射性废物 核素迁移 非均匀介质 分形理论
在线阅读 下载PDF
复杂非均匀介质伪谱法波场数值模拟 被引量:10
7
作者 程冰洁 李小凡 徐天吉 《石油物探》 EI CSCD 2007年第1期16-19,共4页
在地震波场数值正演模拟方法的研究中,计算精度和计算效率是评价方法的有效性及优越性的2个关键问题。从—阶速度—应力弹性波动方程出发,利用伪谱法求解波动方程,对复杂非均匀介质模型中的波场进行了正演模拟,并利用经典的Mamousi速度... 在地震波场数值正演模拟方法的研究中,计算精度和计算效率是评价方法的有效性及优越性的2个关键问题。从—阶速度—应力弹性波动方程出发,利用伪谱法求解波动方程,对复杂非均匀介质模型中的波场进行了正演模拟,并利用经典的Mamousi速度模型验证了该方法所具有的优势及存在的问题。将伪谱法模拟结果与交错网格高阶有限差分法的模拟结果比较可知:对于较为简单的非均匀模型,伪谱法和交错网格高阶有限差分法生成了几乎相同的波场;而当模型非常复杂且存在变化较剧烈的速度间断面时,伪谱法的模拟结果比较差。尽管如此,伪谱法计算速度快,计算效率高,能够直观、高效地反映介质中波场的传播规律。因而仍不失为一种很好的地震波模拟方法。 展开更多
关键词 非均匀介质 伪谱法 数值模拟 波动方程 计算精度 计算效率
在线阅读 下载PDF
非均匀介质散射问题的体积分方程数值解法 被引量:9
8
作者 符力耘 杨慧珠 牟永光 《力学学报》 EI CSCD 北大核心 1998年第4期488-494,共7页
将非均匀介质视为某一均匀背景介质中的扰动,可建立用均匀背景介质格林函数作基本解的体积分方程.给出了配置法求解体积分方程的数值方法,首先解得扰动域内各点以速度扰动为权的波场函数,然后回代计算得到观测面上各接收点的散射波... 将非均匀介质视为某一均匀背景介质中的扰动,可建立用均匀背景介质格林函数作基本解的体积分方程.给出了配置法求解体积分方程的数值方法,首先解得扰动域内各点以速度扰动为权的波场函数,然后回代计算得到观测面上各接收点的散射波场.与边界元法和Born近似法计算结果比较表明该方法具有很高的精度,可得到穿过非均匀介质传播的全波场,包括所有的透射波、反射波、衍射波、首波及多重散射波.最后将体积分方程法用于非均匀介质散射问题的研究和地震模型计算. 展开更多
关键词 非均匀介质 散射问题 体积分方程法 地震模拟
在线阅读 下载PDF
倾斜分层非均匀介质中磁偶极子场的数值分析 被引量:5
9
作者 史晓锋 李铮 +1 位作者 王宝发 肖红兵 《北京航空航天大学学报》 EI CAS CSCD 北大核心 2003年第2期136-139,共4页
研究倾斜分层非均匀介质条件下磁偶极子源在介质空间的电磁响应 ,给出了场表达式和数值递推算法 ,提出了在大倾角条件下利用电阻率曲线划分界面和进行前向预测的方法 ,数值模拟和应用实践证明该方法是有效的 。
关键词 倾斜分层非均匀介质 磁偶极子场 数值分析 电磁波 石油测井 电阻率测量
在线阅读 下载PDF
非均匀介质弹性波动方程的不规则网格有限差分方法 被引量:7
10
作者 孙卫涛 杨慧珠 舒继武 《计算力学学报》 EI CAS CSCD 北大核心 2004年第2期135-142,共8页
从弹性波动方程出发,提出了一种新的空间不规则网格有限差分方法,并用于求解非均匀各向异性介质中的弹性波正演问题。这种方法简单易行,对于复杂几何结构,例如低速层、套管井和非平面界面等,在较细的不规则网格上进行离散,计算时间和占... 从弹性波动方程出发,提出了一种新的空间不规则网格有限差分方法,并用于求解非均匀各向异性介质中的弹性波正演问题。这种方法简单易行,对于复杂几何结构,例如低速层、套管井和非平面界面等,在较细的不规则网格上进行离散,计算时间和占用内存更少。与多重网格差分方法相比,该方法不需要粗、细网格之间的插值,所有网格差分计算在同一次空间迭代中完成。具有复杂几何交界面的模型计算,包括地下透镜体、套管井眼等,在确定弹性常数和密度后,用不规则网格的差分方法更易实现。该方法使用了Higdon吸收边界条件解决人工边界反射问题,引入了新的稳定性条件和网格频散条件,很好地消除了非物理散射波。理论模型的数值计算表明,该方法具有良好的稳定性和计算精度,在模拟非均匀介质弹性波传播时,比相同精度的规则网格有限差分方法计算速度更快。该方法易于推广到非结构网格和三维问题中。 展开更多
关键词 有限差分 不规则网格 非均匀介质 弹性波动方程 固体力学 网格离散化 弹性常数
在线阅读 下载PDF
非均匀介质热蠕变流动的数值求解 被引量:6
11
作者 刘洁 刘启元 宋惠珍 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2006年第4期1029-1036,共8页
针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格-粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场... 针对非均匀介质中热蠕变流动问题,给出了有限单元方法与网格-粒子方法联合求解新技术,即有限单元方法求解欧拉网格节点上的未知量,分布于单元内部作为物质成分标记的粒子反映变形过程.有限元法求解动量方程和连续性方程时引入了速度场和压力场等阶插值的压力场稳定的Petrov-Galerkin方法,求解能量方程时采用了流线迎风Petrov-Galerkin方法,网格-粒子算法中采用双线性插值与有限单元插值函数对应.有限单元计算与网格-粒子计算相对独立,两种方法计算的数据通过有限单元节点传递.同时,实现了三角形单元的算法和程序,解决了复杂结构条件下不规则网格计算的问题.通过经典方腔热对流问题验证了程序,给出了不规则形态块体沉降算例,并分析了数值解的稳定性. 展开更多
关键词 非均匀介质 热蠕变流动 网格-粒子 有限单元方法
在线阅读 下载PDF
合成三维横向非均匀介质远震体波接收函数的Maslov方法 被引量:11
12
作者 陈九辉 刘启元 《地球物理学报》 SCIE EI CAS CSCD 北大核心 1999年第1期84-93,共10页
三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要.本文发展了基于Maslov理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪... 三维横向非均匀介质远震体波接收函数的合成算法对于三维介质远震体波接收函数研究至关重要.本文发展了基于Maslov理论的横向非均匀介质接收函数合成算法,针对远震体波接收函数计算的特点,利用遗传算法完成三维动力学射线追踪,并采用递归算法组织射线的追踪过程,从而在保证追踪精度的同时提高了射线追踪的效率,并可完全避免传统射线追踪过程中的射线编码问题.正演计算结果表明,此法是一种高效、实用的三维横向非均匀介质接收函数计算方法. 展开更多
关键词 Maslov法 接收函数 非均匀介质 远震波 地震勘探
在线阅读 下载PDF
非均匀介质的光学常数 被引量:8
13
作者 徐文兰 褚君浩 《红外与毫米波学报》 SCIE EI CAS CSCD 北大核心 1994年第5期391-395,共5页
提出了一种计算非均匀介质光学常数的新方法,该方法仔细地考虑了光在非均匀介质中的吸收和散射.研究了非均匀介质层的厚度、介质中颗粒粒径、衬底反射性能等对介质光学常数的影响.这种研究将有利于同有关的测量实验结合和比较.
关键词 非均匀介质 光学常数 光吸收 散射 厚度 粒径
在线阅读 下载PDF
垂向非均匀介质中首波特征分析 被引量:6
14
作者 张坚 张海明 陈晓非 《地震学报》 CSCD 北大核心 2002年第6期559-568,共10页
利用合成理论地震图方法研究了存在高速层或低速层的地壳模型及壳-幔过渡带模型中首波等震相的动力学特征,指出首波的特征对于高速层的结构变化比较敏感.当高速层厚度与特征波长相比较小时,地震波的衍射现象明显,这种情况下高速层不能... 利用合成理论地震图方法研究了存在高速层或低速层的地壳模型及壳-幔过渡带模型中首波等震相的动力学特征,指出首波的特征对于高速层的结构变化比较敏感.当高速层厚度与特征波长相比较小时,地震波的衍射现象明显,这种情况下高速层不能屏蔽在其下面一层的上界面传播的首波,且该首波震相的强度随高速层厚度或速度的增加而递减;当高速层厚度与特征波长相当时,高速层底面的反射波震相与首波震相到时接近,会因互相干涉而减弱;对于低速层,首波震相强度较弱且随低速层速度的减小而递减;在壳-幔间断处引入一定厚度和速度变化范围的过渡带,可以得到更加明显的视首波震相,强度随过渡带厚度或速度变化范围的增加而增大. 展开更多
关键词 垂向非均匀介质 理论地震图 首波 高速层 低速层 壳-幔过渡带
在线阅读 下载PDF
非均匀介质地震AVO响应模拟及分析 被引量:3
15
作者 刘财 裴思嘉 +3 位作者 郭智奇 符伟 张宇生 刘喜武 《吉林大学学报(地球科学版)》 EI CAS CSCD 北大核心 2018年第5期1512-1521,共10页
地震AVO(amplitude versus offset)技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率... 地震AVO(amplitude versus offset)技术是一项利用振幅信息研究岩性、检测油气的地震勘探技术。常规方法基于Zoeppritz方程计算模型界面处的反射系数,而实际地下非均匀介质中地震反射特征不仅与入射角度、物性差异有关,还与入射波频率、地层厚度、薄互层结构等因素有关。为此,应用传播矩阵理论充分考虑与这些因素有关的调谐干涉等传播效应,针对实际数据计算非均匀地下介质的高精度合成地震记录,对比Zoeppritz方程、Shuey二项近似方程、Shuey三项近似方程以及传播矩阵算法的模拟效果。研究发现:在小角度入射时Zoeppritz方程、Shuey二项近似方程和Shuey三项近似方程的反射振幅和波形基本一致,大角度入射时Zoeppritz方程与Shuey三项近似方程接近;Zoeppritz算法的模拟结果在小角度入射和浅层情况下与传播矩阵算法差别较小,而在大角度入射和深层情况下与传播矩阵差别较大,说明层间多次波的调谐干涉等传播效应不可忽略。 展开更多
关键词 ZOEPPRITZ方程 Shuey方程 传播矩阵理论 非均匀介质
在线阅读 下载PDF
高斯束射线法在二维非均匀介质复杂构造中的应用 被引量:11
16
作者 吴立明 许云 乌达巴拉 《地球物理学报》 SCIE EI CSCD 北大核心 1995年第A01期144-152,共9页
高斯束射线法是地震勘探中一种非常有效的数值模拟方法,是一般射线方法的概括和推广.它克服了一般射线方法在奇异性区域(如焦散区,临界区和阴影区等)振幅趋于无穷以及对介质近似和模型小细节变化敏感的缺点,高斯束射线法在奇异性... 高斯束射线法是地震勘探中一种非常有效的数值模拟方法,是一般射线方法的概括和推广.它克服了一般射线方法在奇异性区域(如焦散区,临界区和阴影区等)振幅趋于无穷以及对介质近似和模型小细节变化敏感的缺点,高斯束射线法在奇异性区域比一般射线方法得到的结果要好得多.高斯束射线法是沿集中于中心射线邻近的正则性区域求解波动方程的高频近似解,该解主分量的振幅分布沿垂直于中心射线方向是高斯分布,所以称为高斯束射线法.高斯束射线法能对二维甚至三维的横向非均匀介质复杂构造计算合成地震记录,而且从计算的观点看,它避开了耗时的两点射线追踪,所需内存少,计算速度快. 展开更多
关键词 地震勘探 高斯束射线法 非均匀介质 地震波
在线阅读 下载PDF
伪谱法模拟各向异性/非均匀介质中弹性波的传播 被引量:7
17
作者 张光莹 曾新吾 《国防科技大学学报》 EI CAS CSCD 北大核心 2002年第3期18-22,共5页
数值模拟一直是研究弹性波在复杂介质中传播问题的有效手段之一。各向异性问题和非均匀问题日益引起关注。伪谱法是近年来发展较快、建立在弹性波动方程基础上的一种数值模拟方法。它的主要特征是将波动方程中的空间微分变换成频域中的... 数值模拟一直是研究弹性波在复杂介质中传播问题的有效手段之一。各向异性问题和非均匀问题日益引起关注。伪谱法是近年来发展较快、建立在弹性波动方程基础上的一种数值模拟方法。它的主要特征是将波动方程中的空间微分变换成频域中的乘法运算 ,从而只需要较少的空间格点就可以得到较高的计算精度。介绍了伪谱法及其稳定性条件 ;编制了一个适合于小型工作站或者PC机的模拟地震波在二维各向异性均匀 /非均匀介质中传播的数值模拟程序 ,模拟了各向异性 /非均匀介质中弹性波传播理论 。 展开更多
关键词 弹性波 传播 数值模拟 伪谱法 稳定性 各向异性 非均匀介质 复杂介质 地震波
在线阅读 下载PDF
非均匀介质地震波传播交错网格高阶有限差分法模拟 被引量:26
18
作者 裴正林 牟永光 《石油大学学报(自然科学版)》 CSCD 北大核心 2003年第6期17-21,共5页
采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精... 采用常规的二阶声波方程有限差分方法对于非均匀介质进行了数值模拟时 ,其数值模拟精度较低。而采用一阶双曲型标量波动方程 ,则无须对介质的弹性常数进行空间求导。根据Taylor级数展开式 ,推导出了交错网格一阶空间导数的任意偶数阶精度展开式和相应差分系数计算式以及一阶双曲型标量波动方程交错网格任意偶数阶精度差分格式 ,并给出了该差分算法的稳定性条件。用该差分算法对均匀介质模型、非均匀介质模型和Marmousi模型进行了数值模拟试验 ,并与伪谱法进行了对比。结果表明 ,一阶双曲型标量波动方程交错网格高阶差分法的模拟精度与伪谱法的精度非常接近 ,计算效率高 ,且适合于模拟非均匀介质。 展开更多
关键词 非均匀介质 地震波 交错网格 高阶有限差分法 数值模拟 一阶双曲型标量波动方程 油气勘探 地震勘探
在线阅读 下载PDF
垂向非均匀介质中转换波转换点计算 被引量:6
19
作者 陈雨红 魏修成 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2007年第4期1225-1231,共7页
转换点位置的计算是转换波资料处理中的一个关键问题.本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法.研究了通过速度... 转换点位置的计算是转换波资料处理中的一个关键问题.本文提出了分别基于速度随深度线性变化、速度随垂直走时线性变化、慢度随深度线性变化和慢度随垂直走时线性变化四种等效垂向非均匀介质情况下转换点位置的计算方法.研究了通过速度拟合、走时近似和相似系数谱三种方式选择合适的等效速度方法.结合理论模型对非均匀介质转换点计算方法、渐进转换点计算方法、Thomsen近似公式和均匀介质解析计算方法的误差进行了分析,结果表明非均匀介质转换点计算方法能更准确地计算转换点位置. 展开更多
关键词 转换点 转换波 非均匀介质 走时近似 相似系数
在线阅读 下载PDF
横向非均匀介质远震体波接收函数的波场特征 被引量:8
20
作者 陈九辉 刘启元 《地震学报》 CSCD 北大核心 2000年第6期614-621,共8页
系统研究了台站下方存在倾斜界面时 ,远震体波接收函数的波场特征及多方位接收函数对台站下方介质横向非均匀特征的反映能力 .我们利用合成三维横向非均匀介质接收函数的 Maslov方法 ,具体模拟了台站下方存在倾斜界面时的接收函数响应 ... 系统研究了台站下方存在倾斜界面时 ,远震体波接收函数的波场特征及多方位接收函数对台站下方介质横向非均匀特征的反映能力 .我们利用合成三维横向非均匀介质接收函数的 Maslov方法 ,具体模拟了台站下方存在倾斜界面时的接收函数响应 .结果表明 ,当台站下方存在倾向一致的倾斜界面时 ,远震体波接收函数的径向分量和切向分量分别随震源方位角的变化呈现对称和反对称的规律性变化 .利用不同震源方位角接收函数径向分量和切向分量的变化规律 ,可以估计台站下方界面的倾向和倾角 .当台站下方各界面的倾向不同时 ,随方位角的变化 ,接收函数只能直观给出界面的整体倾向 .实际观测数据的分析结果表明 ,对于利用单个台站接收函数研究台站下方介质的横向非均匀特征来说 ,简单倾斜界面可以看作较好的一级近似 . 展开更多
关键词 波场特征 接收函数 横向非均匀介质 远震体波 Maslov法
在线阅读 下载PDF
上一页 1 2 10 下一页 到第
使用帮助 返回顶部