应用声波测井仪对苏里格地区岩石可钻性分级

合集下载

声波测井仪在勘察中的运用

声波测井仪在勘察中的运用

声波测井仪在勘察中的运用通常的声波测井如声速测井和声幅测井,只记录纵波头波的传播时间和第一个波的波幅,而且只是利用了井孔中非常少的波列。

实际上,换能器在井孔中激发出的波列携带着很多的地层信息。

声波全波列测井采用数字记录方式记录了井孔中的全部波列,利用数字信号处理的方法从全波列中提取所感兴趣的信息,用于声波测井资料的地质解释。

1井孔中的声波及其波型成分在钻孔中,由点声源激发的全波列是由多种波列成分组成的,主要包括纵波、横波、伪瑞雷波和斯通利波等(见图1)。

(1)纵波纵波(又称滑行纵波)是由声源发出的以第一临界角入射到井壁后,在井外地层并靠近井壁且以图1全波列波形图地层中的纵波速度沿井壁滑行的波。

这种波在沿井壁传播的同时,又会以第一临界角为折射角折回井中,被接收器接收到。

(2)横波横波又称为滑行横波,它类似于纵波,从射线声学的角度来看,横波头波是由声源发出且以第二临界角入射到井壁后在井外地层并靠近井壁以地层中的横波速度传播的波,这种波在沿井壁传播时又会以第二临界角为折射角折回井中,被接收器接收到。

(3)伪瑞雷波以相速度介于井内流体中的纵波速度和地层中的横波速度传播的无几何衰减的高频散波。

(4)斯通利波以大于且近似等于井内流体中的纵波速度传播的无几何衰减的微频散波。

2声波全波列测井仪系统组成及工作原理2.1系统组成目前在国内工程勘察行业,适用的声波全波列测井仪器主要为北京大地华龙公司生产的XG-Ⅱ长源距全波列测井仪,该仪器是一套双通道高分辨率、数字化的测井仪,具有分时采样、迭加、滤波、信号增强、抑制噪声以及现场实时计算、实时显示实测波形和测试结果等功能。

测井仪系统由主机、井中全波列声系、连接电缆、平面换能器(用于岩芯的波速测试)和数据处理软件组成。

主要技术指标见表1。

井中全波列声系由一个发射探头和两个接收探头组成(见图2),发射探头距接收探头1的距离为1.05m,距接收探头2的距离为1.25m,两接收探头间距为0.20m。

苏里格地区致密储层岩性测井精细识别方法研究

苏里格地区致密储层岩性测井精细识别方法研究

70随着常规油气勘探开发难度的不断增加,非常规油气资源逐渐得到越来越多的关注[1]。

致密砂岩中赋存的油气资源在非常规资源中占据十分重要的地位[2]。

研究区位于鄂尔多斯盆地伊陕斜坡西北缘,苏里格气田西区。

储层的品质与砂岩的类型有着十分重要的联系,而研究区目的层位沉积环境复杂,储层非均质性强,岩性识别困难。

因此准确识别岩性对于寻找有利储层,对指导油气勘探开发,以及后续石油化工生产都有着十分重大的意义。

1 主要岩石类型及其测井响应特征1.1 主要岩石类型的确定研究区复杂的沉积环境导致该区域的岩性较为复杂,不同岩石类型的测井响应特征较为相似,这是研究区岩性识别困难的主要原因。

通过野外剖面观察以及取心井段的岩心观察分析可知,研究区盒8段与山1段主要发育砂、泥岩地层,并且局部有煤层发育。

对研究区取心井进行岩心取样分析,其中盒8段205个砂岩样本,山1段122个共327个砂岩样本,使用气测法测量样品的孔隙度和渗透率。

统计分析样品的孔隙度及渗透率规律,如表1所示。

研究区盒8段与山1段砂岩的孔隙度和渗透率均与砂岩的粒度有着较大的相关性。

随着砂岩粒度的不断增大,孔隙度和渗透率也随之增大。

在砂岩的主要矿物成分差别不大的情况下,砂岩粒度是影响测井响应值的主要原因之一[4]。

通过镜下薄片观察鉴定,研究区盒8段与山1段的砂岩成分主要为石英岩屑砂岩,成分差异小,在测井上响应上区分不大。

