设为首页 | 加入收藏
文献检索:
  • 海上拖缆时移地震采集设计实例
  • 两期海上拖缆地震数据采集时海况等因素的不同,会导致时移地震数据中激发点、接收点和拖缆羽角的差异。特别是在拖缆地震方位角较窄和海上钻井平台位置变化时,这些差异会严重影响时移地震的可重复性。在实施南海M油田时移地震项目过程中,对于前期较差地震观测数据(缆数少、激发航线复杂、拖缆羽角大),二期时移地震采集应采用更多的缆数、较密和较直的激发航线、较小拖缆羽角进行观测;时移地震数据处理中,在满足激发点误差与接收点误差之和不超过100m的条件下,通过时移地震数据重构、相对保持拼接处理和适当的互均衡处理可获得高于99%可重复的时移地震数据,能满足储层沉积解释和剩余油气预测的要求。
  • 面向目标与成像的海上多方位观测系统的评价与优化
  • 针对南海深水区崎岖海底、陡坡带、断裂带和较大向斜、背斜构造下的地震成像模糊问题,开展面向目标照明和成像的多方位采集设计评价和优化研究。主要采用正演模拟、照明分析及偏移成像等多项技术对观测系统做精细评价和优化;综合考虑采集成本、施工难度、数据处理及成像效果等因素,认为沿0°、90°、120°和150°等4个方位的多方位观测系统为最佳;建立了一套完整、有效、适用于海上多方位观测系统评价的方法和流程,有利于实际海上复杂区地震数据采集。
  • 利用近可控震源检波器资料提高单炮相关质量
  • 可控震源地震采集是利用参考信号与母记录进行相关得到地震单炮记录,且两种信号相关性越强则相关结果越好。可控震源激发的参考信号通过震板振动传递给地表,因震板与地表存在耦合关系,导致震板激发振动信号与地表接收振动信号并不完全相等,故利用参考信号与母记录做相关难以得到最佳相关结果。基于此,提出一种利用距可控震源较近(约为15m)检波器记录到的地表振动信号与母记录进行相关的方法;经过实际资料验证,该方法可有效提高相关单炮记录的信噪比和分辨率。
  • 检波器耦合响应定量分析方法及应用
  • 针对检波器与大地耦合响应对地震数据的影响,设计了检波器与大地耦合响应的测定装置和定量分析方法:利用单自由度有阻尼的自由振动系统,定量描述检波器与大地耦合响应过程;通过等效固有频率和等效阻尼比两个参数,实现检波器与大地耦合响应的定量分析;基于QT编程语言,研发了地震检波器耦合响应定量分析软件。检波器耦合试验数据分析结果验证了该定量分析方法的可行性;将集成的定量分析软件应用于现场试验,进一步证实了该方法的实用性和软件的稳定性。
  • 高保真高信噪比地震资料的获取方法
  • 以检波器组合为例,通过射线追踪模拟,揭示了组合基距和炮检距变化与反射波波形特征和信号频谱的关系,讨论了组合基距变化对信号保真度和分辨率的影响。利用波动方程正演,对比分析了组合基距变化对小断层和砂岩尖灭的成像结果,说明常规检波器组合对小尺度地质目标存在平均效应,会降低横向分辨率。结合实际资料分析,进一步说明常规检波器组合方法获得的地震记录是非保真的,对有效信号造成极大的伤害,地震分辨率降低明显。因此常规的组合压噪方法是不科学的,应该先保真原始地震道信号并保证适度能量,然后在数据处理中去噪、提高信噪比。由此提出的一点法地震技术获取的地震资料具有高保真、高信噪比、高分辨率的特点,在苏里格气田开发的多个区块中取得了很好的应用效果。
  • OVT域内插炮检线压制采集脚印
  • 塔里木盆地塔北MS工区奥陶系目的层呈现的严重采集脚印缘于过大的炮检线距和目的层的低信噪比。为此,提出在数据处理阶段压制采集脚印的两种策略:基于常规数据规则化+四维去噪+叠后去噪技术组合,将采集脚印的压制从叠后提到叠前,但难以根除采集脚印;基于OVT域五维插值的内插炮检线技术+OVT域偏移技术,将采集脚印的压制回溯到观测系统设计阶段,即模拟炮检线距减半的现场数据采集,消除了采集脚印产生的根源,用处理技术使得采集阶段形成的脚印最小化,达到"源头控制"的效果。
  • OBC水陆检数据标定因子估算方法
  • 标定因子估算是OBC水陆检数据合并处理的关键。针对常规标定因子扫描估算方法耗时的缺陷,通过计算水中检波器数据和陆地检波器数据的自相关函数及其互相关函数,建立标定因子特征方程;求解标定因子特征方程以确定最佳标定因子,用于水中检波器数据和陆地检波器数据的合并处理。不但可以消除海水鸣震多次波干扰,而且还可吸收陆检数据的低频成分和水检数据的高频成分,有效拓宽海底电缆数据的有效频率带宽,提高地震数据信噪比和分辨率。实际数据处理结果表明了该方法的有效性和实用性。
  • 地震低频信号对关键处理环节的影响分析
  • 首先研究与低频关系最直接的面波压制问题,然后分别讨论了低频信号对反褶积、速度分析、动校正、叠加、叠前偏移等关键处理环节的影响。理论模型测试与实际资料处理表明:1面波压制会损失部分低频有效信号,区域滤波法比带通滤波和F-K滤波的效果要好;2低频缺失较多时反褶积会出现虚假同相轴,使一次波失真;3低频缺失会影响速度谱的能量团聚焦性,甚至会出现假的能量团,继而影响速度场的精度;4缺失低频信息会导致动校拉伸增大,缺失10Hz以下的低频信息会使动校拉伸增大达20%;5叠加后信噪比提高,但低频缺失带来的虚假同相轴比叠前道集明显;6从模型研究的结果看,低频缺失对偏移的影响比人们预期的要小。
  • 一种稳定高效的等效Q值反Q滤波算法及应用
  • 提出了一种基于等效Q值的反Q滤波算法,在补偿高频耗散能量的同时,对速度频散造成的相位畸变也进行了校正。该算法通过光滑的阈值控制增益函数实现稳定性控制,在补偿高频时,引入了一种变频带的计算方法,压制由于频率域的折返效应产生的高频噪声。与一般品质因子Q不同,由于在每一样点的吸收补偿是由该点的Q值唯一决定的,可利用扫描方法直接求取等效Q值。对Q扫描后的叠后反射数据采用对数谱比求频率导数的算法求取等效Q值场,避免了薄层调谐对频谱的不利影响。理论测试与实际资料应用的结果表明,等效Q值反Q滤波算法对含噪地震数据是稳健的,可在保持信噪比的同时提高地震资料的分辨率。
  • 应用同步挤压小波变换去除面波
  • 同步挤压小波变换(SWT)是一种新的连续可逆时频分析方法,具有较高的时间和频率分辨率,利用该变换可去除地震资料中的面波。将时间域地震记录变换到时频域,在时频域准确识别面波区域并切除,再将切除面波后时频域信号反变换到时间域,就得到去除了面波的地震记录。理论模型合成数据和实测资料的处理结果表明:该变换可以有效去除面波而不损失有效反射波,且效果优于连续小波变换、S变换和f-k滤波等方法。
  • 步进迭代法测井与井旁地震资料匹配处理
  • 针对测井资料在地震频带范围内与井旁道难以匹配的问题,提出了步进迭代法对测井资料进行降频处理:1根据测井资料计算初始反射系数,利用自相关法从井旁道提取初始地震子波;2用初始子波作为约束,从井旁道提取新的地震子波;3用初始反射系数作为约束,从井旁道提取新的反射系数;4判断新的地震子波与反射系数的褶积是否与井旁道匹配,如果不匹配则用新的地震子波和反射系数替换初始地震子波和反射系数,并回到步骤2。如此迭代,直到地震子波与反射系数的褶积与井旁道吻合。该方法降低了测井资料的频率,使两种资料在中、低频范围内相互匹配。匹配后的反射系数与高频地震子波褶积,形成拟井间地震记录,用于联合拓频,有效地提高了地震资料分辨率。
  • 横向各向同性介质拟声波一阶速度—应力方程
  • 基于弹性波的基本理论,首先推导了VTI介质拟声波一阶速度—应力方程,通过坐标变换,还得到了TTI介质拟声波一阶速度—应力方程,并从能量守恒的角度证明了上述方程比基于频散关系推导出的拟声波方程更加稳定。针对各向异性参数ε〈δ的模型,对新推导的方程引入横波项,消除了残留的横波产生的数值模拟不稳定现象。通过对均匀二维模型和逆冲模型的正演模拟表明:新推导的方程不仅能很好地维持qP波的运动学特征,而且能够很好地适应任意旋转角度的TTI模型,并且在旋转角变化的TTI介质中数值稳定性要好于基于频散关系推导出的拟声波方程。
  • 基于横纵波速度比值迭代的纵横波联合反演方法
  • 在利用弹性阻抗EI反演纵、横波速度和密度时,若无速度信息,通常设定横纵波速度比k为常数,从而与真实的k值存在误差,将导致反演结果不准确。进一步研究发现,当缺乏大入射角信息时单独基于EI反演纵、横波速度和密度,反演结果易受噪声影响而很不稳定。为此,利用参数分辨矩阵分析了k值对反演结果的影响,发现纵波速度和密度几乎不受k值误差的影响,横波速度对k值的变化十分敏感。基于此提出了非线性k值迭代方法,通过逐次迭代求取准确的k值,相应得到准确的横波速度,进而基于Aki-Richards转换横波反射系数近似方程推导了转换横波弹性阻抗(SEI)方程,提出了基于k值迭代的纵横波联合反演方法。利用层状模型和Marmous 2的部分模型测试了上述方法,结果表明该方法能使k值较快地收敛到真实值,可更加稳定地反演纵、横波速度和密度。
  • 基于子波重构的时—空域高斯束正演方法
  • 近年来高斯波包法和频率域高斯束法均已广泛应用于地震波场模拟中,但高斯波包法只能模拟Gabor子波点源波场,频率域高斯束法虽能模拟任意点源波场,但需要频繁的傅里叶变换,计算效率较低。为此,将基于Gabor分解的子波重构技术应用于时—空域高斯束波场计算中,提出了一种可以模拟任意点源函数透射与反射波场的地震正演方法。采用子波重构技术可将任意震源函数分解为一系列可调制的且具有不同时间延迟的Gabor子波,而调制的Gabor子波波场在非均匀介质中可使用不同角度出射的时—空域高斯束向外传播,因此只需将不同角度、不同时间延迟的时—空域高斯束加权叠加就可得到地下或地表上任意一点记录的波场。光滑非均匀介质和层状模型的数值模拟结果表明:文中方法利用不同时间延迟的Gabor基函数重构震源子波,可模拟任意点源函数产生的地震波场,无需频繁的傅里叶变换,大大提高了计算效率,模拟精度可与有限差分法相比拟;文中方法具有局部的时—空域特性,计算时不需要波场的整体外推,避免了有限差分法中因差分代替微分所引起的数值噪声,可以较准确地描述地震波传播过程中的反射、透射和几何扩散等现象以及介质的非均匀性对波前相位的影响,可以较准确地刻画起伏地表对地震波走时和振幅的影响以及低幅度向斜形成的回转波等特殊波场现象。
  • 一般TI介质中多震相旅行时层析成像——以井间成像为例
  • 传统的各向异性介质中的正、反演研究大多基于弱各向异性假设下的VTI介质,且主要是有关初至波的射线追踪和相应的旅行时反演。然而,随着勘探精度要求的不断提高,一般普适的(含强)TI介质中的多震相射线正、反演方法,则成为急需解决的问题。鉴于此,本文将基于各向同性介质多震相分区多步最短路径射线追踪算法推广到一般各向异性TI介质,结合相、群速度导数的一阶旅行时扰动方程,采用共轭梯度求解带约束的阻尼最小二乘问题,进而提出了一种利用多震相旅行时进行一般TI介质弹性参数反演的方法。井间地震数值实验结果表明,多震相旅行时联合反演可有效提高成像分辨率,验证了方法的有效性和正确性。
  • 多波深度域建模与成像一体化技术
  • 在复杂油气藏勘探中,常规方法在处理多波数据时通常将纵波和转换横波分别做类纵波处理,忽略了地震波的弹性矢量波特征,并且成像精度严重依赖于波场分离的精度。针对这些问题,本文给出了解耦的高斯束弹性波波场反向延拓公式及多波成像公式以实现叠前深度偏移;然后建立了针对弹性多波建模的最小二乘反演方程,以实现基于高斯束角道集的反射层析建模。将深度偏移和反射层析联合在同一处理流程中,前者为后者提供建模角道集和成像质量监控,后者为前者提供成像速度场,两者交替迭代形成深度域建模与成像一体化技术。模型试算和实际资料处理结果验证了一体化技术的正确性和有效性。
  • 逆时深度偏移中的子波拉伸校正
  • 文中提出一种在逆时偏移的波场延拓后、应用成像条件前对地震数据进行子波拉伸校正的方法,在延拓过程中分别提取各网格点震源波场和接收点波场在成像时刻的地震子波,对其进行压缩处理,再应用归一化波场分离互相关成像条件,得到拉伸校正后的成像结果。该方法不改变成像结果的振幅与相位。其中,子波压缩计算需要的各网格点的反射角和局部地层倾角由地震波场的坡印廷矢量(Poynting vector)计算得到。数值算例表明,该方法能够校正由于反射角增大引起的子波拉伸效应,提高逆时偏移的成像品质,为AVA/AVO分析提供高精度的输入道集。
  • 利用双路径积分算法进行高密度偏移速度建模
  • 针对基于双路径积分的高密度偏移速度分析算法求出的偏移速度值偏高及有效速度选取方法不够完善的问题,通过详细分析加权参数对速度建模精度的影响,提出采用拉冬变换求取CRP道集的曲率作为加权参数,有效地解决了拾取速度偏高的问题;在共成像点道集中应用相干阈值约束成像结果,选取有效成像点,显著提高了有效偏移速度的精度,改进了基于双路径积分的高密度自动偏移速度建模方法。该方法的应用不仅获得了高精度的速度模型,而且大幅度降低了人力成本,实现了低信噪比地震资料的高质量偏移成像。模型试算和实际资料处理验证了该方法的有效性。
  • 井震联合识别复合砂体中单一河道——以朝44区块扶余油层为例
  • 以朝阳沟阶地朝44区块扶余油层为例,分析了河流—三角洲复合河道砂体的地质与测井识别标志、地震响应特征以及井震结合精细刻画方法。研究发现:1根据相邻河道的期次及规模差异、复合砂体组合部位的微相、砂体厚度及宽厚比的异常变化,构建出六类测井与地质识别标志;2单期河道砂体接触关系可分为同期对接型、非同期切割型、河间泥及河间砂型四类,分别对应"同相轴能量降低、视厚度变小"、"类断层式波形错动"、"同相轴突然下断"、"能量减弱、视厚度增大"的地震响应特征;3采用"地质认识定位、测井识别定点、属性趋势引导定向、地震剖面定型"的方法,井震结合有效刻画复合砂体中单一河道边界及连通关系。
  • 敏感流体因子定量分析与直接提取
  • 基于叠前反演的流体因子提取是当今油藏地球物理的研究热点,目前的流体因子提取存在两个问题:1流体因子的种类繁多,增加了问题的复杂性和盲目性;2流体因子往往是通过叠前反演的结果间接计算获得,具有较大的累计误差。针对这两个问题,建立了一种敏感流体因子定量分析与直接提取方法,提出了流体因子敏感系数的概念,用于定量分析流体因子区分不同流体的能力,并得到研究工区最敏感的流体因子;在得到敏感流体因子之后,将该流体因子近似方程和叠前反演相结合,建立流体因子直接提取的目标函数,通过求解该目标函数,直接得到了敏感流体因子。实际应用结果表明,该方法的流体因子敏感性定量分析过程合理,流体因子直接提取结果客观准确,益于进一步推广应用。
  • 横波预测技术在苏里格气田储层预测中的应用
  • 横波资料的缺乏严重制约了苏里格气田A区储层岩性、物性及流体识别等工作的开展,以往横波预测结果精度低,识别能力差,已不能满足储层预测需求。为此,基于地质、钻井、测井、地震等资料,针对苏里格气田A区上古生界砂泥岩地层的实际情况,确定了构建Xu-White模型所需的各项岩石物理参数,进行了横波预测研究,结合岩石物理分析、叠前弹性参数反演等技术实现由点到面的储层预测,获得的预测横波曲线与实测横波曲线吻合度较高。结果表明:1使用纵横波速度比、剪切模量以及拉梅模量/剪切模量、体积模量/剪切模量可有效划分砂、泥岩,利用拉梅模量/剪切模量—纵波阻抗、体积模量/剪切模量—纵波阻抗、泊松比—纵波阻抗交会等是识别流体的有效手段;2利用文中方法预测的盒8段储层、气层厚度与实钻数据误差小,储层分布与该区地质规律一致。
  • 砂泥岩薄互层层位精细标定方法及应用——以东营凹陷滨东地区沙四上亚段为例
  • 以滨东地区沙四上亚段(Es_4~s)滩坝砂岩为例,首次提出了在地震相模式指导下,采用"二次标定法则"进行砂泥岩薄互层精细标定的方法。针对研究区滩坝砂岩位于油页岩集中段下部、岩性以砂泥岩薄互层为主、埋藏较深、延伸范围小的特点,在对地震剖面进行极性判断的基础上,以地震相模式为指导,利用"二次标定法则"与"时变子波"相结合的方法对砂泥岩薄互层进行精细标定,最后通过分析目的层段反射系数序列与振幅的对应关系验证精细标定结果的准确性。
  • 复数道动态匹配追踪算法的改进
  • 针对复数道动态匹配追踪算法中采用瞬时频率求取的匹配原子的主频经常是没有物理意义的负数,且很不稳定的问题,在匹配原子主频的求取过程中,提出了基于连续相位的瞬时频率求取方法。该方法摒弃了主值相位,引入连续相位,解决了主值相位存在不连续突变点的问题;而且瞬时频率并不是直接由相位求导获得,而是通过阻尼最小二乘反演的方法求得。实验结果表明,利用连续相位求取的瞬时频率准确度高,且不存在无物理意义的负频率值。模型和实际资料表明改进的匹配追踪算法具有很好的稳定性,能够适应复杂地震信号的分解和重构。
  • 巴布亚前陆盆地东部构造特征及油气圈闭类型
  • 在调研前人研究成果的基础上,以地震资料最为详实的巴布亚湾海域区块为例,从区域构造演化入手,探讨研究区构造特征及油气圈闭类型,得出以下认识:1盆地主要经历了裂谷前基底期、冈瓦纳同生裂谷—裂后期、珊瑚海同生裂谷—裂后期、达拉伊弧后伸展期和前陆盆地期等演化阶段,沉积地层垂向呈"二元结构",即前陆构造层叠加在达拉伊弧后及其下的裂谷—裂后期构造层之上,并呈楔状向南西西快速减薄超覆于东部高地之上;2巴布亚湾海域区块形成了以欧拉褶皱带为前陆冲断带、以莫尔比拗陷为前渊带、以东部高地为前陆隆起和以布莱盆地为隆外拗陷的前陆盆地系统;3初步分析圈闭类型和油气成藏条件后认为,裂谷—裂后和弧后伸展阶段形成优质海相烃源岩,前陆期形成良好的储盖组合;前陆冲断带发育的成排、成带逆冲背斜为首选钻探目标,前渊带的扇三角洲岩性体、浊积水道和前陆隆起带潜山构造及复合圈闭为潜在目标。
  • 海洋可控源电磁与地震一维联合储层参数反演
  • 海洋可控源电磁法能识别高阻油气藏,可大幅提高钻井成功率。以Archie公式和Gaaamann公式为基础,建立油气储层的岩石物理关系,并将其引入目标函数,采用模拟退火结合粒子群算法的优化反演方法进行一维海洋可控源电磁和地震数据联合反演,直接得到储层孔隙度和饱和度参数。数值模拟反演结果表明,该算法不依赖初始值,反演效果好。
  • EMD在大地电磁信号分析中的问题及解决方法
  • 以经验模态分解(Empirical Mode Decomposition,EMD)为核心的希尔伯特—黄变换(Hilbert-Huang Transform,HHT)被越来越多地应用于大地电磁信号处理中,且EMD分解中的曲线拟合、端点效应和筛选停止准则等几个问题处理的好坏直接决定HHT方法应用成败。结合实际大地电磁信号分析的特点,以仿真和实测大地电磁数据为例,采用对比分析的方法,对上述三个主要问题进行了研究,给出了问题的解决方法,结果为大地电磁信号分析中EMD分解方法的选择提供了有价值的参考。
  • 东方地球物理勘探有限责任公司研究院地质研究中心
  • 地质研究中心是在原物探地质研究院基础上发展而成的物探资料解释综合研究单位,参与了国内多个油田的发现,是地震资料解释行业标准的主要起草者。现有员工330人,其中教授级高级工程师4人,高级工程师75人,15年以上丰富工作经验的176人。
  • ABSTRACTS
  • [采集技术]
    海上拖缆时移地震采集设计实例(吴意明;凌云;张永江;郭向宇)
    面向目标与成像的海上多方位观测系统的评价与优化(吕晓春;孙常新;成景旺;李鹏;邹明俊;林泽进)
    利用近可控震源检波器资料提高单炮相关质量(宋智强)
    检波器耦合响应定量分析方法及应用(于富文;魏继东;冯玉苹;徐维秀)
    高保真高信噪比地震资料的获取方法(欧阳永林;宋炜;曾庆才;胡新海;黄家强;姜仁)
    [处理技术]
    OVT域内插炮检线压制采集脚印(段文胜;裴家定;李飞;孙建华;刘正文;宋文杰)
    OBC水陆检数据标定因子估算方法(高轩;弓小平;高少武;王成祥;李建峰;张红英)
    地震低频信号对关键处理环节的影响分析(张军华;张在金;张彬彬;梁鸿贤;傅金荣;郭见乐)
    一种稳定高效的等效Q值反Q滤波算法及应用(吴吉忠;杨晓利;龙洋)
    应用同步挤压小波变换去除面波(刘晗;张建中;黄忠来)
    步进迭代法测井与井旁地震资料匹配处理(宋建国;程亮;孟宪军;魏国华)
    [地震模拟]
    横向各向同性介质拟声波一阶速度—应力方程(徐文才;杨国权;李振春;李河昭)
    基于横纵波速度比值迭代的纵横波联合反演方法(杜启振;霍国栋;陈刚)
    基于子波重构的时—空域高斯束正演方法(印兴耀;王欣;杨继东)
    [偏移成像]
    一般TI介质中多震相旅行时层析成像——以井间成像为例(黄国娇;白超英;钱卫)
    多波深度域建模与成像一体化技术(秦宁;王延光;杨晓东;单联瑜;梁鸿贤;王常波)
    逆时深度偏移中的子波拉伸校正(杨佳佳[1,2];何兵寿[1,2];陈婷[1,2])
    利用双路径积分算法进行高密度偏移速度建模(梁瑶[1,2,3];王真理[1,2];陈雨红;李学良[1,2];段秋梁)
    [综合研究]
    井震联合识别复合砂体中单一河道——以朝44区块扶余油层为例(肖佃师[1,2];张飞飞;卢双舫;陆正元;谷美维;张鲁川)
    敏感流体因子定量分析与直接提取(杨培杰[1,2];董兆丽;刘昌毅;管晓燕;穆星)
    横波预测技术在苏里格气田储层预测中的应用(刘雅杰[1,2];李生杰;王永刚;夏勇辉)
    砂泥岩薄互层层位精细标定方法及应用——以东营凹陷滨东地区沙四上亚段为例(栗宝鹃[1,2];董春梅[1,2];林承焰[1,2];张宪国[1,2];任丽华[1,2];宋亚民)
    复数道动态匹配追踪算法的改进(张繁昌;刘汉卿[1,2];张立强;陈发亮;印兴耀)
    [地震地质]
    巴布亚前陆盆地东部构造特征及油气圈闭类型(李建英[1,2];陈旭;董玉文;王红梅;杨建房;周晓明)
    [非地震]
    海洋可控源电磁与地震一维联合储层参数反演(徐凯军;杜润林;刘展)
    EMD在大地电磁信号分析中的问题及解决方法(蔡剑华;王先春)
    东方地球物理勘探有限责任公司研究院地质研究中心
    ABSTRACTS
    《石油地球物理勘探》封面

    主管单位:中国石油天然气集团公司

    主办单位:东方地球物理勘探有限责任公司

    主  编:钱荣钧

    地  址:河北省涿州市11信箱石油学会

    邮政编码:072751

    电  话:0312-3822282 3822804

    国际标准刊号:issn 1000-7210

    国内统一刊号:cn 13-1095/te

    单  价:42.00

    定  价:294.00


    关于我们 | 网站声明 | 合作伙伴 | 联系方式 | IP查询
    金月芽期刊网 2017 触屏版 电脑版 京ICP备13008804号-2