基于井的碎屑岩储层成岩相横向预测方法

摘要

本发明公开了一种基于井的碎屑岩储层成岩相横向预测方法，该方法以沉积相、岩性、酸碱度、沉积旋回、成岩阶段、断层指数六项地质参数为碎屑岩储层成岩相类型边界约束条件，并结合井上成岩相数据，对成岩相横向展布进行预测，填补了现有技术在成岩相横向预测研究方面工作的不足。该方法确定储层有效性，为油藏评价提供依据。

说明书

技术领域

本发明涉及油气勘探开发技术领域，具体地指一种基于井的碎屑岩储层成岩 相横向预测方法。

背景技术

成岩相是指在成岩环境控制下，由各种成岩作用形成并具有一定几何形态和 特定成岩组构、特定成岩矿物组合的地质体。与沉积相、测井相中相的概念类似， 成岩相所反映的意义是构造、流体、温度、压力等多个因素共同作用下一个储集 层最终状态的集合，是成岩环境下岩石学特征、地球化学特征和岩石物理特征的 总和，其核心内容是现今的矿物成分和组构面貌，主要是表征储集体性质、类型 和优劣的成因性标志，可借以研究储集体形成机理、空间分布与定量评价。预测 有利孔、渗性成岩相是储集层研究和油气勘探的重点。

基于测井资料或岩心数据的成岩相横向预测，只能根据成岩相在二维或三维 空间的变化趋势，大体勾绘出成岩相的横向展布，这样的预测带有很强随意性。 地震资料具有较好的平面连续性，因而对储层成岩相平面上的特征可具有较连续 的反映，然而限于地震资料分辨率影响，在薄层碎屑岩储层以及深层碎屑岩储层 识别上，地震资料具有很大局限性，所以传统储层成岩相预测已不能较好地反映 出储层的发育，影响了薄层碎屑岩储层尤其是深层-超深层储层油气勘探开发。

发明内容

本发明针对现有技术存在的问题，提供了一种基于井的碎屑岩储层成岩相横 向预测方法，该方法以沉积相、岩性、酸碱度、沉积旋回、成岩阶段、断层指数 六项地质参数为碎屑岩储层成岩相类型边界约束条件，并结合井上成岩相数据， 对成岩相横向展布进行预测，填补了现有技术在成岩相横向预测研究这方面工作 的不足。本发明提供了一种碎屑岩储层成岩相横向预测新方法，进而确定储层有 效性，为油藏评价提供依据。

为了实现上述目的，本发明提供的一种基于井的碎屑岩储层成岩相横向预测 方法，包括以下步骤：

1、收集研究区地质资料

2、建立原始资料数据库

原始资料数据库包括碎屑岩沉积相类型数据库F_m、岩性数据库R_n、酸碱 度数据库P_o、沉积旋回数据库C_p、成岩阶段数据库S_q和断层指数数据库 N_r；

3、建立成岩相划分方案

根据研究区孔隙类型、成岩特征以及面孔率，建立成岩相类型D_s划分方 案如下：

实相4、建立成岩相预测模型蚀、弱胶结

(1)储层网格化

储层网格化是将空间上不均匀分布的数据按一定方法(如滑动平均法、克里 格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程； 研究区储层的每个网格可用Wi(X，Y)表示；

(2)确定网格属性

a.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；

b.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；

c.根据研究区酸碱度数据确定网格Wi(X，Y)的酸碱度属性P_o；

d.根据研究区沉积旋回数据确定网格Wi(X，Y)的旋回属性C_p；

e.根据研究区成岩阶段数据确定网格Wi(X，Y)的成岩阶段属性S_q；

f.根据研究区断层指数数据确定网格Wi(X，Y)的断层指数属性N_r；

g.根据上述步骤a、b、c、d、e和f确定每个网格Wi(X，Y)的综合属性 G(F_m，R_n，P_o，C_p，S_q，N_r)，即：

Wi(X，Y)＝F_m+R_n+P_o+C_p+S_q+N_r；

(3)建立成岩相预测模型

a.基于研究区岩石薄片资料、分析测试资料，确定研究区已知井Hj所在网 格Wj(X，Y)的成岩相D_s_j；

b.确定已知井Hj所在网格Wj(X，Y)的综合属性Gj(F_m_j，O_n_j， S_o_j，C_p_j，P_q_j，N_r_j)，建立已知井Hj所在网格Wj(X，Y)成岩相D_s 与综合属性G(F_m_j，O_n_j，S_o_j，C_p_j，P_q_j，N_r_j)的对应关系，即：

D_s_j＝F_m_j+O_n_j+S_o_j+C_p_j+P_q_j+N_r_j；

(4)成岩相预测

a.任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)的综合属性G (F_m，R_n，P_o，C_p，S_q，N_r)；