因此,将砂岩按粒度分为细砂岩、中砂岩与粗砂岩三类。

见表1。

苏里格地区致密储层岩性测井精细识别方法研究马浩星 宋荣彩 李可赛 郭恒玮 成都理工大学 能源学院 四川 成都 610059摘要:苏里格地区致密砂岩储层沉积环境复杂,非均质性强,依靠传统的岩性识别方法难以对储层岩性进行有效区分,影响油气开发生产。

根据砂岩粒度将主要岩石类型划分为煤、泥岩、细砂岩、中砂岩与粗砂岩,在此基础上选取自然伽马、声波时差、补偿密度和补偿中子作为主要特征参数,建立两步剥离的Fisher判别模型,判别准确度可达85.7%,识别准确度高,满足了研究区有利储层识别的需求。

苏5、桃7区块不同粒度碎屑岩测井识别方法

苏5、桃7区块不同粒度碎屑岩测井识别方法

苏5、桃7区块不同粒度碎屑岩测井识别方法罗利;朱心万;常俊;周政英;胡振平【期刊名称】《天然气工业》【年(卷),期】2007(027)012【摘要】苏里格气田的苏5和桃7区块石盒子组及山西组储层为碎屑岩储层,但碎屑岩粒度变化大.在岩石成分差别不大的情况下,岩石粒度是造成电性特征差别大的原因之一,能用测井资料较准确地识别出岩石粒度,再按不同粒度碎屑岩建立储层评价模型,则能提高碎屑岩储层评价精度.以苏5和桃7区块取心井为目标,从常规测井资料中提取指示不同粒度碎屑岩的9个特征参数,经聚类分析后形成包括泥岩在内的粉砂岩、细砂岩、中砂岩、粗砂岩、砾岩共6类碎屑岩样本841个,再经非线性映射变换处理,得出不同粒度碎屑岩点在X-Y平面上大多呈不同半径的圆形分布特征;根据矢量点间距离最小原则对测井曲线进行分层,再输入测井特征参数和不同粒度碎屑岩样本,就能用模式识别方法从测井资料中识别出不同粒度的碎屑岩.3口样本井取心段岩性识别符合率在93%以上;对非样本井进行了岩性识别处理,与岩心分类结果比较,也有着较好的一致性.【总页数】3页(P36-38)【作者】罗利;朱心万;常俊;周政英;胡振平【作者单位】四川石油管理局测井公司;四川石油管理局苏里格项目经理部;四川石油管理局测井公司;四川石油管理局测井公司;四川石油管理局测井公司【正文语种】中文【中图分类】P61【相关文献】1.基于阵列声波测井的海陆过渡相碎屑岩地层裂缝识别方法 [J], 尹帅;丁文龙;赵威;孙圆辉;袁江如;丛森2.碎屑岩储层成岩相测井识别方法综述及研究进展 [J], 赖锦;王贵文;王书南;信毅;吴庆宽;郑懿琼;李鉴伦;苍丹3.中子测井在苏里格气田苏东41-X区块马五5段碳酸盐岩储层测井评价中的应用[J], 兰义飞;陈志华;石林辉;刘莉莉4.声波测井在苏里格气田苏东41-X区块碳酸盐岩储层测井评价中的应用 [J], 张吉;陈志华;吴小宁;冯敏5.电阻率测井在苏里格气田苏东41-X区块马五5段碳酸盐岩储层测井评价中的应用 [J], 陈志华;石林辉;王龙;曹艳因版权原因,仅展示原文概要,查看原文内容请购买。

应用阵列声波测井资料计算大北克深地区岩石力学参数

应用阵列声波测井资料计算大北克深地区岩石力学参数

阵列声波测井信号处理的目的是最终获得声波
慢度。以 8 阵列声波测井仪为例,信号接收探头等
距排列,设定接收探头距发射源距离分别为 ,
, … , ,则任意探头 m 距发射源距离为:
= + ( − 1)∆ (1≤ m≤8) (1)
设各个探头接收到的波形信号分别为 ( ),
( ), … , ( ),任意一个采样点 i 的时刻为:
关关系, =0 表示各波形形态间无任何关联, =1 表
示各波形形态完全一致。
1.3 慢度提取
从偶极横波成像测井全波波形提取纵波、横波
时,首先需选取适当的时窗长度 ,在全波波形上
移动窗长并在合理的慢度范围内移动 ,当 与纵波、
横波的慢度接近时,相关系数 最大,这时相关系数
收稿日期:2018–02–11 作者简介:康利伟,工程师,1980 年生,2012 年毕业于西 北大学地质系矿物学、岩石学、矿床学专业,现从事地质 导向技术研究及测井资料解释工作。
昂贵的弹性力学参数测试费用,能预测和判断深部地层发育情况。
关键词:大北克深地区;阵列声波测井;慢度提取;岩石力学参数
中图分类号:P631.8
文献标识码:A
1 阵列声波测井
1.1 首波检测 阵列声波测井资料处理的首要问题是初至波检
测。初至波检测具有多种方法,目前较好的检测首 波波至方法是变阈值法,这也就是首先为目的波形 预先设定一个阈值,当设备检测到该阈值时,即记 为该波形首波波至时间。这种方法能有效地平滑噪 音并突出信号波至,对于连续波形信号,当新的能 量到来时,振幅包络将产生一个起跳,同时瞬时相 位和瞬时频率也将发生跃变,利用这种起跳或跃变 就可以求得首波的初至点,并可以分析波形信号的 包络变化特征[1-3]。

声波测井在苏里格气田苏东41-X区块碳酸盐岩储层测井评价中的应用

声波测井在苏里格气田苏东41-X区块碳酸盐岩储层测井评价中的应用
摘 要:随着全球油气勘探程度的提高,碳酸盐岩油气藏已经成为一个重要的勘探新领域。多年 的勘探实践表明,碳酸盐岩储层与碎屑岩相比研究难度大。碳酸盐岩储集空间以次生孔隙为主的 特点,决定了其储层非均质性更强,受后生成岩改造的影响更大,从而加大了储层预测的难度。因 此,碳酸盐岩储层测井评价一直是国内外测井解释的难题之一。碳酸盐岩储层测井评价的关键是 处理好碳酸盐岩储层的非均质性的问题,其中的难点主要是储集空间的判断、储层流体性质评价、 储层参数的精确评价。 关键词:碳酸盐岩;储层;测井技术;储层参数
1 苏东 41-X 区块基本概况
苏里格气田勘探面积 4×104km2,主要位于内蒙 古鄂托克旗、鄂托克前旗、乌审旗和陕西的定边境 内。该区内天然气资源丰富,资源评价结果显示该 区具有天然气资源量 3.8×1012m3,占盆地天然气总 资源量的 35%。
苏里格气田苏东 41-X 区块位于长庆靖边气田 北部的乌审旗地区。该区域构造属于鄂尔多斯盆地 陕北斜坡北部中带,行政区属内蒙古自治区鄂尔多 斯市乌审旗,勘探面积约为 600km2。地表为沙漠、 草地,地面海拔一般为 1250~1350m,地形相对平 缓,高差 20m 左右。
碳酸盐岩时差较高,胶结物性质和含量对时差 有影响,如钙质较泥质碳酸盐岩时差低,随钙质增加 时差下降;孔隙度、气体含量增加会使时差增大;泥 岩一般显示高值,但泥岩中含钙质、石膏都会使下 降;页岩介于泥岩和致密碳酸盐岩之间。碳酸盐岩 骨架时差一般 182us/m,泥岩为 328us/m。时差曲线 能很好反映碳酸盐岩的致密性,能区分致密和渗透 性碳酸盐岩。
(1)碳酸岩盐剖面:致密石灰岩和白云岩时差 最低(49us/ft、43us/ft);含泥质时差会增大,泥岩和泥灰 岩时差较高,孔隙性、裂缝性石灰岩和白云岩时差明 显增大,甚至出现周波跳跃。

测井成岩相自动识别及其在鄂尔多斯盆地苏里格地区的应用

测井成岩相自动识别及其在鄂尔多斯盆地苏里格地区的应用

