据岩相转换优势度的沉积相判别与古地貌恢复--以柴达木盆地西部跃进二号油田第三系为例.pdfVIP

据岩相转换优势度的沉积相 判别和古地貌恢复+ ——以柴达木盆地西部跃进二号油田第三系为例 王随继 (中国科学院地理科学与资源研究所，北京100101) E ac．cn mail：wangsj@igsnrr 在碎屑岩沉积相的研究实践中，人们逐渐认识到仅仅依据少数钻井的岩芯资料来判别 沉积相存在着一定的难度，因而强调进行三维沉积剖面的研究[1”]。但是，在缺乏露头剖面 和钻井密度不高的地区，难以查明地下沉积体系的三维分布情况，这就要求采用其他有效途 径，比如二维及三维地球物理勘探、浅层地层的探地雷达勘探等““]。但是如果没有必须的 装备、充足的经费和人员，是难以达到目的的。如果钻井岩芯比较连续、完整，其地下沉积相 垂向上的变化可以利用马尔科夫链方法来判别n]，因为钻井岩芯能够揭示岩性最细微的变 化，且比露头剖面更便于进行岩性垂向转换分析。尽管马尔科夫链分析方法提出很久了，在 各个领域早已得到应用，但在沉积学中迄今所做的马尔科夫链岩相转换分析只是抽取了有 关目标地层的优势岩相转换方向，而对于沉积相的判别没能达到模式化。本文利用马尔科夫 链分析方法对柴达木盆地西部跃进二号油田有关钻井的岩相垂向转换优势度作了分析，并 判别了研究区不同地层的沉积相属性及形成的古地貌背景。 1马尔科夫链分析的基本原理 利用马尔科夫链分析岩相的垂向转换关系时，把岩相由下向上发生转换的实际次数作 为矩阵的元素，由此可以建立一个实测的岩相转换矩阵，记作F。由实测岩性转换矩阵可以 得到3个与之相关的概率矩阵，分别定义为： (1) (2) 肋黝肼 鬻 (3) 式 R。——独立试验矩阵，表示一个给定的岩相转换矩阵随机发生的概率； 尸。——实际转换概率矩阵，表示实测矩阵中岩相实际转换的概率； D，，——差值转换矩阵，为岩相实际转换概率与随机转换概率之差； S，——表示实际转换矩阵中各行元素的和； S，一表示实际转换矩阵中各列元素的和； T一11表示实际转换矩阵中所有元素的总和； z，J——分别表示实际转换矩阵的行号和列号。 矩阵(D，，一P。，一尺。)就是所期望的结果，其中大于零的数值表明它所代表的岩相转换占 ；中国科学院地理科学与资源研究所知识刨新工程资助项目(CXlOC-A0202) ·165· 优势，这些优势转换方向的组台基本能够反映沉积相特征。 2岩相转换的马尔科夫链模型 跃进2号油田位于柴达木盆地西部，下渐新统和中新统～下上新统砂岩体是其主要储 集层，分别称做E。1油藏和N，～N：1油藏。现以其中的342号钻井为例来简述其马尔科夫链 模型的建立过程和结果。参照野外粒度计把岩芯中碎屑岩依据粒度分为泥岩、粉砂岩、极细 砂岩、细砂岩、中砂岩、粗砂岩、极粗砂岩和砾岩8类岩相。根据该钻井岩芯中不同岩相由下 向上实际发生转换的次数进行统计01，得到表l所示的岩相实测转换的频率矩阵。 表1 342号钻井岩芯中岩相实际转换的频率矩阵 岩性 泥岩 粉砂岩极细砂岩细砂岩中砂岩粗砂岩极粗砂岩砾岩 = 泥岩(Y) o 56 4 8 0 0 o 1 69 粉砂岩(s) 55 0 8 8 2 2 1 2 78 极细砂岩(Vf)7 6 0 2