Gi＝(F_m_i，R_n_i，P_o_i，C_p_i，S_q_i，N_r_i)

b.将网格Wi(X，Y)的成岩相综合属性Gi与井Hj所在网格Wj(X，Y) 的综合属性Gj相比较，即：

Gi-Gj＝(F_m_i，R_n_i，P_o_i，C_p_i，S_q_i，N_r_i)-(F_m_j+ R_n_j+P_o_j+C_p_j+S_q_j+N_r_j)；

如果满足Gi-Gj＝0，则未知网格Wi(X，Y)的成岩相D_s_i与井j具有相 同的成岩相D_s_j；

如果Gi-Gj≠0，则按不同属性优先级顺序，

即：沉积相F_m一级、岩性R_n二级、沉积旋回C_p三级、酸碱度P_o四 级、成岩阶段S_q五级、断层指数N_r六级对未知网格Wi(X，Y)的成岩相进 行判识，即：

ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，C_p_i-C_p_j≠0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j≠0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅲ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅳ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0， S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅴ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0， S_q_i-S_q_j＝0，N_r_i-N_r_j≠0；

未知网格Wi(X，Y)的成岩相D_s_i满足条件i时，未知网格Wi(X，Y) 与具有相同属性F_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j 的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅲ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅳ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j、P_q_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅴ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j、P_o_j、S_q_j的网格Wj(X，Y)的成岩相相同；

c.根据每个网格Wi(X，Y)的成岩相确定研究区横向分布。

进一步地，所述步骤1)中，收集的地质资料包括以下内容：

①研究区目标层沉积相研究成果图；

②研究区物源研究成果图；

③研究区断裂分布研究成果图；

④研究区流体史(酸碱度)资料；

⑤研究区目标层成岩阶段研究成果图；

⑥研究区目标层沉积旋回研究成果图；

⑦研究区岩石薄片分析测试数据。

再进一步地，所述步骤2)中，

(1)沉积相类型数据库F_m

根据前人所做的研究，对碎屑岩储层沉积相类型进行概括和分类，其中包括： 冲积扇-旱扇-扇根(F_1)、冲积扇-旱扇-扇中(F_2)、冲积扇-旱扇-扇缘(F_3)、 冲积扇-湿扇-扇根(F_4)、冲积扇-湿扇-扇中(F_5)、冲积扇-湿扇-扇缘(F_6)、 河流相-曲流河-河床亚相-河床滞留(F_7)、河流相-曲流河-河床亚相-边滩(F_8)、 河流相-曲流河-堤岸亚相-天然堤(F_9)、河流相-曲流河-堤岸亚相-决口扇(F_10)、 河流相-曲流河-河漫亚相-河漫滩(F_11)、河流相-曲流河-泛滥盆地-河漫湖泊 (F_12)、河流相-曲流河-河漫沼泽(F_13)、河流相-辫状河-牛轭湖(F_14)、河 流相-辫状河-河床滞留(F_15)、河流相-辫状河-心滩(F_16)、河道(F_17)、河 流相-辫状河-泛滥平原(F_18)、湖泊相-断陷型-湖成三角洲(F_19)、湖泊相- 断陷型-滨湖(F_20)、湖泊相-断陷型-浅湖(F_21)、湖泊相-断陷型-半深湖(F_22)、 湖泊相-断陷型-深湖(F_23)、湖泊相-断陷型-湖湾(F_24)、湖泊相-坳陷型-湖 成三角洲(F_25)、湖泊相-坳陷型-滨湖(F_26)、湖泊相-坳陷型-浅湖(F_27)、 湖泊相-坳陷型-半深湖(F_28)、湖泊相-坳陷型-深湖(F_29)、湖泊相-坳陷型- 湖湾(F_30)、湖泊相-前陆型-湖成三角洲(F_31)、湖泊相-前陆型-滨湖(F_32)、 湖泊相-前陆型-浅湖(F_33)、湖泊相-前陆型-半深湖(F_34)、湖泊相-前陆型- 深湖(F_35)、湖泊相-前陆型-湖湾(F_36)、三角洲相-辫状河三角洲-三角洲平 原-分支(F_37)、三角洲相-辫状河三角洲-三角洲平原-河道(F_38)、三角洲相- 辫状河三角洲-三角洲平原-陆上天然堤(F_39)、三角洲相-辫状河三角洲-三角洲 平原-决口扇(F_40)、三角洲相-辫状河三角洲-三角洲平原-沼泽(F_41)、三角 洲相-辫状河三角洲-三角洲平原-淡水湖泊(F_42)、三角洲相-辫状河三角洲-三 角洲前缘-水下分支河道(F_43)、三角洲相-辫状河三角洲-三角洲前缘-水下天然 堤(F_44)、三角洲相-辫状河三角洲-三角洲前缘-支流间湾(F_45)、三角洲相- 辫状河三角洲-三角洲前缘-分支河口砂坝(F_46)、三角洲相-辫状河三角洲-三角 洲前缘-远砂坝(F_47)、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂(F_48)、 三角洲相-辫状河三角洲-前三角洲(F_49)、三角洲相-扇三角洲-扇三角洲平原- 泥石流(F_50)、三角洲相-扇三角洲-扇三角洲平原-河道充填(F_51)、三角洲 相-扇三角洲-扇三角洲平原-漫滩(F_52)、三角洲相-扇三角洲-扇三角洲前缘-碎 屑流(F_53)、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道(F_54)、三角 洲相-扇三角洲-扇三角洲前缘-支流间湾(F_55)、三角洲相-扇三角洲-扇三角洲 前缘-河口砂坝(F_56)、三角洲相-扇三角洲-扇三角洲前缘-远砂坝(F_57)、三 角洲相-扇三角洲-前扇三角洲(F_58)。

(2)岩性数据库R_n

根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩 性数据库。其中包括：细砾岩R_1、中砾岩R_2、粗砾岩R_3、巨砾岩R_4、石 英砂岩R_5、长石质石英砂岩R_6、岩屑质石英砂岩R_7、长石岩屑质石英砂岩 R_8、长石砂岩R_9、岩屑质长石砂岩R_10、岩屑长石砂岩R_11、岩屑砂岩R_12、 长石质岩屑砂岩R_13、长石岩屑砂岩R_14、黏土R_15、泥岩R_16、页岩R_17。

(3)酸碱度数据库P_o

根据储层孔隙中流体酸碱性大小，建立酸碱度数据库；其中包括：强酸性 P_1：pH＜3；酸性P_2：3≤pH＜5；弱酸性P_3：5≤pH＜7；弱碱性P_4：7≤ pH＜9；碱性P_5：9≤pH＜11；强碱性P_6：11≤pH≤14；

(4)沉积旋回数据库C_p

对碎屑岩储层沉积旋回类型进行概括和分类，其中包括：正旋回C_1、反旋 回C_2、复合旋回C_3；

(5)成岩阶段数据库S_q成岩阶段数据库S_q包括：早成岩阶段A期ⅠA， 即为：S_1；早成岩阶段B期ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3； 中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；

(6)断层指数数据库N_r

根据研究区断层发育情况确定断层指数Ω(断层指数Ω为单位面积断层发育 的条数)，建立断层指数数据库，其中包括：一级指数N_1：Ω＜1；二级指数 N_2：1≤Ω≤3；三级指数N_3：3＜Ω≤5；四级指数N_4：5＜Ω。

本发明的有益效果在于：

本发明基于碎屑岩储层沉积相、岩性、沉积旋回、流体环境以及成岩阶段研 究，建立研究区不同属性类型下“环境属性-成岩相”模型，并以单井成岩相为 约束，对碎屑岩储层成岩相进行预测，确定有效储层空间展布，为油藏评价提供 依据。