(1)利用薄片、铸体薄片和扫描电镜识别取心 岩性和成岩作用。在对取心深度归位的基础上,提取 相应深度测井数据和R值,对井径异常取样点进行 排除后,进行数据标准化,建立成岩相测井数据库。 (2)提取不同岩性成岩相测井信息,针对每种 岩性组成神经网络学习集。训练前随机提取学习集 中70%作为学习样本,剩余30%数据作为检验样 本。对已训练好的神经网络,利用检验样本测试不 同测井组合识别的结果.以衡量不同曲线组合下对 测井成岩相的识别准确性。 (3)利用最佳测井组合。对多口钻井整井进行 测井成岩相识别。总结扩容相、破坏相及保护相空间
2.1识别方法
岩石骨架颗粒是成岩作用改造的主要目标,在相 同流体介质情况下,不同成分颗粒能够产生的成岩作 用类型具有差异性。以碎屑岩为例.其骨架颗粒一般 为石英、长石和岩屑。在石英颗粒含量较高的石英砂 岩中,由于石英的稳定性,因此成岩作用一般以硅质 胶结为主,只有少量石英颗粒边缘发生溶蚀作用。而 在岩屑砂岩、长石砂岩及其过渡岩性中,胶结作用常 表现为硅质胶结、钙质胶结及黏土矿物充填等,溶蚀 作用则可作用于长石、岩屑颗粒及胶结物。这类岩性 成岩相类型具有多类别、多期次的特点。
Abstract:Lithology was used thology influence
on
unit to extracting log feature and identifying diagenetic facies.which could remove the li—
the well-logging informations and control the development types of diagenetic facies in the identification
个神经元的连接权重;P(・)为Parzen窗口函数;c。i 为第一隐层的第k个类别的第i个隐中心矢量;n。为

基于测井资料的岩石可钻性研究

基于测井资料的岩石可钻性研究
声 波 时 差 间 接 确定 地 层 的力 学 性 质 是 石 油 行 业 中 的研 究
声波时 燕I ,- psm
图1 岩 石可 钻性 与声 波 时差 的统 计 关系
重点 ,但 是声波 时差是 岩石 力学性 质 的动 力学 响应参
数 ,岩 石 的 声 波 时差 在很 大 程 度 上 受 到 泥 质 含 量 、矿 物 成 分 、埋 藏 深 度 等 因 素 的影 响 ,因此 必 须 寻 找 一 种新 的 厅法 解 决 这一 问题 。
讲 ,利用测井资料的声 波时差可 以建立地层 的可钻性剖 面。其实现方法如下 :在室 内利用可钻性测定仪确定地 层 岩石 的可钻性级值 ,结合测井资料 中的声波时差进行
信 息系统 工程 I2 l .. 0 1 2 90
9 7

E XCHAN XP I E E 经 验 交 流 GE OFE E ̄ NC d
( 4)岩石埋深 和地质年代 :正常情况下 ,随埋深 的增加 和年代 的久远 ,岩石压实程度不断增加 ,可钻性
变差 ,声波时差 随之下降。
根据 以上分析可知 ,岩石声波时差能够体现岩石的
可钻性 。因而 ,可以运用声波时差曲线进行岩石可钻性
预测 ,但仅用它作为岩石可钻性预测 的唯一变量可能会
造成 同一声波时差曲线对地层岩石可钻性多解的情况 , 因此 ,还必须辅 以其他测井资料 ,即基于测井参 数预测
岩 石 可钻 性 。
机 的 ,它们不能反映层间以及整个井剖面岩石可钻性变
化 情 况 ,不 能 够 建立 起 岩 石 可 钻 性 连续 的变 化 剖 面 另 一
方 面 ,井下 的岩层 与试样 的岩心是有区别 的;此外 ,建
立地区岩石可钻性剖 面需要进行大量 的岩石性质测定 实 验 ,这将花费大量 的人力和资金。基 于以上不足 ,所 以

苏里格气田水平井随钻地质导向技术及应用

苏里格气田水平井随钻地质导向技术及应用

