[0001] 本发明涉及沉积相恢复方法，具体为一种地震相驱动深层碳酸盐岩高精度沉积微相恢复方法，属于石油勘探与开发技术领域。

[0002] 海相碳酸盐岩地层是中国油气勘探的重要方向，优质储层通常发育在滩、礁等优势沉积相带中。四川、塔里木盆地碳酸盐岩地层勘探过程中通过优势沉积相带为指导发现了一大批大中油气田。随着勘探程度的加深，对于沉积相恢复的精细程度有了更高的要求。

[0003] 传统沉积相恢复方法主要依据露头、钻井获取的实际地质信息，通过各种相标志划分沉积相，继而根据多因素综合分析单因素作图法绘制平面沉积相图。该方法虽然对于具有较大局限性，以点窥面，对于各露头、钻井之间广大的区域主要依赖于插值或专家经验，主观性大，难以精确划分相带边界。其次传统平面沉积相绘制成图时间单元较大，对于层序内沉积相演化难以精细刻画。同时随着勘探向盆地深部进军，同一区域深层碳酸盐岩地层钻遇井密度相对于中浅层地层钻遇井密度显著降低，传统沉积相恢复方法对于少井区域难以使用。

[0037] 结合图1，以一实例对本发明作进一步的描述：一种地震相驱动深层碳酸盐岩高精度沉积微相恢复方法，包括以下步骤：

[0038] S1、基础资料收集：收集研究区域的三维地震体，目的层钻测井、露头资料、相关地质概况等资料。

[0039] S2、井‑震标定：开展精细合成记录标定，建立深度域钻测井地质信息同时间域地震信息关系。

[0040] 其中，该实例研究区内地震数据为S地区某三维叠后地震体，区内钻井包括JT‑1井、JG‑1井、YB‑6井、YB‑7井等，研究地层为二叠系茅口组地层。

[0041] S3、建立区域等时层序地层格架：在层序地层学理论指导下，划分三、四级层序、以井‑震结合的思路建立研究区区域等时格架，该过程中须不断校验重复S2、S3步骤以达到最优匹配，然后开展地震层位追踪。

[0042] 其中该实例研究区域茅口组可划分为SQ1、SQ2、SQ3等3个三级层序，可识别出4个层序界面，分别为SB1、SB2、SB3和SB4，划分的三个层序在S地区整体可对比性强(见图2)。SQ1在表现出厚度大不同区域差异小、GR值高的特点，表明SQ1沉积期整体而言水体较深；
SQ2表现出厚度较薄不同区域厚度存在一定差异、GR低的特点；SQ3表现出不同区域厚度差异大，南厚北薄的特点，GR值与之具有对应南低北高的特点。

[0043] S4、测井解释：对于全区测井资料进行标准化，去除因采集的设备、年份差异造成的影响，继而开展测井解释，通过岩心标定，确定颗粒岩门槛值，既而统计相关厚度并分层序计算颗地比(颗粒滩厚度与地层厚度比值)，作为定量沉积相指标。

[0044] 其中该实例中茅口组SQ3共统计18口井、SQ2统计16口井、SQ1统计10口井的单层序颗地比(图3和图4)。

[0045] S5、层序格架内沉积相划分：基于岩心、薄片、测井解释，根据各类相标志以三级层序及以下为单元进行沉积相划分。

[0046] 其中沉积相划分与对比如图5，该地区茅口组沉积期沉积相为碳酸盐岩台地。根据颗地比的不同，将颗地比与沉积微相相应，结合取心将颗地比达到0.6以上对应为高能滩微相，颗地比介于0.3到0.6之间对应为低能滩微相，颗地比小于0.3对应为开阔海\局限海微相。

[0047] S6、地震相划分：以三级层序为基本单元，利用三维地震开展波形聚类地震相划分，同时在典型过井格架线基于专家经验开展相面法分析，两者结合互相验证后，总结典型地震相类型。

[0048] 其中SQ1地震相划分为4类，SQ2地震相划分为6类，SQ3划分为6类，茅口组整体划分为5类，各地震相平面边界清晰(图6到图9)。

[0049] S7、建立地震相‑沉积相转化关系：以S5划分层序格架内沉积相和S6划分地震相进行比对分析，建立典型地震相‑沉积相转化关系，根据沉积相边界确定单个三级层序内沉积相边界。

[0050] 其中以茅口组整体地震相为例，结合典型井地震剖面、波形类型及响应特征、地震相特征、典型岩芯及薄片特征，典型测井相及沉积结构以及单因素沉积指标，划分地震相‑沉积相转化模式(图10)。

[0051] S8、等时层序格架内划分地层切片：基于地震沉积学原理，利用三维地震资料在单层序单元内开展等时地层切片划分，具体地层切片数目应与划分四级及以下测井可识别旋回相对应。

[0052] 其中地层切片以三级层序为单位分别制作，具有较好的精度。根据研究重点的不同，3个三级层序切片数量也不同，对于每个三级层序依据以追踪层位两两之间按照厚度均分的原则内插地层切片，其中SQ1对应地层7分，茅口组底界切片编号为SQ1_1，其上依次编号至SQ1_7；SQ2对应地层7分，SQ2底界切片编号为SQ2_1，其上依次编号至SQ2_7；SQ3对应地层8分，SQ3底界切片编号为SQ3_1，其上依次编号至SQ3_8，茅口组顶界编号为SQ3_9；茅口组内部共计制作23个地层切片。

[0053] S9、优选地震属性：通过三维地震数据体提取地震属性，通过地震属性同井上测井识别单因素沉积指标相关性，地震属性平面敏感性进行优选，提取最能反应地层沉积特征的属性。

[0054] 其中按照三级层序提取了均方根振幅、最大振幅、能量半衰时、平均能量、瞬时Q因子、瞬时相位、瞬时频率、t振幅、道积分属性等进行研究，本实例将三维资料范围内钻穿SQ3的井共计6口，包括YB‑3、YB‑6、YB‑7、YS‑1、YB‑224、JG‑1，逐层提取不同属性不同切片的井点处属性值与井上的颗地比进行相关性分析得到相关性热力图(图11)。

[0055] 其中以颗地比为沉积相定性评价参数的条件下，瞬时Q因子、瞬时频率、道积分属性与其具有良好的相关性，理论上更加能够反映真实的地层情况。其次是振幅类属性最为敏感有着较为良好的相关性，和地层的拟合关系更好，可以为沉积相的研究提供支持。

[0056] 其中为了更清晰的展示相同地层切片不同属性的平面响应特征差异，提取SQ3_8号切片的均方根振幅、最大振幅、能量半衰时、平均能量、瞬时Q因子、瞬时频率、t振幅、道积分属性等进行研究，并绘制对比图件，结合井上信息开展地震属性滩体响应敏感性分析(图12)。

[0057] 其中SQ_8号切片位于SQ3的上部，距离茅口组顶界较近，由于受到茅口组末期暴露剥蚀和上覆吴家坪组泥页岩地层的原因，使得茅口组顶界为一强波峰，对于三维地震资料邻近区域造成影响，使得振幅类属性对于台缘边界的刻画并不清晰(图12中的a，图12中c，图12中g)，但是在台内地区通过和井上对比对于颗粒滩分布有相对良好的反映。能量半衰时属性则完整的刻画茅口组沉积晚期SQ3层序中的台地边缘边界(图12中的f)，这与S7中应用神经网络波形分类计算的SQ3期地震相的边界相一致(图6到图9)。瞬时频率、瞬时Q因子属性也在一定程度上刻画了台地边缘边界(图12中b,图12中e)。道积分属性的相对波阻抗不仅对于台缘边界有比较好的刻画，同时体现台缘地区颗粒滩发育程度的差异，其红色、黄色区域为颗粒滩发育地区，绿色区域则为台内开阔海范围(图12中d)，其范围与地震相范围相一致。

[0058] S10、建立单因素沉积指标‑地震属性关系：将优选的地震属性同单因素沉积指标进行拟合，建立拟合方程，使地震属性能够转换为单因素沉积指标。

[0059] S11、进行属性体等时地层切片：依据S8划分地层切片对于地震属性体进行地震属性体等时地层切片制作。

[0060] 其中按照层序格架分层序制作共计23个地层切片。选取其中具有代表性的地层展现实例中茅口组台地边缘滩体演化，其中SQ1对应地层选择靠近顶底SQ1_3、SQ1_6，SQ2对应地层按照下中上选择SQ2_1、SQ2_4、SQ2_7，SQ3对应地层按照下中上选择SQ3_3、SQ3_5、SQ3_8(图13)，从图中可以看出在在SQ1沉积早期的地层的滩体响应的红色、黄区域较少呈点状零星分布，在晚期则明显增加同时展现出一定的条带性；SQ2沉积期则红色、黄色低阻抗区域广泛分布，呈现明显的条带性，其在研究区的南部最为集中，结合图6到图9可知SQ2期该区域的主要地震相为丘状‑稧状、中弱振幅地震相；SQ3期全区表现出弱波阻抗的响应，北部出现大片的红色、黄色区域，但是在SQ3期由于槽台快速分异，台地边缘以北沉积了大套的深水相泥页岩沉积，其波阻抗是低于颗粒灰岩波阻抗的，应与台内区域的指示颗粒灰岩的滩相沉积有所区别。

[0061] S12、地震相‑沉积相约束下高精度等时沉积微相图绘制：根据前述内容在地震相‑沉积相约束下，依据S10拟合方程将S11制作属性体地层切片转化为平面单因素沉积指标切片，据此绘制等时沉积微相图，并用预留井进行检验。

[0062] 其中选取道积分属性其中8个典型等时地层切片转换为等时沉积微相(图14)。茅口组沉积期不同时期，台地边缘滩体分布范围有所差异，SQ1为零星条带与点状分布，逐步过渡到SQ2、SQ3期大面积连带分布，垂向有较好的连续性，并且SQ3期在元坝地区向北迁移，形成自北向南的斜坡‑高能滩‑低能滩‑开阔海的格局。

[0063] 以上所述仅为本发明实施实例而已，并不用于限制本发明，对本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、同等替换、改进等，均应包含在本发明的权利要求范围之内。