附图说明

图1为东营组二段储层沉积相图；

图2为东营组二段储层沉积旋回图；

图3为东营组二段储层酸碱度图；

图4为东营组二段储层岩性图；

图5为东营组二段储层成岩阶段图；

图6为东营组二段储层断裂分布图；

图7为东营组二段成岩相图.

具体实施方式

为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容， 但本发明的内容不仅仅局限于以下实施例。

实施例1

基于井的碎屑岩储层成岩相横向预测方法，包括以下步骤：

1、收集研究区地质资料，包括以下内容：

①研究区目标层沉积相研究成果图；

②研究区岩性研究成果图；

③研究区断裂分布研究成果图；

④研究区流体史(酸碱度)资料；

⑤研究区目标层成岩阶段研究成果图；

⑥研究区目标层沉积旋回研究成果图；

⑦研究区岩石薄片分析测试数据；

2、建立原始资料数据库

(1)沉积相类型数据库F_m

对碎屑岩储层沉积相类型进行概括和分类，其中包括：冲积扇-旱扇-扇根 (F_1)、冲积扇-旱扇-扇中(F_2)、冲积扇-旱扇-扇缘(F_3)、冲积扇-湿扇-扇 根(F_4)、冲积扇-湿扇-扇中(F_5)、冲积扇-湿扇-扇缘(F_6)、河流相-曲流河 -河床亚相-河床滞留(F_7)、河流相-曲流河-河床亚相-边滩(F_8)、河流相-曲 流河-堤岸亚相-天然堤(F_9)、河流相-曲流河-堤岸亚相-决口扇(F_10)、河流 相-曲流河-河漫亚相-河漫滩(F_11)、河流相-曲流河-泛滥盆地-河漫湖泊(F_12)、 河流相-曲流河-河漫沼泽(F_13)、河流相-辫状河-牛轭湖(F_14)、河流相-辫状 河-河床滞留(F_15)、河流相-辫状河-心滩(F_16)、河道(F_17)、河流相-辫状 河-泛滥平原(F_18)、湖泊相-断陷型-湖成三角洲(F_19)、湖泊相-断陷型-滨湖 (F_20)、湖泊相-断陷型-浅湖(F_21)、湖泊相-断陷型-半深湖(F_22)、湖泊相 -断陷型-深湖(F_23)、湖泊相-断陷型-湖湾(F_24)、湖泊相-坳陷型-湖成三角 洲(F_25)、湖泊相-坳陷型-滨湖(F_26)、湖泊相-坳陷型-浅湖(F_27)、湖泊相 -坳陷型-半深湖(F_28)、湖泊相-坳陷型-深湖(F_29)、湖泊相-坳陷型-湖湾(F_30)、 湖泊相-前陆型-湖成三角洲(F_31)、湖泊相-前陆型-滨湖(F_32)、湖泊相-前陆 型-浅湖(F_33)、湖泊相-前陆型-半深湖(F_34)、湖泊相-前陆型-深湖(F_35)、 湖泊相-前陆型-湖湾(F_36)、三角洲相-辫状河三角洲-三角洲平原-分支(F_37)、 三角洲相-辫状河三角洲-三角洲平原-河道(F_38)、三角洲相-辫状河三角洲- 三角洲平原-陆上天然堤(F_39)、三角洲相-辫状河三角洲-三角洲平原-决口扇(F_40)、三角洲相-辫状河三角洲-三角洲平原-沼泽(F_41)、三角洲相-辫状 河三角洲-三角洲平原-淡水湖泊(F_42)、三角洲相-辫状河三角洲-三角洲前缘- 水下分支河道(F_43)、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤(F_44)、 三角洲相-辫状河三角洲-三角洲前缘-支流间湾(F_45)、三角洲相-辫状河三角 洲-三角洲前缘-分支河口砂坝(F_46)、三角洲相-辫状河三角洲-三角洲前缘- 远砂坝(F_47)、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂(F_48)、三 角洲相-辫状河三角洲-前三角洲(F_49)、三角洲相-扇三角洲-扇三角洲平原- 泥石流(F_50)、三角洲相-扇三角洲-扇三角洲平原-河道充填(F_51)、三角洲 相-扇三角洲-扇三角洲平原-漫滩(F_52)、三角洲相-扇三角洲-扇三角洲前缘- 碎屑流(F_53)、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道(F_54)、三 角洲相-扇三角洲-扇三角洲前缘-支流间湾(F_55)、三角洲相-扇三角洲-扇三角 洲前缘-河口砂坝(F_56)、三角洲相-扇三角洲-扇三角洲前缘-远砂坝(F_57)、 三角洲相-扇三角洲-前扇三角洲(F_58)。

(2)岩性数据库R_n

根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩 性数据库；其中，包括：细砾岩R_1、中砾岩R_2、粗砾岩R_3、巨砾岩R_4、 石英砂岩R_5、长石质石英砂岩R_6、岩屑质石英砂岩R_7、长石岩屑质石英砂 岩R_8、长石砂岩R_9、岩屑质长石砂岩R_10、岩屑长石砂岩R_11、岩屑砂岩 R_12、长石质岩屑砂岩R_13、长石岩屑砂岩R_14、黏土R_15、泥岩R_16、页 岩R_17。

(3)酸碱度数据库P_o

根据储层孔隙中流体酸碱性大小，建立酸碱度数据库；其中包括：强酸性 P_1：pH＜3；酸性P_2：3≤pH＜5；弱酸性P_3：5≤pH＜7；弱碱性P_4：7≤ pH＜9；碱性P_5：9≤pH＜11；强碱性P_6：11≤pH≤14；

(4)沉积旋回数据库C_p

对碎屑岩储层沉积旋回类型进行概括和分类，其中包括：正旋回C_1、反旋 回C_2、复合旋回C_3；

(5)成岩阶段数据库S_q

成岩阶段数据库包括：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期 ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即 为S_4；晚成岩阶段Ⅲ，即为：S_5；

(6)断层指数数据库N_r

根据研究区断层发育情况确定断层指数Ω(断层指数Ω为单位面积断层发育 的条数)，建立断层指数数据库，其中包括：一级指数N_1：Ω＜1；二级指数 N_2：1≤Ω≤3；三级指数N_3：3＜Ω≤5；四级指数N_4：5＜Ω；

3、建立成岩相划分方案

根据研究区孔隙类型、成岩特征以及面孔率，建立成岩相类型D_s划分方 案如下：

4、建立成岩相预测模型

(1)储层网格化

储层网格化是将空间上不均匀分布的数据按一定方法(如滑动平均法、克里 格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程； 研究区储层的每个网格可用Wi(X，Y)表示；

(2)确定网格属性

a.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；

b.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；

c.根据研究区酸碱度数据确定网格Wi(X，Y)的酸碱度属性P_o；

d.根据研究区沉积旋回数据确定网格Wi(X，Y)的旋回属性C_p；

e.根据研究区成岩阶段数据确定网格Wi(X，Y)的成岩阶段属性S_q；

f.根据研究区断层指数数据确定网格Wi(X，Y)的断层指数属性N_r；

g.根据上述步骤a、b、c、d、e和f确定每个网格Wi(X，Y)的综合属性 G(F_m，R_n，P_o，C_p，S_q，N_r)，即：

Wi(X，Y)＝F_m+R_n+P_o+C_p+S_q+N_r；

(3)建立成岩相预测模型

a.基于研究区岩石薄片资料、分析测试资料，确定研究区已知井Hj所在网 格Wj(X，Y)的成岩相D_s_j；

b.确定已知井Hj所在网格Wj(X，Y)的综合属性Gj(F_m_j，O_n_j，S_o_j， C_p_j，P_q_j，N_r_j)，建立已知井Hj所在网格Wj(X，Y)成岩相D_s与综合 属性G(F_m_j，O_n_j，S_o_j，C_p_j，P_q_j，N_r_j)的对应关系，即：

D_s_j＝F_m_j+O_n_j+S_o_j+C_p_j+P_q_j+N_r_j；

(4)成岩相预测

a.任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)的综合属性G (F_m，R_n，P_o，C_p，S_q，N_r)；

Gi＝(F_m_i，R_n_i，P_o_i，C_p_i，S_q_i，N_r_i)

b.将网格Wi(X，Y)的成岩相综合属性Gi与井Hj所在网格Wj(X，Y) 的综合属性Gj相比较，即：

Gi-Gj＝(F_m_i，R_n_i，P_o_i，C_p_i，S_q_i，N_r_i)-(F_m_j+ R_n_j+P_o_j+C_p_j+S_q_j+N_r_j)；

如果满足Gi-Gj＝0，则未知网格Wi(X，Y)的成岩相D_s_i与井j具有相 同的成岩相D_s_j；

如果Gi-Gj≠0，则按不同属性优先级顺序，

即：沉积相F_m一级、岩性R_n二级、沉积旋回C_p三级、酸碱度P_o四 级、成岩阶段S_q五级、断层指数N_r六级对未知网格Wi(X，Y)的成岩相进 行判识，即：

ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，C_p_i-C_p_j≠0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j≠0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅲ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅳ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0， S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅴ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0，S_q_i-S_q_j＝0，N_r_i-N_r_j≠0；

未知网格Wi(X，Y)的成岩相D_s_i满足条件i时，未知网格Wi(X，Y) 与具有相同属性F_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j 的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅲ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅳ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j、P_q_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅴ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j、 C_p_j、P_o_j、S_q_j的网格Wj(X，Y)的成岩相相同；

c.根据每个网格Wi(X，Y)的成岩相确定研究区横向分布。

实施例2

基于井的碎屑岩储层成岩相横向预测方法预测渤中凹陷东营组二段储层，具 体步骤如下：

1、收集研究区地质资料

收集渤中凹陷东营组二段储层的地质资料，并按照原始资料数据库对各种类 型数据进行分类。收集的地质资料包括以下内容：东营组二段储层沉积相图、东 营组二段储层岩性图、东营组二段储层断裂分布图、东营组二段储层酸碱度图、 东营组二段储层成岩阶段图、东营组二段储层沉积旋回图(如图1～6所述)。

2确定储层属性

(1)网格化工区

将渤中凹陷东营组二段碎屑岩储层进行网格化，每个网格用Wi(X，Y)表 示。

(2)确定网格属性

a.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；

b.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；

c.根据研究区酸碱度数据确定网格Wi(X，Y)的酸碱度属性P_o；

d.根据研究区沉积旋回数据确定网格Wi(X，Y)的旋回属性C_p；

e.根据研究区成岩阶段数据确定网格Wi(X，Y)的成岩阶段属性S_q；

f.根据研究区断层指数数据确定网格Wi(X，Y)的断层指数属性N_r；

g.根据上述步骤a、b、c、d、e和f确定每个网格Wi(X，Y)的综合属性 G(F_m，R_n，P_o，C_p，S_q，N_r)，即：

Wi(X，Y)＝F_m+R_n+P_o+C_p+S_q+N_r；

3.确定研究区成岩相

(1)确定单井成岩相

收集研究区岩石薄片分析测试数据，根据岩石颗粒接触关系、溶蚀类型、胶 结类型、原生孔隙、次生孔隙以及胶结含量，结合成岩相划分方案，确定不同井 东营组二段的成岩相类型，即已知井Hi所在网格Wi(X，Y)的成岩相所对应的 综合属性G(F_m，R_n，P_o，C_p，S_q，N_r)。

(2)确定井间成岩相

a.任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)的综合属性G (F_m，R_n，P_o，C_p，S_q，N_r)；即为：

Gi＝(F_m_i，R_n_i，P_o_i，C_p_i，S_q_i，N_r_i)

如网格W30＝(F_43，R_6，ⅠB，C_2，P_2，N1)；

如网格W40＝(F_43，R_10，ⅠB，C_3，P_2，N3)；

如网格W50＝(F_26，R_11，ⅡA，C_1，P_4，N2)；

如网格W60＝(F_45，R_6，ⅠB，C_3，P_2，N2)；

b.将网格Wi(X，Y)的成岩相综合属性Gi与井Hj(CFD16-1-2D、CFD16-1-1、 Hz4、CFD17-3-1、Hz5、CFD18-1-2D、CFD18-2-2D、CFD18-2E-1、BZ13-1-1、 BZ13-1-2、BZ13-1-3、CFD23-1-1、CFD23-3-1、BZ19-4-1、CFD24-1-1、CFD16-2-1、 CFD16-3-1、CFD18-1N-1、BZ19-2-1)所在网格Wj(X，Y)的综合属性Gj相 比较，即：

Gi-Gj＝(F_m_i+O_n_i+S_o_i+C_p_i+P_q_i+N_r_i)

-(F_m_j+O_n_j+S_o_j+C_p_j+P_q_j+N_r_j)

如果满足Gi-Gj＝0，则未知网格Wi(X，Y)的成岩相D_s_i与井j具有相 同的成岩相D_s_j；

如果Gi-Gj≠0，则按不同属性优先级(即为：沉积相一级、岩性二级、旋 回三级、酸碱度四级、成岩阶段五级、断层指数六级)对未知网格Wi(X，Y) 的成岩相进行判识，即：

ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，C_p_i-C_p_j≠0，P_o_i-P_o_j ≠0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j≠0，P_o_i-P_o_j≠ 0，S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅲ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j≠0， S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅳ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0， S_q_i-S_q_j≠0，N_r_i-N_r_j≠0；

ⅴ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，C_p_i-C_p_j＝0，P_o_i-P_o_j＝0， S_q_i-S_q_j＝0，N_r_i-N_r_j≠0；

未知网格Wi(X，Y)的成岩相D_s_i满足条件i时，未知网格Wi(X，Y) 与具有相同属性F_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、O_n_j 的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅲ时，未知网格Wi(X，Y)与具有相同属性F_m_j、O_n_j、 C_p_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅳ时，未知网格Wi(X，Y)与具有相同属性F_m_j、O_n_j、 C_p_j、P_q_j的网格Wj(X，Y)的成岩相相同；

D_s_i满足条件ⅴ时，未知网格Wi(X，Y)与具有相同属性F_m_j、O_n_j、 C_p_j、P_q_j、S_o_j的网格Wj(X，Y)的成岩相相同；

如网格W30＝(F_43，R_6，ⅠB，C_2，P_2，N1)- (F_m_CFD17-3-1+O_n_CFD17-3-1+S_o_CFD17-3-1+C_p_CFD17-3-1+P_q_ CFD17-3-1+N_r_CFD17-3-1)＝0，即网格W30与井CFD17-3-1具有相同的成岩 相，也即：W30＝D_4；

如网格W40＝(F_43，R_10，ⅠB，C_3，P_2，N3)-(F_ CFD18-2E-1+O_n_CFD18-2E-1+S_o_CFD18-2E-1+C_p_CFD18-2E-1+P_q_ CFD18-2E-1+N_r_CFD18-2E-1)＝0，即网格W40与井CFD18-2E-1具有相同的 成岩相，也即：W40＝D_2；

如网格W50＝(F_26，R_11，ⅡA，C_1，P_2，N1)网格属性满足条件 ⅳ，即：F_m-F_m_BZ-13-1-3＝0，O_n-O_n_BZ-13-1-3＝0，C_p-C_p_BZ-13-1-3＝0， P_q-P_q_BZ-13-1-3＝0，S_o-S_o_BZ-13-1-3≠0，N_r-N_r_BZ-13-1-3≠0；则 W80与井BZ-13-1-3具有相同成岩相，也即：W50＝D_4；

如网格W60＝(F_26，R_11，ⅡA，C_1，P_2，N1)网格属性满足条件 ⅴ，即：F_m-F_m_CFD-16-2-1＝0，O_n-O_n_CFD-16-2-1＝0，C_p-C_p_ CFD-16-2-1＝0，P_q-P_q_CFD-16-2-1＝0，S_o-S_o_CFD-16-2-1＝0，N_r-N_r_ CFD-16-2-1≠0；则W60与井CFD-16-2-1具有相同成岩相，也即：W60＝D_2；

c.根据每个网格Wi(X，Y)的成岩相确定研究区成岩相横向分布(如图7 所示)。

其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽 的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据 本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范 围。