苏里格气田水平井随钻地质导向技术及应用万波;周扬;冯春珍;杨大千;孙建伟;李江博【摘要】确保水平井轨迹顺畅、提高储层钻遇率和实施效果已成为水平井开发的技术难点.以苏54区块所实施水平井为例,通过对水平井实施情况进行分析,结合苏里格气田基本地质特征,从水平段地质模型及地质导向等关键环节入手,总结出了适合苏里格气田水平井开发相关随钻地质导向配套技术.通过现场应用,效果明显,也为同类气藏水平井随钻地质导向提供了有益借鉴和指导作用.【期刊名称】《测井技术》【年(卷),期】2016(040)004【总页数】5页(P508-512)【关键词】随钻地质导向;低渗透砂岩气藏;水平井;配套技术;苏里格气田;苏54区块【作者】万波;周扬;冯春珍;杨大千;孙建伟;李江博【作者单位】中国石油集团测井有限公司长庆事业部,陕西西安710200;中国石油集团测井有限公司长庆事业部,陕西西安710200;中国石油集团测井有限公司长庆事业部,陕西西安710200;中国石油集团测井有限公司长庆事业部,陕西西安710200;中国石油集团测井有限公司长庆事业部,陕西西安710200;中国石油集团测井有限公司长庆事业部,陕西西安710200【正文语种】中文【中图分类】P631.840 引言苏54区块位于苏里格气田西区北部,横跨鄂尔多斯盆地伊陕斜坡、伊盟隆起2个构造单元,区域构造为一宽缓区域性西倾大单斜,坡降7~10 m/km,倾角不足1°,尚处于评价建产阶段,面临着沉积相以及砂体分布复杂、开发难度大等特点。

目的层及以上主要地层有本溪组、太原组、山西组、下石盒子组、上石盒子组和石千峰组地层。

主要目的层为二叠系石盒子组盒8段砂岩和山西组山1段、山2段砂岩,气藏埋深3 200~3 500 m,含气面积662.8 km2。

储层为河流相沉积,呈南北向条带状分布,有效储层岩性主要为岩屑砂岩,少量石英砂岩、岩屑石英砂岩。

储集砂体非均质性强、连续性较差[1]。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
AB
fl l
、 f

【 .0 ^ .0 5 0一 ) e (2 00 8 +00 4 g ] g 1)
传 统 的 岩 石 可 钻 性 分 级 是 用 岩 芯 做 力 学 性 质 、 实钻 速 度 、微 钻 速度 或 碎 岩 比功 等 四种 试 验 来进 行 。用 岩 芯 做试 验 是直 接 的方 法 ,可 靠 性 比 较 高 ,但 是 ,取 岩 芯 不仅 费 用 高 ,而 且 不可 能连 续取芯 ,实现 起来 难度 较大 。
些 学 者 在 研 究 岩 石 力 学 时 ,发 现 声 波 测 井 ( 波 时 差 )、密 度 测 井 与 岩 石 的 岩性 、弹 性 模 纵

接 牧 羹
量 有 着 密 切 关 系 ,利 用 测 井 资料 研 究岩 石 力 学 给 岩 石 可 钻 性 分级 。另 外 电阻 率测 井 、 中子 测 井 、 自然 咖 玛 、 自然 电位 测 井 也 与岩 石 可 钻 性 有 着 间 接 的关系 。
一 3一 3
井 仪 ,声 系 如 图 1 示 , 由 两 个 压 电 陶 瓷 发 射 器 所 ( 声源 )T 、T 和 两 个 压 电陶 瓷 接 收 器 R1 1 2 、R2 组成 。T1 到R1 的间距 称 为源 距 。Tl 2 、T 轮流 发 射 声波 ;T1 射后 ,R1 2 发 和R 分别 接 收 并经 首波 到 时 检 测 ,得 到纵 波 到 时t1 t1 l 和 2 。在 图1 ,发 自T1 中 的声 波 ,到达 R1 R 的传 播 路 径 分别 为 AB 和 和 2 CE AB F D ,旅行 时间 ( 时 )分 别为 : 到
来 ,就 是 声波 测 井 。接 受 器 、 声源 统 称 为 声 系 , 根 据 声 系 排列 尺 寸 的不 同 ,声 波 测井 仪 可 分 为 补 偿 测井 仪 ( HC)、长源 距 声波测 井仪 ( S )和 B LS 阵列 声 波测 井 仪 。声 波 在 井 内地 层 中传 播 由于地 层 岩 石 成 分 、结 构 、 孔 隙 中流 体 成 分 的变 化 ,其 波速 、 幅度 甚 至 频 率 都会 发 生变 化 。只 记录 声波 速度 变 化 的称 为 声速测 井 ( AC),记录 声 幅变 化 的称 为 声幅 测 井 。声 波速 度 测 井 中短源 声系 仅 记 录纵 波 ( 即首 波 )传 播 时差 ,长 源距 声系 可 记 录 纵 波 、横 波 、伪瑞 利 波 、斯 通 莱 波 等 各种 波 列 的 传 播 时差 ,所 以又 称 全 波 声 波 测 井 ,而 阵列 声 波 仪 由于 声 系 复杂 既 可记 录 纵 波 声速 ,又 可 记 录全 波 列声速 ,还 可记 录声 幅 。
1 . 偿声 波测井 仪器 的基 本 原理 2补 目前 测 井 主 要 应 用 的 双 发 双 收 补 偿 声 波 测
作者简 介:刘新田 ( 7_),男, 1 0_ 9 陕西西安人,在读研究生,测
井工程师, 从事现场测井和测井资料解释工作。
第5 期
刘新 田等
应用声波测井仪对苏里格地 区岩石可钻性分级
[ 要] 本文分析研 究 了测井资料与岩石硬度 、强度 及其可钻性 的关系;现场实时钻 井钻速与声波 时差的关 系,分析 结果 摘 表 明声波时差 大的岩石可钻性级值 明显 高于声波时差小的岩石 的可钻 性级值 ,机械 钻速也随声波时差的 变化 而变化。文章 对鄂 尔多斯 盆地 苏里格地 区现场 测井资料 进行 了 计分析 ,对计算岩石可钻 性的数 学模 型进行 了修正 ,并与岩芯 实测岩石 统
可 钻 性 进 行 了对 比 分析 ,结 果误 差 在 6 左 右 。 %
[ 关键 词】声波测井;可钻性分级;钻速 ;数 学模型
Hale Waihona Puke 声 波 测 井 是 目前 最 常用 的岩 性 、 孑 隙度 测 井 L 方 法 之 一 ,应 用 在 不 同 钻压 下建 立 声 波 时 差与 钻 速模 型 ,对钻井 工程 与钻 头选取 提供 支持 。
岩 石 可钻 性 表 达 的是 岩 石 硬度 及岩 石 抗 压 强 度 ,在 钻 井 过 程 中的体 现 就 是 岩 石 破 碎 的难 易 程 度 ,岩 石 可 钻 性 是钻 井 工 程 进 行 钻 头选 型 的重 要 依据 。
岩 石 硬 度 是 岩 石 表 面 对 工 具 压 入 的 反 抗 特 性 。衡 量单位 是P 或MP 。 a a 岩 石 抗 压 强 度 是 岩 石整 块 抗 破 碎 的阻 力 , 根 据 理 论 分 析 ,岩 石 抗 压 入 硬度 为 单 向抗 压 强度 的 (+ 1 2n)倍 。试 验 表 明 ,压 入 硬 度 与 单 向抗压 强 度之 比在5 0 间。 ~2 之
3 一 ■ 技术交流 2

2 化 设备 石 与 工 0 油
1 年第1 卷 0 3
应用声波测井仪对苏里格地区岩石可钻性分级
刘新 田 ,李泰武 ,钟吉彬 ,汤宏平
(. 1 中国石油大学 ( 华东)地球 资源与信息学院 , 山东 东营2 7 6 ) 50 1 (. 2长庆油 田勘探开发研究院低渗透油气田勘探开发 国家工程实验室 , 陕西 安7 0 1 ) 1 0 8 (. 3中国石油测井长庆 事业部 , 陕西 西安 7 0 2 ) 10 1
嚣 射 ●
1 声波时差与岩石力学性质 间的关系
11 . 声波测 井原 理
图 1补偿声波测井仪示意图
将 一 受控 声波 振 源 放 入 井 中 ,声源 发 出 的声 波 引起 周 围质 点 的 振动 ,在 地 层 中产 生体 波 即纵 波和 横 波 ,在 井 壁 一 钻 井 液 界 面 上产 生诱 导 的界 面 波 即伪 瑞 利波 和 斯通 莱 波 。这 些波 作 为地 层 信 息 的载 体 ,被井 下 接 受器 接 受 , 送 至地 面 记 录 下
相关文档
最新文档