CN113050164B - 分析沉积序列的方法、装置、服务器和存储介质 - Google Patents
分析沉积序列的方法、装置、服务器和存储介质 Download PDFInfo
- Publication number
- CN113050164B CN113050164B CN201911365712.2A CN201911365712A CN113050164B CN 113050164 B CN113050164 B CN 113050164B CN 201911365712 A CN201911365712 A CN 201911365712A CN 113050164 B CN113050164 B CN 113050164B
- Authority
- CN
- China
- Prior art keywords
- dividing
- phase axis
- indication
- stratum
- deposition
- 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.)
- Active
Links
- 230000008021 deposition Effects 0.000 title claims abstract description 167
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 59
- 238000010183 spectrum analysis Methods 0.000 claims description 6
- 238000000638 solvent extraction Methods 0.000 claims description 5
- 208000035126 Facies Diseases 0.000 abstract description 13
- 238000011161 development Methods 0.000 abstract description 4
- 238000005553 drilling Methods 0.000 abstract description 4
- 230000002349 favourable effect Effects 0.000 abstract description 2
- 238000005755 formation reaction Methods 0.000 description 54
- 238000010586 diagram Methods 0.000 description 10
- 238000012545 processing Methods 0.000 description 7
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- 230000009286 beneficial effect Effects 0.000 description 4
- 230000003628 erosive effect Effects 0.000 description 4
- 239000013049 sediment Substances 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000000654 additive Substances 0.000 description 2
- 230000000996 additive effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000002679 ablation Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- -1 gravel Substances 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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/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
技术领域
本公开涉及地震勘探领域，特别涉及一种分析沉积序列的方法、装置、服务器和存储介质。
背景技术
水流中所夹带的岩石、砂砾、泥土等在河床和海湾等低洼地带沉淀、淤积，这样沉下来的物质形成沉积序列。沉积序列包括堆叠在一起的多个地质层位，每个层位对应于在预定义的地质时间段(从几千年到几百万年)内沉积的物质。分析沉积序列，有利于深入了解地质特征，提高钻探成功率和扩大勘探开发成果。
发明内容
本公开实施例提供了一种分析沉积序列的方法、装置、服务器和存储介质，可以利用地震剖面形成沉积序列。所述技术方案如下：
第一方面，本公开实施例提供了一种分析沉积序列的方法，所述方法包括：
获取地震剖面图像；
在所述地震剖面图像上确定地层接触关系的指示特征；
根据所述指示特征，在竖直方向上划分地层单元；
根据所述指示特征，在水平方向上划分沉积相带；
根据所述地层单元和所述沉积相带，生成沉积序列图像。
可选地，所述获取地震剖面图像，包括：
对地震数据进行频谱分析，得到所述地震数据的频率分布范围；
去除所述频率分布范围内的低频部分，得到所述地震剖面图像，所述低频部分的下限值为所述频率分布范围的下限值，所述低频部分的上限值与下限值之差为所述频率分布范围的上限值与下限值之差的1/5～1/3。
可选地，在所述地震剖面图像上确定地层接触关系的指示特征，包括以下步骤中的至少一种：
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的上倾方向迁移，确定所述指示特征为上超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的下倾方向迁移，确定所述指示特征为下超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴没有相对所述第一同相轴迁移，确定所述指示特征为平行；
若所述地震剖面图像中同相轴的两端均延伸到所述地震剖面图像的边缘，确定所述指示特征为削截。
可选地，所述根据所述指示特征，在竖直方向上划分地层单元，包括：
将所述指示特征相同的相邻地层划分到一个地层单元；
将所述指示特征不同的相邻地层划分到不同的地层单元。
可选地，所述根据所述指示特征，在水平方向上划分沉积相带，包括：
将所述指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
将没有所述指示特征且靠近断层线的地层划分到扇根亚相；
将所述指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
可选地，所述方法还包括：
根据所述指示特征，确定所述地层单元的沉积方向。
第二方面，本公开实施例提供了一种分析沉积序列的装置，所述装置包括：
图像获取模块，用于获取地震剖面图像；
关系确定模块，用于在所述地震剖面图像上确定地层接触关系的指示特征；
单元划分模块，用于根据所述指示特征，在竖直方向上划分地层单元；
相带划分模块，用于根据所述指示特征，在水平方向上划分沉积相带；
图像建立模块，用于根据所述地层单元和所述沉积相带，生成沉积序列图像。
可选地，所述图像获取模块包括：
分析子模块，用于对地震数据进行频谱分析，得到所述地震数据的频率分布范围；
去除子模块，用于去除所述频率分布范围内的低频部分，得到所述地震剖面图像，所述低频部分的下限值为所述频率分布范围的下限值，所述低频部分的上限值与下限值之差为所述频率分布范围的上限值与下限值之差的1/5～1/3。
可选地，所述关系确定模块用于，
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的上倾方向迁移，确定所述指示特征为上超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的下倾方向迁移，确定所述指示特征为下超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴没有相对所述第一同相轴迁移，确定所述指示特征为平行；
若所述地震剖面图像中同相轴的两端均延伸到所述地震剖面图像的边缘，确定所述指示特征为削截。
可选地，所述单元划分模块用于，
将所述指示特征相同的相邻地层划分到一个地层单元；
将所述指示特征不同的相邻地层划分到不同的地层单元。
可选地，所述相带划分模块包括：
扇中划分子模块，用于将所述指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
扇根划分子模块，用于将没有所述指示特征且靠近断层线的地层划分到扇根亚相；
扇端划分子模块，用于将所述指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
可选地，所述装置还包括：
方向确定模块，用于根据所述指示特征，确定所述地层单元的沉积方向。
第三方面，本公开实施例提供了一种服务器，所述服务器包括：
处理器；
被配置为存储所述处理器的可执行指令的存储器；
其中，所述处理器被配置为执行第一方面提供的方法。
第四方面，本公开实施例提供了一种计算机可读存储介质，当所述计算机可读存储介质中的指令由服务器的处理器执行时，使得所述服务器能够执行第一方面提供的方法。
本公开实施例提供的技术方案带来的有益效果是：
通过在地震剖面图像上确定出地层接触关系的指示特征，进而根据指示特征，在竖直方向上进行地层单元的划分，在水平方向上进行沉积相带的划分，再根据划分出的地层单元和沉积相带，即可生成沉积序列图像，在地震剖面的基础上实现沉积序列的分析，有利于深入了解地质特征，提高钻探成功率和扩大勘探开发成果。
附图说明
为了更清楚地说明本公开实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本公开的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本公开实施例提供的一种分析沉积序列的方法的流程图；
图2是本公开实施例提供的上超的示意图；
图3是本公开实施例提供的下超的示意图；
图4是本公开实施例提供的平行的示意图；
图5是本公开实施例提供的削截的示意图；
图6是本公开实施例提供的一种分析沉积序列的方法的流程图；
图7是本公开实施例提供的某一地区地震数据的频谱图；
图8是本公开实施例提供的频谱范围截取的示意图；
图9是本公开实施例提供的频谱范围截取前后地震剖面图像的对比图；
图10是本公开实施例提供的地震剖面图像上地层之间的接触关系的示意图；
图11是本公开实施例提供的地震剖面图像上沉积相带的示意图；
图12是本公开实施例提供的沉积序列图像的示意图；
图13是本公开实施例提供的一种分析沉积序列的装置的结构示意图；
图14是本公开实施例提供的一种服务器的结构示意图。
具体实施方式
为使本公开的目的、技术方案和优点更加清楚，下面将结合附图对本公开实施方式作进一步地详细描述。
本公开实施例提供了一种分析沉积序列的方法。图1为本公开实施例提供的一种分析沉积序列的方法的流程图。参见图1，该方法包括：
步骤101：获取地震剖面图像。
在本实施例中，地震剖面图像是一张标示测线的地震资料图。采用地震勘探技术在地表以人工方法激发地震波向地下传播，地震波在传播过程中遇到介质性质不同的岩层分界面会反生反射与折射，反射或折射的地震波被地表或井中的检波器接收并记录下来，经过一定的处理即可形成地震剖面图像。在地震剖面图像中，黑线为各道反射波振动相位相同的极值(波峰或波谷)的连线，即同相轴。
步骤102：在地震剖面图像上确定地层接触关系的指示特征。
在本实施例中，地层接触关系的指示特征可以包括上超、下超、平行和削截中的至少一种。
图2为本公开实施例提供的上超的示意图。参见图2，最下方的同相轴是沉积面a，沉积面a呈倾斜状，沉积面a的左侧向上倾斜，左侧的延伸方向为沉积面a的上倾方向；沉积面a的右侧向下倾斜，右侧的延伸方向为沉积面a的下倾方向。地层b、地层c、地层d、地层e依次层叠在沉积面a上，并且沿沉积面的上倾方向逐层超覆尖灭，指示特征为上超，表示水域不断扩大，向地势较高处迁移。
图3为本公开实施例提供的下超的示意图。参见图3，最下方的同相轴是沉积面a，沉积面a呈倾斜状，沉积面a的左侧向上倾斜，左侧的延伸方向为沉积面a的上倾方向；沉积面a的右侧向下倾斜，右侧的延伸方向为沉积面a的下倾方向。地层b、地层c、地层d、地层e依次层叠在沉积面a上，并且沿沉积面的下倾方向作底部超覆，指示特征为下超，表示水域不断缩小，向地势较低处迁移。
图4为本公开实施例提供的平行的示意图。参见图4，最下方的同相轴是沉积面a，沉积面a呈倾斜状，沉积面a的左侧向上倾斜，左侧的延伸方向为沉积面a的上倾方向；沉积面a的右侧向下倾斜，右侧的延伸方向为沉积面a的下倾方向。地层b、地层c、地层d、地层e依次层叠在沉积面a上，并且地层b、地层c、地层d、地层e近似平行，指示特征为平行，表示水平面相对静止。
图5本公开实施例提供的削截的示意图。参见图5，最下方的同相轴是沉积面a，沉积面a呈倾斜状，沉积面a的左侧向上倾斜，左侧的延伸方向为沉积面a的上倾方向；沉积面a的右侧向下倾斜，右侧的延伸方向为沉积面a的下倾方向。地层b、地层c、地层d、地层e依次层叠在沉积面a上，并且地层b、地层c、地层d、地层e的顶部在同一平面上，指示特征为削截，是由侵蚀作用或者构造作用引起的顶部反射终止，表示一种侵蚀作用，说明在下伏地层沉积之后，经过强烈的构造运动或者强烈的切割侵蚀。
在实际应用中，在介质性质不同的岩层分界面，地震波会反生反射与折射，在地震剖面图像上表现为强反射的黑线，这些黑线称为同相轴。相邻两条同相轴之间的位置关系，可以确定地层接触关系的指示特征。
步骤103：根据地层接触关系的指示特征，在竖直方向上划分地层单元。
在实际应用中，地层接触关系的指示特征，可以表示地层的沉积方向，从而划分出地层单元。地层单元为沉积序列中的地质层位，每个层位对应于在预定义的地质时间段(从几千年到几百万年)内沉积的物质。
步骤104：根据地层接触关系的指示特征，在水平方向上划分沉积相带。
在实际应用中，地层接触关系的指示特征表明了其地震相特征，利用地震相特征和沉积相特征之间的对应关系，可以根据地层接触关系的指示特征，划分沉积相带。
步骤105：根据地层单元和沉积相带，生成沉积序列图像。
在本实施例中，沉积序列图像是表示地层沉积顺序的图像。基于在地震剖面图像上划分出的地层单元和沉积相带，可以绘出相应的沉积序列图像。
本公开实施例通过在地震剖面图像上确定出地层接触关系的指示特征，进而根据指示特征，在竖直方向上进行地层单元的划分，在水平方向上进行沉积相带的划分，再根据划分出的地层单元和沉积相带，即可生成沉积序列图像，在地震剖面的基础上实现沉积序列的分析，有利于深入了解地质特征，提高钻探成功率和扩大勘探开发成果。
本公开实施例提供了一种分析沉积序列的方法，图6为本公开实施例提供的一种分析沉积序列的方法的流程图。参见图6，该方法包括：
步骤201：对地震数据进行频谱分析，得到地震数据的频率分布范围。
在实际应用中，地震数据是采用地震勘探技术在地表以人工方法激发地震波向地下传播，地震波在传播过程中遇到介质性质不同的岩层分界面反生反射与折射，并被地表或井中的检波器接收和记录下来的数据。地震数据承载的是地震波的相关信息，具有一定的频率范围。
图7为本公开实施例提供的某一地区地震数据的频谱图。参见图7，示例性地，地震数据的频率分布范围约为5Hz～55Hz。
步骤202：去除频率分布范围内的低频部分，得到地震剖面图像。
在实际应用中，可以采用中高通滤波器滤除地震数据中的低频部分。
可选地，低频部分的下限值为频率分布范围的下限值，低频部分的上限值与下限值之差可以为频率分布范围的上限值与下限值之差的1/5～1/3。通过滤除较低频率的数据，可以筛选出主要信息，有利于提高地震剖面图像的分辨率，方便从中找到需要的信息。
图8为本公开实施例提供的截取频率分布范围的示意图。参见图8，以虚线为界，去除从频率分布范围的起点至虚线之间的频段，去除的频率范围约为5Hz～20Hz；留下从虚线至频率分布范围终点之间的频段，留下的频率范围约为20Hz～55Hz。
在实际应用中，地震数据的频率分布范围可以按照频率高低，划分为高频、中频和低频，或者高频和低频。低频段的频率低，携带的信息多且杂，导致图像分辨率较低。从地震数据中去除低频部分形成地震剖面图像，可以提高地震剖面图像的分辨率，方便地层接触关系指示特征的识别。
图9为本公开实施例提供的截取前后地震剖面图像的对比图。参见图9，将截取之后的地震剖面图像与截取之前的地震剖面图像相比，虚线框所围的区域中的同相轴明显更加清楚、连续、分辨率高，有利于后续识别的便利性和准确性。
步骤203：在地震剖面图像上确定地层接触关系的指示特征。
在实际应用中，可以由人在地震剖面图像上识别同相轴，并根据同相轴之间的位置关系，分析出地层接触关系的指示特征；也可以利用图像处理技术从地震剖面图像中提取图像特征进行处理，得到地层接触关系的指示特征。
可选地，在步骤203之前，该方法还可以包括：
在地震剖面图像上识别断层线；
以断层线为分界线，将地震剖面图像划分为至少两个区域；
在目标区域的地震剖面图像中，确定地层接触关系的指示特征，目标区域为至少两个区域中的任一个。
通过识别断层线，减小识别区域，并且可以辅助确定地层接触关系的指示特征、以及划分沉积相带。
可选地，该步骤203可以包括以下步骤中的至少一种：
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴相对第一同相轴朝沉积面的上倾方向迁移，确定指示特征为上超；
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴相对第一同相轴朝沉积面的下倾方向迁移，确定指示特征为下超；
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴没有相对第一同相轴迁移，确定指示特征为平行；
若地震剖面图像中同相轴的两端均延伸到地震剖面图像的边缘，确定指示特征为削截。
其中，沉积面的上倾方向为沉积面向上倾斜的方向，如图2至图5中沉积面a左侧的延伸方向；沉积面的下倾方向为沉积面向下倾斜的方向，如图2至图5中沉积面a右侧的延伸方向。
结合各类指示特征的定义，利用地震剖面图像上相邻两条同相轴之间的位置关系，将地震剖面图像上的同相轴与地层接触关系的指示特征对应，以便从地震剖面图像上确定海岸迁移规律。示例性地，指示特征为上超，则表示海平面上升，向陆地迁移；指示特征为下超，则表示海平面下降，向盆地中央方向迁移；指示特征为平行，则表示海平面相对静止。另外，利用削截特别代表的侵蚀作用造成层序的顶部反射终止，有利于准确划分层序。
图10为本公开实施例提供的地震剖面图像上地层接触关系的指示特征的示意图。参见图10，示例性地，S为断层线，作为沉积面；沉积面S的左侧向上倾斜，左侧的延伸方向为沉积面S的上倾方向；沉积面S的右侧向下倾斜，右侧的延伸方向为沉积面S的下倾方向。
同相轴A1、同相轴A2、同相轴A3、同相轴A4依次形成在沉积面S上，同相轴A2相对同相轴A1朝沉积面S的上倾方向迁移，同相轴A3相对同相轴A2朝沉积面S的上倾方向迁移，同相轴A4相对同相轴A3朝沉积面S的上倾方向迁移，地层接触关系的指示特征为上超。
同相轴B1、同相轴B2、同相轴B3、同相轴B4依次形成在沉积面S上，同相轴B2相对同相轴B1朝沉积面S的上倾方向迁移，同相轴B3相对同相轴B2朝沉积面S的上倾方向迁移，同相轴B4相对同相轴B3朝沉积面S的上倾方向迁移，地层接触关系的指示特征为上超。
同相轴C1、同相轴C2依次形成在沉积面S上，同相轴C2相对同相轴C1朝沉积面S的下倾方向迁移，地层接触关系的指示特征为下超。
同相轴D1、同相轴D2、同相轴D3依次形成在沉积面S上，同相轴D2相对同相轴D1朝沉积面S的下倾方向迁移，同相轴D3相对同相轴D2朝沉积面S的下倾方向迁移，地层接触关系的指示特征为下超。
同相轴E1、同相轴E2依次形成在沉积面S上，同相轴E2相对同相轴E1朝沉积面S的下倾方向迁移，地层接触关系的指示特征为下超。
同相轴F1、同相轴F2、同相轴F3依次形成在沉积面S上，同相轴F1、同相轴F2、同相轴F3相互之间没有相对迁移，地层接触关系的指示特征为平行。
同相轴K的两端均延伸到去除图像的边缘，地层接触关系的指示特征为削截，可以作为地层单元的分界线。
步骤204：根据地层接触关系的指示特征，在竖直方向上划分地层单元。
可选地，该步骤204可以包括：
将指示特征相同的相邻地层划分到一个地层单元；
将指示特征不同的相邻地层划分到不同的地层单元。
地层接触关系的指示特征相同，地层的沉积方向相同，通常属于同一个地质时间段，因此划分为同一个地层单元，方便统一研究。地层接触关系的指示特征不同，地层的沉积方向改变，表示地质时段出现跨度，因此划分为不同的地层单元，方便进行区分。
示例性地，如图10所示，同相轴A1、A2、A3、A4对应地层和同相轴B1、B2、B3、B4对应地层的指示特征为上超，划分到一个地层单元L；同相轴C1、C2对应地层、同相轴D1、D2、D3对应地层和同相轴E1、E2的指示特征为下超，划分到一个地层单元M；同相轴F1、F2、F3对应地层的指示特征为平行，划分到一个地层单元N。另外，同相轴K对应地层的指示特征为削截，作为地层单元M和地层单元N之间的分界线。
可选地，该方法还可以包括：
根据地层单元内地层接触关系的指示特征，确定地层单元的沉积方向。
由于地层接触关系的指示特征相同的相邻地层划分到一个地层单元，因此地层单元内地层接触关系的指示特征相同。根据地层接触关系的指示特征，可以确定地层单元的沉积方向，有利于了解沉积情况，分析沉积序列的演化过程。
在实际应用中，根据地层单元内地层接触关系的指示特征，确定地层单元的沉积方向，可以包括：
若指示特征为上超，确定地层单元的沉积方向是从地势低处向地势高处；
若指示特征为下超，确定地层单元的沉积方向是从地势高处向地势低处；
若指示特征为平行，确定地层单元的沉积方向是垂直向上。
示例性地，如图10所示，地层单元L中，地层接触关系的指示特征为上超，表示湖平面上升，属于从湖盆中心向物源方向的退积型沉积；而且地层单元L中有同相轴A1、A2、A3、A4形成有一个上超，同相轴B1、B2、B3、B4形成有一个上超，表明地层单元L沉积时经历过两次退积型沉积。由于后面一次上超的起点比前面一次上超的起点靠近湖盆中心，因此可以推测出两次退积型沉积之间有一次快速的进积型沉积，估计是因为沉积时间短、厚度薄等原因而无法在地震剖面图像上识别出来。
地层单元M中，地层接触关系的指示特征为下超，表示物源供给增强，属于从物源向湖盆中心方向的进积型沉积；而且地层单元M中同相轴C1、C2形成有一个下超，同相轴D1、D2、D3形成有一个上超，同相轴E1、E2形成有一个下超，表明地层单元M沉积时经历过三次进积型沉积。由于地层单元M顶部的同相轴K形成有削截，表示第三次进积型沉积之后，地层抬升遭受剥蚀。
地层单元N中，地层接触关系的指示特征为平行，表明地层单元N沉积时湖盆沉积较为稳定，表现为向上垂向加积。
步骤205：根据地层接触关系的指示特征，在水平方向上划分沉积相带。
在本实施例中，沉积相带可以包括扇根亚相、扇中亚相和扇端亚相。
在实际应用中，沉积物从物源向湖盆中心移动并在移动过程中逐渐沉积下来，重量最大的石块等物质首先在靠近物源的位置沉积下来，在地震剖面图像上表现为反射杂乱、有不规则空白反射，划分为扇根亚相；然后是重量中等的砂石等物质在物源和湖盆中心的中间位置沉积下来，在地震剖面图像上表现为有明显退积、进积、加积特征的连续反射，划分为扇中亚相；最后是重量较轻的沙粒等物质在靠近湖盆中心的位置沉积下来，在地震剖面图像上表现为连续强反射，划分为扇端亚相。因此，可以直接在地震剖面图像上划分出沉积相带。
但是考虑到扇中亚相在地震剖面图像上表现为有明显退积、进积、加积特征的连续反射，可以识别出至少两种指示特征。为了提高划分的准确度，通常会采用地层接触关系的指示特征，辅助沉积相带的划分。
可选地，该步骤205可以包括：
将地层接触关系的指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
将没有地层接触关系的指示特征且靠近断层线的地层划分到扇根亚相；
将地层接触关系的指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
利用指示特征的数量和断层线配合，可以准确划分沉积相带。
图11为本公开实施例提供的地震剖面图像上沉积相带的示意图。参见图11，示例性地，靠近断层线的区域X有不规则空白反射，不能识别地层接触关系的指示特征，划分为扇根亚相；远离断层线的区域Z有连续强反射，地层接触关系的指示特征只有平行，划分为扇端亚相；区域X和区域Z之间的区域Y有明显退积、进积、加积特征的连续反射，地层接触关系的指示特征至少有上超、下超和平行三种，划分为扇中亚相。
步骤206：根据指示特征，生成沉积序列图像。
图12为本公开实施例提供的沉积序列图像的示意图。参见图12，竖直方向上从下到上依次为地层单元L、地层单元M和地层单元N，水平方向上从左到右依次为扇根亚相X、扇中亚相Y和扇端亚相Z，扇中亚相Y中的深色区域表示砂体。
本公开实施例提供了一种析沉积序列的装置，适用于图1或图6所示的分析沉积序列的方法。图13为本公开实施例提供的一种分析沉积序列的装置的结构示意图。参见图13，该装置包括：
图像获取模块301，用于获取地震剖面图像；
关系确定模块302，用于在地震剖面图像上确定地层接触关系的指示特征；
单元划分模块303，用于根据指示特征，在竖直方向上划分地层单元；
相带划分模块304，用于根据指示特征，在水平方向上划分沉积相带；
图像建立模块305，用于根据地层单元和沉积相带，生成沉积序列图像。
可选地，图像获取模块301可以包括：
分析子模块，用于对地震数据进行频谱分析，得到地震数据的频率分布范围；
去除子模块，用于去除频率分布范围内的低频部分，得到地震剖面图像，低频部分的下限值为频率分布范围的下限值，低频部分的上限值与下限值之差为频率分布范围的上限值与下限值之差的1/5～1/3。
可选地，关系确定模块302可以用于，
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴相对第一同相轴朝沉积面的上倾方向迁移，确定指示特征为上超；
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴相对第一同相轴朝沉积面的下倾方向迁移，确定指示特征为下超；
若地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，第二同相轴没有相对第一同相轴迁移，确定指示特征为平行；
若地震剖面图像中同相轴的两端均延伸到地震剖面图像的边缘，确定指示特征为削截。
可选地，单元划分模块303可以用于，
将指示特征相同的相邻地层划分到一个地层单元；
将指示特征不同的相邻地层划分到不同的地层单元。
可选地，相带划分模块304可以包括：
扇中划分子模块，用于将指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
扇根划分子模块，用于将没有指示特征且靠近断层线的地层划分到扇根亚相；
扇端划分子模块，用于将指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
可选地，该装置还可以包括：
方向确定模块，用于根据指示特征，确定地层单元的沉积方向。
图14是根据一示例性实施例示出的一种服务器的结构示意图。参见图14，所述服务器700包括中央处理单元(CPU)701、包括随机存取存储器(RAM)702和只读存储器(ROM)703的系统存储器704，以及连接系统存储器704和中央处理单元701的系统总线705。所述服务器700还包括帮助计算机内的各个器件之间传输信息的基本输入/输出系统(I/O系统)706，和用于存储操作系统713、应用程序714和其他程序模块715的大容量存储设备707。
所述基本输入/输出系统706包括有用于显示信息的显示器708和用于用户输入信息的诸如鼠标、键盘之类的输入设备709。其中所述显示器708和输入设备709都通过连接到系统总线705的输入输出控制器710连接到中央处理单元701。所述基本输入/输出系统706还可以包括输入输出控制器710以用于接收和处理来自键盘、鼠标、或电子触控笔等多个其他设备的输入。类似地，输入输出控制器710还提供输出到显示屏、打印机或其他类型的输出设备。
所述大容量存储设备707通过连接到系统总线705的大容量存储控制器(未示出)连接到中央处理单元701。所述大容量存储设备707及其相关联的计算机可读介质为服务器700提供非易失性存储。也就是说，所述大容量存储设备707可以包括诸如硬盘或者CD-ROM驱动器之类的计算机可读介质(未示出)。
不失一般性，所述计算机可读介质可以包括计算机存储介质和通信介质。计算机存储介质包括以用于存储诸如计算机可读指令、数据结构、程序模块或其他数据等信息的任何方法或技术实现的易失性和非易失性、可移动和不可移动介质。计算机存储介质包括RAM、ROM、EPROM、EEPROM、闪存或其他固态存储其技术，CD-ROM、DVD或其他光学存储、磁带盒、磁带、磁盘存储或其他磁性存储设备。当然，本领域技术人员可知所述计算机存储介质不局限于上述几种。上述的系统存储器704和大容量存储设备707可以统称为存储器。
根据本发明的各种实施例，所述服务器700还可以通过诸如因特网等网络连接到网络上的远程计算机运行。也即服务器700可以通过连接在所述系统总线705上的网络接口单元711连接到网络712，或者说，也可以使用网络接口单元711来连接到其他类型的网络或远程计算机系统(未示出)。
所述存储器还包括一个或者一个以上的程序，所述一个或者一个以上程序存储于存储器中，中央处理器701通过执行该一个或一个以上程序来实现图1或图6所示的分析沉积序列的方法。
在示例性实施例中，还提供了一种包括指令的非临时性计算机可读存储介质，例如包括指令的存储器，上述指令可由服务器的处理器执行以完成本发明各个实施例所示的长输管道的监控方法。例如，所述非临时性计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
需要说明的是：上述实施例提供的分析沉积序列的装置在沉积序列时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，上述实施例提供的分析沉积序列的装置与分析沉积序列的方法实施例属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成，也可以通过程序来指令相关的硬件完成，所述的程序可以存储于一种计算机可读存储介质中，上述提到的存储介质可以是只读存储器，磁盘或光盘等。
以上所述仅为本公开的可选实施例，并不用以限制本公开，凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。
Claims (8)
1.一种分析沉积序列的方法，其特征在于，所述方法包括：
获取地震剖面图像；
在所述地震剖面图像上识别断层线；
以所述断层线为分界线，将所述地震剖面图像划分为至少两个区域；
在目标区域的所述地震剖面图像上确定地层接触关系的指示特征，目标区域为至少两个区域中的任一个；
根据所述指示特征，在竖直方向上划分地层单元；
根据所述指示特征，在水平方向上划分沉积相带；
根据所述地层单元和所述沉积相带，生成沉积序列图像；
其中，在目标区域的所述地震剖面图像上确定地层接触关系的指示特征，包括以下步骤中的至少一种：
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的上倾方向迁移，确定所述指示特征为上超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的下倾方向迁移，确定所述指示特征为下超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴没有相对所述第一同相轴迁移，确定所述指示特征为平行；
若所述地震剖面图像中同相轴的两端均延伸到所述地震剖面图像的边缘，确定所述指示特征为削截；
其中，所述根据所述指示特征，在竖直方向上划分地层单元，包括：
将所述指示特征相同的相邻地层划分到一个地层单元；
将所述指示特征不同的相邻地层划分到不同的地层单元；
其中，所述根据所述指示特征，在水平方向上划分沉积相带，包括：
将所述指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
将没有所述指示特征且靠近断层线的地层划分到扇根亚相；
将所述指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
2.根据权利要求1所述的方法，其特征在于，所述获取地震剖面图像，包括：
对地震数据进行频谱分析，得到所述地震数据的频率分布范围；
去除所述频率分布范围内的低频部分，得到所述地震剖面图像，所述低频部分的下限值为所述频率分布范围的下限值，所述低频部分的上限值与下限值之差为所述频率分布范围的上限值与下限值之差的1/5～1/3。
3.根据权利要求1或2所述的方法，其特征在于，所述方法还包括：
根据所述指示特征，确定所述地层单元的沉积方向。
4.一种分析沉积序列的装置，其特征在于，所述装置包括：
图像获取模块，用于获取地震剖面图像；
关系确定模块，用于在所述地震剖面图像上识别断层线；以所述断层线为分界线，将所述地震剖面图像划分为至少两个区域；在目标区域的所述地震剖面图像上确定地层接触关系的指示特征，目标区域为至少两个区域中的任一个；
单元划分模块，用于根据所述指示特征，在竖直方向上划分地层单元；
相带划分模块，用于根据所述指示特征，在水平方向上划分沉积相带；
图像建立模块，用于根据所述地层单元和所述沉积相带，生成沉积序列图像；
其中，所述关系确定模块用于，
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的上倾方向迁移，确定所述指示特征为上超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴相对所述第一同相轴朝所述沉积面的下倾方向迁移，确定所述指示特征为下超；
若所述地震剖面图像中第一同相轴和第二同相轴依次形成在沉积面上，所述第二同相轴没有相对所述第一同相轴迁移，确定所述指示特征为平行；
若所述地震剖面图像中同相轴的两端均延伸到所述地震剖面图像的边缘，确定所述指示特征为削截；
其中，所述单元划分模块用于，
将所述指示特征相同的相邻地层划分到一个地层单元；
将所述指示特征不同的相邻地层划分到不同的地层单元；
其中，所述相带划分模块包括：
扇中划分子模块，用于将所述指示特征在竖直方向上有两种以上的地层划分到扇中亚相；
扇根划分子模块，用于将没有所述指示特征且靠近断层线的地层划分到扇根亚相；
扇端划分子模块，用于将所述指示特征在竖直方向上只有平行且远离断层线的地层划分到扇端亚相。
5.根据权利要求4所述的装置，其特征在于，所述图像获取模块包括：
分析子模块，用于对地震数据进行频谱分析，得到所述地震数据的频率分布范围；
去除子模块，用于去除所述频率分布范围内的低频部分，得到所述地震剖面图像，所述低频部分的下限值为所述频率分布范围的下限值，所述低频部分的上限值与下限值之差为所述频率分布范围的上限值与下限值之差的1/5～1/3。
6.根据权利要求4或5所述的装置，其特征在于，所述装置还包括：
方向确定模块，用于根据所述指示特征，确定所述地层单元的沉积方向。
7.一种服务器，其特征在于，所述服务器包括：
处理器；
被配置为存储所述处理器的可执行指令的存储器；
其中，所述处理器被配置为执行权利要求1至3任一项所述的方法。
8.一种计算机可读存储介质，其特征在于，当所述计算机可读存储介质中的指令由服务器的处理器执行时，使得所述服务器能够执行权利要求1至3任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911365712.2A CN113050164B (zh) | 2019-12-26 | 2019-12-26 | 分析沉积序列的方法、装置、服务器和存储介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911365712.2A CN113050164B (zh) | 2019-12-26 | 2019-12-26 | 分析沉积序列的方法、装置、服务器和存储介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113050164A CN113050164A (zh) | 2021-06-29 |
| CN113050164B true CN113050164B (zh) | 2024-01-30 |
Family
ID=76505261
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911365712.2A Active CN113050164B (zh) | 2019-12-26 | 2019-12-26 | 分析沉积序列的方法、装置、服务器和存储介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113050164B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117233843B (zh) * | 2023-11-13 | 2024-02-02 | 福瑞升(成都)科技有限公司 | 一种多模式地层接触关系Wheeler变换方法、装置及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103913769A (zh) * | 2014-02-24 | 2014-07-09 | 中海石油(中国)有限公司上海分公司 | 一种下降体系域的快速识别方法 |
| CN104977613A (zh) * | 2015-07-01 | 2015-10-14 | 中国石油天然气股份有限公司 | 基于多信息的碳酸盐岩岩相古地理重构方法及装置 |
| CN105653834A (zh) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | 基于层序地层原理的断裂瞬时活动定量计算方法 |
| CN109782344A (zh) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | 沉积层序边界识别方法及装置 |
-
2019
- 2019-12-26 CN CN201911365712.2A patent/CN113050164B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103913769A (zh) * | 2014-02-24 | 2014-07-09 | 中海石油(中国)有限公司上海分公司 | 一种下降体系域的快速识别方法 |
| CN105653834A (zh) * | 2014-11-10 | 2016-06-08 | 中国石油化工股份有限公司 | 基于层序地层原理的断裂瞬时活动定量计算方法 |
| CN104977613A (zh) * | 2015-07-01 | 2015-10-14 | 中国石油天然气股份有限公司 | 基于多信息的碳酸盐岩岩相古地理重构方法及装置 |
| CN109782344A (zh) * | 2018-12-13 | 2019-05-21 | 中国石油天然气集团有限公司 | 沉积层序边界识别方法及装置 |
Non-Patent Citations (1)
| Title |
|---|
| 孙建孟 等.《地球物理资料综合应用》.石油大学出版社,2001,296-327. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113050164A (zh) | 2021-06-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Galice | The continental margin off Galicia and Portugal: acoustical stratigraphy, dredge stratigraphy and structural evolution | |
| CN104765065B (zh) | 一种确定沉积地层古地貌的方法及装置 | |
| CN109870719B (zh) | 一种碳酸盐岩致密薄储层的井位布设方法、装置及系统 | |
| CN111025384B (zh) | 基于波形分类交汇融合的储层预测方法和装置 | |
| Von Huene | Structure of the Outer Convergent Margin Off Kodiak Island, Alaska, from Multichannel Seismic Records: Convergent Margins | |
| CN109143399B (zh) | 一种识别碳酸盐岩层序界面的方法 | |
| US11221425B1 (en) | Generating a model for seismic velocities in a subsurface region using inversion with lateral variations | |
| CN102109611A (zh) | 一种利用低频地震属性预测油藏优质储层的方法 | |
| Di Giulio et al. | Seismic response of a deep continental basin including velocity inversion: the Sulmona intramontane basin (Central Apennines, Italy) | |
| Anderson et al. | Facies, architecture, and origin of a reservoir-scale sand-rich succession within submarine canyon fill: insights from Wagon Caves Rock (Paleocene), Santa Lucia Range, California, USA | |
| CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
| CN113050164B (zh) | 分析沉积序列的方法、装置、服务器和存储介质 | |
| Curry Jr et al. | South Glenrock, a Wyoming stratigraphic oil field | |
| Chen et al. | Occurrence of submarine canyons, sediment waves and mass movements along the northern continental slope of the South China Sea | |
| Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
| CN109782344B (zh) | 沉积层序边界识别方法及装置 | |
| CN110727027B (zh) | 多期次河道砂的精细刻画方法 | |
| Ghalayini et al. | Using polygonal layer-bound faults as tools to delimit clastic reservoirs in the Levant Basin offshore Lebanon | |
| CN115903038A (zh) | 一种深水沉积岩性圈闭识别和定量刻画方法 | |
| CN113534253B (zh) | 页岩气三维地震甜点区优选方法和装置 | |
| Parchkoohi et al. | Automatic detection of karstic sinkholes in seismic 3D images using circular Hough transform | |
| Bersezio et al. | Trends in bed thickness and facies of turbiditic sandstone bodies: unravelling the effects of basin confinement, depositional processes, and modes of sediment supply | |
| CN113376692A (zh) | 致密砂岩气水平井压裂改造方案优化方法及装置 | |
| Bellatreche et al. | Detailed seismic lithofluid distribution using Bayesian stochastic inversion for a thinly bedded reservoir: A case study over Huntington UK Central North Sea | |
| Duvbiama et al. | 3D Static Modelling of an offshore field in the Niger-delta |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |