分析河流的沉积微相平面展布与沉积演化的方法和系统与流程

1.本发明属于油气勘探、开发领域，涉及一种分析河流的沉积微相平面展布与沉积演化的方法和系统。


背景技术：

2.经典层序地层学在我国油气勘探方面取得了巨大成就，但通常其应用于低频三级层序及以上级别层序的地层，在开发阶段成功应用的实例较少。近年来，四级及其以下高频层序，即经典层序地层学中的准层序组级别的高频层序，在油气藏开发阶段被越来越重视。精细的高频层序划分与对比往往依赖于地质露头的详细观察，而在缺乏地质露头对比或钻井取心厚度有限的条件下，基于测井曲线进行四级及其以下级别的高频层序的划分方法相对较少。而且所应用的测井曲线往往集中于自然伽玛(gr)、声波时差(ac)等一条或两条，存在多解性大、描述方法也都过于笼统、未对层序界面的测井曲线响应及变化特征进行详细分析的问题，因而实现全井段高频层序划分的实例较少。
3.近年来，碎屑岩微相的分析研究工作日益受到关注，并得以迅速发展，已成为储层研究中不可缺少的内容。许多学者对沉积微相的研究也做出了重要的论述。“微相”一词最早由brown(1943)提出，fl
ü
gel(2004)进一步将“沉积微相”定义为可在薄片、揭片和磨光片中区分出所有古生物学和沉积学特征的综合。沉积微相的研究在方法和技术上，通常强调野外工作和室内薄片研究的结合，利用野外露头、岩心、常规薄片、扫描电镜、阴极发光、x衍射、矿物学与地球化学等来分析沉积微相。在盆地或区带勘探的阶段应用广泛，即通过沉积微相的研究，分析沉积环境，从而寻找有利的勘探区带。相对而言，油气藏开发阶段，在四级或五级高频层序格架的约束下，更加精细的开展沉积微相的纵向和横向对比及沉积演化研究的实例较少。
4.碎屑岩地层中的高效油气勘探和开发，离不开沉积微相的指导。传统的平面沉积微相的研究一般是统计单井上各层的地层厚度、砂岩厚度和含砂率等数据，考虑各井的各层段的测井曲线形态，以井点为控制点，勾绘工区内的等值线，进而在平面上分析砂体的展布状况以及确定工区的沉积微相展布。目前，利用地震资料进行沉积相的研究主要包括反演与地震属性的分析。反演方法目前常用多井约束进行反演，但由于过多考虑了井信息，其预测性受到一定限制。地震属性方法在目前应用较广，更遵循地震信息。
5.对于一个油气田来说，在勘探阶段，井距较大，单凭井资料研究的精度是远远不够的。而在开发阶段，虽然井网密度较大，但仅靠单井信息是达不到开发调整要求，需要借助地震信息辅助刻画沉积微相的展布分析，以此达到勘探、开发调整的要求。但由于地震分辨率有限的原因，一般认为不能直接利用地震资料进行高频层序单元的沉积相刻画。在垂向上，陆相油藏属于薄互层结构，单层厚度往往低于地震分辨率——四分之一主波长，无法达到真正薄层意义上的地震成像。
6.针对复杂的辫状河油气藏在勘探开发阶段所遇到沉积微相的多变、纵向上多期的切割叠加、沉积微相的平面展布复杂、沉积演化的规律难以总结等难题，有必要提供一种方
案来解决上述问题，达到油气藏精细描述的目的。


技术实现要素：

