• / 130
  • 下载费用:10 下载币  

岩石物理

关 键 词:
地质 储层 沉积 地化 层序地层
资源描述:
第二章地震岩石物理基础岩石骨架孔隙 /流体干燥岩石 (孔隙空 )饱和岩石 (孔隙充满流体 ) • on of • to on on 石物理研究目的• 地球物理勘探所获得的物理量与地下储层参数的对应关系通常是非常复杂的,它们之间的关系通常是通过间接途径获得的,这二者的对应关系的确定程度必然影响到地球物理探测结果的储层解释能力及其效果。。岩石物理研究目的• 对于地震勘探来说确定地震波特点的因素除激发、接收条件外主要受岩石的弹性模量、密度和吸收等特性的影响。而这些特性又与岩石成分、孔隙度、埋深、孔隙内流体性质、压力、岩层的不均匀性以及其它的地质特性密切相关,了解岩石物性与地震波特性的关系可以更好地研究储层性质(孔隙度、渗透率)及其状态(饱和度、孔隙压力等)。岩石物理研究目的• 随着地震勘探技术的发展,现在地震数据通常还用来分析岩性、孔隙度、孔隙流体和饱和度等储层特性。随着地震勘探技术研究的深入,近年来关于岩石及其流体性质与波速关系的实验室研究方面取得了很大的成果岩石物理研究目的• 以 些研究结果意味着储层参数可能可以由地震波中得到,由于岩石物理是连接地震数据和储层属性与参数的桥梁,它在近年来新发展起来的一些新方法、新技术中起了重要的作用,如四维地震、地震岩性识别、 “亮点”技术和球物理的未来——精确描述储层岩石特性• 地球物理的未来在于岩石物理性质描述• 为了更好地服务于油气勘探业务,地球物理的最终目标将是地下岩石所有重要性质的一个完整而精确的再现。• 要达到这一理想目标,我们不仅要知道那里有油气沉积,还要知道沉积有多大,产出有多快,开发该油藏所需井的最佳数量是多少,以及钻井和生产时会遇到何种条件。重要岩石性质有以下几点:• 岩性是什么 ——是储层岩石性质、非储层岩石性质还是盖层性质?• 岩石的几何形状是什么 ——它们有多深,储层是在圈闭中吗?• 储层岩石的孔隙度是多少• 储层岩石的渗透率是多少? — 储层的流体类型是什么?• 岩石中的孔隙压力是多少?• 大多数情况下,这些问题的答案是包含在地震数据中的。• 地球物理进步是显著的。在大公司大型研究与开发预算的推动下,已经看到共深度点成像、三维和多分量数据采集、多种震源、地球物理工作站和虚拟现实系统。数据准备•解释将永远成为地震数据处理的关键部分。为了获得可解释的数据集,我们可以考虑三级数据准备:• — 第一级是普遍地产生叠加记录剖面数据(静校正、噪声抑制、时差、偏移)。• 目的是获得清晰可见的数据集。数据准备•— 第二级是提取在叠加(或未叠加)剖面中不能轻易辨别的各种属性。这些属性包括:希尔伯特变换(瞬时相位、频率和振幅)、 λμ地震道和远地震道叠加、阻抗和速度。数据准备的第二级产生了使熟练解释人员能用来指导和约束解释的数据格式。• — 第三级是从地震数据中提取基本岩石特性(岩性、孔隙度、渗透率、流体类型和孔隙压力)。• 这一级数据是以岩石为基础的综合的关键部分;以这样一种方式提交数据是为了对所有参与发现和开发油气的学科都有用。• 早期以及正在进行的是将各种第二级地震属性与基本岩石特性联系在一起的尝试一般是通过统计分析实现的。• 当你离开用来导出统计关系的控制点时,统计关系会变得不可靠。• 根据定义,确定性方法是建立在物理定律基础上的,因此应在整个数据集中都是相对不变的。• 这两种方法都有它们的用途,统计技术在完全确定性方法还没有开发的地方尤其有用。地震岩石物理学, 二十年来,地震岩石物理学已在地球物理界的某些方面进行了实践,它的影响正在增加,它正在演化成为以岩石为基础的综合方法的一个整体部分。地震岩石物理学•在地震岩石物理学中,为了获得基本岩石性质,要测量两个独立的参数并绘制交会图。最普遍应用中,地震岩石物理可从地震数据中提取对剪切和压缩阻抗的估计并绘制交会图,以获得岩性、孔隙度和流体信息。• 如果地球物理的长远未来在于 岩石 ,那么通向未来的桥梁就是以岩石为基础的 综合 。地震数据是唯一的,它们可提供全三维的数据体,不仅能够产生有价值的信息,而且还能提供与这些信息关联的框架。• 在地质学家谈论岩相的地方,岩石物理学家谈论的是孔隙度,而工程师想知道渗透率,只提供阻抗的岩石物理方案确实不能提供真正有价值的信息或真正有用的框架。• 地震岩石物理学是一门相对较新的学科,目的在于获得地震数据以外的有关岩石物理类型的信息。• 但更重要的是,地震岩石物理的具体目的应是从地震数据中获得像岩性、孔隙度和流体成分这样的信息,并且不考虑对油藏工程师和地质学家没有帮助的特性,如 λ 地球物理的未来是明确的,它就在岩石中。测井分析提供了有用的类比,测井告诉我们,像岩性、孔隙度和流体成分这样有用的岩石性质可从一系列测量中获得。• 地震和地震岩石物理的未来在于把地震数据转变成像岩性、孔隙度和流体这类性质的创新方法 ——这是以岩石为基础的综合方法的关键元素,并且是通向地球物理未来的桥梁 。 ( & 岩石的组成 ,包括组成岩石的矿物成分,岩石内部的孔隙度,岩石的饱和状态和孔隙流体的性质等;2. 岩石内部的结构 ,包括矿物颗粒的大小、形状及胶结情况,岩石内部的裂隙和其它不连续界面等;3. 岩石所处的热力学环境 ,包括温度、压力和地应力场等。影响弹性波速因素弹性常数了解弹性参数对于我们准确理解岩石物性有很大的帮助,同时它是我们研究储层和识别油气的重要参数。普遍使用的弹性参数主要有:弹性模量、密度、纵横波速度、泊松比等。弹性模量是岩石在外力作用下发生的伸缩、剪切和体积变化的特征参数,是联系应力 —反映在外力作用下,岩石发生的伸缩变化剪切模量 :在外力作用下,岩石外形发生的剪切位移;体积模量 :在外力作用下,岩石体积发生的变化,反映岩石的可压缩性;2、纵波速度和横波速度 横波速度 23/4  2  (1)/( ,平面波模量 M , ,M ,E 2 3 43 K 221 )1(    E   )1(2   3 43 M )21(2 E )23(43()(2  )(2 23(2 23     2k 221  松比与岩性及流体成分的关系(横坐标为纵波速度)只利用纵波速度,难以区分砂岩和泥岩。而综合泊松比和纵波速度,情况就不一样。s s P/VS be P/ 2, = 0– P/ = P/ 2, = 1/3 ( P/ , = 0)1as a of et us s to ,  . be by 3422  23123222222 in () (). No we us or 波速度和横波速度的关系泥岩: c=b=c=b=c=b=c=b=c=b=  993)993)还提出用抛物关系来回归不同岩石的密度(或横波速度)与纵波速度之间的关系:  2) 1:不同岩石的密度与纵波速度的关系系数 石灰岩 白云岩 砂岩 页岩 膏岩A 灰岩 白云岩 砂岩 页岩 煤A - - 不同岩石的横波速度与纵波速度的关系is as of to is of of to in   v v T  v 隙度分为绝对孔隙度和有效孔隙度,一般在沉积岩中不连通的孔隙其体积一般不超过孔隙总体积的百分之几,可是当 φ值较低时(φ<10%),此值可能增大。纵横波速度与孔隙度的关系岩石埋深增大时,孔隙度一般变小。此外,孔隙度的变化不仅沿垂线方向,而且沿岩层层理方向都很大。处于原生埋藏状态的沉积岩,孔隙里始终充填着液体,如矿化水,有时还有石油或天然气。• 在其它因素相同时,孔隙度大的岩石波速低。这是因为岩石孔隙度的变化会有效地改变岩石密度,而且与孔隙度随应力变化而变化的变化率相关。这种变化率有赖于孔隙形状。研究指出,孔隙总体积对速度的影响是主要的,孔隙大小和形态及矿物颗粒大小和形状对速度的影响要小得多。• 在相同应力条件下,圆形孔隙形变比狭长孔隙形变小得多。• 当孔隙度相同时,骨架颗粒细、孔隙小的岩石比颗粒粗、孔隙大的岩石的速度高约百分之几。• 考虑到这一点,在理论研究时,作为第一级近似,可以假定固结岩石的固体物质成分和性质是均匀分布的,其骨架则为宏观均匀的各性同性介质。• 随着孔隙度的不同,同一种岩石的速度值可以在很大范围内变化。另一方面,固相成分不同的岩石也可能具有相同的波速。1956) 时间平均方程mf  11此方程最适合的是中等孔隙度的砂岩( 10%<φ<25%)。时间平均公式在声波测井中曾经被广泛应用于计算岩层的孔隙度,但是人们很快发现此公式的不足之处,特别是这个公式过高地估计了波在含粘土砂岩中的速度。1980)提出了一个非线性经验公式用于描述孔隙度和波速的关系,即:fm   )1( 2该公式可适用于更大的孔隙范围,包括高孔隙度的非固结沉积物。但是这个公式也忽略了其他参数的作用。1984)把时间平均公式修改为纯经验公式,用于描述速度和岩石孔隙度的关系,即1其中 A、 个公式也过分简化了速度与岩石物性的关系,特别是忽略了粘土含量对砂岩横波速度的作用。986)等人对此作过综合研究。通过对 80块不同泥质含量和孔隙度的砂岩样品,在不同的压力下测试它们的纵波速度和横波速度。研究了纵横波速度比υp/υm。泥质含量、含水饱和度、弹性模量等与速度的关系。研究表明,纵横波速度 υ何拟合泥质砂岩的纵横波速度 υ面公式给出了围压为 40隙压力为 p p p ps s s B C B C N    对于非纯砂岩,实验室测量和测井统计表明:孔隙度、泥质百分比含量 数 多作者的统计结果见表p s 岩样测试系统实验设备主要包括三轴高压釜、声波测试系统、压力控制系统、流体饱和系统、计算机伺服控制系统。轴高压釜2、 算机伺服控制系统饱和岩样测试系统加载路径储层中流体对波速的影响储层中流体对地震波性质影响的研究是利用地震资料进行储层描述和储层监测的最重要内容之一,虽然弹性波在含流体的孔隙固体中传播问题早在 1956年 立了粘滞流体饱和的孔隙固体中弹性波的传播理论。但不幸的是在传统的地震勘探中往往忽略了它们的影响由于孔隙中流体对波速有非常明显的影响。在储层地球物理的研究中,含流体介质的地震波特性的研究具有特别重要的意义。气的速度• 在储层条件下 , 炭氢类气体和非炭氢类气体的速度通常在 200—400m/s , 他们随着温度增加而增加 , 在低压范围内随压力增加而减少 , 随密度或他们的分子量增大而减小 , 与其它流体相比较 ,含气储层的地震波速度的研究相对简单 。水的速度• 研究含水储层的地震波速度的变化比较复杂,图给出了速度随温度和压力变化的曲线在常压下约 73℃ 时速度最大,大于 730时随温度增加而有所减小。压力越大,与温度的关系越弱,其最大值对应的温度越高。图 2 - 3 水中速度与温度和压力的关系 ( L aw s H 1963 ) 轻烃速度• 与烃类气相 , 烃类液体的速度随温度增加 , 其速度近似线性减小 ( 图 2, 样本的拟合曲线为:V=和 和 它与给出的烃有关 , 这种线性变化在许多液体中都可观测到 , 例如 基苯 、 乙酸甲酯 、 乙酸乙酯等 。 在轻质原油中速度有类似关系 。图 2 - 4 正构烷烃 ( a) 和 1 - 烯烃 ( b) 的速度 - 温度关系 ( W N 1987) • 在轻质原油中速度有类似关系。图 2 - 5 34 I 轻质原油在不同压力下的速度 - 温度关系( W an g ) 重烃速度•重烃的速度与温度压力关系与轻烃有些不同 ,图 200 其关系曲线不是线性下降 , 其关系曲线与原油中所含的固体成分有关 , 例如沥青 、 腊等 , 当它们溶化时使得速度降低 。图 2 - 6 100A P I 原油在不同压力处的温度与速度关系( W 1988 ) • 图 2110 220温度关系曲线很明显,随着温度的增加速度明显下降,其主要原因是由于温度增加其可压缩性也增加。图 2 - 7 蒸汽驱油开采的沥青 (1 10 A P I ) 和 220A P I 原油的温度 - 速度关系曲线 含气石油• 图 2度、压力的关系曲线。速度随压力线性减小,当压力达到饱和点(在 72℃ 时,压力 3000磅/英寸 度下降与理论有些差异,可能是泡点太小而影响速度,然后饱和气溢出,再进行测试,速度与压力具有线性关系但是明显偏高,这意味着原油中的气体使速度下降(约百分之十五)。图 2 - 8 两种不同温度下活油的压力 - 速度关系 温度对速度的影响• 1984年 , 某些岩样当温度从25℃ 上升到 150℃ 时速度下降达 40%而这种变化是以前所未能预见到的 , 针对此种情况开展了许多系统性的研究 , 主要结论如下:• ⑴ 完全重油饱和未固结砂岩对温度变化非常敏感,比对压力变化敏感得多。 1987)实验的两种样品,在温度由 25℃ 增大到 150℃时(有效压力保持在 100纵波速度降低 22%~ 40%。图 重 油饱和度 100%)试验结果。由图可见,在30300 10100定有效压力下,随着温度升高纵波速度急剧下降图 2 - 1 1 100 % 重油饱和末固结砂岩在两种有效压力下速度 - 温度关系 ( T a , 1987 ) • 图 00%盐水饱和时情况,除样品孔隙空间里重油变为盐水外,其它样品特征及外部条件都与图 们看到,同样的未固结砂岩,水饱和时的纵波速度几乎与温度没有关系。图 2 - 12 100 % 盐水饱和末固结砂岩在两种有效压力下速度 - 温度关系 ( T a , 1987 ) • ⑵ 含重油未固结砂岩纵波速度随温度升高而降低的幅度与含油饱和有关 。 含油饱和度越高 , 降低幅度越大 。• 图 图 其它试验条件也相同。我们看到,此时纵波速度依然随温度升高而降低,但在同一温度范围内降低幅度仅为 20%,即介于重油饱和与盐水饱和情况之间。图 2 - 13 重油和盐水各 50% 的末固结砂岩在两种有效压力下速度 - 温度关系 ( T a , 1987 ) • 含沥青未固结砂岩实验结果 (图 2实验有效压力 15条曲线对应含 5%和 11%沥青(重量比),另两条曲线分别含沥青(重量比)为 20%,由图 2波速度随温度升高而急剧降低(约降低 30%)。同时,我们看到,沥青含量在 20%之间变化时,对速度降低的幅度影响很小。但当部分含水时,则速度较高而随温度下降轻慢图 2 - 14 含气、水的未固结砂岩及其混合沥青时的温度 - 速度关系 ( et • 含重油或沥青的固结良好砂岩中 , 纵波速度也随温度升高而降低 , 这与未固结砂岩类似 , 但速度降低幅度较小 。温度在 20℃ ~ 120℃ 之间变化时 , 纵波速度降低约 10%~ 15%。• 图 200岩 (孔隙度为 19%)的纵波和横波速度与温度关系。可以看到,随温度升高,纵波速度降低的幅度比横波大许多。图 2 - 15 含 100A P I 重油的 B 岩 ( 孔隙度为 19%) 的纵波和 横波速度与温度关系 • 下表概述了碳氢化合物和含碳氢化合物,水、气岩石的温度对 的来说,在碳氢化合物饱和的岩石中每升高 100℃ 根据不同的碳氢化合物,纵波速度下降 9~ 40%,含重油的固结岩石下降范围约为 9~ 22%。含碳氢化合物的未固结岩石下降 15~ 40%,相对比含水和含气岩石的速度下降相对地要小得多;上述研究成果对监测热采是非常有用的。温度从22℃到122℃时 绝对值(米/秒) 百分比 %纯 烃重 油含气砂含水砂含纯烃砂含重油砂含气砂岩含水砂岩含纯烃砂岩含重油砂岩350- 450270- 80080110280- 320300- 700130- 160200- 260200- 400300- 90020- 3015- 405715- 1715- 303- 65- 76- 89- 22孔隙介质理论当我们把注意力转向多孔流体饱和的岩石时,情况变得相当复杂。自从 1928年起已经出现了近 20种理论。最早的要数1928~1929年的 最早最实用的要数 1951年的 1956年出现了最著名的 、弹性模量2、纵波速度和横波速度 3/42  K 3、(1)/(K ,平面波模量 M , ,M ,E 2 3 43 K 221 )1(    E   )1(2   3 43 M )21(2 E )23(43()(2  )(2 23(2 23     2k 221  石基质、孔隙流体系统、干岩石框架 (骨架 )和饱和岩石本身 。岩石骨架孔隙 /流体干燥岩石 (孔隙空 )饱和岩石 (孔隙充满流体 )该形式下的多孔介质的特征推动了弹性理论的扩展-多孔弹性理论诞生• 多孔弹性理论描述了在正常压力和温度下介质的性态,同时考虑了孔隙中充填的流体的作用。• 在岩石四部分组成中,人们用抽空岩石骨架和湿润岩石两项替换了上面所说的干岩石骨架和饱和岩石两项。• 当孔隙中含有流体时,介质变为 双相介质• 研究双相介质理论就是研究组成骨架及流体的各个成分的物性对岩石整体物性的作用和贡献。• 因为各种理论纷繁复杂,总的说来,大致有三种理论:有效介质理论、接触理论和自适应理论。• 岩石本质上是非均质的,其性质随测量尺度而变化,当研究尺度与介质中的“夹杂”尺度相比大得多时才可以按均匀介质处理。• 。• 为了描述多孔流体饱和的岩石性质,常用有效介质理论,即假设有这样一种单相的介质,其性质与多相介质的宏观平均相同,这种假设的单相介质就称为该多相介质的 “ 有效介质 ” 。• 有效介质理论认为岩石总体的物性参数是由各成分各自的物性参数综合而成的,称为有效物性参数。• 如 与各成分的压缩系数及体积模量之间的关系为:ni 1 i为各成分的体积百分比• 接触理论 研究颗粒物质的有效弹性,它适用于非固结储层,用于估计孔隙度和深度对速度的影响。常见的接触理论有地震勘探中应用的范围有限。• 自适应理论 是对波动方程作了自适应假设后导出的,如 个方程假设岩石是均匀各向同性的,孔隙是全部连通的,孔隙中充满无摩擦的流体,固体 动时流体与固体的相对运动可以忽略不计,波动频率是低频。事实上, 于自适应理论的还有 951)提出骨架的特性可以用某种方法测量,并且他为任何一种性质已知的饱和流体,都推导出了相应岩石的表达式。在这种推导中,他假定流休和固体间的任何相对运动与饱和岩层自身的运动相比可以忽略不计,在低频时,从直观上看是正确的。他指出,骨架中的任何各向异性必然加在饱和的岩层上。为了简单起见,我们讨论的骨架将是由固有的各向同性的弹性固体组成的,就其平均特性来说,骨架是各向同性固体假设流体和固体一起运动,介质密度是固体骨架和孔隙流体密度的按体积加权平均:)1()1(**mf介质总密度孔隙度孔隙流体密度 为岩石骨架 (基质 )密度利用 盐水密度烃类密度岩石干燥时密度假定在任意改变骨架的切变模量时,流体与固体相互不发生影响。则有:* d*d介质总剪切模量, 骨架剪切模量])1([)1(22*介质总体积模量, 基质 (颗粒 )体积模量骨架体积模量, 孔隙流体体积模量*由此可以得到流体饱和介质中纵波和横波的速度分别为*******)3/4K• 据统计,有关岩石物理研究的论文, 1/3以上都涉及 见其重要性。 与孔隙几何形状无关,惯性和散射可忽略不计,外力引起的孔隙压力可在比单个颗粒或孔隙大得多的尺度上实现平衡。(1)观测岩石(骨架和基质)宏观上是均匀分布的 。这个假设是多孔介质中波传播的普遍理论,它保证了波长大于颗粒和孔隙的尺寸。对于大多数的岩石,从实验室频率到地震频率的波一般能符合这个假设。 1975)曾将 ( 2)所有的孔隙是连通的 。这个假设条件意味着岩石高孔隙度和高渗透率,岩石中不存在孤立或连通性差的孔隙。这个假设保证在半个波周期的时间框架内,波传播引发的孔隙流体流动得到充分均衡。对于 假设无限波长(频率为零)的前提下,无论孔隙的连通性如何,大多数的岩石都能满足这个假设。对于测井和实验室用的高频率波,大多数岩石不能够符合这个假设。因此测井和实验室测量的速度通常高于 (3)孔隙中充满无摩擦流体(液体、气体或者混合物)该条件意味着饱和流体的粘度是零。这个假设的目的在于在此确保孔隙流体流动的充分均衡。这个假设也跟频率和波长有关,如果频率为零,任何粘度的流体都将在半波周期内(无限时间)达到均衡。而所有的波都有限定的波长,因此利用 (4)研究的岩石 排出液体)此条件意味着对于实验室岩石样品来说,岩石 岩石表面上没有流体可以流进流出。(5)孔隙流体不会对岩石骨架产生软化或硬化作用。这个假设的目的在于排除岩石基质和孔隙流体之间的任何化学 /物理相互作用的影响。实际上,孔隙流体不可避免地与岩石固体基质发生相互作用而改变表面能量。当岩石为流体所饱和时,流体会削弱或强化岩石基质。• 利用测量的骨架岩石的速度,可以计算出骨架体积模量和剪切模量)34( 22 2 骨架模量与干燥模量是不同的。常为水)条件下测量的。不可还原流体是岩石骨架的一部分,而不属于孔隙空间。实验室中过于干燥的岩样会产生错误的 1分别是水、油、气的体积模量wK oK S 别是水、油、气的饱和度混合流体的体积密度可以由下式计 括通常在实验室测得的干燥骨架体积模量和剪切模量、孔隙度、颗粒密度以及流体体积模量(不可压缩率)。如果没有实验室数据可供使用,这些参数也可通过测井资料或经验关系式估算出来。例如,孔隙度能从中子或声波测井中求得,等等。孔隙流体泥 岩泥 岩1 21 31 51 71 81 91 110含水100%含气100%含油100%含气95%含水5%含水100%含水95%油50%含气50%砂 岩 孔隙度 20%道号含油95%含气5%模型基本参数流体模型层位 1层位 2含水100%含气100%含油100%含气95%含水5%含水100%含水95%油50%含气50%含油95%含气5%多相岩石弹性模量模型• 石骨架应该由单一矿物构成,为一相。孔隙流体应该由气体或液体构成,为另一相。这与地震勘探的实际有差异。作为油气主要储层的砂岩,除主要成岩矿物石英之外,总含有一定量的粘土及其他矿物,岩石骨架本身就是多相体。天然气混有一定量的水份。原油也有一定的挥发成份。孔隙流体有可能是气、油、水构成的多相体。理论研究和实验室测定都证明,流体中只要含有少量气体,就足以使流体的体积压缩模量大大降低。这个影响比较简单。粘土、云母等矿物的影响就比较复杂。这是 理论上,若想计算多相体的的模量,需要知道• (1)每一相体积百分比• (2)每相的模量• (3)各相的几何形状和空间分布、排列情况,但是在一般情况下难以知道。如果我们只有每一相的体积百分比和其模量,那么,我们所能做得至多是估计其等效模量和速度的 上下限给出等效模量的上下限之后,集合体的模量值落在上下限之间,颗粒间“刚性”排列使多相集合体模量趋于上限,而“柔性”排列使集合体模量趋于下限。其确切值则取决于岩石的微观组构。有许多方法可以估算多相体的等效体积模量的上下限,能给出较窄范围的估算方法应是最好的。上限下限0 物质 2体积百分比 100%有效弹性模量边界示意图 已知岩石由 中第 如每种成分的弹性模量 ,剪切模量,杨氏模量 知, 910)和 929)提供了岩石弹性模量的两个估算公式。• 有效弹性模量的 • 此也称为等应变平均。1• 有效弹性模量的 Reus
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:岩石物理
链接地址:http://www.oilwenku.com/p-58688.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开