期刊文献+
共找到18篇文章
< 1 >
每页显示 20 50 100
SegNet-based first-break picking via seismic waveform classification directly from shot gathers with sparsely distributed traces 被引量:1
1
作者 San-Yi Yuan Yue Zhao +2 位作者 Tao Xie Jie Qi Shang-Xu Wang 《Petroleum Science》 SCIE CAS CSCD 2022年第1期162-179,共18页
Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to a... Manually picking regularly and densely distributed first breaks(FBs)are critical for shallow velocitymodel building in seismic data processing.However,it is time consuming.We employ the fullyconvolutional Seg Net to address this issue and present a fast automatic seismic waveform classification method to pick densely-sampled FBs directly from common-shot gathers with sparsely distributed traces.Through feeding a large number of representative shot gathers with missing traces and the corresponding binary labels segmented by manually interpreted fully-sampled FBs,we can obtain a welltrained Seg Net model.When any unseen gather including the one with irregular trace spacing is inputted,the Seg Net can output the probability distribution of different categories for waveform classification.Then FBs can be picked by locating the boundaries between one class on post-FBs data and the other on pre-FBs background.Two land datasets with each over 2000 shots are adopted to illustrate that one well-trained 25-layer Seg Net can favorably classify waveform and further pick fully-sampled FBs verified by the manually-derived ones,even when the proportion of randomly missing traces reaches50%,21 traces are missing consecutively,or traces are missing regularly. 展开更多
关键词 First-break picking Deep learning irregular seismic data Waveform classification
下载PDF
南海北部中生界地震成像制约因素分析
2
作者 曾婷 郑金云 +4 位作者 张青林 李洪博 陶文芳 李杰 马艺璇 《石油物探》 CSCD 北大核心 2024年第2期318-335,共18页
南海北部中生界油气勘探积累了多个年度、多次采集、多种处理方案的地震资料,但地震成像品质总体上仍难以满足目前中生界油气勘探研究需求,对制约地震成像的影响因素认识不一,缺乏系统对比论证。通过对地质背景及沉积特征的分析,以现有... 南海北部中生界油气勘探积累了多个年度、多次采集、多种处理方案的地震资料,但地震成像品质总体上仍难以满足目前中生界油气勘探研究需求,对制约地震成像的影响因素认识不一,缺乏系统对比论证。通过对地质背景及沉积特征的分析,以现有二维地震资料为基础,建立典型地质结构模型并正演模拟,对制约中生界地震成像品质的地质、地球物理因素进行剖析对比和论证,从定性和定量的角度明确了不同地质条件下造成现有地震资料品质差的主控因素。研究表明,制约中生界地震成像的主要因素是海底及中新生界之间的两个强反射界面、中生界地层内部波阻抗差异小所导致的低信噪比;绕射波、多次波等各种干扰波异常发育使有效信号难以识别和成像;崎岖海底造成照明能量散射等。南海北部中生界制约地震成像因素的明确为后续中生界地震数据采集、处理提供了有效指导及技术支持,并进一步提出通过适当加大震源容量、提高覆盖次数、深沉放枪缆、高维度,获得更多深层有效反射波能量的采集建议;通过自适应鬼波压制、多域提高信噪比处理及基于频率谐振能量补偿等技术解决弱反射问题,宽频组合解决多次波发育问题,基于地质、折射波增强速度规律认识,沿层高密度速度拾取提高速度分析精度等处理建议。 展开更多
关键词 南海北部中生界 地震资料品质 强反射界面 波阻抗差异 崎岖海底 地震照明能量
下载PDF
一种基于双调和样条插值的完全非规则地震数据重建 被引量:1
3
作者 许昌昊 曹静杰 +1 位作者 高康富 蔡志成 《河北地质大学学报》 2023年第2期44-51,共8页
在地震勘探过程中,由于地形、环境等影响,采集到的地震数据常常是不规则且不完整的,这会对后续地震数据的处理造成不良影响,因此地震数据重建是一个常用的处理手段。目前已经发展了多种地震数据重建方法,大部分的地震数据的重建工作都... 在地震勘探过程中,由于地形、环境等影响,采集到的地震数据常常是不规则且不完整的,这会对后续地震数据的处理造成不良影响,因此地震数据重建是一个常用的处理手段。目前已经发展了多种地震数据重建方法,大部分的地震数据的重建工作都是建立在一个规则网格上进行,而对于完全非规则地震数据的重建研究较少。本文基于双调和样条插值的方法提出了一种非规则地震数据重建方法,通过求取一个能够通过n个点的双调和样条函数,将残缺的且非规则的地震数据插值到规则的网格中,从而得到完整的地震数据。本文进一步提出了分块插值的思想,在插值之前将每一时间或频率切片进行分块,对每一块分别进行插值。此方法改善了插值对内存需求过大的弊端,并极大的提高了计算效率。通过合成地震数据进行测试,实验证明了提出的双调和样条插值方法的插值效果。 展开更多
关键词 非规则地震数据 插值 双调和样条 重建
下载PDF
不规则地震道数据规则化重建方法研究 被引量:17
4
作者 高建军 陈小宏 +1 位作者 王芳芳 马剑 《地球物理学进展》 CSCD 北大核心 2011年第3期983-991,共9页
不规则地震数据会对地震多道处理技术的正确运行产生不良影响,降低地震资料处理质量.本文依据不规则地震数据的表现特征将其划分为四种类型并针对第三类不规则地震道数据采用抗泄露Fourier变换方法进行规则化重建.不规则采样数据会破坏F... 不规则地震数据会对地震多道处理技术的正确运行产生不良影响,降低地震资料处理质量.本文依据不规则地震数据的表现特征将其划分为四种类型并针对第三类不规则地震道数据采用抗泄露Fourier变换方法进行规则化重建.不规则采样数据会破坏Fourier基函数的正交性并产生频谱泄漏现象.抗泄露Fourier变换方法通过递归相减来压制谱泄露,实现Fourier基函数的重正交,进而估计出相对精确的Fourier系数,达到规则化重建的目的.抗泄露Fourier变换方法具有一定的抗噪声和抗假频重建能力,这使它优于传统的非均匀Fourier变换重建方法.理论模型和实际资料处理表明用抗泄露Fourier变换重建方法来解决第三类不规则地震道数据的规则化重建问题是有效可行的. 展开更多
关键词 不规则地震数据 地震数据规则化 地震数据重建 抗泄露Fourier变换
下载PDF
一种边缘保持的地震数据插值方法 被引量:9
5
作者 陆艳洪 陆文凯 翟正军 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2012年第3期991-997,共7页
在地震数据处理中,地震数据插值方法常常用来解决地震数据空间采样率低和不规则的问题.本文提出了一种基于边缘保持滤波器的地震数据插值方法.在该方法中,对于一个1D信号,逐点滑动一个处理窗口,将信号分成多个信号片段.对于某一个待恢... 在地震数据处理中,地震数据插值方法常常用来解决地震数据空间采样率低和不规则的问题.本文提出了一种基于边缘保持滤波器的地震数据插值方法.在该方法中,对于一个1D信号,逐点滑动一个处理窗口,将信号分成多个信号片段.对于某一个待恢复的缺失采样点,存在多个包含(或邻近)该采样点的信号片段可以用来估计它.采用多项式来拟合这些信号片段,并选择拟合误差最小的信号片段估计此缺失采样点,达到边缘保持的目的.对于2D地震信号,先沿不同方向扫描抽取1D信号,然后采用上述1D边缘保持插值算法分别进行处理,得到沿不同方向的插值结果.对于任一待插值采样点,选取对应拟合误差最小的方向的插值结果作为最后输出的2D数据的插值结果.理论模型和实际资料的处理结果表明,所提方法具有保边缘、抗假频及能够进行不规则数据重建等特点,既能有效的实现不规则地震数据的重建,又能很好的保持原有数据的边缘特征. 展开更多
关键词 插值 不规则地震数据 边缘保持 最小拟合误差 抗假频
下载PDF
不规则地震数据的抗假频重建方法 被引量:6
6
作者 高建军 陈小宏 +2 位作者 李景叶 张南南 刘志鹏 《石油地球物理勘探》 EI CSCD 北大核心 2010年第3期326-331,共6页
不规则采样地震数据会对地震数据多道处理的效果造成严重影响。本文将单步预测滤波拓展到多步预测滤波,基于多步自回归预测滤波方法,实现对不规则地震数据的分步、抗假频重建方法。文中首先运用带限Fourier重建方法对不规则地震数据进... 不规则采样地震数据会对地震数据多道处理的效果造成严重影响。本文将单步预测滤波拓展到多步预测滤波,基于多步自回归预测滤波方法,实现对不规则地震数据的分步、抗假频重建方法。文中首先运用带限Fourier重建方法对不规则地震数据进行规则化重建,得到无假频低频数据的频率谱;接着使用多步自回归算子从重建数据的无假频低频谱中提取整个有效频带的预测滤波因子;然后利用已知道数据和预测滤波因子重建得到完整数据的频率谱,最终实现不规则地震数据的抗假频重建。理论模型和实际资料的处理结果验证了本文方法的有效性。 展开更多
关键词 不规则地震数据 带限Fourier重建 多步自回归算法 地震数据重建
下载PDF
三维不规则地震数据重建方法 被引量:12
7
作者 高建军 陈小宏 李景叶 《石油地球物理勘探》 EI CSCD 北大核心 2011年第1期40-47,164+169,共8页
不规则采样地震数据会对多道处理技术产生不良影响,降低地震资料的处理质量。本文针对有限带宽三维不规则地震数据,将二维空间非均匀Fourier变换理论和贝叶斯参数反演方法相结合,进行反演重建。首先,采用分频重建策略,对每一个时间频率... 不规则采样地震数据会对多道处理技术产生不良影响,降低地震资料的处理质量。本文针对有限带宽三维不规则地震数据,将二维空间非均匀Fourier变换理论和贝叶斯参数反演方法相结合,进行反演重建。首先,采用分频重建策略,对每一个时间频率依据最小视速度确定出待重建数据的空间频率带宽,从不规则地震数据中估计出重建数据的空间Fourier谱。然后,将不规则地震数据重建视为谱重建的地球物理反演问题,运用贝叶斯参数反演理论估计出空间Fourier谱。在反演求解时,采用Delaunay三角网格剖分方法来确定不规则采样点的权值。此外,为避免复数矩阵求逆,使用预条件共轭梯度算法来保证求解的稳定性和收敛速度。理论模型和实际资料处理结果验证了本文方法的有效性和实用性。 展开更多
关键词 不规则地震数据 非均匀傅里叶变换 地震数据重建 贝叶斯反演
下载PDF
基于POCS方法指数阈值模型的不规则地震数据重建(英文) 被引量:16
8
作者 高建军 陈小宏 +2 位作者 李景叶 刘国昌 马剑 《Applied Geophysics》 SCIE CSCD 2010年第3期229-238,292,293,共12页
不规则地震数据会对地震多道处理技术的正确运行造成不良影响,降低地震资料的处理质量。本文将广泛用于图形图像重建的凸集投影方法应用到地震数据重建领域,实现规则样不规则道缺失数据的插值重建。对于整道缺失地震数据,将POCS迭代重... 不规则地震数据会对地震多道处理技术的正确运行造成不良影响,降低地震资料的处理质量。本文将广泛用于图形图像重建的凸集投影方法应用到地震数据重建领域,实现规则样不规则道缺失数据的插值重建。对于整道缺失地震数据,将POCS迭代重建过程由时间域转移到频率域实现,避免每次迭代都对时间做正反Fourier变换,节约了计算量。在迭代过程中,阈值参数的选择方式对重建效率有重要影响。本文设计了两种阈值集合模型进行重建试验,试验结果表明:在相同重建效果下,指数型阈值集合模型可以有效减少迭代次数,提高重建效率。此外,分析了POCS重建方法的抗噪性能和抗假频性能。最后,理论模型和实际资料处理效果验证了本文重建方法的正确性和有效性。 展开更多
关键词 不规则地震数据 地震数据重建 凸集投影 阈值模型
下载PDF
非规则测区三维地震偏移的边界吸收层法 被引量:6
9
作者 王玉学 张关泉 周瑞芬 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2002年第5期700-706,共7页
对不规则测区的地震资料直接进行三维偏移不仅可减少计算量 ,而且可提高偏移剖面的质量 .本文在单程波方程中增加了一项简单的阻尼因子 ,其作用如同在实际不规则测区的外部有一个薄吸收层 .因为波穿越这个薄层时 ,波场值迅速衰减 ,所以... 对不规则测区的地震资料直接进行三维偏移不仅可减少计算量 ,而且可提高偏移剖面的质量 .本文在单程波方程中增加了一项简单的阻尼因子 ,其作用如同在实际不规则测区的外部有一个薄吸收层 .因为波穿越这个薄层时 ,波场值迅速衰减 ,所以偏移时把吸收层外边界处的波场值取为零所产生的反射非常微弱 .又因为不规则测区零边界条件偏移程序极易编制 ,所以就可以实现对不规则测区的地震资料直接进行三维偏移 .通过对实际地震资料的处理 。 展开更多
关键词 边界吸收层法 不规则测区 偏移 吸收层 地震资料
下载PDF
OVT域地震数据规则化技术及应用 被引量:18
10
作者 李博 《石油物探》 EI CSCD 北大核心 2019年第1期53-62,共10页
宽方位矢量偏移距(OVT)地震处理方法对于复杂地质体的成像照明有良好的效果,但OVT域数据很难实现理想的规则采集,因此其关键处理技术之一就是规则化技术。通过分析OVT域的数据排列特征,提出了一种基于非规则傅里叶变换的最小平方反演数... 宽方位矢量偏移距(OVT)地震处理方法对于复杂地质体的成像照明有良好的效果,但OVT域数据很难实现理想的规则采集,因此其关键处理技术之一就是规则化技术。通过分析OVT域的数据排列特征,提出了一种基于非规则傅里叶变换的最小平方反演数据重建技术及实施流程,引入迭代非规则傅里叶变换加权范数正则化约束,使得重建结果在有限频宽内保持数据信号不受损失,同时在能量谱约束下,利用由低频信号估计的加权函数压制高频信号的假频问题。经过模型资料抽稀后插值的数据与原始数据残差对比验证了反演精度和可行性;将该技术应用于实际资料的联片处理,结果表明该技术可消除非规则采集对OVT域处理的影响,在保护微幅构造和断层成像方面有良好的效果,证实了方法的有效性和实用性。 展开更多
关键词 偏移距矢量面元 地震数据规则化 非规则采集 非规则傅里叶变换 最小平方反演
下载PDF
基于CRP不规则二维地震采集方法研究(英文) 被引量:1
11
作者 赵虎 尹成 +2 位作者 何光明 陈爱萍 敬龙江 《Applied Geophysics》 SCIE CSCD 2015年第1期73-78,122,共7页
中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集... 中国西部复杂山地是地震勘探的难点区域,地表和地下都复杂,从而导致地表激发困难,地震资料干扰严重、信噪比低。地下构造形态复杂,导致反射点分散严重,前人尝试过弯线采集,但很少考虑地下实际构造特征。本文章提出基于CRP的不规则采集设计方法,使地表避开难以激发的区域,地下考虑构造特征及CRP点分布位置,优化地表接收线位置,通过地表有目的布置接收点位置,最大程度的保证地下CRP点不分散,从而提高地震数据的信噪比.并以四川盆地某工地的实际资料通过模型数据验证了方法的合理性与可行性。表明该方法能够解决复杂地区地震采集存在的部分难题,为复杂地区的地震勘探采集提供一个新的思路。 展开更多
关键词 观测系统设计 不规则观测系统 采集设计 地震采集 CRP分析
下载PDF
加吸收层的三维45°上行波方程的隐式差分格式 被引量:1
12
作者 冯凤萍 周瑞芬 刘畅 《大庆石油学院学报》 CAS 北大核心 2005年第4期98-100,共3页
给出了附加吸收层的三维45°上行波方程,构造了该方程的一种隐式差分格式.利用该方程可实现对不规则测区直接进行三维偏移.实际地震资料处理的结果说明了该方法的有效性.
关键词 不规则测区 偏移 吸收层 地震资料
下载PDF
曲波域不规则地震数据提高分辨率技术研究
13
作者 孔雪 王德营 +2 位作者 张瑞香 朱学娟 胡秋媛 《CT理论与应用研究(中英文)》 2017年第6期707-713,共7页
受野外观测条件的限制,采集的地震数据体通常不规则,并缺失一部分数据道。传统的单道提高分辨率方法无法兼顾横向地震信息,处理结果存在空间一致性问题。为此,本文提出在曲波域内进行不规则地震数据,通过曲波变换实现对地震数据的稀疏表... 受野外观测条件的限制,采集的地震数据体通常不规则,并缺失一部分数据道。传统的单道提高分辨率方法无法兼顾横向地震信息,处理结果存在空间一致性问题。为此,本文提出在曲波域内进行不规则地震数据,通过曲波变换实现对地震数据的稀疏表征,将提高分辨率问题转化为曲波域1-范数约束的稀疏促进求解,得到规则化的高分辨率地震数据体。该方法避免传统单道提高分辨率方法存在的局限性,在提高分辨率的同时,能够恢复缺失的地震数据、压制随机噪声,进而提高地震数据的完备性,模型和实际资料试算,验证了该方法的正确性、有效性和适用性。 展开更多
关键词 提高分辨率 曲波域 不规则地震数据 地震数据规则化
下载PDF
三维FKK滤波技术在叠前去噪中的应用 被引量:4
14
作者 陈习峰 薛永安 +1 位作者 俞华 管文华 《复杂油气藏》 2013年第4期34-38,共5页
针对陆上三维地震数据普遍存在的采集不均匀、空间采样率不足及振幅非一致性,不能满足FKK滤波技术前提条件的问题,对现有的FKK滤波配套技术方法进行认真分析,指出了其存在的不足,并给出空间假频门槛分析法,最后提出了三维地震数据规则... 针对陆上三维地震数据普遍存在的采集不均匀、空间采样率不足及振幅非一致性,不能满足FKK滤波技术前提条件的问题,对现有的FKK滤波配套技术方法进行认真分析,指出了其存在的不足,并给出空间假频门槛分析法,最后提出了三维地震数据规则化及能量一致性处理的FKK滤波配套处理方法,实际应用表明该配套处理方法是有效的。 展开更多
关键词 地震数据处理 滤波噪音 空间假频 空间不规则采样 门槛分析
下载PDF
基于压缩傅里叶算法的五维插值技术在改善振幅一致性中的应用 被引量:5
15
作者 齐鹏 胡玮 +2 位作者 管文华 汤国松 刘立民 《石油物探》 EI CSCD 北大核心 2018年第2期213-221,共9页
地震数据采集时受地表因素影响常常变观,使得炮检点分布不规则,影响了振幅特征的一致性,不利于后续的叠前属性分析及储层预测。为此提出了基于压缩傅里叶插值算法的数据规则化技术,在插值位置对不同空间带宽波数进行扫描,优选出满足均... 地震数据采集时受地表因素影响常常变观,使得炮检点分布不规则,影响了振幅特征的一致性,不利于后续的叠前属性分析及储层预测。为此提出了基于压缩傅里叶插值算法的数据规则化技术,在插值位置对不同空间带宽波数进行扫描,优选出满足均方根误差条件的最大空间波数,用于计算插值权系数。为了提高插值结果的振幅一致性,插值时考虑了纵向、横向、炮检距、方位角和时间等五个维度的信息,并增加了局部倾角校正处理以提高插值的抗假频能力。利用该技术对某工区不规则地震数据进行了五维插值规则化处理,处理后的地震数据不仅具有规则的炮检点分布,而且提高了覆盖次数的均匀性。对比分析了五维插值前后地震数据的振幅特征,结果表明,基于压缩傅里叶算法的五维插值技术可以很好地降低不规则炮检点分布对振幅特征的影响,提高了振幅的一致性。 展开更多
关键词 五维插值技术 不规则观测系统 数据规则化 压缩傅里叶插值算法 振幅一致性
下载PDF
基于压缩感知的非规则地震勘探观测系统设计与数据重建 被引量:15
16
作者 周松 吕尧 +3 位作者 吕公河 舒国旭 石太昆 霍守东 《石油物探》 EI CSCD 北大核心 2017年第5期617-625,共9页
提出了基于压缩感知的非规则观测系统优化设计和地震数据重建方法。基于贪心序贯策略,利用逐个选取采样点位置的方法来降低感知矩阵的最大互相关值,从而完成非规则观测系统的优化;利用L_0正则化和L_1正则化混合迭代的方法求解欠定采样矩... 提出了基于压缩感知的非规则观测系统优化设计和地震数据重建方法。基于贪心序贯策略,利用逐个选取采样点位置的方法来降低感知矩阵的最大互相关值,从而完成非规则观测系统的优化;利用L_0正则化和L_1正则化混合迭代的方法求解欠定采样矩阵,重建地震记录,获得良好的效果。进行了TFT工区实际观测系统的设计,鉴于接收电缆长度的限制,添加最大道间距和最小道间距作为观测系统设计过程中的约束,得到了符合实际勘探要求的非规则观测系统的优化设计。该方法为三维高密度地震勘探的大规模开展提供了有效的方案。 展开更多
关键词 压缩感知 稀疏 非规则 观测系统设计 数据重建
下载PDF
3D高阶抛物Radon变换在不规则地震数据保幅重建中的应用 被引量:9
17
作者 唐欢欢 毛伟建 詹毅 《地球物理学报》 SCIE EI CAS CSCD 北大核心 2020年第9期3452-3464,共13页
3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换(NHOPRT)地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而... 3D地震数据不规则采样缺失重建是地震勘探数据处理流程中的重要问题.本文提出了一种基于具有保幅特性的非均匀高阶抛物Radon变换(NHOPRT)地震数据重建方法.在最小二乘反演方程中引入Delaunay三角网格剖分来计算空间不规则加权系数,从而获得最接近完整规则数据的高阶抛物Radon变换域系数.在用SVD求解反演方程过程中,利用高阶抛物Radon变换算子在频率域为指数函数,具有线性可分解特性,将二维空间的高阶抛物Radon变换算子分解为两个独立的一维空间变换算子,减小了变换算子的矩阵大小,从而很大程度地提高了计算效率.理论模型和实际地震数据重建测试证明了本文方法的有效性以及实用性. 展开更多
关键词 不规则地震数据 3D高阶Radon变换 指数函数分解 保幅重建
下载PDF
基于奇异值谱约束的叠前平面波最小二乘逆时偏移方法(英文)
18
作者 李闯 黄建平 +1 位作者 李振春 王蓉蓉 《Applied Geophysics》 SCIE CSCD 2017年第1期73-86,190,共15页
最小二乘偏移能够压制不规则地震数据偏移时产生的偏移假象,但存在计算量巨大、收敛速度缓慢等问题。为了快速地压制偏移噪音,发展了基于奇异值谱约束的叠前平面波最小二乘逆时偏移方法(Plane-wave prestack least-squares reverse time... 最小二乘偏移能够压制不规则地震数据偏移时产生的偏移假象,但存在计算量巨大、收敛速度缓慢等问题。为了快速地压制偏移噪音,发展了基于奇异值谱约束的叠前平面波最小二乘逆时偏移方法(Plane-wave prestack least-squares reverse time migration,PLSRTM)。一方面,通过奇异值谱分析(Singular spectrum analysis,SSA)方法对角度域共成像点道集(Take-Off Angle Domain Common-Image Gathers,TADCIGs)进行预处理,以压制偏移假象,同时引入随机奇异值分解减小SSA的计算量。另一方面,由于偏移速度不准确时叠加成像结果不能有效地压制偏移噪音,将相邻角度平面波偏移结果的差异作为正则化约束加入误差泛函,以改善这一问题。对Marmousi模型数据的成像测试结果表明该方法能够快速压制平面波道集及不规则地震数据偏移产生的偏移假象,改善PLSRTM的成像质量;当偏移速度不准确时,该方法能够得到偏移噪音更少、构造更加连续的成像结果。 展开更多
关键词 最小二乘偏移 平面波 不规则地震数据 奇异值谱分析 共成像点道集
下载PDF
上一页 1 下一页 到第
使用帮助 返回顶部