CN111913232A - 地层平衡剖面恢复方法和装置 - Google Patents
地层平衡剖面恢复方法和装置 Download PDFInfo
- Publication number
- CN111913232A CN111913232A CN202010771313.2A CN202010771313A CN111913232A CN 111913232 A CN111913232 A CN 111913232A CN 202010771313 A CN202010771313 A CN 202010771313A CN 111913232 A CN111913232 A CN 111913232A
- Authority
- CN
- China
- Prior art keywords
- stratum
- period
- geological structure
- rectangle
- rectangles
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本申请提供一种地层平衡剖面恢复方法和装置，该方法包括：通过获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，在当前时期地质结构剖面中设置矩形阵列，矩形阵列包括(N‑1)*M的矩形，矩形阵列的(N‑1)行矩形分别设置在最接近地表的地层之外的N‑1层地层，每列的矩形位于同一垂向方向，从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N‑1个不同时期的平衡恢复后的地质结构剖面，因此，能更加准确地对复杂构造变形地质结构进行平衡剖面的恢复。
Description
技术领域
本申请涉及地质勘探技术领域，尤其涉及一种地层平衡剖面恢复方法和装置。
背景技术
平衡剖面技术是一种遵循几何守恒原则而建立的地质剖面正演与恢复方法，是构造演化定量分析的有效手段。平衡剖面的恢复，是借助平衡剖面技术，利用现今的地质剖面来还原各时期地层形成时的地质构造特征，对地质构造变形的时间、变形的强度和变形的机理进行分析研究，为深刻认识地质构造发育史，分析油气的形成时期、运移时期以及成藏时期，进而寻找可能的含油气区域，对后续油气的勘探开发提供有效的指导。
目前，平衡剖面恢复的方法主要有两种：一种是基于层长守恒原理的层位拉平恢复方法，该方法恢复的前提是假定地层的层面变形前后长度不发生改变，地层受构造运动的影响只发生了褶皱变形、断裂，并没有发生层间的岩石流动，恢复时通过计算变形后地层的长度，确定变形前原始地层层面长度，进而对原始地层恢复；另一种是基于面积守恒原理的断块调整恢复方法，主要适于于地质结构相对简单的地区，在恢复过程中，以断层为边界，剖面被划分为若干断块，作为恢复的基本地质单元，对每一个断块进行平移、旋转，使层位逐段拉平。
上述各方法不能准确地对复杂构造变形地质结构进行平衡剖面的恢复。
发明内容
本申请提供一种地层平衡剖面恢复方法和装置，用以解决不能准确地对复杂构造变形地质结构进行平衡剖面的恢复的问题。
第一方面，本申请提供一种地层平衡剖面恢复方法，包括：
获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数；
在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向；
从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面；
其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
可选的，每层地层内的各矩形的顶边与所述地层的上地层线贴齐，且底边与所述地层的下地层线贴齐；
所述矩形的高度与所在地层的厚度有关，所述矩形的宽度为所述矩形的高度的K倍，所述K为大于等于2的整数。
可选的，相邻列的矩形间的距离与所述相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
可选的，所述从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面，包括：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，所述第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层；
将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，所述第j+1个时期的平衡恢复后的地质结构剖面包括所述N-(j+1)层不同属性的地层；
其中，j为大于等于1且小于等于N-2的整数。
可选的，所述在所述当前时期地质结构剖面中设置矩形阵列，包括：
若所述N-1层地层中包括遭受剥蚀的地层，则根据所述当前时期地质结构剖面中所述遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定所述遭受剥蚀的地层的原始沉积厚度；
根据所述遭受剥蚀的地层的原始沉积厚度，为所述遭受剥蚀的地层设置剥蚀段的矩形。
可选的，还包括：
若所述N-1层地层中包括被断裂错断的地层，则在垂向位移被断裂错断的地层的过程中，将所述被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；
其中，所述被断裂错断的地层中逆断层内的矩形向地底方向移动，所述被断裂错断的地层中正断层内的矩形向地表方向移动。
可选的，还包括：
获取所述当前时期地质结构剖面中目标地层的当前长度，所述目标地层为遭受剥蚀的地层或被断裂错断的地层；
根据所述目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中所述目标地层的长度等于所述当前长度。
第二方面，本申请提供一种地层平衡剖面恢复装置，包括：
获取模块，用于获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数；
设置模块，用于在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向；
处理模块，用于从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面；
其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
可选的，每层地层内的各矩形的顶边与所述地层的上地层线贴齐，且底边与所述地层的下地层线贴齐；
所述矩形的高度与所在地层的厚度有关，所述矩形的宽度为所述矩形的高度的K倍，所述K为大于等于2的整数。
可选的，相邻列的矩形间的距离与所述相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
可选的，所述处理模块，具体用于：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，所述第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层；
将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，所述第j+1个时期的平衡恢复后的地质结构剖面包括所述N-(j+1)层不同属性的地层；
其中，j为大于等于1且小于等于N-2的整数。
可选的，所述设置模块，具体用于：
若所述N-1层地层中包括遭受剥蚀的地层，则根据所述当前时期地质结构剖面中所述遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定所述遭受剥蚀的地层的原始沉积厚度；
根据所述遭受剥蚀的地层的原始沉积厚度，为所述遭受剥蚀的地层设置剥蚀段的矩形。
可选的，所述处理模块，还用于：
若所述N-1层地层中包括被断裂错断的地层，则在垂向位移被断裂错断的地层的过程中，将所述被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；
其中，所述被断裂错断的地层中逆断层内的矩形向地底方向移动，所述被断裂错断的地层中正断层内的矩形向地表方向移动。
可选的，所述处理模块，还用于：
获取所述当前时期地质结构剖面中目标地层的当前长度，所述目标地层为遭受剥蚀的地层或被断裂错断的地层；
根据所述目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中所述目标地层的长度等于所述当前长度。
第三方面，本申请提供一种地层平衡剖面恢复装置，包括：存储器和处理器；
所述存储器用于存储程序指令；
所述处理器用于调用所述存储器中的程序指令执行如本申请第一方面所述的地层平衡剖面恢复方法。
第四方面，本申请提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序指令，所述计算机程序指令被执行时，实现如本申请第一方面所述的地层平衡剖面恢复方法。
本申请提供的地层平衡剖面恢复方法和装置，通过获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数，在当前时期地质结构剖面中设置矩形阵列，矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向，从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面，其中，在垂向位移每层地层时，其它各层地层均随之垂向位移，通过上述方式，能够更加准确地对复杂构造变形地质结构进行平衡剖面的恢复，可对先存构造特征进行复原，认识构造发育史，解析构造活动的时期、活动的强度，有助于分析和判断油气的运移、聚集的时期，对于多期叠合盆地研究、超深层油气资源的勘探、开发工作均起到重要的促进作用。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本申请一实施例提供的地层平衡剖面恢复方法的流程图；
图2为本申请一实施例提供的水平地层平衡剖面恢复过程示意图；
图3为本申请一实施例提供的矩形阵列在复杂构造区域设置示意图；
图4为本申请另一实施例提供的剥蚀区地层平衡剖面恢复过程示意图；
图5为本申请另一实施例提供的断层活动引起的复杂地层平衡剖面恢复过程示意图；
图6为本申请另一实施例提供的矩形阵列在测地层倾角上的应用示意图；
图7为本申请一实施例提供的地层平衡剖面恢复方法的应用效果图；
图8a和图8b为本申请另一实施例提供的地层平衡剖面恢复方法的应用效果图；
图9为本申请一实施例提供的地层平衡剖面恢复装置的结构示意图；
图10为本申请另一实施例提供的地层平衡剖面恢复装置的结构示意图；
图11为本申请另一实施例提供的地层平衡剖面恢复装置的结构示意图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
平衡剖面的恢复，是借助平衡剖面技术，利用现今的地质剖面来还原各时期地层形成时的地质构造特征，对地质构造变形的时间、变形的强度和变形的机理进行分析研究，为深刻认识地质构造发育史，分析油气的形成时期、运移时期以及成藏时期，进而寻找可能的含油气区域，对后续油气的勘探开发提供有效的指导。目前，平衡剖面恢复的方法主要有两种：一种是基于层长守恒原理的层位拉平恢复方法，另一种是基于面积守恒原理的断块调整恢复方法。上述各方法不能准确地对复杂构造变形地质结构进行平衡剖面的恢复。
本申请提供的地层平衡剖面恢复方法和装置，通过获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数，在当前时期地质结构剖面中设置矩形阵列，矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向，从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面，其中，在垂向位移每层地层时，其它各层地层均随之垂向位移，通过上述方式，能够更加准确地对复杂构造变形地质结构进行平衡剖面的恢复，可对先存构造特征进行复原，认识构造发育史，解析构造活动的时期、活动的强度，有助于分析和判断油气的运移、聚集的时期，对于多期叠合盆地研究、超深层油气资源的勘探、开发工作均起到重要的促进作用。
图1为本申请一实施例提供的地层平衡剖面恢复方法的流程图，本实施例的方法可以应用于电子设备中，该电子设备可以是终端设备、服务器等，终端设备又例如可以是手机、平板电脑、笔记本电脑、台式电脑等。如图1所示，本实施例的方法可以包括：
S101、获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分。
本实施例中，综合钻井、非地震资料、野外露头资料，获取地层的埋深、地层展布特征、地层岩性等信息对地震剖面进行解释，解释的内容包括：在地震剖面上划分地层层序、确定地层的形成时期以及断层，通过解释好的地震剖面可以直接观察到复杂构造变形地质结构包括多期叠合盆地在内的各套地层在地下的真实形态。
具体地，地层展布特征是指一套地层在形成的时侯，在平面上所覆盖的区域范围和在不同区域的厚度会有变化，这种分布范围和厚度变化总称为地层展布特征。例如：新地层形成的时候，地层的下方老地层是一个隆起的样貌，即四周低，中间高，四周在水下，中间露出水面，那么新地层形成的时候，新地层中间部位沉积的少，就厚度小，四周沉积的多，就厚度大。
具体地，在地震剖面上解释地层的形成时期的方法为：地层是由岩石组成的，岩石中又含有多种矿物，在地震剖面上解释地层的形成时期时，需要结合油田钻井取出的岩芯(即井里面取的岩石，一个圆柱形的岩石柱子)，根据取出岩石的种类、岩石中存在的某些可用来确定岩石形成时间的矿物，来确定该岩石的形成时间。不同时期形成的岩石有区别，这就为确定地层的形成时期提供了依据，地震剖面纵向上的深度和油井取芯层位的深度可以相互对应，相同深度下，将岩芯与地震剖面对应上，就可以确定地震反射界面的岩石种类与形成时期，进而对地震界面做出标定，解释出地震剖面上岩层形成时间，即地震剖面上地层的形成时期。
具体地，对地震剖面进行解释的步骤包括如下三步：
第一步，通过钻井资料标定主要地层单元所对应的地震反射界面，通过地震反射界面的解释划分出上下地层的界限，即地层线；
第二步，对解释的地层线赋予不同的属性进行区分，相同的地层线应用相同的属性，属性具体包括线条颜色、粗细、线型；
第三步，对解释好的地震剖面进行时深转换，即根据钻井取的岩心、野外岩石露头、地震数据等资料，计算出时间与深度的换算关系，将时间域地震剖面转化为深度域的地质结构剖面。
通过以上三步，就可以获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数，例如：图2为本申请一实施例提供的水平地层平衡剖面恢复过程示意图，如图2的
所示，地质结构剖面包含A、B、C、D、E、F、G、H、I、J、K、L共12层不同属性的地层，不同的地层间用不同属性的地层线划分。
S102、在当前时期地质结构剖面中设置矩形阵列，矩形阵列包括(N-1)*M的矩形，矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向。
本实施例中，当前时期地质结构剖面包括N层不同属性的地层，在除最接近地表的地层之外的N-1层地层中设置矩形阵列，矩形阵列包括(N-1)*M的矩形，即包括N-1行、M列的矩形，M为大于等于1的整数，N-1行对应N-1层地层，每一层地层对应每一列都包含一个矩形，且每列的矩形位于同一垂向方向。例如：如图2的
为在N-1层地层设置的一列矩形，该列的矩形位于同一垂向方向，图2的
为在N-1层地层设置的M列矩形，每一列的矩形位于同一垂向方向，构成了(N-1)*M的矩形阵列。同一地层中矩形的颜色属性一致。
可选的，每层地层内的各矩形的顶边与地层的上地层线贴齐，且底边与地层的下地层线贴齐。
本实施例中，每层地层内的包含多个矩形，矩形的属性与地层的层位以及地层的厚度有关，每层地层内的矩形的顶边与地层的上地层线贴齐，且底边与地层的下地层线贴齐。例如：在水平地层段，如图2的
所示，矩形顶面与地层线贴齐竖直放置；图3为本申请一实施例提供的矩形阵列在复杂构造区域设置示意图，如图3所示，地层为倾斜地层段、构造变形强烈段，如图3的
和图3的
矩形阵列根据地层产状进行适度旋转调整，使其上下边与地层线顶底贴合。
可选的，矩形的高度与所在地层的厚度有关，矩形的宽度为矩形的高度的K倍，K为大于等于2的整数。
本实施例中，每层地层内的包含多个矩形，矩形的属性与地层的层位以及地层的厚度有关，矩形高度取决于地层的厚度变化，单一矩形的宽度一般为高度的2-3倍。
可选的，相邻列的矩形间的距离与相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
本实施例中，根据地质结构剖面中地层的长度、厚度和变形强度在地质结构剖面中设置不同密度的矩形阵列为M列，矩形阵列的列密度决定了恢复的精度。根据剖面在水平延伸方向(剖面走向)及变形特征设置矩形阵列的列密度，在地层厚度均匀段每隔15±2.5km设置矩形阵列；在地层厚度渐变的古隆起、凹陷、断块发育段每隔10±2.5km设置矩形阵列；在构造变形强烈的逆冲断裂、同沉积构造变形、盐构造、复杂褶皱等段每隔5±2.5km设置矩形阵列。
可选的，若N-1层地层中包括遭受剥蚀的地层，则根据当前时期地质结构剖面中遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定遭受剥蚀的地层的原始沉积厚度。
本实施例中，图4为本申请另一实施例提供的剥蚀区地层平衡剖面恢复过程示意图，如图4所示，该剥蚀地层在设置矩形阵列的时候，先要根据地层厚度变化与地层的接触关系判断抬升时期。判断抬升时期即判断遭受剥蚀的地层的年代，根据被抬升地层的年代确定抬升开始与结束的时间，例如：侏罗纪形成的地层发生变形，比其更晚形成的白垩纪地层未变形，则可判断为地层抬升年代为侏罗纪，即地层抬升时期为侏罗纪。然后测量地质结构剖面上未剥蚀段的矩形高度，将其求和取平均值，所得到的数值代表同一层位未剥蚀段的矩形平均高度，将其视为遭受剥蚀段地层的原始沉积厚度。
可选的，根据遭受剥蚀的地层的原始沉积厚度，为遭受剥蚀的地层设置剥蚀段的矩形。
本实施例中，在上述实施例的基础上，在确定了遭受剥蚀的地层的原始沉积厚度后，按地层由老到新依次用原始沉积厚度对应的矩形对剥蚀区的矩形进行同底替换，若局部地层完全剥蚀，也按照此方式由老至新进行剥蚀量的恢复(如图4的
)，通过这种方式来完成遭受剥蚀的地层剥蚀量的恢复。
S103、从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面。
本实施例中，地质结构剖面有N-1层地层需要做平衡恢复，逐层恢复时按照从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，将位移后具有相同矩形属性的矩形底边进行连线，该连线即为进行平衡恢复后的地层底面，并为该连线附上同位移之前的属性，获得该层地层对应时期的平衡恢复后的地质结构剖面，逐层恢复后，可以得到N-1层地层对应的N-1个不同时期的平衡恢复后的地质结构剖面。
可选的，在垂向位移每层地层时，其它各层地层均随之垂向位移。
本实施例中，通过将各地层同步垂向位移，来保持每层地层内的矩形阵列垂向相对关系不变。
其中，上述目标区域深度域的当前时期地质结构剖面可以是用户向执行本方法实施例的电子设备输入的，或者，是其它设备向执行本方法实施例的电子设备发送的。
在图1所示实施例的基础上，S103的具体实现方式为：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至第2层地层的矩形的顶边与地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层。
将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且第2层地层下方的其它各层地层的矩形也随之垂向位移，直至第2层地层的矩形的顶边与地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，第j+1个时期的平衡恢复后的地质结构剖面包括N-(j+1)层不同属性的地层。j为大于等于1且小于等于N-2的整数。
假设当前时期地质结构剖面包括4层不同属性的地层，且除最接近地表的第1层地层之外的3层地层中设置了矩形阵列，矩形阵列包括3*3的矩形，将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至第2层地层的矩形的顶边与地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，第1个时期的平衡恢复后的地质结构剖面包括3层不同属性的地层；将第1个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且第2层地层下方的其它各层地层的矩形也随之垂向位移，直至第2层地层的矩形的顶边与地表贴齐，获得最接近当前时期的第2个时期的平衡恢复后的地质结构剖面，第2个时期的平衡恢复后的地质结构剖面包括2层不同属性的地层；将第2个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且第2层地层下方的其它各层地层的矩形也随之垂向位移，直至第2层地层的矩形的顶边与地表贴齐，获得最接近当前时期的第3个时期的平衡恢复后的地质结构剖面，第3个时期的平衡恢复后的地质结构剖面包括1层不同属性的地层。通过这种方式，获得了当前时期之前3个不同时期的平衡恢复后的地质结构剖面。
进行位移恢复时，矩形阵列顶层的矩形决定了阵列重排的位置，该矩形的顶、底边对齐方式由地质构造背景而定。
可选的，地质构造为水平地层，如图2所示，该水平地层的具体位移恢复步骤包括如下三步：
第一步，恢复地层原始变形特征，即：将要恢复的地层A(如图2的
)的层内矩形垂向平移至地表，即该层矩形顶边与0海拔线贴齐，所有的矩形阵列保持垂向相对关系不变，跟随层A进行垂向位移(如图2的
)；
第二步，将位移后具有相同矩形属性的矩形底边进行连线，该连线即为进行构造恢复后的地层底面(如图2的
)，并为该线附上属性，连线的属性同位移前的连线的属性；
第三步，按照此方法，依次逐层恢复各时期的地层变形特征，在恢复下一层的地层时，移除上一层矩形(如图2的
-图2的
)，恢复至最古老的一层沉积岩结束。
可选的，地质构造为倾斜剥蚀地层，图4为本申请一实施例提供的倾斜剥蚀地层平衡剖面恢复过程示意图，如图4所示，该倾斜剥蚀地层的具体位移恢复步骤包括如下三步：
第一步：插入矩阵(如图4的
)，根据地层厚度变化与地层的接触关系判断抬升时期；
第二步，补充剥蚀量：剥蚀量的恢复按地层由老到新依次进行，恢复时，通过测量剖面上未剥蚀段的矩形高度，将其求和取平均值，所得到的数值代表同一层位未剥蚀段的矩形平均高度，将其视为被剥蚀段地层的原始沉积厚度，用该矩形对剥蚀区的矩形进行同底替换。若局部地层完全剥蚀，也按照此方式由老至新进行剥蚀量的恢复(如图4的
)；
第三步，恢复抬升剥蚀前的构造形态：恢复剥蚀量后的剖面进行阵列的垂向平移，使阵列顶部与0海拔线对齐，将相同属性矩阵底面进行连线，获得地质剖面即为构造抬升剥蚀前的原始地层剖面。
可选的，地质构造为断层活动引起的复杂地层，断层引起的复杂地层主要包括：逆冲推覆构造、伸展构造，图5为本申请另一实施例提供的断层活动引起的复杂地层平衡剖面恢复过程示意图，如图5所示，地层中包括被断裂错断的地层，首先根据地震资料，判断断层的活动时期，判断方法为分析被断裂错断的地层年代，被错断的最晚地层形成的年代，即为断裂的活动时期，然后恢复被断层错断的地层，将设置好的矩形阵列平移至地表，地表即为0海拔线，在垂向位移被断裂错断的地层的过程中，将被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；其中，被断裂错断的地层中逆断层内的矩形向地底方向移动(如图5的
-图5的
)，被断裂错断的地层中正断层内的矩形向地表方向移动(如图5的
-图5的
)，获得断裂前的原始地层剖面。
可选的，在上述实施例的基础上，获取当前时期地质结构剖面中目标地层的当前长度，目标地层为遭受剥蚀的地层或被断裂错断的地层；根据目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中目标地层的长度等于当前长度。
本实施例中，遭受剥蚀的地层或被断裂错断的地层在获得原始地层剖面后，需要做层长校正，例如：被断裂错断的地层，在错断断层活动前后，它的地层长度是不会发生改变的，改变的只是被错开的地层垂向上的错断位置(如图5的
)，恢复时若直接把地层平移至地表，地层的长度会发生误差，不符合平衡剖面恢复的层长守恒的原则，正断层若不校正，会使得地层的长度相较于未变形前变长，不准确(如图5的
)；逆断层若不校正，会使的地层的长度较未变形前变短，不准确(如图5的
)。具体校正步骤包括如下四步：
第一步，利用矩形阵列中矩形的调整角度，计算变形期前后地层的倾角变化。
图6为本申请另一实施例提供的矩形阵列在测地层倾角上的应用示意图，如图6的
所示，依据地层的变形程度，设置不同密度的矩形阵列为n列，这样地层线就被分为长度不等的n段，矩形阵列的每一段的长度为mi，若某一段内矩形倾斜角为θi，该段地层倾角为αi，则：
θi＝αi 公式一
某一区域的地层倾角平均值为：
第二步，对于遭受剥蚀的地层，可以根据对地层剥蚀量的恢复，计算该时期地层抬升角度、沉积-剥蚀界线。
如图4的
所示，若地层的最大抬升铅直高度为h，抬升变形段地层的水平长度为r，则地层抬升角度(即该变形段地层αi的平均值)：
α＝arctan(h/r) 公式三
第三步，计算剖面构造变形后的地层长度。
依据地层的变形程度，将地层线分为长度不等的n段，每一段的长度为mi，该处的地层倾角为αi。
变形后的地层长度：
第四步，根据剖面构造变形后的地层长度，即目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，即：用M1表示变形后的地层长度，即目标地层的当前长度为M1，用M2表示目标地层变形前的长度，则对于遭受剥蚀的地层，使得变形前的地层长度M2等于M1的长度，就完成了对剖面的层长校正(如图4的
-图4的
)；对于被断裂错断的地层，使得变形前的地层长度M2等于M1的长度，就完成了对剖面的层长校正(如图5的
和图5的
)，根据相同属性的矩形做出恢复后的地层线，删除断层线。
可选的，地层剖面完成平衡恢复后，需要给地层剖面中的每层地层填充颜色，每层地层用不同的颜色区分。
可选的，在上述实施例公式的基础上，还可以根据如下公式五计算剖面构造变形前的地层长度：
可选的，在上述实施例公式的基础上，还可以根据如下公式六和公式七计算剖面构造变形前后的伸展率、缩短率。
剖面伸展率：
R伸展＝(M2-M1)/M1*100％ 公式六
剖面缩短率：
R缩短＝(M1-M2)/M1*100％ 公式七
根据上述实施例中的公式，可以得到如表1所示的一地质结构剖面地层变形的定量分析数据示例，根据地层倾角、地层的抬升角度、地层的伸展率或缩短率可以定量表征地层的变形强度，通过对地层变形的定量分析数据可为后续油气勘探工作的盆地模拟、油藏模拟提供数据支撑，为认识构造发育史、分析油气成藏、聚集规律提供基础数据。
表1一地质结构剖面地层变形的定量分析数据示例
本实施例提供的地层平衡剖面恢复方法，通过获取目标区域深度域的当前时期地质结构剖面，当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数，在当前时期地质结构剖面中设置矩形阵列，矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向，从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面，其中，在垂向位移每层地层时，其它各层地层均随之垂向位移，通过上述方式，能够更加准确地对复杂构造变形地质结构进行平衡剖面的恢复，可对先存构造特征进行复原，认识构造发育史，解析构造活动的时期、活动的强度，有助于分析和判断油气的运移、聚集的时期，对于多期叠合盆地研究、超深层油气资源的勘探、开发工作均起到重要的促进作用。
本申请提供的地层平衡剖面恢复方法，能够更加准确地对复杂构造变形地质结构进行平衡剖面的恢复。塔里木盆地是由古生界的克拉通盆地与中新生界前陆盆地组成的一个大型多期叠合盆地，属于复杂构造变形地质结构。受到多期构造运动的影响，盆地内构造变形特征复杂，盆地由多个构造单元组成，包含塔北、塔中、塔东南等隆起构造单元；库车、满加尔、阿瓦提、塔西南、塔东南等坳陷构造单元。塔里木盆地复杂的构造变形与地层接触关系，给盆地平衡剖面的恢复带来了较大的困难。应用本方法，可对塔里木盆地不同尺度的剖面进行平衡剖面恢复，并对不同时期的构造变形做出定量分析。应用本方法可以达到以下效果：(1)针对百公里级的大尺度剖面，本方法可以有效恢复各层的厚度变化，准确表达各地质时期内地层的变形特征，保留局部构造形态，并且可定量分析地层变形量。(2)针对十公里的较小尺度剖面，特别是复杂构造变形区，可以通过高密度的矩形阵列恢复，保留层内变形细节，对变形前地层进行准确恢复，并对地层的倾角等产状进行测量。该方法在塔里木盆地区域基础地质研究，库车坳陷，塔北—塔中古隆起构造演化研究，构造对油气的控制作用研究中均发挥了重要的作用。图7为本申请一实施例提供的地层平衡剖面恢复方法的应用效果图，如图7所示，通过本申请提供的地层平衡剖面恢复方法，依次得到了古近系沉积前、白垩系沉积前、三叠-侏罗系沉积前、二叠系沉积前、石炭系沉积前、泥盆系沉积前、志留系沉积前共7个不同时期的地层平衡剖面，说明本申请提供的地层平衡剖面恢复方法在实际应用中，可以精确恢复复杂构造变形区的平衡剖面。图8a和图8b为本申请另一实施例提供的地层平衡剖面恢复方法的应用效果图，如图8a和图8b所示，是一幅横贯塔里木盆地东西向长达五百公里的大剖面，通过本申请提供的地层平衡剖面恢复方法，依次得到了古近系沉积前、白垩系沉积前(三叠纪剥蚀期)、三叠纪-侏罗纪沉积期、三叠系沉积前、二叠系沉积前、石炭系沉积前、泥盆系沉积前、志留系沉积前(晚奥陶世剥蚀期)、上奥陶统沉积期、上奥陶统沉积前、中下奥陶统沉积前、下寒武统沉积前、中寒武统沉积前、下寒武统沉积前共14个不同时期的地层平衡剖面，说明本申请提供的地层平衡剖面恢复方法在对大尺度的平衡剖面恢复中也同样的适用。
图9为本申请一实施例提供的地层平衡剖面恢复装置的结构示意图。如图9所示，本实施例的地层平衡剖面恢复装置900包括：获取模块901、设置模块902和处理模块903。
获取模块901，用于获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数。
设置模块902，用于在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向。
处理模块903，用于从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面。
其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
本实施例的装置，可以用于执行图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
在上述任一所示实施例的基础上，每层地层内的各矩形的顶边与所述地层的上地层线贴齐，且底边与所述地层的下地层线贴齐；所述矩形的高度与所在地层的厚度有关，所述矩形的宽度为所述矩形的高度的K倍，所述K为大于等于2的整数。
在上述任一所示实施例的基础上，相邻列的矩形间的距离与所述相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
在上述任一所示实施例的基础上，所述处理模块903，具体用于：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，所述第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层；将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，所述第j+1个时期的平衡恢复后的地质结构剖面包括所述N-(j+1)层不同属性的地层；其中，j为大于等于1且小于等于N-2的整数。
在上述任一所示实施例的基础上，所述设置模块902，具体用于：
若所述N-1层地层中包括遭受剥蚀的地层，则根据所述当前时期地质结构剖面中所述遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定所述遭受剥蚀的地层的原始沉积厚度；根据所述遭受剥蚀的地层的原始沉积厚度，为所述遭受剥蚀的地层设置剥蚀段的矩形。
在上述任一所示实施例的基础上，所述处理模块903，还用于：
若所述N-1层地层中包括被断裂错断的地层，则在垂向位移被断裂错断的地层的过程中，将所述被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；其中，所述被断裂错断的地层中逆断层内的矩形向地底方向移动，所述被断裂错断的地层中正断层内的矩形向地表方向移动。
在上述任一所示实施例的基础上，所述处理模块903，还用于：
获取所述当前时期地质结构剖面中目标地层的当前长度，所述目标地层为遭受剥蚀的地层或被断裂错断的地层；根据所述目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中所述目标地层的长度等于所述当前长度。
本实施例的装置，可以用于执行上述任一所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
图10为本申请另一实施例提供的地层平衡剖面恢复装置的结构示意图。如图10所示，本实施例的装置1000包括：存储器1001和处理器1002。其中，处理器1001、存储器1002通过总线连接。
所述存储器1001用于存储程序指令。
所述处理器1002用于调用所述存储器中的程序指令执行：
获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数；在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向；从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面；其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
可选的，每层地层内的各矩形的顶边与所述地层的上地层线贴齐，且底边与所述地层的下地层线贴齐；所述矩形的高度与所在地层的厚度有关，所述矩形的宽度为所述矩形的高度的K倍，所述K为大于等于2的整数。
可选的，相邻列的矩形间的距离与所述相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
可选的，所述处理器1002，具体用于：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，所述第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层；将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，所述第j+1个时期的平衡恢复后的地质结构剖面包括所述N-(j+1)层不同属性的地层；其中，j为大于等于1且小于等于N-2的整数。
可选的，所述处理器1002，具体用于：
若所述N-1层地层中包括遭受剥蚀的地层，则根据所述当前时期地质结构剖面中所述遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定所述遭受剥蚀的地层的原始沉积厚度；根据所述遭受剥蚀的地层的原始沉积厚度，为所述遭受剥蚀的地层设置剥蚀段的矩形。
可选的，所述处理器1002，还用于：
若所述N-1层地层中包括被断裂错断的地层，则在垂向位移被断裂错断的地层的过程中，将所述被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；其中，所述被断裂错断的地层中逆断层内的矩形向地底方向移动，所述被断裂错断的地层中正断层内的矩形向地表方向移动。
可选的，所述处理器1002，还用于：
获取所述当前时期地质结构剖面中目标地层的当前长度，所述目标地层为遭受剥蚀的地层或被断裂错断的地层；根据所述目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中所述目标地层的长度等于所述当前长度。
图11为本申请另一实施例提供的地层平衡剖面恢复装置的结构示意图，如图11所示，例如，地层平衡剖面恢复装置1100可以被提供为一服务器或计算机。参照图11，地层平衡剖面恢复装置1100包括处理组件1101，其进一步包括一个或多个处理器，以及由存储器1102所代表的存储器资源，用于存储可由处理组件1101执行的指令，例如应用程序。存储器1102中存储的应用程序可以包括一个或一个以上的每一个对应于一组指令的模块。此外，处理组件1101被配置为程序指令，以执行上述任一方法实施例。
地层平衡剖面恢复装置1100还可以包括一个电源组件1103被配置为执行装置1100的电源管理，一个有线或无线网络接口1104被配置为将装置1100连接到网络，和一个输入输出(I/O)接口1105。装置1100可以操作基于存储在存储器112的操作系统，例如Windows ServerTM，Mac OS XTM，UnixTM,LinuxTM，FreeBSDTM或类似。
本申请还提供一种计算机可读存储介质，所述计算机可读存储介质中存储有计算机程序指令，当处理器执行所述计算机程序指令时，实现如上地层平衡剖面恢复方法。
上述的计算机可读存储介质，上述可读存储介质可以是由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(SRAM)，电可擦除可编程只读存储器(EEPROM)，可擦除可编程只读存储器(EPROM)，可编程只读存储器(PROM)，只读存储器(ROM)，磁存储器，快闪存储器，磁盘或光盘。可读存储介质可以是通用或专用计算机能够存取的任何可用介质。
一种示例性的可读存储介质耦合至处理器，从而使处理器能够从该可读存储介质读取信息，且可向该可读存储介质写入信息。当然，可读存储介质也可以是处理器的组成部分。处理器和可读存储介质可以位于专用集成电路(Application Specific IntegratedCircuits，简称：ASIC)中。当然，处理器和可读存储介质也可以作为分立组件存在于地层平衡剖面恢复装置中。
本领域普通技术人员可以理解：实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时，执行包括上述各方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上各实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述各实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的范围。
Claims (10)
1.一种地层平衡剖面恢复方法，其特征在于，包括：
获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数；
在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向；
从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面；
其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
2.根据权利要求1所述的方法，其特征在于，每层地层内的各矩形的顶边与所述地层的上地层线贴齐，且底边与所述地层的下地层线贴齐；
所述矩形的高度与所在地层的厚度有关，所述矩形的宽度为所述矩形的高度的K倍，所述K为大于等于2的整数。
3.根据权利要求2所述的方法，其特征在于，相邻列的矩形间的距离与所述相邻列的矩形所在各地层段的长度、厚度和变形强度中至少一项有关。
4.根据权利要求1所述的方法，其特征在于，所述从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面，包括：
将当前时期地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层中的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第1个时期的平衡恢复后的地质结构剖面，所述第1个时期的平衡恢复后的地质结构剖面包括N-1层不同属性的地层；
将第j个时期的平衡恢复后的地质结构剖面中从地表到地底的顺序第2层地层的矩形，朝地表的方向垂向位移并且所述第2层地层下方的其它各层地层的矩形也随之垂向位移，直至所述第2层地层的矩形的顶边与所述地表贴齐，获得最接近当前时期的第j+1个时期的平衡恢复后的地质结构剖面，所述第j+1个时期的平衡恢复后的地质结构剖面包括所述N-(j+1)层不同属性的地层；
其中，j为大于等于1且小于等于N-2的整数。
5.根据权利要求1所述的方法，其特征在于，所述在所述当前时期地质结构剖面中设置矩形阵列，包括：
若所述N-1层地层中包括遭受剥蚀的地层，则根据所述当前时期地质结构剖面中所述遭受剥蚀的地层中未剥蚀段的各矩形的高度，确定所述遭受剥蚀的地层的原始沉积厚度；
根据所述遭受剥蚀的地层的原始沉积厚度，为所述遭受剥蚀的地层设置剥蚀段的矩形。
6.根据权利要求1所述的方法，其特征在于，还包括：
若所述N-1层地层中包括被断裂错断的地层，则在垂向位移被断裂错断的地层的过程中，将所述被断裂错断的地层中的矩形沿断裂错断面平行的方向位移；
其中，所述被断裂错断的地层中逆断层内的矩形向地底方向移动，所述被断裂错断的地层中正断层内的矩形向地表方向移动。
7.根据权利要求5或6所述的方法，其特征在于，还包括：
获取所述当前时期地质结构剖面中目标地层的当前长度，所述目标地层为遭受剥蚀的地层或被断裂错断的地层；
根据所述目标地层的当前长度，对垂向位移后的目标地层的长度进行校正，使得在N-1个不同时期的平衡恢复后的地质结构剖面中所述目标地层的长度等于所述当前长度。
8.一种地层平衡剖面恢复装置，其特征在于，包括：
获取模块，用于获取目标区域深度域的当前时期地质结构剖面，所述当前时期地质结构剖面包括N层不同属性的地层，相邻地层间以地层线区分，N为大于等于2的整数；
设置模块，用于在所述当前时期地质结构剖面中设置矩形阵列，所述矩形阵列包括(N-1)*M的矩形，M为大于等于1的整数，所述矩形阵列的(N-1)行矩形分别设置在最接近地表的地层之外的N-1层地层，每列的矩形位于同一垂向方向；
处理模块，用于从地表到地底的顺序，依次将每层地层的矩形朝地表的方向垂向位移，直至该层地层的矩形的顶边与所述地表贴齐，获得该层地层对应时期的平衡恢复后的地质结构剖面，以获得当前时期之前N-1个不同时期的平衡恢复后的地质结构剖面；
其中，在垂向位移每层地层时，其它各层地层均随之垂向位移。
9.一种地层平衡剖面恢复装置，其特征在于，包括：存储器和处理器；
所述存储器用于存储程序指令；
所述处理器用于调用所述存储器中的程序指令执行如权利要求1-7任一项所述的地层平衡剖面恢复方法。
10.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机程序指令，所述计算机程序指令被执行时，实现如权利要求1至7任一项所述的地层平衡剖面恢复方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010771313.2A CN111913232A (zh) | 2020-08-04 | 2020-08-04 | 地层平衡剖面恢复方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010771313.2A CN111913232A (zh) | 2020-08-04 | 2020-08-04 | 地层平衡剖面恢复方法和装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111913232A true CN111913232A (zh) | 2020-11-10 |
Family
ID=73287796
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010771313.2A Pending CN111913232A (zh) | 2020-08-04 | 2020-08-04 | 地层平衡剖面恢复方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111913232A (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4821242A (en) * | 1986-06-26 | 1989-04-11 | Hennington Willard M | Depositional reconstruction for petroleum location |
| CN105388526A (zh) * | 2015-10-28 | 2016-03-09 | 中国石油化工股份有限公司 | 一种陆相湖盆的古地质结构恢复方法 |
| CN106483579A (zh) * | 2015-08-28 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于获取断裂活动情况的方法 |
| CN106842293A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 基于三角洲前缘古地形特征的平衡剖面恢复方法 |
| CN107015290A (zh) * | 2017-03-13 | 2017-08-04 | 西北大学 | 一种改造型断陷盆地原始沉积面貌恢复的方法 |
| CN107742016A (zh) * | 2017-09-29 | 2018-02-27 | 中国石油化工股份有限公司 | 一种适用于上升盘地层缺失型断层的活动性计算方法 |
-
2020
- 2020-08-04 CN CN202010771313.2A patent/CN111913232A/zh active Pending
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4821242A (en) * | 1986-06-26 | 1989-04-11 | Hennington Willard M | Depositional reconstruction for petroleum location |
| CN106483579A (zh) * | 2015-08-28 | 2017-03-08 | 中国石油化工股份有限公司 | 一种用于获取断裂活动情况的方法 |
| CN105388526A (zh) * | 2015-10-28 | 2016-03-09 | 中国石油化工股份有限公司 | 一种陆相湖盆的古地质结构恢复方法 |
| CN106842293A (zh) * | 2015-12-04 | 2017-06-13 | 中国石油化工股份有限公司 | 基于三角洲前缘古地形特征的平衡剖面恢复方法 |
| CN107015290A (zh) * | 2017-03-13 | 2017-08-04 | 西北大学 | 一种改造型断陷盆地原始沉积面貌恢复的方法 |
| CN107742016A (zh) * | 2017-09-29 | 2018-02-27 | 中国石油化工股份有限公司 | 一种适用于上升盘地层缺失型断层的活动性计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 孙永壮 等: "《古潜山构造分析与储层裂缝预测 以济阳坳陷桩海地区为例》", 31 October 2006, 中国地质大学出版社 * |
| 王亚东 等: "平衡剖面反映的柴西新生代变形对青藏高原隆升的响应", 《兰州大学学报（自然科学版）》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mayolle et al. | Scaling of fault damage zones in carbonate rocks | |
| Serpelloni et al. | Vertical GPS ground motion rates in the Euro‐Mediterranean region: New evidence of velocity gradients at different spatial scales along the Nubia‐Eurasia plate boundary | |
| CN110795878B (zh) | 一种隧道涌水量预测方法 | |
| Anderson et al. | The Quaternary marine terraces of Santa Cruz, California: Evidence for coseismic uplift on two faults | |
| CN106646613B (zh) | 深度域多尺度井控建模与成像联合处理方法 | |
| CN105386756B (zh) | 一种应用应变量计算脆性地层孔隙度的方法 | |
| CN112147684B (zh) | 表征同沉积断层活动强度的方法和装置 | |
| CN105114063A (zh) | 一种斜井双侧向测井曲线校正方法 | |
| CN115324560A (zh) | 利用地应力场模拟确定压裂诱发油气套管变形位置的方法 | |
| CN105929452A (zh) | 基于地震数据预测地下裂缝空间展布的方法及装置 | |
| Novopashina et al. | Methodical approach to isolation of seismic activity migration episodes of the northeastern Baikal rift system (Russia) | |
| CN104564069A (zh) | 一种基于方格网法的地面动态沉陷预测与复垦方法 | |
| Little | Transpressive ductile flow and oblique ramping of lower crust in a two‐sided orogen: Insight from quartz grain‐shape fabrics near the Alpine fault, New Zealand | |
| CN115993103B (zh) | 一种采空区体积确定方法和装置 | |
| CN111913232A (zh) | 地层平衡剖面恢复方法和装置 | |
| Griffith et al. | Mechanical validation of the three-dimensional intersection geometry between the Puente Hills blind-thrust system and the Whittier fault, Los Angeles, California | |
| Konstantinou | Accurate relocation of seismicity along the North Aegean Trough and its relation to active tectonics | |
| CN111638552A (zh) | 一种古地貌恢复方法 | |
| CN114200538B (zh) | 地应力方向确定方法、装置、存储介质和计算机设备 | |
| CN114491774A (zh) | 一种深部背斜构造及地层结构三维数值模型构建方法 | |
| CN110618455B (zh) | 一种沉积盆地构造隆升的定量评价方法 | |
| CN113934806A (zh) | 一种土壤地球化学测量数据精细化处理方法 | |
| Zimmaro | Seismic response of the farneto del principe dam in Italy using hazard-consistent and site-specific ground motions | |
| CN114114404B (zh) | 一种用于挤压盆地的古地貌恢复方法 | |
| CN114076982B (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20201110 |
|
| RJ01 | Rejection of invention patent application after publication |