CN114994748A - 古地貌恢复方法和装置 - Google Patents
古地貌恢复方法和装置 Download PDFInfo
- Publication number
- CN114994748A CN114994748A CN202110226941.7A CN202110226941A CN114994748A CN 114994748 A CN114994748 A CN 114994748A CN 202110226941 A CN202110226941 A CN 202110226941A CN 114994748 A CN114994748 A CN 114994748A
- Authority
- CN
- China
- Prior art keywords
- restoration
- thickness
- ancient landform
- layer
- stratum
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 107
- 238000011084 recovery Methods 0.000 claims abstract description 44
- 230000002159 abnormal effect Effects 0.000 claims abstract description 18
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 17
- 239000013049 sediment Substances 0.000 claims abstract description 10
- 230000008021 deposition Effects 0.000 claims description 57
- 238000004458 analytical method Methods 0.000 claims description 16
- 230000002547 anomalous effect Effects 0.000 claims description 15
- 238000005553 drilling Methods 0.000 claims description 7
- 238000000638 solvent extraction Methods 0.000 claims 1
- 239000010410 layer Substances 0.000 description 119
- 238000009826 distribution Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- 230000008859 change Effects 0.000 description 6
- 230000003628 erosive effect Effects 0.000 description 5
- 230000009466 transformation Effects 0.000 description 5
- 230000015556 catabolic process Effects 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 230000002349 favourable effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 229910000514 dolomite Inorganic materials 0.000 description 2
- 239000010459 dolomite Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 206010027336 Menstruation delayed Diseases 0.000 description 1
- 235000010678 Paulownia tomentosa Nutrition 0.000 description 1
- 240000002834 Paulownia tomentosa Species 0.000 description 1
- 238000002679 ablation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000013505 freshwater Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000001502 supplementing effect Effects 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种古地貌恢复方法和装置,属于油气勘探技术领域。所述方法包括:建立目标区域的层序地层框架,对层序地层框架进行精细层序划分;分别利用印模厚度法、残余地层厚度法对目标区域的目的层进行古地貌恢复;当所述目的层含有异常沉积体时,若印模厚度法对应的古地貌恢复趋势与残余地层厚度法对应的古地貌恢复趋势一致或不一致时,分别通过不同的公式确定目的层的相对古地貌恢复值;将相对古地貌恢复值与印模厚度法对应的初始古地貌恢复值相加,确定目的层的最终古地貌恢复值,进行古地貌恢复。本发明可对目的层的古地貌进行准确恢复。
Description
技术领域
本发明涉及油气勘探技术领域,特别涉及一种古地貌恢复方法和装置。
背景技术
古地貌是风化、剥蚀、岩溶作用与多种地质条件综合作用的结果,不同地貌形态对储层的发育有较大的影响。因此,有必要提供一种古地貌恢复方法,以准确有效地划分储层有利相带。
目前,主要采用残余厚度法、印模法进行古地貌恢复。其中,残余厚度法是用侵蚀面至其下伏的水平基准面的残余厚度来反应古地貌;印模法是用侵蚀面至其上覆的水平基准面的厚度来反应古地貌。
但是上述两种方法均没有考虑内部层序是否发生变化以及地层剥蚀的影响,不能准确地对古地貌进行恢复。
发明内容
本发明实施例提供了一种古地貌恢复方法和装置,能够克服上述技术问题。所述技术方案如下:
一方面,提供了一种古地貌恢复方法,所述古地貌恢复方法包括:
建立目标区域的层序地层框架,对所述层序地层框架进行精细层序划分;
分别利用印模厚度法和残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复;
当所述目的层含有异常沉积体时,若所述印模厚度法对应的初始古地貌恢复趋势与所述残余地层厚度法对应的初始古地貌恢复趋势一致,则采用如下公式确定所述目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若不一致时,则采用如下公式获取所述目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H2+Hq2)/4
式中:
H—所述目的层的相对古地貌恢复值,m;
H1max—所述目的层的上覆地层沉积最大厚度值,m;
H1—所述目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—所述目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
将所述相对古地貌恢复值与所述印模厚度法对应的初始古地貌恢复值相加,确定所述目的层的最终古地貌恢复值,进行古地貌恢复。
可选的,所述建立目的层的层序地层框架,包括:
根据所述目标区域的地震剖面资料,建立所述目标区域的三维地层模型;
基于所述三维地层模型,拾取所述目标区域的所有地层面,并提取等时层序地层信息,建立所述层序地层框架。
可选的,所述对所述层序地层框架进行精细层序划分,包括:
采用小波分析法以及时频分析方法,对所述层序地层框架进行精细层序划分。
可选的,所述利用印模厚度法对所述目标区域的目的层进行初始古地貌恢复,包括:
选取所述目的层所对应的上覆地层顶界,并获取所述上覆地层顶界与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
可选的,所述利用残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复,包括:
确定所述目的层的下伏基准面,并获取所述下伏基准面与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
可选的,所述异常沉积体包括前积体、退积体中的至少一种。
可选的,所述方法还包括:
利用钻井资料对所述最终古地貌恢复值进行校准。
另一方面,还提供了一种古地貌恢复装置,所述古地貌恢复装置包括:
划分模块,用于建立目标区域的层序地层框架,对所述层序地层框架进行精细层序划分;
初始古地貌恢复模块,用于分别利用印模厚度法和残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复;
第一确定模块,用于当所述目的层含有异常沉积体时,若所述印模厚度法对应的初始古地貌恢复趋势与所述残余地层厚度法对应的初始古地貌恢复趋势一致,则采用如下公式确定所述目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若不一致时,则采用如下公式获取所述目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H2+Hq2)/4
式中:
H—所述目的层的相对古地貌恢复值,m;
H1max—所述目的层的上覆地层沉积最大厚度值,m;
H1—所述目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—所述目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
第二确定模块,用于将所述相对古地貌恢复值与所述印模厚度法对应的初始古地貌恢复值相加,确定所述目的层的最终古地貌恢复值,进行古地貌恢复。
可选的,所述划分模块,用于:
根据所述目标区域的地震剖面资料,建立所述目标区域的三维地层模型;
基于所述三维地层模型,拾取所述目标区域的所有地层面,并提取等时层序地层信息,建立所述层序地层框架。
可选的,所述划分模块,用于:
采用小波分析法以及时频分析方法,对所述层序地层框架进行精细层序划分。
本发明实施例提供的技术方案带来的有益效果是:
本发明实施例所提供的古地貌恢复方法,通过识别目标区域目的层内是否含有异常沉积体,来考虑内部层序是否发生变化以及地层剥蚀的影响,并同时结合印模厚度法和残余地层厚度法,对目的层的古地貌进行准确恢复,弥补了现有技术所存在的技术问题,使得恢复结果精度高,降低了结果的多解性,实现过程工作量小,可推广性强。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种目的层古地貌的求取原理示意图;
图2是本发明实施例提供的一种利用印膜厚度法形成的灯四段沉积前古地貌的示意图;
图3是本发明实施例提供的一种利用残余地层厚度法形成的灯四段沉积前古地貌的示意图;
图4是本发明实施例提供的一种利用残余地层厚度法形成的灯三段沉积前古地貌的示意图;
图5是本发明实施例提供的一种利用残余地层厚度法形成的灯一、灯二段沉积前古地貌的示意图;
图6是本发明实施例提供的一种利用残余地层厚度法形成的灯一、灯二、灯三段沉积前古地貌的示意图;
图7是本发明实施例提供的一种灯四段内部的沉积体的厚度分布示意图;
图8是本发明实施例提供的一种灯四段沉积前的最终古地貌恢复示意图;
图9是本发明实施例提供的一种古地貌恢复装置的结构示意图;
图10是本发明实施例提供的一种古地貌恢复装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种古地貌恢复方法,该古地貌恢复方法包括:
步骤S1、建立目标区域的层序地层框架,对层序地层框架进行精细层序划分。
步骤S2、分别利用印模厚度法、残余地层厚度法对目标区域的目的层进行古地貌恢复。
步骤S3、当所述目的层含有异常沉积体时,若印模厚度法对应的古地貌恢复趋势与残余地层厚度法对应的古地貌恢复趋势一致时,采用如下公式确定目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若印模厚度法对应的古地貌恢复趋势与残余地层厚度法对应的古地貌恢复趋势不一致时,采用如下公式确定目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H3+Hq2)/4
式中:
H—目的层的相对古地貌恢复值,m;
H1max—目的层的上覆地层沉积最大厚度值,m;
H1—目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
步骤S4、将相对古地貌恢复值与印模厚度法对应的初始古地貌恢复值相加,确定目的层的最终古地貌恢复值,进行古地貌恢复。
其中,上述各个参数之间的关系,可以参见附图1。需要说明的是,利用印模厚度法对目的层的古地貌恢复趋势与利用残余地层厚度法对目的层的古地貌恢复趋势一致是指,针对于目的层同一位置要么均是高势区,要么均是低势区。
可见,本发明实施例所提供的古地貌恢复方法,通过识别目标区域目的层内是否含有异常沉积体,来考虑内部层序是否发生变化以及地层剥蚀的影响,并同时结合印模厚度法和残余地层厚度法,对目的层的古地貌进行准确恢复,弥补了现有技术所存在的技术问题,使得恢复结果精度高,降低了结果的多解性,实现过程工作量小,可推广性强。
下面就本发明实施例所提供的古地貌恢复方法的各个步骤进行描述:
针对于步骤S1,主要是通过层序地层框架的高分辨率层序地层学理论对目的层进行精细层序划分,以对目的层的层序结构进行精确描述,达到能准确判断目的层中是否具有异常沉积体的目的。
需要说明的是,所谓“高分辨率”的实质是指对不同级次地层基准面旋回进行划分和等时对比的高精度时间辨率,即高分辨率的时间—地层单元划分。
其中,可通过如下方法可建立目标区域的层序地层框架:
步骤101、根据目标区域的地震剖面资料,建立目标区域的三维地层模型。
具体为,获取目标区域的地震资料,并根据地震资料中的地震剖面,查找目标区域的已有层位以及断层,进而通过全局等时层序构建方法建立目标区域的三维地层模型。
步骤102、基于三维地层模型,拾取目标区域的所有层序界面以及异常体的等时层序地层信息,从而能够精准建立层序地层框架。
其中,层序地层框架还可以称为等时地层框架。
该步骤可以具体为:在三维地层模型的全三维数据体中拾取目的层的所有层序界面,以及在三维地层模型的地震数据体中提取等时层序地层信息(包括解释不整合面,断层和主要的层位信息),对其进行曲线参数化处理,建立参数空间域(UVT)变换,关联三维地层模型空间域中每一个(x,y,z)点到上述空间域(UVT)中的一个(u,v,t)点。
其中,全三维数据体是包含层序地层信息的地震数据体,通过全局等时层序构建方法,拾取目标区域层序界面信息后,建立三位地层模型,便可以提取目标层序信息。
其中,地震数据体包含振幅、频率和波形等地震波反射的基本信息。
之后,为了提高层序划分的精度,采用小波分析法以及时频分析方法,对所述目标区域进行精细层序划分。
具体为:根据目标区域的地震数据以及地震解释数据,对目标区域的地震剖面进行时频分析,并根据目标区域的测井资料,进行小波变换,待小波变换之后,将小波变换的结果对时频分析划分层序的结果进行标定和修正,实现层序的精细划分。
其中,地震资料(即包括地震数据与地震解释数据)在横向上具有良好的连续性,能反映沉积特征的特点,但纵向分辨率差;测井资料具有较强的纵向分辨率,但是与地震剖面的吻合程度差。基于地震资料与测井资料的特点,本发明实施例采用小波变换与时频分析相结合的方法进行层序划分,不仅可提高与地震剖面的吻合程度,也达到了对层序进行精细划分的目的。
针对于步骤S2,分别利用印模厚度法、残余地层厚度法对目标区域的目的层进行初始古地貌恢复。
其中,利用印模厚度法对目标区域的目的层进行初始古地貌恢复,包括:
步骤S21、选取目的层所对应的上覆地层顶界,并获取上覆地层顶界与目的层顶界之间的厚度(即附图1所示的H1)。
其中,上覆地层顶界的选取需要满足以下几个条件:能够全区对比追踪、等时性强且沉积前没有大的构造运动。
步骤S22、根据厚度,对目的层进行古地貌形态描述。
其中,通过上覆地层顶界与目的层顶界之间厚度的镜像关系来表征古地貌形态,厚度小的区域代表古地貌高势区,厚度大的区域代表古地貌低势区。
另外,利用残余地层厚度法对目标区域的目的层进行初始古地貌恢复,包括:
步骤S23、确定目的层的下伏基准面,并获取下伏基准面与目的层顶界之间的厚度(即附图1所示的H2、H3)。
其中,下伏基准面可以为目的层底界未遭受大量剥蚀且未经历剧烈的构造运动的可全区对比追踪的界面。
步骤S24、根据厚度,对目的层进行古地貌形态描述。
具体为,通过下伏基准面与目的层顶界之间的厚度,来表征古地貌形态,厚度小的区域代表古地貌高势区,厚度大的区域代表古地貌低势区。
其中,目的层的下伏地层厚度变化对目的层沉积前古地貌存在一定的继承性。因此,可结合上一步印模法古地貌恢复结果来分析其与下伏地层厚度之间趋势的差异和吻合之处。
针对于步骤S3而言,异常沉积体包括前积体、退积体中的至少一种。其中,尽管在一个时期的古地貌上局部前积或退积形成的沉积体在后期会因地质构造运动而使得地形发生变化直至演化成现今的构造,但是其上的沉积体依然会保留有其在沉积时期基本一致的分布或形态。换句话说,在现今地震剖面上发现的异常沉积体(参见附图1),如前积体,它的分布范围和厚度变化受该沉积时期的古地貌影响较大,因此有必要考虑这些异常沉积体对目的层的古地貌恢复的影响。
针对于步骤S4而言,可具体为,利用对目的层的古地貌影响较大的下伏地层厚度和异常沉积体厚度变化趋势对目的层的初始古地貌进行约束,即直接将印模厚度法对应的初始古地貌恢复值加上相对古地貌恢复值,确定目的层的最终古地貌恢复值。
进一步地,为了保证目的层的古地貌恢复值的精度,本发明实施例中,还利用钻井资料对目的层的最终古地貌恢复值进行校准。具体为,利用钻井资料对目的层的古地貌进行恢复,然后将该恢复结果与本发明实施例所获取的恢复结果进行比对,若这两者的恢复趋势不一致的话,重新执行步骤S1至步骤S3,直至这两者恢复趋势一致。
综上所述,本发明实施例所提供古地貌恢复方法,一方面,同时考虑了上覆地层厚度和下伏地层厚度对目的层沉积的影响;另一方面将目的层的地震地层接触关系、厚度变化、沉积间断等地质信息从层序地层框架中加以量化,对局部不整合沉积体系进行雕刻成像,将其厚度变化趋势对前期结果进行约束,避免了仅仅利用残余地层厚度法和印模厚度法进行古地貌恢复无法兼顾地层内沉积变化的缺点。
本发明实施例以高石梯-磨溪地区灯影组为例,并利用上述所提供的古地貌恢复方法进行古地貌恢复。
其中,该地区的深层(埋藏深度5000~6000m)古老白云岩岩溶且存在地层剥蚀,且该地区目的层局部还存在前积体的异常沉积体。该高石梯-磨溪地区根据目前的地层划分方案,震旦系灯影组可划分为4段,自下而上分别为:灯影组一段、灯影组二段、灯影组三段、灯影组四段。其中,灯影组四段(以下简称灯四段)储集层是丘滩相白云岩受桐湾Ⅱ幕大气淡水岩溶改造而成,且相控岩溶作用是优质储集层形成的关键,故将灯四段作为目的层,进行古地貌恢复。
首先,建立高石梯-磨溪地区的三维地层模型,并拾取区域的所有地层层面信息以及等时层序地层信息,建立层序地层框架,对层序地层框架进行精细层序划分。
其次,利用印模厚度法进行初始古地貌恢复,其恢复结构参见附图2。其中,首先优选合理的上覆地层顶界,具体为,在震旦系灯影组上覆的寒武系筇竹寺组—沧浪铺组沉积时期为一个完整的海侵—海退旋回,主要为补偿沉积,对震旦系灯影组受剥蚀古地貌基本填平补齐,并且筇竹寺组—沧浪铺组沉积晚期乐山—龙女寺古隆起区构造运动相对稳定;而震旦系灯影组沉积时期随着海平面的下降,台地边界整体向东迁移,虽然盆地西部遭受了强烈剥蚀,并且不同位置剥蚀量有差异,但总的来看表层剥蚀厚度基本在0~50m之间,而灯四段残余地层厚度在270~340m之间,因此,震旦系灯四段底—沧浪铺组顶(即龙王庙底)的印模厚度能基本反映灯四段沉积前的古地貌特征,但是精确度不够。加之该地区内高品质地震资料三维连片面覆盖,因此,选取沧浪铺组顶界作为上覆地层顶界,与灯四段底界的印模厚度趋势来表征灯四段沉积前的古地貌是可行的。
然后,形成灯三段、灯一、二、三段和灯一、二段的残余地层厚度图(参见附图3、附图4、附图5、附图6)。由于灯四段下伏地层灯一、二段沉积较稳定,后期遭受剥蚀量较小,地层厚度较大,普遍在400m以上;而灯三段是离灯四段最近的下伏地层,沉积较稳定,后期遭受剥蚀量较小,该层厚度小,普遍在20-80m之间。同时,灯三段与灯四段之间剥蚀程度低,灯四段沉积受其下伏地层古地貌影响较大。因此,灯四段下伏地层厚度变化对灯四段沉积前古地貌存在一定的继承性。结合上一步印模法古地貌恢复结果来分析其与下伏地层厚度之间趋势的差异和吻合之处。
其次,在三维地层模型的基础上识别出灯四段内部存在的前积沉积层,其识别结果如附图6所示,刻画其分布范围并形成前积层厚度图。这些前积沉积体在磨溪地区灯四段内部普遍存在由西向东的超覆现象(参见附图1),主要是在沉积期受海水进退和古地貌影响形成的进积层,前积沉积层的结构、分布范围和厚度变化直接反应了沉积期的地貌差异。
另外,研究该地区古地貌对古岩溶分布的控制作用较明显,岩溶斜坡为最有利于古岩溶发育的部位,这些部位水动力较强,交替活跃,致使水平和垂直形态的岩溶普遍发育,形成叠加的多套水平溶洞层相。岩心分析显示,灯四段大量发育溶蚀孔洞,在镜下薄片上也见到粒间溶蚀孔洞,该地区在古地貌斜坡区域取心段溶洞所占比例大同时洞密度多大于30个/m。本区钻井资料较丰富,完钻井已超过70口,在大量实测数据面前,可以反向刻画岩溶古地貌发育的基本规律。
最后,在以上条件建立之后,当印模法与残余地层厚度法古地貌恢复结果趋势变化一致时,则用公式一恢复古地貌;在印模厚度与灯四段残余地层厚度变化趋势背离时,则重点参考异常沉积体厚度变化趋势和实测岩溶分布规律,运用公式二对符合这一规律的结果进行加权。最后以完钻井古地貌恢复值(参见表1)对恢复结果进行标定,获得最终古地貌恢复值(参见附图8)。
表1该区完钻井古地貌恢复值统计表
井名 | 上覆沉积层厚度值(m) | 古地貌恢复值(m) |
MX8 | 353.13 | 246.87 |
MX9 | 358.28 | 241.72 |
MX10 | 464.84 | 135.16 |
GS18 | 427.04 | 172.96 |
GS19 | 341.60 | 258.40 |
GS21 | 231.01 | 368.99 |
GS101 | 325.33 | 274.67 |
GS103 | 331.58 | 268.42 |
GS105 | 334.94 | 265.06 |
GS109 | 442.89 | 157.11 |
GS111 | 296.76 | 303.24 |
… | … | … |
其中,表1中第二列是钻井古地貌的厚度值,第三列为最终古地貌恢复值。
可见,本发明实施例提供的方法,便于地质研究,将灯四沉积前古地貌划分为古地貌高地、古地貌陡坡、古地貌缓坡和古地貌洼地(参见附图8),分析灯四沉积前后古地貌与灯四段岩溶储集层发育分布的关系,预测出古地貌缓坡为岩溶储层发育优势区,解决了深层古老碳酸盐岩岩溶储层有利储集区划分的一大难题。同时该成果可指导高石梯-磨溪地区的开发井位部署工作,操作人员在有利于岩溶储层发育的古地貌分布区(古地貌缓坡区)部署4口开发井,经过实钻均位于古地貌缓坡区,同时经过测试其无阻流量均超过100×104m3/d,实现了开发井高产的目标。因此通过该本发明实施例得到的灯四段沉积前后古地貌成果,为研究该区储层分布及岩溶分布创造了有利条件,同时也为研究区气田井位部署和增储上产提供了重要的参考依据。
本实施例还提供了一种古地貌恢复装置,如图9所示,所述古地貌恢复装置包括:
划分模块910,用于建立目标区域的层序地层框架,对所述层序地层框架进行精细层序划分;
初始古地貌恢复模块920,用于分别利用印模厚度法和残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复;
第一确定模块930,用于当所述目的层含有异常沉积体时,若所述印模厚度法对应的初始古地貌恢复趋势与所述残余地层厚度法对应的初始古地貌恢复趋势一致,则采用如下公式确定所述目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若不一致时,则采用如下公式获取所述目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H2+Hq2)/4
式中:
H—所述目的层的相对古地貌恢复值,m;
H1max—所述目的层的上覆地层沉积最大厚度值,m;
H1—所述目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—所述目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
第二确定模块940,用于将所述相对古地貌恢复值与所述印模厚度法对应的初始古地貌恢复值相加,确定所述目的层的最终古地貌恢复值,进行古地貌恢复。
可选的,所述划分模块910,用于:
根据所述目标区域的地震剖面资料,建立所述目标区域的三维地层模型;
基于所述三维地层模型,拾取所述目标区域的所有地层面,并提取等时层序地层信息,建立所述层序地层框架。
可选的,所述划分模块910,用于:
采用小波分析法以及时频分析方法,对所述层序地层框架进行精细层序划分。
可选的,所述初始古地貌恢复模块920,用于:
选取所述目的层所对应的上覆地层顶界,并获取所述上覆地层顶界与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
可选的,所述初始古地貌恢复模块920,用于:
确定所述目的层的下伏基准面,并获取所述下伏基准面与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
可选的,所述异常沉积体包括前积体、退积体中的至少一种。
可选的,如图10所示,所述装置还包括校准模块950,用于利用钻井资料对所述最终古地貌恢复值进行校准。
需要说明的是:上述实施例提供的古地貌恢复装置在古地貌恢复时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将装置的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的古地貌恢复装置与古地貌恢复方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
上述所有可选技术方案,可以采用任意结合形成本公开的可选实施例,在此不再一一赘述。
以上所述仅为本发明的一种实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种古地貌恢复方法,其特征在于,所述古地貌恢复方法包括:
建立目标区域的层序地层框架,对所述层序地层框架进行精细层序划分;
分别利用印模厚度法和残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复;
当所述目的层含有异常沉积体时,若所述印模厚度法对应的初始古地貌恢复趋势与所述残余地层厚度法对应的初始古地貌恢复趋势一致,则采用如下公式确定所述目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若不一致时,则采用如下公式获取所述目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H2+Hq2)/4
式中:
H—所述目的层的相对古地貌恢复值,m;
H1max—所述目的层的上覆地层沉积最大厚度值,m;
H1—所述目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—所述目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
将所述相对古地貌恢复值与所述印模厚度法对应的初始古地貌恢复值相加,确定所述目的层的最终古地貌恢复值,进行古地貌恢复。
2.根据权利要求1所述的古地貌恢复方法,其特征在于,所述建立目的层的层序地层框架,包括:
根据所述目标区域的地震剖面资料,建立所述目标区域的三维地层模型;
基于所述三维地层模型,拾取所述目标区域的所有地层面,并提取等时层序地层信息,建立所述层序地层框架。
3.根据权利要求1所述的古地貌恢复方法,其特征在于,所述对所述层序地层框架进行精细层序划分,包括:
采用小波分析法以及时频分析方法,对所述层序地层框架进行精细层序划分。
4.根据权利要求1所述的古地貌恢复方法,其特征在于,所述利用印模厚度法对所述目标区域的目的层进行初始古地貌恢复,包括:
选取所述目的层所对应的上覆地层顶界,并获取所述上覆地层顶界与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
5.根据权利要求1所述的古地貌恢复方法,其特征在于,所述利用残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复,包括:
确定所述目的层的下伏基准面,并获取所述下伏基准面与所述目的层顶界之间的厚度;
根据所述厚度,对所述目的层进行初始古地貌恢复。
6.根据权利要求1所述的古地貌恢复方法,其特征在于,所述异常沉积体包括前积体和退积体中的至少一种。
7.根据权利要求1所述的古地貌恢复方法,其特征在于,所述方法还包括:
利用钻井资料对所述最终古地貌恢复值进行校准。
8.一种古地貌恢复装置,其特征在于,所述古地貌恢复装置包括:
划分模块,用于建立目标区域的层序地层框架,对所述层序地层框架进行精细层序划分;
初始古地貌恢复模块,用于分别利用印模厚度法和残余地层厚度法对所述目标区域的目的层进行初始古地貌恢复;
第一确定模块,用于当所述目的层含有异常沉积体时,若所述印模厚度法对应的初始古地貌恢复趋势与所述残余地层厚度法对应的初始古地貌恢复趋势一致,则采用如下公式确定所述目的层的相对古地貌恢复值:
H=(H1max-H1+H2+Hq2)/3,
若不一致时,则采用如下公式获取所述目的层的相对古地貌恢复值:
H=(H1max-H1+Hq1+H2+Hq2)/4
式中:
H—所述目的层的相对古地貌恢复值,m;
H1max—所述目的层的上覆地层沉积最大厚度值,m;
H1—所述目的层中待恢复位置的上覆地层沉积厚度值,m;
H2—所述目的层的下伏地层沉积厚度,m;
Hq1—目的层上覆异常沉积体的沉积厚度,m;
Hq2—目的层下伏异常沉积体的沉积厚度,m;
第二确定模块,用于将所述相对古地貌恢复值与所述印模厚度法对应的初始古地貌恢复值相加,确定所述目的层的最终古地貌恢复值,进行古地貌恢复。
9.根据权利要求8所述的古地貌恢复装置,其特征在于,所述划分模块,用于:
根据所述目标区域的地震剖面资料,建立所述目标区域的三维地层模型;
基于所述三维地层模型,拾取所述目标区域的所有地层面,并提取等时层序地层信息,建立所述层序地层框架。
10.根据权利要求8所述的古地貌恢复装置,其特征在于,所述划分模块,用于:
采用小波分析法以及时频分析方法,对所述层序地层框架进行精细层序划分。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110226941.7A CN114994748A (zh) | 2021-03-01 | 2021-03-01 | 古地貌恢复方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110226941.7A CN114994748A (zh) | 2021-03-01 | 2021-03-01 | 古地貌恢复方法和装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114994748A true CN114994748A (zh) | 2022-09-02 |
Family
ID=83018370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110226941.7A Pending CN114994748A (zh) | 2021-03-01 | 2021-03-01 | 古地貌恢复方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114994748A (zh) |
-
2021
- 2021-03-01 CN CN202110226941.7A patent/CN114994748A/zh active Pending
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104765065B (zh) | 一种确定沉积地层古地貌的方法及装置 | |
CN111766630B (zh) | 一种盆地古地貌的恢复方法 | |
CN107589470B (zh) | 小区域范围内岩溶古地貌恢复方法 | |
O'Byrne et al. | Sequence, parasequence, and intraparasequence architecture of the Grassy Member, Blackhawk Formation, Book Cliffs, Utah, USA | |
CN112394404B (zh) | 一种渐进式储层精细表征方法 | |
CN112327355B (zh) | 基于源-汇系统的古地貌恢复方法 | |
Back et al. | Growth faults above shale–Seismic-scale outcrop analogues from the Makran foreland, SW Pakistan | |
CN111624651B (zh) | 基于古地貌约束的储层预测方法及装置 | |
CN107167853A (zh) | 一种精细恢复岩溶古地貌的方法 | |
CN104076394A (zh) | 一种碳酸盐岩潜山洞穴分层解释方法 | |
CN112130210A (zh) | 基于古今地质要素的井位优选方法及装置 | |
Cabello et al. | Modelling facies belt distribution in fan deltas coupling sequence stratigraphy and geostatistics: The Eocene Sant Llorenç del Munt example (Ebro foreland basin, NE Spain) | |
CN114994748A (zh) | 古地貌恢复方法和装置 | |
CN111638552A (zh) | 一种古地貌恢复方法 | |
CN114578431A (zh) | 一种古地貌的岩性恢复方法 | |
CN112052591A (zh) | 一种储层构型约束下的隔夹层精细刻画及嵌入式建模方法 | |
CN113885096A (zh) | 高频层序划分与小层对比方法、装置、电子设备及介质 | |
Maynard et al. | Appraisal of a complex, platform carbonate, Bukit Tua discovery, Ketapang PSC, East Java Basin, Indonesia | |
CN114076982B (zh) | 一种基于波形特征差异的岩溶古地貌恢复方法及装置 | |
CN112099085B (zh) | 一种基于转换波的碳酸盐岩微古地貌雕刻及储层预测方法 | |
CN118050819A (zh) | 一种岩溶古地貌图的确定方法 | |
CN117741802A (zh) | 一种基于岩性和相序的混合沉积古地貌恢复方法 | |
Zhang et al. | Identification of the subaerial unconformity, genesis, and distribution in Lower Triassic Baikouquan Formation, Mahu Sag, Junggar Basin, Western China | |
CN115542421A (zh) | 一种不连续分布剥蚀残余灰岩的岩溶古地貌恢复方法 | |
CN114114404A (zh) | 一种用于挤压盆地的古地貌恢复方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |