CN113204855B - 地层缩短速率确定方法、系统、设备及可读介质 - Google Patents

地层缩短速率确定方法、系统、设备及可读介质 Download PDF

Info

Publication number
CN113204855B
CN113204855B CN202010078505.5A CN202010078505A CN113204855B CN 113204855 B CN113204855 B CN 113204855B CN 202010078505 A CN202010078505 A CN 202010078505A CN 113204855 B CN113204855 B CN 113204855B
Authority
CN
China
Prior art keywords
stratum
formation
fault
determining
line
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
Application number
CN202010078505.5A
Other languages
English (en)
Other versions
CN113204855A (zh
Inventor
吴海
柳少波
卓勤功
鲁雪松
邓亚
于志超
桂丽黎
李伟强
王龙
白东来
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010078505.5A priority Critical patent/CN113204855B/zh
Publication of CN113204855A publication Critical patent/CN113204855A/zh
Application granted granted Critical
Publication of CN113204855B publication Critical patent/CN113204855B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Remote Sensing (AREA)
  • Computer Graphics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种地层缩短速率确定方法、系统、设备及可读介质,所述方法包括:通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理;根据数值化处理后的地层和断层线确定目标区域的地层和断块;根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率,本发明可使地层缩短速率的计算过程简便、快速和准确。

Description

地层缩短速率确定方法、系统、设备及可读介质
技术领域
本发明地质勘探技术领域,尤其涉及一种地层缩短速率确定方法、系统、设备及可读介质。
背景技术
地层横向缩短速率是指地层在横向挤压应力作用下,会发生缩短、错断、褶皱和隆升,地层缩短长度随时间的变化速度就是地层缩短率。挤压背景下断层相关褶皱地层横向缩短速率可以推测出挤压背景下主要构造地质活动或主要造山活动时期,推测可供油气聚集的圈闭形成的主要时期,研究逆冲断层系统下伏构造的几何形态、地震的复发周期和潜在的地震危险性等。研究发现传统的计算和测量地层缩短速率的手段具有操作方法复杂、成本高、监测时间长等不利特点。
发明内容
本发明的一个目的在于提供一种地层缩短速率确定方法,使地层缩短速率的计算过程简便、快速和准确。本发明的另一个目的在于提供一种地层缩短速率确定系统。本发明的再一个目的在于提供一种服计算机设备。本发明的还一个目的在于提供一种可读介质。
为了达到以上目的,本发明一方面公开了一种地层缩短速率确定方法,包括:
通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理;
根据数值化处理后的地层和断层线确定目标区域的地层和断块;
根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。
优选的,所述通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理具体包括:
获取目标区域的二维构造解释剖面图的图像文件;
将所述图像文件输入Midland Valley Move软件对图像文件中的地层和断层线进行数值化处理。
优选的,所述根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率具体包括:
确定每个地层的地层线长度;
确定所有相邻的两个地层的地层线长度差值和地质年龄差值;
根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率。
优选的,所述确定每个地层的地层线长度具体包括:
若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度;
若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
本发明还公开了一种地层缩短速率确定系统,包括:
剖面图处理模块,用于通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理;
地层划分模块,用于根据数值化处理后的地层和断层线确定目标区域的地层和断块;
缩短速率确定模块,用于根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。
优选的,所述剖面图处理模块具体用于获取目标区域的二维构造解释剖面图的图像文件,将所述图像文件输入Midland Valley Move软件对图像文件中的地层和断层线进行数值化处理。
优选的,所述缩短速率确定模块具体用于确定每个地层的地层线长度,确定所有相邻的两个地层的地层线长度差值和地质年龄差值,根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率。
优选的,所述缩短速率确定模块进一步用于若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度,若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
本发明还公开了一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述方法。
本发明还公开了一种计算机可读介质,其上存储有计算机程序,该程序被处理器执行时实现如上所述方法。
本发明收集目标区域的二维构造解释剖面图,通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理,并根据数值化处理后的地层和断层线划分出目标区域的地层和断块,从而可根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。本发明的地层缩短速率确定方法对设备要求不高,通过简单的设备即可实现地层缩短速率的定量计算,整个操作流程便捷,不需要对地层二维构造解释剖面图进行复杂的处理,计算步骤少,定量计算速度快,成本低,实用性强。并且本发明对于一般挤压背景下形成的断层相关褶皱地层均可定量计算,适用性较好。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1示出本发明一种地层缩短速率确定方法一个具体实施例的流程图之一;
图2示出本发明一种地层缩短速率确定方法一个具体实施例的流程图之二;
图3示出本发明一种地层缩短速率确定方法一个具体实施例的流程图之三;
图4示出本发明一种地层缩短速率确定方法一个具体实施例的流程图之四;
图5示出本发明一种地层缩短速率确定方法一个具体例子中的区块剖面图;
图6示出本发明一种地层缩短速率确定方法一个具体例子中数值化处理后的剖面图;
图7示出本发明一种地层缩短速率确定方法一个具体例子中断块划分后的剖面图;
图8示出本发明一种地层缩短速率确定方法一个具体例子中地层线拉平后的示意图;
图9示出本发明一种地层缩短速率确定方法一个具体例子中地层缩短速率的演化图;
图10示出本发明一种地层缩短速率确定系统一个具体实施例的结构图。
图11示出适于用来实现本发明实施例的计算机设备的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,关于定量计算地层缩短速率的方法主要是平衡地质剖面法、GPS法和古地磁法。其中,平衡地质剖面法需要将各个地质历史时期的剖面进行恢复,然后计算缩短量。其中,基于平衡地质剖面的恢复计算地层褶皱变形量的方法背后的原理主要4种:面积守恒法、深度(厚度)-隆升面积法、面积-深度-应变分析法和轴面分析法,运用这几种原理各自具有不同的适应性条件,其在一定程度上具有缺陷和不确定性,同时在应用过程中需要满足许多假设和前提,应用过程较为复杂,剖面恢复周期长。
GPS法主要是通过卫星监测技术来监测现代地表地层的位置随时间的移动速度来计算地层缩短速率,该方法具有很高的精度,但只能监测现今地表的地层移动量,其监测年代和地层范围有限,同时长时间的卫星监测成本较高。对于定量计算地震解释剖面中地层的缩短速率该方法不适用。
古地磁法主要是通过研究地层古地磁的变化来推测地层的活动、地层的变形甚至地层的倒转,但一般情况下古地磁法多会结合其它方法(如平衡地质剖面、地质模型分析等手段)来综合研究地层的变形速率,对于定量计算地震解释剖面中地层的缩短速率该方法具有明显的限制性。总体而言,传统的方法具有其优点,也有其明显的劣势,研究过程通常较为复杂。
根据本发明的一个方面,本实施例公开了一种地层缩短速率确定方法。如图1所示,本实施例中,所述方法包括:
S100:通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理。
S200:根据数值化处理后的地层和断层线确定目标区域的地层和断块。可按照地层的年限分别对不同的地层进行命名,例如,剖面自上到下依次命名地层为Hi(i=0,1,2,3……m),其中,m为地层的数目,取正整数。进一步的,可将剖面根据断层的切割特点按照从左至右、从上至下的顺序依次给断块进行命名,断块名称分别为A,B,C,D……N。
S300:根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。
本发明收集目标区域的二维构造解释剖面图,通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理,并根据数值化处理后的地层和断层线划分出目标区域的地层和断块,从而可根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。本发明的地层缩短速率确定方法对设备要求不高,通过简单的设备即可实现地层缩短速率的定量计算,整个操作流程便捷,不需要对地层二维构造解释剖面图进行复杂的处理,计算步骤少,定量计算速度快,成本低,实用性强。并且本发明对于一般挤压背景下形成的断层相关褶皱地层均可定量计算,适用性较好。
在优选的实施方式中,如图2所示,所述S100具体可包括:
S110:获取目标区域的二维构造解释剖面图的图像文件。
具体的,二维构造解释剖面图的图像文件可设置为位图格式,并具有长度比例尺、剖面深度和长度标线等信息。根据图片文件的比例尺或者构造解释剖面标线确定剖面的长、宽(深)数值,可分别记作为X0和Z0。实际应用时,可向地质建模软件导入二维地震解释剖面,设置剖面坐标尺度参数为X(长度)=X0,Y(厚度)=0,Z(宽度或深度)=Z0,由于二维构造解释剖面图是一个二维平面的图片,认为其没有厚度,所以设置Y=0。
S120:将所述图像文件输入Midland Valley Move软件对图像文件中的地层和断层线进行数值化处理。本实施例中,选择Midland Valley Move软件对图像文件进行数值化处理。在其他实施方式中,也可以通过其他地质建模软件或编写程序对图像文件进行数值化处理,本发明对此并不作限定。
可以理解的是,在该优选的实施方式中,通过地质建模软件Midland Valley Move的建模模块(Modeling Building)可对目标区域的二维构造解释剖面图的图像文件的地层和断层线进行数值化处理,建立剖面地层格架。通过建模模块将地层线和断层线进行数值化处理,建立剖面地层格架,不同的地层线运用不同的颜色区分表示。
进一步的,S200中,根据数值化处理后的图像文件上的地层和断层线,根据现有的地质研究方法可进一步确定出目标区域的地层和断块。可以理解的是,确定目标区域的地层和断块时,用户可以对数值化后的图像文件进行人为分析得到目标区域的地层和断块,也可通过根据目前的地质研究方法设置自动分析程序实现目标区域的地层和断块的划分。优选的,当通过地质建模软件得到断块后,可通过软件对各个断块进行编号并记录断块信息,以通过各个断块信息进行地层缩短速率的计算。
在优选的实施方式中,如图3所示,所述S300具体可包括:
S310:确定每个地层的地层线长度SUMLi。其中,地层i的地层线长度为SUMLi
S320:确定所有相邻的两个地层的地层线长度差值ΔLi和地质年龄差值ΔTi
S330:根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率Vi,即Vi=ΔLi/ΔTi
在优选的实施方式中,如图4所示,所述S310具体可包括:
S311:若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度。
S312:若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
具体的,计算得到各个地层界线的长度值Li,若有断层切割,各个断块地层界线的长度值为Lij,j为断块数量,j=A,B,C,D……N。当地层没有被断层切割时,该地质时期地层线总长度为SUMLi=Li;当地层被断层切割时,总共剖分出j个断块,地层线总长度为SUMLi=LiA+LiB+……+Lij。则相邻两个地层的地层界线长度差值ΔLi=SUMLi-SUMLi-1,对应的地层界线地质年龄差值ΔTi=Ti-Ti-1。第Hi套地层在Ti-1到Ti的时间段内横向缩短速率为Vi=ΔLi/ΔTi。
本发明通过将断层相关褶皱地层图件数值化后,将各个地层界线进行拉平处理。测量拉平后的各个断块地层长度值,得到各地层的总长度值。通过计算相邻地层界线长度差值与对应地层地质年龄差值的比值得到横向地层缩短速率。本发明借助Move软件实现位图曲线数值化并测量曲线拉平后长度,进而定量计算地层的缩短速率,通过各个时期地层缩短速率大小判断地层变形强弱,为构造变形、油气圈闭主要形成时期、地质活动和造山运动时期的确定提供方法。
综上,本发明具有以下优点:1.对设备要求不高。仅需要收集带有比例尺或深度、长度信息的位图和Move软件。图件的分辨率不需要太高要求,仅需要看清地层和断层线即可。2.操作流程便捷。整个方法操作流程简捷,不需要对图件进行太过复杂的处理,无需二维面积的计算,将二维的面积计算转化为对线条长度的计算。3.适用性较好。一般对于挤压背景下形成的断层相关褶皱地层均可适用。4.计算方法简单。整个计算方法简便直观,计算步骤少,定量计算速度快,成本低,实用性强。
下面通过一个具体例子来对本发明作进一步的说明。以一个区块挤压背景下断层相关褶皱地层模型来计算各地质时期地层的缩短速率。收集目标区域带有深度、长度标线的二维构造解释剖面图件,如图5所示,图件为PNG格式。根据剖面图件的深度、长度标线读取长度数值为:X0=8000m;深度数值为:Z0=4000m。将位图图件导入Move软件,设置图件的坐标参数为X(长度)=8000,Y(厚度)=0,Z(宽度或深度)=4000,单位:m(米)。利用Move软件建模(Modeling Building)模块将地层线(Horizon)和断层线(Fault)进行数值化处理,区块数值化后的地层格架如图6所示。
根据区块地质剖面特点,划分地层和断块。该区块总共有5条地层界线(包括顶界和底界),将地层从新到老(剖面上为从上到下)地层界线名依次为H0,H1,H2,H3,H4。各个地层对应的地质年代为T0=0Ma,T1=3Ma,T2=23.3Ma,T3=35.8Ma,T4=65Ma。将地质剖面根据断层切割情况进行断块剖分,区块地层可划分为两个断块,分别为断块A和断块B,如图7所示。
运用Move软件中的二维剖面分析(2D Section Analysis)模块进行分析计算,勾选分析整个剖面,同时勾选断层分析选项框中的显示端点(Display End Points)、地层线图层(Horizon Labels)和(地层线长度图层),让拉平后的地层线显示长度和断点,计算拉平后的地层线如图8所示。分别读取各个地层线的长度L0=8000m,L1=8000m,L2A=2475.1m,L2B=6207.2m,L3A=4233.4m,L3B=5018.2m,L4为本剖面基底,未读取长度,数据详情见表1。
表1
Figure BDA0002379375260000071
Figure BDA0002379375260000081
其中,各地层线总长度为:SUML0=L0=8000m,SUML1=L1=8000m,SUML2=L2A+L2B=8682.3m,SUML3=L3A+L3B=9251.5m。相邻地层缩短量为:H1与H0长度差值为:ΔL1=SUML1-SUML0=0m,H2与H1长度差值为:ΔL2=SUML1-SUML0=682.3m,H3与H2长度差值为:ΔL3=SUML1-SUML0=569.2m。相邻地层界线时间差为:H1与H0地质时间差值为:ΔT1=T1-T0=3Ma,H2与H1地质时间差值为:ΔT2=T2-T1=20.3Ma,H3与H2地质时间差值为:ΔT3=T3-T2=12.5Ma。地层缩短速率为:T1到T0时间段内地层缩短速率为:V1=ΔL1/ΔT1=0m/Ma,T2到T1时间段内地层缩短速率为:V2=ΔL2/ΔT2=33.611m/Ma,T3到T2时间段内地层缩短速率为:V3=ΔL3/ΔT3=45.536m/Ma。根据得到的地层缩短速率可得到区块断层相关褶皱地质历史时期地层缩短速率演化图,如图9所示,数据详情见表2。
表2
Figure BDA0002379375260000082
基于相同原理,本实施例还公开了一种地层缩短速率确定系统。如图10所示,所述系统包括剖面图处理模块11、地层划分模块12和缩短速率确定模块13。
其中,剖面图处理模块11用于通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理。
地层划分模块12用于根据数值化处理后的地层和断层线确定目标区域的地层和断块。
缩短速率确定模块13用于根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率。
在优选的实施方式中,所述剖面图处理模块11具体用于获取目标区域的二维构造解释剖面图的图像文件,将所述图像文件输入Midland Valley Move软件对图像文件中的地层和断层线进行数值化处理。
在优选的实施方式中,所述缩短速率确定模块13具体用于确定每个地层的地层线长度,确定所有相邻的两个地层的地层线长度差值和地质年龄差值,根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率。
在优选的实施方式中,所述缩短速率确定模块13进一步用于若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度,若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
由于该系统解决问题的原理与以上方法类似,因此本系统的实施可以参见方法的实施,在此不再赘述。
上述实施例阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。一种典型的实现设备为计算机设备,具体的,计算机设备例如可以为个人计算机、膝上型计算机、蜂窝电话、相机电话、智能电话、个人数字助理、媒体播放器、导航设备、电子邮件设备、游戏控制台、平板计算机、可穿戴设备或者这些设备中的任何设备的组合。
在一个典型的实例中计算机设备具体包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现如上所述方法。
下面参考图11,其示出了适于用来实现本申请实施例的计算机设备600的结构示意图。
如图11所示,计算机设备600包括中央处理单元(CPU)601,其可以根据存储在只读存储器(ROM)602中的程序或者从存储部分608加载到随机访问存储器(RAM))603中的程序而执行各种适当的工作和处理。在RAM603中,还存储有系统600操作所需的各种程序和数据。CPU601、ROM602、以及RAM603通过总线604彼此相连。输入/输出(I/O)接口605也连接至总线604。
以下部件连接至I/O接口605:包括键盘、鼠标等的输入部分606;包括诸如阴极射线管(CRT)、液晶反馈器(LCD)等以及扬声器等的输出部分607;包括硬盘等的存储部分608;以及包括诸如LAN卡,调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器610上,以便于从其上读出的计算机程序根据需要被安装如存储部分608。
特别地,根据本发明的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本发明的实施例包括一种计算机程序产品,其包括有形地包含在机器可读介质上的计算机程序,所述计算机程序包括用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分609从网络上被下载和安装,和/或从可拆卸介质611被安装。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
为了描述的方便,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
以上所述仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (6)

1.一种地层缩短速率确定方法,其特征在于,包括:
通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理;
根据数值化处理后的地层和断层线确定目标区域的地层和断块;
根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率;
其中,所述根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率具体包括:
确定每个地层的地层线长度;
确定所有相邻的两个地层的地层线长度差值和地质年龄差值;
根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率;
其中,所述确定每个地层的地层线长度具体包括:
若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度;
若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
2.根据权利要求1所述的地层缩短速率确定方法,其特征在于,所述通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理具体包括:
获取目标区域的二维构造解释剖面图的图像文件;
将所述图像文件输入Midland Valley Move软件对图像文件中的地层和断层线进行数值化处理。
3.一种地层缩短速率确定系统,其特征在于,包括:
剖面图处理模块,用于通过地质建模软件对目标区域的二维构造解释剖面图的地层和断层线进行数值化处理;
地层划分模块,用于根据数值化处理后的地层和断层线确定目标区域的地层和断块;
缩短速率确定模块,用于根据相邻两个地层的地层线长度和地质年龄确定地层缩短速率;
其中,所述缩短速率确定模块具体用于确定每个地层的地层线长度,确定所有相邻的两个地层的地层线长度差值和地质年龄差值,根据相邻的两个地层的地层线长度差值和地质年龄差值的比值确定地层缩短速率;
其中,所述缩短速率确定模块进一步用于若一个地层不包括断块,计算所述一个地层的地层线总长度得到所述一个地层的地层线长度,若一个地层包括至少一个断块,计算每个断块的地层线总长度,计算所有断块的地层线总长度的和得到所述一个地层的地层线长度。
4.根据权利要求3所述的地层缩短速率确定系统,其特征在于,所述剖面图处理模块具体用于获取目标区域的二维构造解释剖面图的图像文件,将所述图像文件输入MidlandValley Move软件对图像文件中的地层和断层线进行数值化处理。
5.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现如权利要求1或2所述方法。
6.一种计算机可读介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1或2所述方法。
CN202010078505.5A 2020-02-03 2020-02-03 地层缩短速率确定方法、系统、设备及可读介质 Active CN113204855B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010078505.5A CN113204855B (zh) 2020-02-03 2020-02-03 地层缩短速率确定方法、系统、设备及可读介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010078505.5A CN113204855B (zh) 2020-02-03 2020-02-03 地层缩短速率确定方法、系统、设备及可读介质

Publications (2)

Publication Number Publication Date
CN113204855A CN113204855A (zh) 2021-08-03
CN113204855B true CN113204855B (zh) 2022-11-01

Family

ID=77024842

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010078505.5A Active CN113204855B (zh) 2020-02-03 2020-02-03 地层缩短速率确定方法、系统、设备及可读介质

Country Status (1)

Country Link
CN (1) CN113204855B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108254784A (zh) * 2017-12-11 2018-07-06 中国石油天然气股份有限公司 一种基于二维地震数据的断层识别方法、装置及系统
CA3063929A1 (en) * 2017-08-18 2019-02-21 Landmark Graphics Corporation Fault detection based on seismic data interpretation
CN109444960A (zh) * 2018-11-26 2019-03-08 中国石油天然气股份有限公司大港油田分公司 用于断层构造的判定方法
CN109597964A (zh) * 2018-11-05 2019-04-09 中国石油天然气集团有限公司 理论油藏渗流参数值确定方法、系统、设备、可读介质
CN109975872A (zh) * 2019-04-26 2019-07-05 河南理工大学 一种利用Croel DRAW软件直接编制平衡剖面的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA3063929A1 (en) * 2017-08-18 2019-02-21 Landmark Graphics Corporation Fault detection based on seismic data interpretation
CN108254784A (zh) * 2017-12-11 2018-07-06 中国石油天然气股份有限公司 一种基于二维地震数据的断层识别方法、装置及系统
CN109597964A (zh) * 2018-11-05 2019-04-09 中国石油天然气集团有限公司 理论油藏渗流参数值确定方法、系统、设备、可读介质
CN109444960A (zh) * 2018-11-26 2019-03-08 中国石油天然气股份有限公司大港油田分公司 用于断层构造的判定方法
CN109975872A (zh) * 2019-04-26 2019-07-05 河南理工大学 一种利用Croel DRAW软件直接编制平衡剖面的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
川东北高陡复杂构造区三维地震精细构造解释技术;徐敏等;《石油物探》;20150325(第02期);全文 *
活断层构造建模与可视化;刘玉娟等;《地壳构造与地壳应力文集》;20070515(第00期);全文 *

Also Published As

Publication number Publication date
CN113204855A (zh) 2021-08-03

Similar Documents

Publication Publication Date Title
US10145227B2 (en) Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves
US8599643B2 (en) Joint structural dip removal
US10007015B2 (en) Methods, systems and devices for predicting reservoir properties
CA2913827C (en) Methods, systems and devices for predicting reservoir properties
US9383465B2 (en) Quantitative analysis of time-lapse seismic data
CN105804730B (zh) 使用历史井数据的资源识别的方法和系统
US9157319B2 (en) Stratigraphic modeling using production data density profiles
CN102576370B (zh) 空隙度分析的系统和方法
CN113552621B (zh) 页岩气地应力确定方法和装置
Wald et al. A new strategy for developing Vs30 maps
Uzielli et al. Risk analysis for the Ancona landslide—I: characterization of landslide kinematics
Bigi et al. Discrete fracture network of the Latemar carbonate platform
CN116663762A (zh) 一种城市规划地下空间勘察测绘方法及系统
CN103628866B (zh) 获取地层的地层压力系数的方法
CN113204855B (zh) 地层缩短速率确定方法、系统、设备及可读介质
CN104199091B (zh) 基于方位地震数据的纹理属性方法
Zou et al. Identification of Paleoearthquakes and Coseismic Slips on a Normal Fault Using High-Precision Quantitative Morphology: Application to the Jiaocheng Fault in the Shanxi Rift, China
Ziesch et al. Subseismic pathway prediction by three-dimensional structural restoration and strain analysis based on seismic interpretation
CN107797148B (zh) 一种基于三维地质建模的航磁异常场分离方法及系统
US20210405234A1 (en) System and method for subsurface structural interpretation
Jarosiński et al. Contemporary state of stress in a stable plate interior (northern Poland): The integration of satellite geodesy, borehole and seismological data
Suckale et al. Probabilistic seismic hazard model for Vanuatu
CN107015289A (zh) 圈闭地质资源量确定方法及装置
Stewart Scale dependence of strike and dip in sedimentary basins: Implications for field measurements and integrating subsurface datasets
CN113267810B (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
GR01 Patent grant
GR01 Patent grant