7.本发明所要解决的技术问题之一是需要提供一种分析河流的沉积微相平面展布和沉积演化的方法，该方法将基于高频层序划分与利用地震沿层切片分析相结合，能够准确描述以高频层序为约束和划分单元的沉积微相展布，并进行河道的迁移摆动演化分析。
8.为了解决上述技术问题，本技术的实施例首先提供了一种分析河流的沉积微相平面展布与沉积演化的方法，该方法包括：步骤一，对目的层的岩心样品进行分析，确定目的层的岩性和沉积环境；步骤二，利用岩心样品的分析结果对选取的测井曲线进行组合和标定，获取岩性敏感曲线及其岩性敏感变化特征；步骤三，利用岩性敏感曲线的岩性敏感变化特征，对标志层和高频层序界面进行识别以及对测井相和地震相进行标定；步骤四，利用所述步骤三的识别结果和标定结果，通过井震结合建立高频层序格架并以此高频层序格架来划分目的层的目的小层，在测井相和地震相的综合响应特征的分析约束下，分析以目的小层为单元的沉积微相平面展布以及主河道的沉积演化过程。
9.在一个实施例中，在所述步骤二中，执行如下步骤：选择多种测井曲线，根据所述岩心样品的岩性分析结果，对这些测井曲线进行组合并标定，确定不同岩性的各种测井曲线的响应特征，并从测井曲线组合中寻找对岩性敏感的测井曲线作为岩性敏感曲线，以及将它们随岩性变化的变化趋势作为岩性敏感变化特征。
10.在一个实施例中，所述测井曲线包括自然伽马测井曲线、自然电位测井曲线、井径测井曲线、深侧向电阻率测井曲线、浅侧向电阻率测井曲线和三孔隙度测井曲线；在进行组合时，将自然伽马测井曲线、自然电位测井曲线和井径测井曲线至于同一测井曲线道，将三孔隙度测井曲线和电阻率曲线分别至于另外两个测井曲线道。
11.在一个实施例中，在所述步骤二中，还包括：根据所述岩心样品的岩性分析结果和测井曲线组合的岩性敏感变化特征选取致密砂岩段，利用实测的岩心孔隙度对计算得到的三条孔隙度曲线进行校正，将校正后的三孔隙度计算曲线在致密砂岩段处重合作为岩性敏感曲线之一。
12.在一个实施例中，在所述步骤三中，通过如下步骤对标志层进行识别：在测井曲线组合的岩性敏感变化特征及重合后的三孔隙度计算曲线的基础上，结合区域钻井地层、岩性特征资料，利用测井曲线突变特征所反映的岩性突变，对目的层已有钻井的标志层进行测井识别和级次划分。
13.在一个实施例中，在所述步骤三中，通过如下步骤对高频层序界面进行识别：利用标志层划分结果，通过地震合成记录来标定地震剖面，并建立标志层的岩性界面、测井界面和地震界面彼此对应的综合响应特征模式；在标志层的综合响应特征模式的地震响应特征的约束下，开展标志层的井间横向对比；以标志层的井间横向对比结果作为约束，总结高频沉积旋回的岩性变化趋势以及不同高频沉积旋回接触面上下的岩性突变差异，识别高频层序界面；利用所述高频层序界面，标定对岩性敏感的测井曲线与重合后的三孔隙度计算曲线，总结高频层序界面的测井响应特征，进而完成全井段的高频层序界面识别。
14.在一个实施例中，在所述步骤三中，通过如下步骤对测井相和地震相进行标定：对目的层利用岩石相标定测井曲线得到测井相，利用测井相标定地震剖面响应特征得到地震
相，进而总结岩石相、测井相和地震相彼此对应的综合响应模式。
15.在一个实施例中，在所述步骤四中，包括：根据识别的高频层序界面建立高频层序格架，并依据砂体分布和实际开发需要，划分目的层的目的小层，并在地震剖面和标志层对比剖面的约束下，开展目的小层的对比。
16.在一个实施例中，在所述步骤四中，还包括：以岩石相、测井相和地震相彼此对应的综合响应模式为沉积微相分析的参考依据，以高频层序格架为约束，完成目的层的单井沉积微相划分和连井沉积微相的对比；以单井沉积微相划分、井间沉积微相对比为基础，依据井震标定结果，计算每个目的小层厚度在地震剖面上所对应的时窗，利用沿层地震切片技术和该时窗提取能够反映各个小层沉积微相的地震相属性，实现目的小层的平面沉积微相展布和河道带主流线的绘制；将目的小层的平面沉积微相的河道主流线在平面上叠置，分析总结主河道的迁移摆动特征得到河流沉积演化过程。
17.根据本发明的另一方面，还提供了一种分析河流的沉积微相平面展布与沉积演化的系统，该系统执行如上所述的方法。
18.与现有技术相比，上述方案中的一个或多个实施例可以具有如下优点或有益效果：
19.本技术实施例提出了一种分析河流的沉积微相平面展布和沉积演化的方法，该方法基于河流相的高频层序格架约束，通过岩心及薄片的沉积相标志分析，开展单井和连井沉积微相分析，以例如四级高频层序为作图单元，建立基于井震结合与地震沿层切片技术的小层沉积微相平面展布与主河道沉积演化。该方法能够准确描述以高频层序为约束和划分单元的沉积微相展布，并进行河道的迁移摆动演化分析。
20.本发明的其它特征和优点将在随后的说明书中阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明的技术方案而了解。本发明的目的和其他优点可通过在说明书、权利要求书以及附图中所特别指出的结构和/或流程来实现和获得。
附图说明
21.附图用来提供对本技术的技术方案或现有技术的进一步理解，并且构成说明书的一部分。其中，表达本技术实施例的附图与本技术的实施例一起用于解释本技术的技术方案，但并不构成对本技术技术方案的限制。
22.图1为本技术实施例的分析河流的沉积微相平面展布与沉积演化的方法的流程示意图。
23.图2为本发明实施例的目的层段的地层岩心的岩性、沉积构造的特征照片。
24.图3为本发明实施例的目的层段的测井曲线的岩性识别的交会图。
25.图4为本发明实施例的利用目的层段a1井的实测岩心孔隙度对三孔隙度校正后，将其在致密砂岩段处重合的结果图。
26.图5为本发明实施例的利用目的层段的标志层岩心界面来标定测井标志层界面的结果图。
27.图6为本发明实施例的将目的层段上古地层的标志层岩心界面、测井界面、地震界面相互标定的结果及响应特征图。
28.图7为本发明实施例的通过目的层段的井震结合得到的不同级别标志层的地震界
面识别结果在地震剖面上的追踪情况的示意图。
29.图8为本发明实施例的使目的层段在地震剖面约束下的标志层在连井剖面上的对比结果图。
30.图9为本发明实施例的利用目的层段a1的井岩心高频层序的界面识别对测井曲线进行标定的结果图。
31.图10为本发明实施例的目的层段的高频层序划分及小层对比结果图。
32.图11为本发明实施例的目的层段通过岩心标定测井和测井标定地震，建立沉积微相、测井相和地震相的综合响应模式的示意图。
33.图12为本发明实施例的使目的层段在研究区高频层序格架和地震剖面的约束下，连井沉积微相对比的结果图。
34.图13为本发明实施例的研究区的目的层段按照h12小层对应的时间切片的地震属性及其对应的沉积微相平面展布和河道主流线的刻画示意图。
35.图14为本发明实施例的研究区目的层段h11小层、h12小层和h13小层的辫状河道主流线的迁移摆动的平面演化过程示意图。
36.图15为本技术实施例的分析河流的沉积微相平面展布与沉积演化的系统的功能框图。
具体实施方式
37.以下将结合附图及实施例来详细说明本发明的实施方式，借此对本发明如何应用技术手段来解决技术问题，并达成相应技术效果的实现过程能充分理解并据以实施。本技术实施例以及实施例中的各个特征，在不相冲突前提下可以相互结合，所形成的技术方案均在本发明的保护范围之内。
38.另外，附图的流程图示出的某些步骤可以在诸如一组计算机可执行指令的计算机系统中执行。并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
39.本发明实施例提供了一种细致、可操作性强和符合沉积学原理的分析河流的沉积微相平面展布和河道迁移规律、沉积演化的方法。该方法针对缺乏野外露头对比的复杂碎屑岩地层的分布区，基于已有的钻井资料，建立了高频层序格架，并以此为约束。具体来说，通过对主要目的层段岩心、岩屑及其薄片的详细观察，分析主要的沉积环境，进而总结得到沉积微相类型。在高频层序的约束及沉积认识的指导下，开展单井沉积微相的分析和连井微相的对比，并以高频层序为分析单元，在地震解释标志层层位上选取合适的沿层时窗切片。通过提取能够反映沉积微相的敏感属性进行井震结合，实现沉积微相的平面展布。最后，利用河流(例如辫状河)河道主流线的迁移摆动分布趋势，来总结河流平面摆动沉积演化过程。
40.需要说明的是，在本实施例中以辫状河为例来进行说明，但并不限定与辫状河，还可以是曲流河，以及辫状河三角洲等。该实施例解决了如下技术问题：由于勘探开发阶段辫状河岩性气藏目的层缺乏露头对比、钻井少、岩性复杂，难以准确有效的开展基于高频层序格架约束的小层对比和精细描述以小层为单元的辫状河道砂体展布，以及河道迁移沉积演化过程分析难度大。
41.实施例一
42.图1为本技术实施例的分析河流的沉积微相平面展布与沉积演化的方法的流程示意图。概括来说，在本发明实施例中，利用不同级别标志层的岩心界面特征详细标定测井曲线，利用测井曲线的组合特征，对标志层、高频层序界面进行识别、对测井相和地震相进行标定。在全井段识别高频层序界面之后，利用地震标志层界面的识别在横向上进行对比约束，开展小层对比、分析小层内沉积微相和辫状河道迁移摆动过程。
43.下面参考图1来具体说明该方法涉及的各个步骤。
44.首先，在步骤s10(步骤一)中，对目的层岩心样品进行分析，确定目的层的岩性与沉积环境。
45.在本步骤中，通过详细观察目的层岩心样品，确定目的层主要岩性与大致沉积环境。具体地，通过对取心段的岩心、岩屑及其薄片进行详细观察来确定该目的层段的岩性和沉积环境。
46.图2为本发明实施例的目的层段处地层岩心的岩性、沉积构造的特征照片，该图像显示了杭锦旗太原组-下石盒子组气藏研究区目的层段(厚约250m-300m)的岩心的岩性和沉积构造。在图2中，从左侧开始说明这些岩心样品的岩性和沉积环境：锦94井：层位在盒1段，沉积环境为杂色砾岩河道滞留沉积；锦69井，层位在盒1段，岩性为下部为泥质粉砂岩、上部为粗砂岩，沉积环境为冲刷面；锦72井，层位在山1段，岩性为灰色含砾粗砂岩，沉积环境为粒序层理；锦75井，岩性为粗砂岩、沉积环境为槽状交错层理；锦93井，层位在盒2段，岩性为含砾粗砂岩，沉积环境为板状交错层理；锦21井，层位在盒一段，岩性为粗砂岩，沉积环境为块状层理；锦21井，为辫状河道非典型二元结构，下部河床沉积厚度远大于上部堤岸、河漫；锦107井，层位在盒3段，岩性为褐色泥岩；锦111井，层位在盒1段，岩性为花斑状泥岩；锦75井，层位在山1段，岩性为浅灰绿色泥岩；锦116井，层位在山1段，岩性为泥质粉砂岩，沉积环境为透镜状层理；锦21井，层位在山1段，河漫浅灰色粉砂岩与灰绿色、灰色泥岩交互沉淀。
47.如图2所示，通过对目的层段的取心进行详细的岩心观察与薄片镜下的鉴定结果显示，目的层以含砾粗砂岩、粗砂岩、中砂岩为主，层理类型多样，为高能(水动力条件比较强的)辫状河道的典型岩相发育。主要目的层盒1段的河道类型多以向上变浅的米级沉积旋回为主，呈非典型二元结构，下部河床沉积厚度远大于上部堤岸、河漫沉积，反映了典型的辫状河道沉积(见图2中间部分)。山西组-盒1段的泛滥平原沉积总体以深灰色、灰绿色、灰色等(弱)还原色为主，盒2+3段褐色等氧化色逐渐增多(图2)，以向上变浅的河漫浅灰色的粉砂岩与灰绿色、灰色的泥岩交互沉积为主(图2的最右侧)。
48.结合岩石相分析结果和前人对本区域的沉积研究结果，可以得到本区太原组-山西组-下石盒子组具有近物源、相变快的特点，主要属于(辫状河)三角洲-辫状河沉积体系，可划分辫状河三角洲、辫状河、曲流河-辫状河等沉积相。目的层h1段为辫状河沉积，进一步划分为辫状河河道心滩、辫状河河道充填和泛滥平原等微相沉积。
49.接下来，在步骤s20(步骤二)中，利用岩心样品的分析结果对选取的测井曲线进行组合和标定，获取岩性敏感曲线及其岩性敏感变化特征。在本例中，岩性敏感曲线包括对岩性敏感的测井曲线和在致密砂岩段处重合的通过计算得到的三孔隙度曲线(可简称“三孔隙度计算曲线”)。
50.具体地，选择多种测井曲线，根据步骤s10中岩心样品的岩性分析结果，对这些常规测井曲线进行组合并标定，确定不同岩性的各种测井曲线的响应特征，然后，从测井曲线组合中寻找对岩性敏感的测井曲线(作为岩性敏感曲线)，以及岩性敏感测井曲线随岩性变化的变化趋势，进而得到测井曲线组合的岩性敏感变化特征。
51.在本例中，优选常规测井曲线包括自然伽马(gr)测井曲线、自然电位(sp)测井曲线、井径(cal)测井曲线、深侧向电阻率(lld)测井曲线、浅侧向电阻率(lls)测井曲线和三孔隙度测井曲线(包括密度(den)、声波时差(ac)和中子孔隙度(cnl))这八条常规测井曲线。在进行组合时，将cal、gr、sp曲线至于同一测井曲线道，此处的组合是因为它们可以展现出基本的岩性特征。将三孔隙度测井曲线和电阻率曲线(lld和lls)分别至于另外两个测井曲线道。
52.此外，发明人经过研究发现，通过岩性鉴定结果来标定测井曲线时，针对单井岩心分析，在识别不同的岩性方面，自然伽马曲线和中子曲线在识别岩性方面的效果比较好，可以被认为是对岩性敏感的测井曲线，因此，本发明提出的岩性识别图版主要是中子曲线(cnl)与自然伽玛曲线(gr)结合的交会图，如图3所示的取心段不同岩性的测井曲线交会图的示例图。从图3中可以看出，砾岩、砂砾岩、含砾粗砂岩、粗砂岩的测井响应重叠性很高，基本可以归为一类，因此，将上述例子中的锦72地区岩性分为三大类，分别是粗岩性类(砾岩、砂砾岩、含砾粗砂岩、粗砂岩)、中细砂岩类和泥岩。自然伽玛和中子曲线随着这三类岩性的粒度变细，测井值呈由低到高变化。
53.需要说明的是，仅是利用选取的测井曲线组合可以识别大部分的岩性，但对于某些特殊情况来说，岩性识别的精确度不高，因此在本例中，提出的上述交会图可以很好地准确反应岩性特征，除交会图之外的其他测井曲线具有辅助的岩性识别作用，它们的结合能够更好地提高岩性识别的精确度。另外，本步骤所选择的测井曲线组合、交会图仅为一个优选示例，根据实际的应用情况可以选取其他的曲线，本发明实施例并不限定其他合理的组合方式。
54.另外，如图1所示，在该步骤s20中，还包括如下步骤：根据步骤s10的岩心样品的岩性分析结果和测井曲线组合的岩性敏感变化特征(具体指反应岩性和岩性界面的特征)，选取致密砂岩段，利用实测岩心孔隙度对计算得到的三条孔隙度曲线进行校正，将校正后的三孔隙度计算曲线在致密砂岩段处重合作为岩性敏感曲线之一。
55.具体地，利用声波和密度的孔隙度解释模型，分别计算声波孔隙度和密度孔隙度，并将计算的这两个孔隙度与中子孔隙度放置于同一测井曲线道，形成孔隙度测井曲线组合作为计算的三孔隙度曲线。然后利用实测岩心孔隙度(依据岩心样本实测出来的三种孔隙度)对计算得到的三条孔隙度曲线进行校正，然后将校正后的三孔隙度计算曲线在致密砂岩段处重合。
56.需要说明的是，该步骤的设置是为了进一步提高岩性区分和岩性界面的区分的精确度，在三孔隙度重合时反应出当前层为致密砂岩，由重合到分开的状态，则反映了另一种岩性。
57.如图4a1井在3027m-3034m处，岩心样品孔隙度值介于1.2％-8.7％，平均4.07％，为典型的致密砂岩段。密度骨架平均值和声波骨架平均值分别取值约为2.63g/cm3和180μs/m，通过计算密度孔隙度和声波孔隙度，并与中子孔隙度放置于同一测井曲线道，形成计
算的孔隙度测井曲线组合，利用实测岩心孔隙度对三孔隙度(密度孔隙度、声波孔隙度与中子孔隙度)测井曲线进行校正，然后将校正后的三孔隙度计算曲线在该层段处大致重合。
58.在步骤s30(步骤三)中，利用岩性敏感曲线的岩性敏感变化特征，对标志层和高频层序界面进行识别以及对测井相和地震相进行标定。
59.在该步骤中包括标志层识别步骤s310、高频层序界面识别步骤s320～s340和测井相、地震相标定步骤s350。
60.在步骤s310中，在步骤s20的测井曲线组合的岩性敏感的特征识别结果及重合后的三孔隙度计算曲线的基础上，结合区域钻井地层、岩性特征资料，利用测井曲线突变特征所反映的岩性突变，对研究区目的层已有钻井的标志层进行测井识别和级次划分。
61.在进行级次划分时，按照利用岩性和测井曲线方式进行标志层识别的难易程度，大致将标志层划分为如下四种不同的级次。
62.1)全区分布稳定的明显标志层界面
63.太原组(t)/下古生界奥陶系(o)界面：奥陶系碳酸盐岩顶面为一区域岩溶风化壳不整合面界面，界面上下岩性与常规曲线(如自然伽玛测井曲线、电阻率测井曲线等)以及计算的三孔隙度曲线呈突变特征，标志层界面特征明显(如图5a所示)。
64.石千峰组(p2sh)/上石盒子组(p1s)界面：上石盒子组的低阻、厚层状湖相泥岩为标志层。标志层界面为高阻、中厚层状河流相砂岩与低阻、厚层状湖相泥岩的分界面。计算的三孔隙度曲线、电阻率测井曲线、自然电位测井曲线为岩性特征曲线，易于识别。
65.2)全区分布较稳定、易识别的标志层界面
66.山西组(s)/太原组(t)和山二段(s2)/山一段(s1)。
67.太原组顶部煤层和山一段顶部煤层为标志层。标志层地层界面为河流-三角洲碎屑岩(底部以砂岩为主)与沼泽相厚煤层的界面。在测井曲线上，煤层呈高阻、高声波时差、高中子、低密度“三高一低”的测井曲线响应特征，与上覆地层电性差异明显，为易识别区域标志层(如图5b所示)。
68.3)大部分地区分布较稳定、较易识别的标志层界面
69.下石盒子组(h)/山西组(s)界面：下石盒子组的骆驼脖子砂岩通常为标志层，标志层界面为下石盒子组的骆驼脖子砂岩与山西组的顶面平原沼泽灰黑色泥岩的界面。在界面的上下部分，自然伽玛(gr)、深侧向电阻率(lld)和计算的三孔隙度曲线突变特征也较明显，为较易识别的标志层(如图5c所示)。
70.山1-1小层(s1-1)/山1-2小层(s1-2)小层界面：山1-1小层顶部煤层为标志层。标志层界面为山1-2小层河流-三角洲碎屑岩(底部以砂岩为主)与山1-1小层顶部沼泽相煤层的界面。在测井曲线上，计算的三孔隙度曲线、深侧向电阻率(lld)曲线突变特征也较明显，为较易识别的标志层。
71.4)局部分布较稳定、可识别的标志层界面
72.上石盒子组(p2s)/下石盒子组(h)界面：上石盒子组的厚层状低阻泥岩为标志层。上石盒子组的岩性主要为暗紫色、灰紫色湖湘泥岩夹杂色砂岩；在测井曲线上，自然电位(sp)明显偏正，电阻率偏低。与下伏h3段泥岩或砂岩具有一定的区别，因而可形成一个标志层界面，但在研究区部分地区，界面差异特征不明显。
73.盒二段(h2)/盒一段(h1)界面：盒一段的顶部低阻泥岩为一局部标志层。研究区大
部分地区盒一段的顶部通常发育低阻泥岩层，与上覆盒二段相对高阻的砂岩层，可形成较明显的岩性界面(如图5d所示)。
74.盒一段内部h1-3/h1-2界面：h1-2顶部通常发育低阻泥岩层，为一局部标志层，与上覆h1-3段相对高阻的砂岩层具有较明显的区别，为可识别的局部标志层界面。
75.在步骤s320中，利用标志层划分结果，通过地震合成记录来标定地震剖面，并建立标志层的岩性界面、测井界面和地震界面彼此对应的综合响应特征模式(可记为“岩性界面-测井界面-地震界面的综合响应特征模式”)。
76.具体地，执行井震结合，通过合成地震记录标定的方式，将标志层测井界面标定地震界面，总结标志层的岩性界面-测井界面-地震界面的综合响应特征模式。
77.如图6所示，通过井震结合和标志层界面横向追踪综合分析，研究区有三个标志层界面在地震剖面上可以很好地被识别和区域追踪，分别是全区分布稳定的标志层界面，即太原组/下古生界界面、石千峰组/上石盒子组界面，和研究区大部分地区稳定、较易识别的下石盒子组/山西组的标志层界面，它们分别对应的地震解释层位为t9bc、t8、t9d(参见图7)。有两个局部分布较稳定、可识别的标志层界面在地震剖面上大致可以识别和追踪，分别对应的地震解释层位为t9e、t9f(参见图7)。另外，还有一个大部分地区较稳定、较易识别的标志层在地震剖面上也可局部追踪，对应的地震解释层位为t9c1(参见图7)。
78.在步骤s330中，在标志层的综合响应特征模式的地震响应特征的约束下，开展标志层的井间横向对比。
79.在以上不同级别的标志层岩性-测井-地震识别结果的逐级约束下，利用地震剖面进行标志层井间约束对比，为高频层序和主要目的小层的进一步精细划分提供依据。
80.首先将区域稳定分布的标志层界面在连井剖面上进行对比，然后进行全区较稳定易识别的标志层开展连井对比，以此类推，采用标志层界面逐级约束的方法，开展对比。在以上两级标志层岩性-测井界面对比格架的约束下，以井震标定的结果作为井间横向对比重要的参考依据，横向上的厚度变化参考地震剖面(如图7所示)，开展下一级别标志层的井间对比，避免对比串层。即在对比过程中，利用地震剖面横向上厚度变化的趋势，对测井曲线标志层相对不明显的下石盒子组/山西组、s1-1/s1-2标志层界面和局部分布较稳定可识别的标志层界面上石盒子组/下石盒子组、盒二段/盒一段、h1-3/h1-2小层界面进行约束。图8即为通过井震紧密结合约束后的不同级次标志层界面的对比结果。
81.在步骤s340中，以步骤s330识别的标志层井间横向对比结果作为约束，总结高频沉积旋回的岩性变化趋势以及不同高频沉积旋回接触面上下的岩性突变差异，识别高频层序界面。然后，利用岩心识别的高频层序界面，进一步标定对岩性敏感的测井曲线与重合后的三孔隙度计算曲线，总结高频层序界面的测井响应特征，进而完成全井段的高频层序界面识别。
82.如在步骤s330的低频层序界面约束下，进一步详细观察岩心，对高频层序界面进行识别，本区四级高频层序界面一般为砂岩/泥岩的岩性突变界面、河道冲刷面等，测井曲线(如伽马曲线)和三孔隙度计算曲线的响应突变特征较明显。如以上岩性岩相界面或河道冲刷面等高频层序的界面上下表现为由高伽马突变为低伽马，三孔隙度计算曲线往往由重合状态突变为分开(密度孔隙度易受井径变化影响，结果用于参考)；在每个高频层序内，伽马测井曲线随着岩石中泥质含量增加逐渐升高，三孔隙度计算曲线总体呈重合和分开的特
征，由此构成一个明显向上变细沉积的高频层序(如图9所示，(1)表示冲刷面滞留砾岩与粉细砂岩的界面；(2)表示砂岩与粉砂质泥岩的界面；(3)表示河道冲刷面砂砾岩与下伏泥岩的分界面；(4)表示粗砂岩与粉砂质泥岩的界面；(5)表示河道冲刷面含砾砂岩与下伏泥岩的界面)。
83.在步骤s350中，在步骤s10、s20、s310分析基础上，进一步对主要目的层利用岩石相详细标定测井曲线得到测井相，利用测井相标定地震剖面响应特征得到地震相，进而总结岩石相、测井相和地震相彼此对应的综合响应模式。
84.通过步骤s10对研究区主要目的层盒1段的岩石相分析，该研究区为辫状河，主要包含以下三种微相：辫状河河道心滩、河道充填和泛滥平原。进一步通过岩石相详细标定测井曲线，使岩心—测井紧密结合，具体岩石相和测井相特征如下：
85.河道心滩微相：主要发育于垂向相序的中下部，可见较典型冲刷面及各种类型交错层理(含砾)粗砂岩相等，测井曲线gr呈高幅齿化箱形，三孔隙度计算曲线多呈略有分开或大致重合的特征。
86.河道充填微相：主要发育于垂向相序的中上部，可见中小交错层理砂岩相、砂纹层理细砂岩等，测井曲线gr呈幅齿化箱形和钟形；三孔隙度计算曲线中，声波时差孔隙度和密度孔隙度大致重合，而中子孔隙度略有升高，而与其他两条孔隙度曲线分离，或者三孔隙度曲线呈低值的大致重合特征。
87.泛滥平原微相：主要发育于垂向相序的上部，岩性以泥岩相为主，夹粉砂岩，偶见泥岩变形层理，测井曲线gr呈低幅平直、齿化，曲线异常幅度低；三孔隙度计算曲线的中子孔隙度曲线较声波时差曲线和密度孔隙度曲线呈明显升高，可与其他两条孔隙度曲线可形成明显的包络。
88.这三种微相在垂向上自下而上依次叠置，总体上构成了沉积粒度由粗变细的沉积相序旋回，与向上变浅的四级高频层序类型基本一致。
89.通过以上钻井岩心的岩石沉积构造相对应的测井响应特征，进一步根据测井曲线的幅度、形态、接触关系、组合特征的研究，建立十里加汗地区山西组-下石盒子组主要沉积相的测井相模式(如图11)。利用测井相分析结果，在步骤s320的利用合成地震记录详细标定地震剖面的响应特征的基础上，依据地质-测井-地震的对应关系(模式)对地震相进行了地质属性分析，结合区域沉积背景，建立地震反射结构、振幅属性与相关沉积微相类型的转化模式，进一步总结与岩石相-测井相对应的地震相特征。
90.例如，杭锦旗地区下石盒子组砂岩的测井曲线主要为箱型，少量钟型、指型等。在垂直物源方向，盒1段砂体在地震剖面具有中短轴低频-中强波谷状反射特征。通常情况下，河道砂体因强烈侵蚀下伏泥岩沉积，因而由具有较强的波阻抗反射界面所形成，地震剖面通常呈强振幅亮点反射，主要分布于辫状河河道心滩微相。在顺物源方向，砂体连续性好，t9d之上的亮点反射相对更连续。中弱振幅和弱振幅反射这种地震反射经钻井资料(合成地震记录)均证实分别为河道充填、泛滥平原沉积更为发育。综合以上分析认为，河道心滩微相通常为中强振幅波谷反射，河道充填微相为中弱振幅波谷反射，泛滥平原微相为典型的弱振幅波谷反射，据此，进一步总结了岩石相-测井相-地震相的综合响应模式(如图11所示)。
91.在步骤s40(步骤四)中，利用所述步骤s30的识别结果和标定结果，通过井震结合
建立高频层序格架并以此高频层序格架来划分目的层的目的小层，在测井相和地震相的综合响应特征的分析约束下，分析以目的小层为单元的沉积微相平面展布以及主河道的沉积演化过程。
92.具体地，该步骤包括如下子步骤。
93.在步骤s410中，在辫状河沉积特征及规律的指导下，根据步骤s340识别的高频层序界面建立高频层序格架，并依据砂体分布和实际开发需要，划分目的层的目的小层，并在步骤s330的地震剖面和标志层对比剖面的约束下，开展目的小层的对比。
94.如研究区主要为辫状河-泛滥平原沉积，平原化特征明显，各高频旋回或各小层厚度变化小，因此，高频旋回间横向对比具有一定的稳定性。依据高频层序在地层中的分布和泥岩隔夹层在地层中的分布，将主要目的层太原组-下石盒子组划分为13个小层，并在标志层对比剖面和地震剖面的约束下，开展小层对比(如图10所示的结果)。
95.在步骤s420中，以步骤s350的岩石相、测井相和地震相彼此对应的综合响应模式为沉积微相分析的参考依据，以高频层序格架为约束，完成目的层的单井沉积微相划分和连井沉积微相的对比。
96.具体地，利用步骤s350的综合响应模式下的测井相开展主要目的层单井沉积相划分，将基于步骤s410的高频层序识别的小层划分对比结果作为小层沉积微相分析的基本格架，并以连井地震剖面作为参考和约束条件，开展主要目的层的井间沉积微相对比。
97.例如图12所示，利用研究区的岩相组合所反映的测井相特征，对研究区主要目的层的下石盒子组进行单井相划分，并基于步骤s410的高频层序识别的小层划分对比结果，作为小层沉积微相分析的基本格架，并以单井相对应的井旁地震相模式及其变化趋势为约束，开展主干连井剖面沉积微相对比，从而实现地震相识别与沉积微相的刻画，然后逐条剖面向无井区或少井区，进行追踪、对比和闭合(如图12所示)。
98.在步骤s430中，以步骤s420获得的单井沉积微相划分、连井沉积微相对比为基础，依据步骤s320的井震标定结果，计算每个目的小层厚度在地震剖面上所对应的时窗，利用沿层地震切片技术和该时窗提取能够反映各个目的小层沉积微相的地震相属性，实现小层平面沉积微相展布和河道带主流线的绘制。
99.以单井沉积微相划分、连井沉积微相对比为基础，平面上以各目的小层的地震振幅属性作为沉积微相边界的约束，井震结合，相互标定，编制各目的小层平面沉积微相相图。例如，对于目的层h1段，利用该时窗沿t9d层(h1段的底部)向上提取能够反映各个目的小层沉积微相的地震敏感振幅属性切片，井震结合刻画平面沉积微相的分布特征，并绘制河道带主流线。
100.例如，如图13所示，锦72井区h11小层、h12小层和h13小层平面上自西向东，大致可以划分为6条复合河道带。进一步在复合河道带内，将砂体粒度粗、厚度大、最大波谷振幅响应强的心滩微相，在河道带内顺河道方向用线条依次相连，绘制辫状河道的主流线。
101.在步骤s440中，将目的小层的平面沉积相的河道主流线在平面上叠置，分析总结辫状主河道的迁移摆动特征，总结河流沉积演化过程。
102.如图14所示，自上而下依次将h11小层、h12小层、h13小层的辫状河主流线在平面上叠置，从三个小层的辫状河道主流线的平面分布特征可以看出，河道弯曲度较小，平面分布具有分叉合并，迁移摆动的特征。相对而言，h11小层和h12小层的河道主流线呈小幅的迁
移摆动特征，表明河道总体具有继承性发育的特点，河流的展布形态可能一定程度上受沉积时期古地貌等条件的限制。h13小层辫状河道主流线与h11小层、h12小层河道主流线相比，弯曲度略有增加，摆动特征相对明显，反映了h13小层沉积时期，影响辫状河发育的古地貌、地形坡度等沉积环境有了一定的变化。
103.本发明实施例通过岩心、薄片观察鉴定出样品的岩性，以高频层序地层学、沉积学和岩石学地质理论为指导，将高频层序格架作为约束，地球物理作为剖面沉积相和平面沉积相的参考，利用地球物理沿层切片技术与高频层序相结合，形成了一种辫状河沉积微相与河道迁移摆动沉积演化分析的方法。该方法解决了缺乏露头对比、钻井少、岩性复杂的辫状河碎屑岩地层，难以准确刻画河道边界，以及河道迁移摆动过程沉积演化分析难度大的问题。该方法可操作性强，符合层序地层学、沉积学和地球物理学原理，经中石化已有辫状河岩性气藏应用证实，可以准确描述以高频层序为约束和划分单元的沉积微相展布，并进行辫状河道的迁移摆动演化分析，为辫状河岩性气藏的沉积认识、储层地球物理预测正演模型建立、气藏综合评价等奠定了良好的基础。
104.实施例二
105.图15为本技术实施例的分析河流的沉积微相平面展布与沉积演化的系统的功能结构框图。下面参考图15来说明该系统的功能模块。
106.该系统1500通过利用岩心-测井曲线相互标定，标志层分级约束，井震结合建立高频层序格架，划分小层；在目的小层内部，利用岩心观察总结沉积微相标志，建立岩石相-测井相-地震相的综合响应模式，随之开展单井沉积微相和连井沉积微相对比分析，以目的小层为作图单元，利用地震沿层切片技术提取地震相等属性，综合开展沉积微相平面展布和辫状河道迁移摆动的沉积演化分析。
107.如图15所示，该系统包括第一模块1500a、第二模块1500b、第三模块1500c和第四模块1500d，这些模块通过合理设计能够实现实施例一的对应步骤，即执行步骤s10、s20、s30和s40，因此不再对这些模块的功能进行赘述。
108.应该理解的是，本发明所公开的实施例不限于这里所公开的特定结构、处理步骤，而应当延伸到相关领域的普通技术人员所理解的这些特征的等同替代。还应当理解的是，在此使用的术语仅用于描述特定实施例的目的，而并不意味着限制。
109.说明书中提到的“一个实施例”或“实施例”意指结合实施例描述的特定特征、结构或特性包括在本发明的至少一个实施例中。因此，说明书通篇各个地方出现的短语“一个实施例”或“实施例”并不一定均指同一个实施例。
110.虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所公开的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。