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

三维地质建模技术的作用

关 键 词:
三维 地质 建模 技术 作用
资源描述:
三维地质建模技术在开发阶段油藏描述中的地位与作用吴键中国石油勘探开发研究院件使用 "试用版本创建 三维地质建模(3D 就是运用计算机技术,在三维环境下,将空间信息管理、地质解译、空间分析和预测、地学统计、实体内容分析以及图形可视化等工具结合起来,并用于地质分析的技术-拿大) 1993。• 三维地质建模技术就是综合地质学、地球物理学、地质统计学、计算机技术对地质体进行综合研究的一种技术方法。三维地质建模的基本概念件使用 "试用版本创建 根据油田勘探和开发的不同阶段、资料的丰富程度以及研究的任务的不同把储层地质模型分为概念模型、静态模型和预测模型-裘怿楠• 地质模型大致包括构造模型、相(沉积相、岩相等)模型、储层物性参数模型、流体参数模型等多种模型。• 三维地质建模以井资料为基础,通过地质统计学方法进行井间参数的预测和插值。三维地质建模的基本概念件使用 "试用版本创建 地质统计学是三维地质建模的主要基础,是法国巴黎国立高等矿业学院马特隆教授(G.1962年创立的。他将传统统计学理论与区域化变量的概念相结合,发展出一套以变差函数为工具研究矿产特征区域分布的数学技术。• 目前在地质统计学上,多点地质统计是最新的研究领域。• 自上世纪80年代后期随着计算机技术的飞速发展逐步开发出应用于三维地质建模的计算机软件。美国前已经有多种大型建模软件被广泛的应用于三维地质建模工作,其中比较某种的有 目前应用最为广泛,发展前景也最为良好的是件使用 "试用版本创建 给领导汇报出漂亮图!• 好看不好用!• 给数模搭架子!• 地质建模就是计算个模型!• …………………件使用 "试用版本创建 三维地质建模经常被称为开发阶段油藏描述的核心,但目前的建模成果主要还是集中在研究成果的集成和显示方面,对其它相关的研究工作并没有起到指导性的作用。• 研究水平的整体提高离不开新技术的应用,三维地质模技术是开发阶段油藏地质研究中最为现实、最具可操作性和系统性的新技术,理应在开发阶段油藏描述工作中起到核心作用。• 油藏描述水平的提高需要在具体的技术操作层面上提出可实现的的技术方法。件使用 "试用版本创建 近二、三十年,勘探阶段的油藏研究通过地震数据的充分应用取得了突飞猛进的进展,相比之下开发阶段油藏描述的进步并不明显。• 三维地质建模技术可以在开发阶段地质研究工作中起到地震勘探在勘探阶段的同等作用。• 提高三维地质建模技术的应用水平就是提高开发阶段油藏精细描述的整体水平。• 三维地质建模技术应用水平的提高依赖于彻底改变建模工作基本思路,将其从简单的“模型计算”提升为油藏地质研究的基础平台。件使用 "试用版本创建 充分利用三维地质建模技术特有的技术优势,将三维地质建模工作转变为综合地质研究的一种工具和手段,将建模工作延伸到基础地质研究和储层综合评价中,可以解决常规方法无法解释的问题;• 将三维地质模型转变为三维地质数据体,以其为基础开展更为全面、丰富的综合研究,体现出建模工作的特殊价值。件使用 "试用版本创建 在中国的许多沉积盆地,由于构造和沉积都很复杂,地层横向变化大,电测曲线可对比性差,缺少足够的对比标志地层,对比工作经常存在一定的误差,尤其是小层和断点的对比;• 实际工作证明,这种偏差很难通过单纯的地层对比本身来解决。件使用 "试用版本创建 地层对比是地质研究工作中最基础的工作,往往也是最难做准的工作。• 地层对比的准确性和可靠性直接影响了油藏研究的可靠性。当断层模型与断点数据吻合时,储层得到很好的描述当断层模型与断点数据不能吻合时,储层的描述将会是完全错误的。油层断点断层井位油层断点断层井位井 井件使用 "试用版本创建 传统的地层对比往往难以及时发现在构造形态和深度上与邻井的差异。• 三维可视化技术十分清楚的显示地层对比中的问题。• 通过重新进行对比可以使问题得到及时、准确的校正。小层对比方案调整:地质模型三维可视化质量控制地层对比方案调整前 地层对比方案调整后齐10件使用 "试用版本创建 在局部区域内,地层的构造形态虽然会存在一定的起伏,但幅度有限;• 在连续沉积的条件下构造形态具有一定的继承性。如果相邻井间存在过大的深度差,往往代表了分层方案存在误差;• 在三维可视化的质量控制下,对分层方案进行重新对比和调整,使对比方案和小层的构造形态得到合理的描述,同时合理的微构造得到保留。小层对比方案调整:构造趋势分析与地层对比结合根据构造趋势分析调整地层分层方案实例图调整前莲花油层Ⅱ23小层顶面三维形态 调整后莲花油层Ⅱ23小层顶面三维形态件使用 "试用版本创建 通过将地震解释、地层对比、三维地质建模技术相结合,为断层精细解释、钻井断点位置对比等问题提供了新的解决方案。齐40区块三维地震640区块断层三维展布图三维地质建模在断层落实中的应用件使用 "试用版本创建 所有过井断层都要在钻井上对比有相应的断点• 所有钻井上对比出的断点都要在地震上有所反映,否则就属于不落实• 所有断点要分布在一个光滑的三维面上,并且要与地震资料相一致。小层对比方案调整:构造模型与钻井断点的对比5整后齐5件使用 "试用版本创建 通过地层对比与断层模型的交互校正,可以使二者的精度都得到大幅的提高27312774?? 27312731277429断层 件使用 "试用版本创建 造继承性、微构造分布等多种因素,因此对储层的构造特征有更准确的描述件使用 "试用版本创建 造趋势的分析,可以在无井控制的层内准确描述微构造件使用 "试用版本创建 三维可视化断层解释质量控制。• 断层面在地震资料分辨率的尺度内为具有一定光滑度、连续性的三维曲面。• 利用三维地质建模软件的三维可视化技术可以对断层的地震解释质量进行直观的分析。件使用 "试用版本创建 三维构造模型相当于在三维空间内进行了一次精细构造落实,所得到的构造图更为合理。• 由于构造模型是构造形态的三维描述,因此可以直接描述逆断层的形成。件使用 "试用版本创建 在复杂断块区,由于断层的准确位置难以确定,开发井往往不敢过于接近断层。• 以高精度的构造模型为基础,可以更充分的开发断层高部位的剩余油。010203040日产液油(t)0204060801007- 26 8- 16 8- 23 8- 30 9)件使用 "试用版本创建 地震反演技术通常被当做储层预测的重要手段;• 当地震资料分辨率较低时,基于模型的反演往往会存在假象,而忠实于地震资料的反演又往往达不到地质研究的要求。件使用 "试用版本创建 基于模型的地震反演需要利用测井资料建立初始模型,而目前通用的反演软件在计算初始模型时都是简单的数学插值,并不考虑原始地震数据的分布。• 三维地质模型的模拟计算过程中可以通过协克里金方法将地震数据做为第二变量对插值进行趋势控制。• 以基于地震的地震反演成果为基础计算出的地质模型可以更好的反映出含油单砂体的分布。件使用 "试用版本创建 以模型为基础提取的砂体与钻井资料更为接近。• 在地质模型的基础上可以、快捷、直接、直观的定义出油藏的分布。基于地震反演提取的砂体 基于地质模型提取的砂体件使用 "试用版本创建 测井曲线往往反映不出单砂体内部的相带变化。• 渤海湾盆地下第三系平均沉积速率为3左右,10 因此多数情况下,测井曲线反映的是沉积体长期叠合的产物。件使用 "试用版本创建 测井曲线反映的单砂体往往不是单一的相带。• 沉积微相图反映的往往是沉积体长期演变的综合结果,无法反映出河道的变迁和摆动。2钻井 砂岩等值线微相分界线水流方向图 例244224422222砂河漫滩砂泛滥平原堤岸砂泛滥平 原 河漫滩砂河道河道砂河道砂泛滥 平原堤 岸砂漫滩砂件使用 "试用版本创建 三维沉积微相模型,可以更细致的反映出河道单砂体的空间形态和发育史,是二维沉积微相研究的进一步扩展和延伸。• 以三维地质模型为基础研究剩余油分布和开发方案调整较常规方法更为细致,更为有效。a bc 件使用 "试用版本创建 三维沉积微相模型应该反映出砂体的三维形态、河道的摆动与变迁。• 二维相图的照搬将失去“三维”的意义和作用。件使用 "试用版本创建 火成岩体是一种特殊的岩性体,通常具有孔隙、裂缝二种储集空间。• 结合地震解释结果建立体控构造模型可以直接、准确的描述火成岩体的空间形态。件使用 "试用版本创建 目前针对裂缝型储层的研究方法主要有应力场构造面曲率分析和地震属性研究二个主要方向。• 地层的曲率与地层的裂缝孔隙度为正相关;单层厚度与裂缝密度为负相关。• 在碳酸盐岩纯裂缝储层内,波阻抗、地震振幅往往与裂缝发育程度相关,但在孔隙、裂缝双重介质储层中情况更为复杂。• 方差体可以反映大尺度的裂缝,地层倾角、倾角偏差体可以反映小尺度的裂缝发育情况。地层倾角体 地层倾角偏差体相干体件使用 "试用版本创建 目前裂缝型三维地质建模工作主要有两类方法,一类是直接利用测井、地震资料建立等效的裂缝孔隙度、裂缝渗透率模型,另一类是建立离散裂缝网络模型。• 直接建立等效的裂缝孔隙度、渗透率模型可以保证与测井资料的完全吻合,但缺少对裂缝的直接描述。• 离散裂缝网络模型(近几年新发展起来的一项技术,可以对裂缝的形态和分布进行直接的描述,然后再通过渗流力学理论,从离散裂缝网络模型直接计算出等效的裂缝孔隙度、裂缝渗透率等模型,其主要缺点是计算出的孔隙度、渗透率等储层物性参数与测井解释直接获得的结果不完全一致。• 二种不同的方法适用于不同的地质情况、基础资料条件及不同的研究目的。件使用 "试用版本创建 合理选择地震属性参数是震控建模的关键。• 相对波阻抗较好的反映了基质孔隙度的总体变化,但与裂缝孔隙度没有明显的相关性。这说明裂缝即使对地震波有所影响,也会被基质孔隙的影响所“淹没”。• 方差、倾角、倾角偏差三个构造参数与裂缝孔隙度都有部分的相关性,可以利用神经网络技术将多种参数综合到一起。• 在砂岩、碳酸盐岩等层状岩石中裂缝的密度与单层厚度成反比,与岩层的曲率关系不明显,但在成层差、横向变化大的厚层火成岩内裂缝密度曲线与地震地层倾角体和倾角偏差之间具有较好的相关性,反映裂缝密度的发育同样受到地层的构造变形的控制。件使用 "试用版本创建 以地震数据做条件约束,利用协克里金方法对模拟计算进行条件约束可以建立起比较合理、可靠的模型系统模型神经网络孔隙度体裂缝孔隙度模型随机模型+协克里金件使用 "试用版本创建 在三维地质模型的基础上可以开展进一步的储层分布研究工作。件使用 "试用版本创建 在致密砂岩储层地区,可以利用地层倾角体研究压裂缝的分布,并进一步研究水窜、气窜的方向。件使用 "试用版本创建 充分利用三维地质建模技术特有的技术优势,将三维地质建模工作转变为综合地质研究的一种工具和手段,将建模工作延伸到基础地质研究和储层综合评价中,可以解决常规方法无法解释的问题;• 将三维地质模型转变为三维地质数据体,以其为基础开展更为全面、丰富的综合研究,体现出建模工作的特殊价值。件使用 "试用版本创建 高精度的三维地质模型只是一个中间成果,模型计算完成后还可以其为基础开展进一步的研究工作。• 在三维地质模型的基础上可以快速的编制大量的基础地质图件,极大的提高地质研究的工作效率。提供丰富的基础地质图件件使用 "试用版本创建 三维地质模型是建立在测井数据基础上的三维插值计算,并进行了地质统计学分析。• 与传统二维编图相比,数据空间分布更为合理,对砂体分布的描述更为细致。9100 80 100 120 140 160 180 200 220 240 260 280 300 320 34040 60 80 100 120 140 160 180 200 220 240 260 280 300 320 340映出砂体纵向上的叠置关系件使用 "试用版本创建 从三维地质模型中可以方便的提取出各种剖面图,可以快捷的进行地质分析。• 尤其是过水平井的剖面图,是常规研究方法难以实现的,为水平井的设计和效果评价提供了方便。件使用 "试用版本创建 以三维地质模型为基础,可以进一步计算出储层非均性研究的成果图件i∑=( )i 2−=异系数求取公式:渗透率均值:渗透率均方差值:渗透率变异系数:• 目前的油藏描述研究主要以含油单砂体为主要目标。• 在油藏开发后期,含油单砂体内部的储层非均质性对开发效果有很大的影响。• 单砂体内部的储层非均质性研究目前还是空白,主要难点是工作量巨大。• 利用三维地质建模技术研究储层内部的非均质性可以极大的提高研究工作的效率,整体提高剩余油分布研究的水平和精度研究含油单砂体内部储层非均质性件使用 "试用版本创建 三维地质模型相当于一个三维地质数据体,可以利用软件的过滤、提取等功能从中方便的提取出油层内部高渗区或泥质夹层,开展进一步的研究。渗透率模型高渗带三维分布火成岩体基质孔隙度高值区分布图火成岩体裂缝密度高值区分布图件使用 "试用版本创建 储层内部的高渗条带往往会引起注水、注气的突进。• 通过建立三维地质模型,可以方便的研究高渗条件的分布。1190e v e l o p e r s ' l i c e n s 503320050334005033600503380050340005034200503440050332005033400503360050338005034000503420050344000 50 100件使用 "试用版本创建 00412 200502 200504 200506 200508 200510 200512 200602 200604日产液(油)吨0102030405060708090100含水%日产液日产油含水高104)件使用 "试用版本创建 将地质建模与油藏工程相结合,可以进一步优化油藏开发方案。从变化找地质依据 开发方案优化静态参数模型初始地质模型油水分布规律油藏数值模拟动态历史拟合模型耦合实现地质油藏一体化修改后(局部网格渗透率2倍)修改前(初始模型)件使用 "试用版本创建 利用不同时间的测井数据建立研究开发不同阶段油藏的变化。• 对测井数据要求较高。不同含水期渗透率模型件使用 "试用版本创建 地质建模基本流程• 地质统计学的几个基本概念• 模型验证方法• 相控&震控• 模型的地质解释件使用 "试用版本创建 地质建模的基本流程为:基础数据准备-构造模型-属性模型-成果输出;• 不存在一个完全固定的工作流程,要根据地质要求和基础数据情况进行调整。• 合理的工作流程要保障建模成果是真正的地质模型,而不是数学计算模型或者成果展示模型。• 要注重地质研究与建模工作的交互:构造建模与地层对比和构造落实的交互;沉积相研究和油藏预测与属性建模的交互;建模成果与油藏研究的交互。件使用 "试用版本创建 确定性建模(:对井间未知区给出确定性的预测结果。• 随机建模(应用随机模拟方法,对井间未知区给出多种可能的预测结果。确定性建模 随机建模件使用 "试用版本创建 实性建模):• 即根据待估点周围的若干已知信息,应用变差函数的性质,对估点的未知值作出最优(估计方差最小)、无偏(估计值的均值与观测值的均值相等)的估计。• 克里金插值法为光滑内插方法,为减小估计方差而对真实观测数据的离散性进行了平滑处理,虽然可以得到由于光滑而更美观的等值线图或三维图,但一些有意义的异常带也可能被光滑作用而“光滑”掉了。• 克里金插值为局部估计方法,对估计值的整体空间相关性考虑不够,它保证了数据的估计局部最优,却不能保证数据的总体最优,因为克里金估值值的方差比原始数据的方差要小。因此,当井点较少且分布不均时可能会出现较大的估计误差,特别是在井点之外的无井区误差可能更大。件使用 "试用版本创建 随机模拟是以随机函数理论为基础的。随机函数由一个区域化变量的分布函数和协方差函数(或变差函数)来表征。 随机模拟的基本思想是从一个随机函数Z(u)中抽取多个可能的实现,即人工合成反映Z(u)空间分布的可供选择的、等概率的高分辨率实现,记为{Z(l)(u),l=1,…,L, 代表变量Z(u)在非均质场 模拟方法首先考虑的是模拟值的全局空间相关性,其次才是局部估计值的精确程度。• 随机模拟结果在在光滑趋势上加上系统的“随机噪音”,这一“随机噪音” 能更好地表现真实曲线的波动情况。• 随机模拟产生许多可选的模型,各种模型之间的差别正是空间不确定性的反映。件使用 "试用版本创建 在三维地质建模工作中,随机地质统计学占有重要的位置。• 所谓随机模拟就是指以已知的信息为基础,人工合成可选的、等概率的和高精度的反映变量Z(u)空间分布的模型。其中,可选是指所合成的模型是多个的,每个模型都是对原始数据的反映;• 等概率是指模拟参数的统计特征与现有样品的统计特征或参数的理论分布是一致的;• 高精度是指所合成的模型能够反映参数的细微变化。各个模型之间的差别就反映了由于资料不足而引起的建模中的不确定性。• 随机建模就是通过提供多个可选的、等概率的和高精度的储层结构或属性的空间分布图象,比较各个图象上的差异性和变化性来恢复变量的细微变化。模型 种类 方法 算法示性点过程法马尔科夫随机域法条件模拟 截断高斯法离散性 基于目标 两点直方图法模型 指示模拟法模拟退火法非条件模拟 布尔模拟模拟退火法多变量联合条件模拟序贯高斯模拟连续性 基于象元 条件模拟 序贯指示模拟模型 分形随机函数法马尔科夫随机域法向带法模拟技术方法简表件使用 "试用版本创建 按照变量的特征,将随机模拟分为离散性和连续性随机模拟• 从模拟实现的过程看,可以分为基于目标(模拟和基于象元(模拟• 从变量类型角度看,各种基于指示克里金的模拟(序贯高斯模拟、基于率场模拟)、截断高斯模拟、布尔模拟适用于离散型变量(模拟;• 序贯高斯模拟、序贯高斯协同模拟、分形随机函数法、矩阵方法、频率域方法)适用于连续型变量(模拟;• 模拟退火模拟(包括浴法、协同模拟退火法、并行模拟退火法)既适用于离散类型变量,也适用于连续型变量。离散和连续要构造目标函数能综合各种信息计算量大,不易收敛随机模拟方法变量类型适用条件 优点 缺点布尔类型离散型 可以重复而易描述的形状原理简单,计算量小很难忠实于具体位置的信息,不能反映砂体内部非均质性高斯类型连续和离散变量必须是正态或多元正态分布计 算速度快,数学上具有一致性很难考虑间接信息,要求变量服从正态分布。指示类型离散和连续没有具体的要求能综合各种信息,是最灵活的随机模拟方法计算量大,需要推断很多协方差函数。不能忠实于试井资料模拟退火类型分形类型连续性变量变量具有分形特征快速和经验性强难考虑间接信息离散性变量 连续性变量件使用 "试用版本创建 协克里金是一种多变量估计技术,通过研究主变量及次级变量的空间相关关系,将次级变量的信息整合到估计结果中,以弥补主变量数据不足的缺点。• 在简单克里金回归方法中,未知值是通过 克里金面u)可被看作是某一拟合()维的u)与Z1(u+h)在距离hα== 1, … 协克里金方法是该种回归方法的拓展,它包含不同于如,如果除了,有u’α)在未采样点处,直接变量的协克里金估计为:• 克里金与协克里金的区别在于前者为单个协变差函数,而后者为四个:h) = 1(u), Z1(u+h)﹜h) = 2(u), Z2(u+h)﹜h) = 1(u), Z2(u+h)﹜h) = 2(u), Z1(u+h)﹜Z*(u)=Z(11∑λ u)ααα=1Z*(u)=Z(11∑λ u)ααα=1α=1件使用 "试用版本创建 用的是将协克里金法结合到序贯高斯模拟计算中。具体步骤为:1. 将原始数据做正态变换;2. 确定一个随机路径,访问所有的待模拟的节点;3. 在每个节点处保留一定数量的邻域条件数据(包括原始的转换数据和先前模拟的网格节点值),应用协克里金以及变差函数模型整合次级变量数据,建立高斯条件累积分布函数(4. 从将该值加载到已有的数据值;5. 返回步骤3,直到访问完所有的点;6. 将数据变换到原始空间;7. 重新确定一条路径进行另一实现。件使用 "试用版本创建 ) ( ) ()( )[ ]221 +=x)在x, x+:件使用 "试用版本创建 选择模型计算参数时要坚持地质合理性优先于数学合理性的原则。• 已知点间距(井距)对变差函数有明显的影响。• 经验值:平面变程3向变程1;块金值0。变程:600m 变程:800m(统计变程) 变程:1200算结果保持稳定纵向变程大于1,可能引起砂体窜层,与地质认识不符块金值不为零造成砂体分布分散,还容易造成油藏数值模型计算的不收敛。件使用 "试用版本创建 多点地质统计学是相对于传统的两点地质统计学而言的;• 传统的地质统计学以变差函数做为研究地质变量空间相关性的重要工具。然而,变差函数只能把握空间上两点之间的相关性,亦即在二阶平稳或本征假设的前提下空间上任意两点之间的相关性,因而难于表征复杂的空间结构和再现复杂目标的几何形态(如弯曲河道)。• 多点地质统计学中,应用“训练图像”代替变差函数表达地质变量的空间结构性。• 多点地质统计学随机模拟方法(如传统的地质统计学随机模拟方法(如序贯指示模拟本质差别在于未取样点处条件概率分布函数的求取方法不同。前者应用多点数据样板扫描训练图像以构建搜索树并从搜索树中求取条件概率分布函数,而后者通过变差函数分析并应用克里金方法求取条件概率分布函数。• 正是这一差别,使多点地质统计学克服了传统二点统计学难于表达复杂空间结构性和再现目标几何形态的不足。件使用 "试用版本创建 虑一种属性S(如沉积相),可取K 个状态(如不同相类型),即{ Sk,k=1,2,⋯K},则一个以小为n 的“数据事件”以下两部分组成:1)由个向量{ha,a=1,2,⋯ }确定的几何形态(数据构形),亦称为数据样板(记为∏n;2) 一个五点构形的数据事件,由一个中心点和四个向量及数值组成。多点统计则可表述为一个数据事件 {S (=,...n}出现的概率,即数据事件中S(分别处于态时的概率,也可表述为:• 在实际建模过程中,上述多点统计或概率难于通过稀疏的井资料来获取,而需要借助于训练图像。• 训练图像为能够表述实际储层结构、几何形态及其分布模式的数字化图像。对于沉积相建模而言,训练图像相当于定量的相模式,它不必忠实于实际储层内的井信息,而只反映一种先验的地质概念件使用 "试用版本创建 001) 应用一种动态数据结构即“搜索树”一次性存储训练图像的条件概率分布,提出了多点统计随机模拟的)建立训练图像。2)准备建模数据,将实测的井数据标注在最近的网格节点上。3)应用自定义的与数据搜索邻域相联系的数据样板δ构建搜索树。4)确定一个访问未取样节点的随机路径。在每一个未取样点得条件数据置于一个以n’表示条件数据的个数,条件数据事件。从搜索树中检索C(和Ck()并求取)从模拟值加入到原来的条件数据集中,作为后续模拟的条件数据。6)沿随机路径访问下一个节点,并重复3)、4)步骤。如此循环下去,直到所有节点都被模拟到为止,从而产生一个随机模拟实现。7)改变随机路径,产生另一随机模拟实现。件使用 "试用版本创建 究区训练图像 手工编制的研究区沉积微相图研究区多点统计随机建摸的一个实现 研究区基于变差函数的序贯指示模拟的一个实现件使用 "试用版本创建 地下储层本身是确定的,但是在现有数据不完善的条件下人们对它的认识存在一些不确定的因素,因此地质建模的结果具有不确定性。• 随机模拟计算的特点之一是计算结果不是单一的模型,而是一系列的实现(即针对同一地区,应用同一套资料、同一建模方法可以计算出多个储层地质模型。• 随机模拟的多次实现常被用来评价地质建模的不确定性。多次实现与模型不确定性评价岩相模型的四个实现,平均井距小于100件使用 "试
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

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

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