CN114996795A - 一种厘定地下走滑断裂滑动方向与走滑距离的方法 - Google Patents
一种厘定地下走滑断裂滑动方向与走滑距离的方法 Download PDFInfo
- Publication number
- CN114996795A CN114996795A CN202210396893.0A CN202210396893A CN114996795A CN 114996795 A CN114996795 A CN 114996795A CN 202210396893 A CN202210396893 A CN 202210396893A CN 114996795 A CN114996795 A CN 114996795A
- Authority
- CN
- China
- Prior art keywords
- sliding
- fracture
- distance
- sliding fracture
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000011084 recovery Methods 0.000 claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 25
- 238000002679 ablation Methods 0.000 claims description 15
- 238000011160 research Methods 0.000 claims description 13
- 238000004590 computer program Methods 0.000 claims description 12
- 238000013519 translation Methods 0.000 claims description 10
- 230000008859 change Effects 0.000 claims description 9
- 230000011218 segmentation Effects 0.000 claims description 9
- 230000003628 erosive effect Effects 0.000 claims description 5
- 230000004044 response Effects 0.000 claims description 4
- 238000009941 weaving Methods 0.000 claims description 4
- 230000002829 reductive effect Effects 0.000 claims description 3
- 239000010410 layer Substances 0.000 description 145
- UONOETXJSWQNOL-UHFFFAOYSA-N tungsten carbide Chemical compound [W+]#[C-] UONOETXJSWQNOL-UHFFFAOYSA-N 0.000 description 10
- 238000010586 diagram Methods 0.000 description 7
- 239000000758 substrate Substances 0.000 description 6
- 230000000694 effects Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 230000008021 deposition Effects 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 241000272814 Anser sp. Species 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 230000003247 decreasing effect Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000012791 sliding layer Substances 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000006071 cream Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000002441 reversible effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明公开了一种厘定地下走滑断裂滑动方向与走滑距离的方法、装置、计算机设备和存储介质，所述方法包括以下步骤：获取走滑断裂的断穿层位以及走滑断裂的走向，在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线，确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复，按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原始厚度‑距离测线一侧末端的距离”折线图，基于每组平行于走滑断裂走向的测线的折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏移量估算走滑距离。上述方法可广泛应用于地下走滑断裂滑动方向与走滑距离的厘定，且仅需常规测井数据与地震数据，估算成本低廉，具有较高的可操作性。
Description
技术领域
本申请涉及地下走滑断裂滑动技术领域，特别是涉及一种厘定地下走滑断 裂滑动方向与走滑距离的方法、装置、计算机设备和存储介质。
背景技术
地下走滑断裂系统的构造特征对古海相盆地的储层形成和油气聚集具有 非常重要的影响(Deng S,Li H,Zhang Z,et al.Structural characterization ofintracratonic strike-slip faults in the central Tarim Basin[J].AAPG Bulletin,2018,103(1):109-137)。目前对地下走滑 断裂的研究主要集中于走滑断裂几何学特征的刻画，系统研究走滑断裂运动学 特征的文献很少。不同于沿倾向活动且可根据断层两盘地层配置关系判断滑动 方向的正断层或逆断层，走滑断裂以较为显著的沿走向滑动为特征，对其滑动 方向及走滑距离的判定是当前主要技术难点之一。目前判断地下走滑断裂滑动方向及走滑距离的方法有以下三种：1)依据走滑断裂上覆雁列正断层的展布 方式判定走滑断裂的滑动方向。雁列正断层组的排列角定义为雁列正断层的走 向与下伏剪切断裂走向之间的锐夹角(Martin,E.S.,2016,The distribution and characterizationofstrike-slip faults on Enceladus: Geophysical Research Letters,v.43,no.6,p.2456–2464, doi:10.1002/2016GL067805)，雁列正断层组排列角位于下伏剪切断裂的一侧 盘时，该侧盘的滑动方向与排列角指向一致。2)两个相邻的走滑断裂形成的 断层组按照其不同的排列方式(左阶或右阶)及断层组中心地层的变形特征(下 沉或抬升)可判定其滑动方向(Cembrano J,G Gonz á lez,Arancibia G, et al.Fault zone developmentand strain partitioning in an extensional strike-slip duplex:A case studyfrom the Mesozoic Atacama fault system, Northern Chile[J].Tectonophysics,2005,400(1-4):105-125)。左旋走滑 运动可使左阶排列断层中心的地层下沉，使右阶排列断层中心的地层抬升，而 右旋走滑运动可使左阶排列断层中心的地层抬升，使右阶排列断层中心的地层 下沉。3)依据水平时间切片小波特征的偏移标记判定断层面的切割及走滑断 裂的滑动方向和走滑距离(Deng S,Li H,Zhang Z,et al.Structuralcharacterization of intracratonic strike-slip faults in the central TarimBasin[J].AAPG Bulletin,2018,103(1):109-137)。当在水平时间切 片上观察到由于受早期倾斜和走滑变形影响的前积层而略微倾斜的上覆岩层 时，可通过恢复小波特征的原始几何形态判定走滑断裂的滑动方向和滑动距 离。
前人虽然已使用该三种方法对盆地地下走滑断裂滑动方向与走滑距离的 判识做出了尝试，但方法本身的运用具有如下三点局限性：1)当以厚层膏盐、 泥岩等为代表的塑性应力释放层(例如四川盆地东南部中-下寒武统高台组、 龙王庙组膏盐滑脱层)位于深层走滑断裂上方时，走滑断裂在某些时期的剪切 形变能量被塑性层吸收，难以在其上覆地层中形成雁列式正断层组，无法依据 浅层变形特征判定下部剪切断裂滑动方向。2)晚期强烈的构造活动可能使得 两个相邻的走滑断裂形成的断层组中心地层变形特征不能反映走滑断裂的运 动学性质。3)以水平时间切片小波特征偏移标记为例的涉及线性特征(如地 层标志或者古河道)的断错体贯穿点很难获得，这使得地下走滑断裂滑动方向 与走滑距离的判定具有一定误差，往往是推测性的。
发明内容
基于此，针对上述技术问题，提供一种厘定地下走滑断裂滑动方向与走滑 距离的方法、装置、计算机设备和存储介质。
第一方面，一种厘定地下走滑断裂滑动方向与走滑距离的方法，所述方法 包括：
S1、获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿 层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂 的延伸长度；
S2、在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线；
S3、确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
S4、按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原 始厚度-距离测线一侧末端的距离”折线图；
S5、基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线 一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏 移量估算走滑距离。
上述方案中，可选地，所述S1包括：
若走滑断裂断穿了某地层，则确定某地层为走滑断裂的断穿层位，若沿走 滑断裂的走向其方位角恒定，则认为走滑断裂走向不变，若沿走滑断裂的走向 其方位角变化，则根据方位角的变化特征将走滑断裂分为数段。
上述方案中，进一步可选地，所述S2包括：
基于走滑断裂的走向确定选取测线的走向，基于走滑断裂的延伸长度确定 选取测线的延伸长度；其中，所选测线的长度与走滑断裂的延伸长度相等，基 于走滑断裂的延伸长度及响应用户第一输入，确定各测线间的垂直距离及测线 数量，所述用户第一输入根据研究精度要求设定；具体为，当走滑断裂的延伸 长度较大时，所选测线之间的垂直距离应适当增大，当走滑断裂的延伸长度较 小时，所选测线之间的垂直距离应适当减小，当研究精度要求较高时，所选测 线之间的垂直距离应适当增大；当研究精度要求较低时，所选测线之间的垂直 距离应适当增大。
将距离走滑断裂的垂直距离相等且分别位于走滑断裂两侧的测线编为一 组，并按照距离走滑断裂的垂直距离从小到大对每一组进行编号。
上述方案中，进一步可选地，所述基于走滑断裂的走向确定选取测线的走 向具体为：若走滑断裂的走向恒定，则所选测线的走向恒定，且与走滑断裂相 同；若走滑断裂的走向变化且具有分段特征，则所选测线的走向变化且与走滑 断裂具有相同分段特征。
上述方案中，进一步可选地，所述S3包括：
获取走滑断裂断穿层位的地震同相轴清晰程度，响应用户的第二输入，选 取目的层的顶面或底面；具体的，地震同相轴清晰程度包括物理清晰程度与地 质清晰程度，所述用户第二输入为物理清晰程度与地质清晰程度都较高的断穿 层位。对选取的目的层顶面和底面开展层位追踪，对走滑断裂两侧局部范围剥 蚀量恢复。
上述方案中，进一步可选地，所述S4包括：按照同一方向获取已选取测 线的残余地层厚度；其中，所述残余地层厚度为选取位置目的层的顶面深度与 底面深度差的绝对值，将残余地层厚度与选取位置的地层剥蚀量相加，按照同 一方向获取已选取测线的原始地层厚度，以获取测线地层厚度方向的起点为原 点，编制“目的层原始厚度-距离测线一侧末端的距离”折线图；其中，同一 组测线的“目的层原始厚度-距离测线一侧末端的距离”折线置于同一折线图 中进行比较。
上述方案中，进一步可选地，所述S5包括：通过每一组中分别位于走滑 断裂两侧测线的“目的层原始厚度-距离测线一侧末端的距离”折线之间的平 移后是否重合性判定走滑断裂的滑动方向；其中，若位于走滑断裂某一侧测线 的“目的层原始厚度-距离测线一侧末端的距离”折线经平移后可与另一侧测 线的“目的层原始厚度-距离测线一侧末端的距离”折线重合，则相反方向为 走滑断裂的滑动方向，基于每一组中“目的层原始厚度-距离测线一侧末端的 距离”折线的几何特征厘定特征点位置；根据每一组中“目的层原始厚度-距 离测线一侧末端的距离”折线特征点的偏移量估算断裂的走滑距；其中，特征 点的选取为“目的层原始厚度-距离测线一侧末端的距离”折线上的最大值点、 最小值点、极大值点、极小值点中的一种或多种。
第二方面，一种厘定地下走滑断裂滑动方向与走滑距离的装置，所述装置 包括：
获取单元：用于获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑 断裂的断穿层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确 定走滑断裂的延伸长度；在走滑断裂两侧等间隔选取数组平行于走滑断裂走向 的测线；确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
编制单元：用于按照同一方向获取已选取测线的目的层原始厚度，并编制 “目的层原始厚度-距离测线一侧末端的距离”折线图；
估算单元：用于基于每组平行于走滑断裂走向的测线的“目的层原始厚度 -距离测线一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特 征点位置偏移量估算走滑距离。
第三方面，一种计算机设备，包括存储器和处理器，所述存储器存储有计 算机程序，所述处理器执行所述计算机程序时实现以下步骤：
S1、获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿 层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂 的延伸长度；
S2、在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线；
S3、确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
S4、按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原 始厚度-距离测线一侧末端的距离”折线图；
S5、基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线 一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏 移量估算走滑距离。
第四方面，一种计算机可读存储介质，其上存储有计算机程序，所述计算 机程序被处理器执行时实现以下步骤：
S1、获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿 层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂 的延伸长度；
S2、在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线；
S3、确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
S4、按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原 始厚度-距离测线一侧末端的距离”折线图；
S5、基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线 一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏 移量估算走滑距离。
本发明至少具有以下有益效果：
本发明基于对现有技术问题的进一步分析和研究，通过获取走滑断裂的断 穿层位以及走滑断裂的走向，在走滑断裂两侧等间隔选取数组平行于走滑断裂 走向的测线，确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢 复，按照同一方向获取已选取测线的目的层原始厚度，编制“目的层原始厚度 -距离测线一侧末端的距离”折线图，基于每组平行于走滑断裂走向的测线的 “目的层原始厚度-距离测线一侧末端的距离”折线图判定断裂滑动方向，基 于地层厚度曲线特征点位置偏移量估算走滑距离。此方法可广泛应用于地下走 滑断裂滑动方向与走滑距离的厘定，且仅需常规测井数据与地震数据，估算成 本低廉，具有较高的可操作性。
附图说明
图1为本发明一个实施例提供的厘定地下走滑断裂滑动方向与走滑距离 的方法的应用环境示意图；
图2是右行走滑断裂两侧地层厚度变化的几何学与运动学模型。
图3是根据本发明的实施例的走滑断裂F1平面展布及部分测线位置图。
图4是根据本发明的实施例的测线1的地层-断裂-构造综合解释剖面图。
图5是根据本发明的实施例的测线2的地层-断裂-构造综合解释剖面图。
图6是根据本发明的实施例的测线3的地层-断裂-构造综合解释剖面图。
图7是根据本发明的实施例的位于走滑断裂F1两侧且平行于其走向的测 线平面分布图。
图8是根据本发明的实施例的选定测线目的层顶面与底面追踪及厚度计 算模式图。
图9、图10、图11、图12、图13是根据本发明的实施例的选定测线目的 层原始厚度曲线图及特征点位置。
图14为一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实 施例，对本申请进行进一步详细说明。应当理解，此处描述的具体实施例仅仅 用以解释本申请，并不用于限定本申请。
本申请提供的厘定地下走滑断裂滑动方向与走滑距离的方法，可
在一个实施例中，如图1所示，提供了一种厘定地下走滑断裂滑动方向与 走滑距离的方法，包括以下步骤：
S1、获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿 层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂 的延伸长度。
具体的，若走滑断裂断穿了某地层，则确定某地层为走滑断裂的断穿层位， 若沿走滑断裂的走向其方位角恒定，则认为走滑断裂走向不变，若沿走滑断裂 的走向其方位角变化，则根据方位角的变化特征将走滑断裂分为数段。
在一个实施例中，所述步骤S1包括：
S11、查明走滑断裂的断穿层位；
S12、明确走滑断裂的走向；
S13、估算走滑断裂的延伸长度。
在一个实施例中，所述步骤S11中，若走滑断裂断穿了某些地层，则认为 该地层沉积时期的构造运动对走滑断裂的形成或演化具有影响或控制作用，作 为走滑断裂的形成期或主要活动期；若走滑断裂未断穿部分地层，则认为该部 分地层沉积时期的构造运动对走滑断裂的形成与演化的影响较低，不作为走滑 断裂的形成期或主要活动期。
在一个实施例中，所述步骤S12中，若沿走滑断裂的走向其方位角恒定， 则认为走滑断裂走向不变；若沿走滑断裂的走向其方位角变化，则应根据方位 角的变化特征将走滑断裂分为数段，每段的走向应大致相同。
S2、在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线。
具体的，基于走滑断裂的走向确定选取测线的走向具体为：若走滑断裂的 走向恒定，则所选测线的走向恒定，且与走滑断裂相同；若走滑断裂的走向变 化且具有分段特征，则所选测线的走向变化且与走滑断裂具有相同分段特征。
基于走滑断裂的走向确定选取测线的走向，基于走滑断裂的延伸长度确定 选取测线的延伸长度；其中，所选测线的长度与走滑断裂的延伸长度相等，基 于走滑断裂的延伸长度及响应用户第一输入，确定各测线间的垂直距离及测线 数量，所述用户第一输入根据研究精度要求设定；具体为，当走滑断裂的延伸 长度较大时，所选测线之间的垂直距离应适当增大，当走滑断裂的延伸长度较 小时，所选测线之间的垂直距离应适当减小，当研究精度要求较高时，所选测 线之间的垂直距离应适当增大；当研究精度要求较低时，所选测线之间的垂直 距离应适当增大。
将距离走滑断裂的垂直距离相等且分别位于走滑断裂两侧的测线编为一 组，并按照距离走滑断裂的垂直距离从小到大对每一组进行编号。
在一个实施例中，所述步骤S2包括：
S21、基于走滑断裂的走向明确选取测线的走向。
S22、基于走滑断裂的延伸长度明确选取测线的延伸长度。
S23、基于走滑断裂的延伸长度及研究精度要求，明确各测线间的垂直距 离及测线数量。
S24、将距离走滑断裂的垂直距离相等且分别位于走滑断裂两侧的测线编 为一组，并按照距离走滑断裂的垂直距离从小到大对每一组进行编号。
在一个实施例中，步骤S21中，若走滑断裂的走向恒定，则所选测线的走 向恒定，且与走滑断裂相同；若走滑断裂的走向变化且具有分段特征，则所选 测线的走向变化且与走滑断裂具有相同分段特征。
在一个实施例中，步骤S22中，所选测线的长度应与走滑断裂的延伸长度 大致相等。
在一个实施例中，步骤S23中，当走滑断裂的延伸长度较大时，所选测线 之间的垂直距离应适当增大；当走滑断裂的延伸长度较小时，所选测线之间的 垂直距离应适当减小；当研究精度要求较高时，所选测线之间的垂直距离应适 当增大；当研究精度要求较低时，所选测线之间的垂直距离应适当增大。
S3、确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复。
具体的，获取走滑断裂断穿层位的地震同相轴清晰程度，响应用户的第二 输入，选取目的层的顶面或底面；具体的，地震同相轴清晰程度包括物理清晰 程度与地质清晰程度，所述用户第二输入为物理清晰程度与地质清晰程度都较 高的断穿层位。对选取的目的层顶面和底面开展层位追踪，对走滑断裂两侧局 部范围剥蚀量恢复。
在一个实施例中，所述步骤S3包括：
S31、评估走滑断裂断穿层位的地震同相轴清晰程度；
S32、选取目的层并对目的层顶面与底面开展层位追踪；
S33、利用适当方法开展走滑断裂两侧局部范围剥蚀量恢复。
在一个实施例中，所述步骤S31中，地震同相轴清晰程度应分为地球物理 方面的清晰程度与地质方面的清晰程度，地球物理方面的清晰程度即同相轴清 晰可辨，平面连续性较强，具有较好的可追踪性，地质方面的清晰程度即无地 层重复等复杂地质现象，以减少层位追踪过程中的不确定因素。地震同相轴清 晰程度较高的地层界面应作为目的层的顶面或底面。
在一个实施例中，所述步骤S33中，目的层若涉及不整合面(既包括目的 层内的界面，也包括目的层的顶面，但不包括目的层的底面)，则应根据不整 合面所代表的地层剥蚀厚度大小及平面变化规律判定是否需要开展剥蚀量恢 复。若地层剥蚀厚度较小，且剥蚀量平面分布稳定，则无需开展剥蚀量恢复； 若地层剥蚀厚度较小，但剥蚀量平面变化较为显著，则仍需开展剥蚀量恢复； 若地层厚度剥蚀量较大，则需开展剥蚀量恢复。进一步的，目的层应尽可能选 取不包含不整合面或不整合面剥蚀厚度较小且分布稳定的地层。
S4、按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原 始厚度-距离测线一侧末端的距离”折线图。
具体的，按照同一方向获取已选取测线的残余地层厚度；其中，所述残余 地层厚度为选取位置目的层的顶面深度与底面深度差的绝对值，将残余地层厚 度与选取位置的地层剥蚀量相加，按照同一方向获取已选取测线的原始地层厚 度，以获取测线地层厚度方向的起点为原点，编制“目的层原始厚度-距离测 线一侧末端的距离”折线图；其中，同一组测线的“目的层原始厚度-距离测 线一侧末端的距离”折线置于同一折线图中进行比较。
在一个实施例中，所述步骤S4包括
S41、按照同一方向获取已选取测线的残余地层厚度(等同于现今地层厚 度)；
S42、将残余地层厚度与相应位置的地层剥蚀量(若有)相加，按照同一 方向获取已选取测线的原始地层厚度；
S43、以获取测线地层厚度方向的起点为原点，编制“目的层原始厚度- 距离测线一侧末端的距离”折线图。
在一个实施例中，所述步骤S41中，残余地层厚度为某一位置目的层的顶 面深度与底面深度差的绝对值。
在一个实施例中，所述步骤S43中，同一组测线的“目的层原始厚度-距 离测线一侧末端的距离”折线应置于同一折线图中进行比较。
S5、基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线 一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏 移量估算走滑距离。
具体的，通过每一组中分别位于走滑断裂两侧测线的“目的层原始厚度- 距离测线一侧末端的距离”折线之间的平移后是否重合性判定走滑断裂的滑动 方向；其中，若位于走滑断裂某一侧测线的“目的层原始厚度-距离测线一侧 末端的距离”折线经平移后可与另一侧测线的“目的层原始厚度-距离测线一 侧末端的距离”折线重合，则相反方向为走滑断裂的滑动方向，基于每一组中 “目的层原始厚度-距离测线一侧末端的距离”折线的几何特征厘定特征点位 置；根据每一组中“目的层原始厚度-距离测线一侧末端的距离”折线特征点 的偏移量估算断裂的走滑距；其中，特征点的选取为“目的层原始厚度-距离 测线一侧末端的距离”折线上的最大值点、最小值点、极大值点、极小值点中 的一种或多种。
在一个实施例中，所述步骤S5包括：
S51、通过每一组中分别位于走滑断裂两侧测线的“目的层原始厚度-距离 测线一侧末端的距离”折线之间的平移相似性判定走滑断裂的滑动方向；
S52、基于每一组中“目的层原始厚度-距离测线一侧末端的距离”折线的 几何特征厘定特征点位置；
S53、基于每一组中“目的层原始厚度-距离测线一侧末端的距离”折线特 征点的偏移量估算断裂的走滑距。
在一个实施例中，所述步骤S51中，若位于走滑断裂某一侧测线的“目的 层原始厚度-距离测线一侧末端的距离”折线经平移后可与另一侧测线的“目 的层原始厚度-距离测线一侧末端的距离”折线大致重合，则认为前者所在的 走滑断裂一盘沿相同方向平移相同的距离后可与后者所在的走滑断裂一盘对 接，即相反方向为走滑断裂的滑动方向。
在一个实施例中，所述步骤S52中，特征点的选取应为“目的层原始厚度-距离测线一侧末端的距离”折线上的最大值点、最小值点、极大值点、极小 值点中的一种或多种，且由于走滑断裂的走滑距往往在中部最大，故所选标志 点应尽量位于折线的中部，不应位于折线的端部。
在一个实施例中，所述步骤S53中，当走滑断裂一侧测线的“目的层原始 厚度-距离测线一侧末端的距离”折线中部的特征点沿某方向平移一定距离后 可与位于走滑断裂另一侧测线的“目的层原始厚度-距离测线一侧末端的距离” 折线中部特征点位置重合，则认为前者所在的走滑断裂一盘沿相同方向平移相 同的距离后可与后者所在的走滑断裂一盘对接，即相反方向为走滑断裂的滑动 方向，该平移距离为断裂的走滑距。计算出各组的走滑距后，应对各组走滑距 的计算结果进行可靠性评估，主要基于同组中两条测线的目的层原始厚度曲线 形态相似程度和走滑距的横向变化特征：当同组中两条测线的目的层原始厚度 曲线形态相似程度较高且走滑距呈现断裂中央最高，断裂两端最低且为0的特 征时，计算出的走滑距具有较高的可靠程度；当同组中两条测线的目的层原始 厚度曲线形态相似程度较低或走滑距不满足断裂中央最高，断裂两端最低且为 0的特征时，计算出的走滑距可靠程度中等；当同组中两条测线的目的层原始 厚度曲线形态相似程度较低且走滑距不满足断裂中央最高，断裂两端最低且为 0的特征时，计算出的走滑距可靠程度较低。
上述方法主要基于沉积地层被走滑断裂错开后，其分别位于走滑断裂两侧 的地层厚度既存在相似性也存在差异性这一特征，相似性在于走滑断裂两侧地 层原始厚度的局部高值和局部低值仍可大致对应，差异性在于走滑断裂两侧地 层原始厚度沿断裂走向变化曲线的极大值点、极小值点等特殊点位置产生偏 移，且该偏移的方向反映走滑断裂的滑动方向，该偏移的距离反映走滑断裂的 走滑距离。附图2(a)中，Fm.B代表西薄东厚的目的层(Formation B)，用较 为稀疏的黑点填充，Fm.A与Fm.C分别为目的层的上覆地层和下伏地层。右行 走滑断层FAULT使得南侧断盘(BlockS，简写为B.S)目的层Fm.B的较厚处 与北侧断盘(BlockN，简写为B.N)目的层Fm.B的较薄处对应，如附图2(b) 所示，剖面xx’(位置见附图2(a))为垂直走滑断裂FAULT走向的近南北向剖 面，表明走滑断裂两侧地层厚度突变，且B.S断盘目的层Fm.B比B.N断盘目 的层Fm.B厚，反映该情况下与走滑断裂走向垂直的剖面的目的层Fm.B厚度变 化特征。附图2(c)所示为走滑断裂两侧目的层Fm.B厚度等值线图(采用时间 厚度，单位为毫秒，ms)，B.N断盘内目的层Fm.B厚度从西至东单调递增，由10ms递增至50ms；南侧目的层厚度亦单调递增，且由20ms递增至60ms，走 滑断裂两侧目的层厚度变化趋势相同，但相同厚度的等值线在B.N断盘和B.S 断盘内呈现明显的平移特征，B.N断盘目的层厚度等值线向左平移L.S.后可与 B.S断盘目的层厚度等值线大致重合，表明在地层Fm.C沉积过程中或之后的 某个地质历史时期形成了右行走滑断裂，使得原本南北走向的目的层Fm.B厚 度等值线被东西向错开且具备L.S.的右行走滑距从而呈现现今状态，故走滑 断裂两侧目的层Fm.B厚度的变化特征可反映走滑断裂的运动学特征，但当目的层遭受剥蚀后其厚度变化不能反映由走滑断裂导致的初始走滑距，故应采用 目的层的原始厚度。本发明通过获取与走滑断裂平行且分别位于走滑断裂两侧 的测线的目的层原始厚度，基于其曲线特征与特征点位置偏移量定量表征走滑 断裂的水平位移量，从而明确走滑断裂的滑动方向与走滑距离，对走滑断裂形 成与演化模式的建立具有较强的理论指导意义。
在一个具体的实施例中，该实例判定了四川盆地南部NE向走滑断裂F1 的滑动方向与走滑距离。走滑断裂F1在志留系底界面(记为TS界面)相干体 属性切片的响应如附图3(a)所示，可见走滑断裂F1为近NE向断裂，且方位 角稳定，无明显分段特征，总长约为7km.其北部受限于NEE向基底走滑断裂 F3，南侧在深部与NEE向基底走滑断裂F2相互交切。为明确走滑断裂F1的断 穿层位，选取三条近NS向剖面开展地层-断裂-构造精细解释，自西到东分别 为测线1、测线2、测线3，具体位置见附图3(b).
附图4为测线1的地层-断裂-构造精细解释剖面图，可见走滑断裂F2和 F3为近直立的基底走滑断裂，F1为F3的分支断裂。附图4中，走滑断裂F1 深部发育至上震旦统灯影组底界面(记为TZ2dy)以下，但未直插基底，而是 终止于走滑断裂F3，其浅部断穿上寒武统洗象池群底界面(记为TC3x)，但未 断穿至下奥陶统底界面(记为TO)；走滑断裂F2深部直插基底，浅部断穿下 寒武统龙王庙组底界面(记为TC1l)，但未断穿上寒武统洗象池群底界面(记 为TC3x)；走滑断裂F3深部直插基底，浅部断穿下奥陶统底界面(记为TO) 但未断穿志留系底界面(记为TS)。
附图5为测线2的地层-断裂-构造精细解释剖面图，可见走滑断裂F1深 部发育至TZ2dy界面附近并终止于走滑断裂F3，浅部断穿TS界面但并未断穿 二叠系底界面(记为TP1)；走滑断裂F2深部直插基底，浅部断穿TS界面但并 未断穿TP1界面；走滑断裂F3深部直插基底，浅部断穿TS界面但并未断穿TP1界面。
附图6为测线3的地层-断裂-构造精细解释剖面图，可见走滑断裂F1深 部发育至下寒武统筇竹寺组底界面(记为TC1q)附近并终止于走滑断裂F3， 浅部断穿TS界面但并未断穿TP1界面；走滑断裂F2深部直插基底，浅部断穿 TS界面但并未断穿TP1界面；走滑断裂F3深部直插基底，浅部断穿TS界面但 并未断穿TP1界面。
综合多条地震测线的地层-断裂-构造精细解释，明确了走滑断裂F1始终 断穿的地层底面为TC1q，顶面为TC3x.
附图7(a)所示为基于走滑断裂F1的平面展布特征及测线选取，根据三维 地震资料精度及研究需要，在走滑断裂F1南北两侧分别选取了5条与走滑断 裂近平行的测线，且测线长度与走滑断裂F1延伸长度大致相同。附图7(b)所 示为近南北向拉长的走滑断裂F1及选取测线平面图，走滑断裂F1北部测线由 南至北分别命名为N1、N2、N3、N4、N5，南部测线由北至南分别命名为S1、 S2、S3、S4、S5，相邻测线之间的垂直距离均相等且为△X＝60m，每条测线的 西侧末端节点命名为“测线号-A”，东侧末端节点命名为“测线号-B”，例如测 线N1的西侧末端节点为N1-A，东侧末端节点为N1-B.按照分别位于走滑断裂 两侧且与走滑断裂的距离相等的测线编为一组这一原则，将上述10条测线编 为5组，第一组为测线N1和测线S1，第二组为测线N2和测线S2，第三组为 测线N3和测线S3，第四组为测线N4和测线S4，第五组为测线N5和测线S5， 5组测线距离走滑断裂的距离由小到大。
在走滑断裂F1必定断穿的所有层位中，TC1q界面与TC3x界面的地震同相 轴清晰度较好且区域上可追踪，故选取TC1q界面为目的层的底面，选取TC3x 界面为目的层的顶面，目的层为TC1q界面与TC3x界面之间的地层。对上述5 组共10条测线的TC1q界面与TC3x界面的追踪如附图8所示，其中附图8(a) 和附图8(b)分别为测线N1和测线S1即第1组两条测线的层位追踪结果，附 图8(c)和附图8(d)分别为测线N2和测线S2即第2组两条测线的层位追踪结 果,附图8(e)和附图8(f)分别为测线N3和测线S3即第3组两条测线的层位 追踪结果，附图8(g)和附图8(h)分别为测线N4和测线S4即第4组两条测线 的层位追踪结果，附图8(i)和附图8(j)分别为测线N5和测线S5即第5组两 条测线的层位追踪结果。以测线N1的目的层顶面和底面追踪结果(附图8(a)) 为例说明目的层厚度读取方法：附图8(a)中共有4个竖直且带有双向箭头的 线段，以最左侧为例，线段上端为该处的目的层顶面TC3x，线段下端为该处的 目的层底面TC1q，线段上端的时间深度为Z1＝1849ms，线段下端的时间深度为 Z2＝1979ms，那么该处目的层的时间地层厚度为Z2-Z1＝130ms.同样的，所有 已选测线各处的目的层厚度均可获取，该厚度为残余地层厚度。TC1q界面与 TC3x界面之间的地层沉积时期并未遭受较为严重的剥蚀，其间也无较大规模的 不整合面，故认为本具体的实施例中，目的层残余厚度等于目的层原始厚度。
按照从西到东的顺序获取上述5组共10条测线的目的层原始厚度如附图 9、图10、图11、图12、图13所示。附图9为第1组测线的“目的层原始厚 度-距离测线一侧末端的距离”折线图，折线图的横轴为距离测线西侧末端的 距离，0处即为测线的西侧末端节点A，每条测线长为7km，故横坐标为7km 的点为每条测线的东侧末端节点B.蓝线为N1测线的目的层原始厚度曲线， 黄线为S1测线的目的层原始厚度曲线。可以较为明显地看出两条曲线的形态 较为类似但将蓝线向右平移一定距离后可与黄线重合，所以走滑断裂F1北侧 N1测线目的层厚度变化曲线可看作由走滑断裂F1南侧S1测线目的层厚度变 化曲线向左平移获得，故依据第1组测线的目的层原始厚度曲线变化特征可判 定断裂F1为左行走滑断裂。蓝线与黄线在曲线中央即3km～4km处均拥有一个 非常明显的特征点，如附图9中两条曲线上的红点所示，红点代表两条曲线在 中央部位的极大值点，即走滑断裂中央部位两侧目的层的厚度最大值所在位 置，且两条曲线上的红点均为目的层厚度快速增大至较慢速减小的拐点，拥有 极为类似的数学特征，故可作为第1组测线的特征点。容易获取蓝线中央红点 的横坐标为3.26km，黄线中央红点的横坐标为3.45km，蓝线上的红点相较黄 线上的红点向左偏移了0.19km，既与两条曲线整体的偏移特征相符合，同时 也定量计算出了特征点的偏移量，该偏移量代表由第1组测线获取的走滑断裂F1左行走滑距为0.19km.
附图10为第2组测线的目的层原始厚度曲线图，与第1组测线的目的层 原始厚度曲线图相比，第2组测线的原始地层厚度曲线平移特征也较为明显， 且也反映断裂F1的左行走滑特征。北侧测线目的层原始厚度曲线特征点横坐 标为3.19km，南侧测线目的层原始厚度曲线特征点横坐标为3.52km，故由第 2组测线获取的走滑断裂F1左行走滑距为0.33km.附图11为第3组测线的目 的层原始厚度曲线图，与第2组测线的目的层原始厚度曲线图相比，第3组测 线的原始地层厚度曲线平移特征较为明显且有平移量更大的趋势，反映断裂 F1的左行走滑特征。北侧测线目的层原始厚度曲线特征点横坐标为3.12km， 南侧测线目的层原始厚度曲线特征点横坐标为3.55km，故由第2组测线获取 的走滑断裂F1左行走滑距为0.43km.附图12为第4组测线的目的层原始厚 度曲线图，与第3组测线的目的层原始厚度曲线图相比，第4组测线的原始地 层厚度曲线平移特征更为明显且曲线平移量突然增大，反映断裂F1的左行走 滑特征。北侧测线目的层原始厚度曲线特征点横坐标为2.94km，南侧测线目 的层原始厚度曲线特征点横坐标为3.79km，故由第4组测线获取的走滑断裂 F1左行走滑距为0.85km.附图13为第5组测线的目的层原始厚度曲线图，与 第4组测线的目的层原始厚度曲线图相比，第5组测线的原始地层厚度曲线平 移特征更为明显且曲线平移量持续增大，反映断裂F1的左行走滑特征。北侧 测线目的层原始厚度曲线特征点横坐标为2.71km，南侧测线目的层原始厚度 曲线特征点横坐标为3.79km，故由第5组测线获取的走滑断裂F1左行走滑距 为1.08km.因理想状态下走滑断裂的水平位移量应为中部最大，两端最小且 为0，故可据此原则判定由上述5组测线获取的走滑断裂F1走滑距的可靠程 度：第4组与第5组测线的目的层原始厚度曲线偏移量呈现断裂中部最大，两 端最小且为0的特征，可靠程度较高；第1组、第2组、第3组的目的层原始 厚度曲线偏移量在走滑断裂F1中央最大，但在走滑断裂两端并不能严格对应， 走滑距不为0，故可靠程度为中等。综上所述，在本具体的实施例中，走滑断 裂F1为左行走滑断裂，且左行走滑距为0.85km～1.09km，各组测线的特征点 横坐标及走滑距计算结果如表1所示。
表1各组测线的特征点横坐标及走滑距计算结果
在一个实施例中，提供了一种厘定地下走滑断裂滑动方向与走滑距离的装 置，包括以下程序单元：获取单元、编制单元和估算单元，其中：
获取单元：用于获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑 断裂的断穿层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确 定走滑断裂的延伸长度；在走滑断裂两侧等间隔选取数组平行于走滑断裂走向 的测线；确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
编制单元：用于按照同一方向获取已选取测线的目的层原始厚度，并编制 “目的层原始厚度-距离测线一侧末端的距离”折线图；
估算单元：用于基于每组平行于走滑断裂走向的测线的“目的层原始厚度 -距离测线一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特 征点位置偏移量估算走滑距离。
关于厘定地下走滑断裂滑动方向与走滑距离的装置的具体限定可以参见 上文中对于厘定地下走滑断裂滑动方向与走滑距离的方法的限定，在此不再赘 述。上述厘定地下走滑断裂滑动方向与走滑距离的装置中的各个模块可全部或 部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立 于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器 中，以便于处理器调用执行以上各个模块对应的操作。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其 内部结构图可以如图14所示。该计算机设备包括通过系统总线连接的处理器、 存储器、通信接口、显示屏和输入装置。其中，该计算机设备的处理器用于提 供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。 该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存 储介质中的操作系统和计算机程序的运行提供环境。该计算机程序被处理器执 行时以实现一种厘定地下走滑断裂滑动方向与走滑距离的方法。该计算机设备 的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可 以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球 或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图14中示出的结构，仅仅是与本申请方案相 关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限 定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些 部件，或者具有不同的部件布置。
在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器 中存储有计算机程序，涉及上述实施例方法中的全部或部分流程。
在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程 序，涉及上述实施例方法中的全部或部分流程。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程， 是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于 一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述 各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、 存储、数据库或其它介质的任何引用，均可包括非易失性和易失性存储器中的 至少一种。非易失性存储器可包括只读存储器(Read-Only Memory，ROM)、磁 带、软盘、闪存或光存储器等。易失性存储器可包括随机存取存储器(Random Access Memory，RAM)或外部高速缓冲存储器。作为说明而非局限，RAM可以是多种形式，比如静态随机存取存储器(Static Random Access Memory，SRAM) 或动态随机存取存储器(Dynamic Random Access Memory，DRAM)等。
以上实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述 实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特 征的组合不存在矛盾，都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详 细，但并不能因此而理解为对发明专利范围的限制。应当指出的是，对于本领 域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形 和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所 附权利要求为准。
Claims (10)
1.一种厘定地下走滑断裂滑动方向与走滑距离的方法，其特征在于，所述方法包括：
S1、获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂的延伸长度；
S2、在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线；
S3、确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
S4、按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原始厚度-距离测线一侧末端的距离”折线图；
S5、基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏移量估算走滑距离。
2.根据权利要求1所述的方法，其特征在于，所述获取走滑断裂的断穿层位以及走滑断裂的走向具体包括：
若走滑断裂断穿了某地层，则确定某地层为走滑断裂的断穿层位；
若沿走滑断裂的走向其方位角恒定，则认为走滑断裂走向不变；若沿走滑断裂的走向其方位角变化，则根据方位角的变化特征将走滑断裂分为数段。
3.根据权利要求1所述的方法，其特征在于，所述S2具体包括：
基于走滑断裂的走向确定选取测线的走向；
基于走滑断裂的延伸长度确定选取测线的延伸长度；其中，所选测线的长度与走滑断裂的延伸长度相等；
基于走滑断裂的延伸长度及响应用户第一输入，确定各测线间的垂直距离及测线数量；所述用户第一输入根据研究精度要求设定；具体为，当走滑断裂的延伸长度较大时，所选测线之间的垂直距离应适当增大；当走滑断裂的延伸长度较小时，所选测线之间的垂直距离应适当减小；当研究精度要求较高时，所选测线之间的垂直距离应适当增大；当研究精度要求较低时，所选测线之间的垂直距离应适当增大；
将距离走滑断裂的垂直距离相等且分别位于走滑断裂两侧的测线编为一组，并按照距离走滑断裂的垂直距离从小到大对每一组进行编号。
4.根据权利要求3所述的方法，其特征在于，所述基于走滑断裂的走向确定选取测线的走向具体为：若走滑断裂的走向恒定，则所选测线的走向恒定，且与走滑断裂相同；若走滑断裂的走向变化且具有分段特征，则所选测线的走向变化且与走滑断裂具有相同分段特征。
5.根据权利要求1所述的方法，其特征在于，所述S3具体包括：
获取走滑断裂断穿层位的地震同相轴清晰程度，响应用户的第二输入，选取目的层的顶面或底面；具体的，地震同相轴清晰程度包括物理清晰程度与地质清晰程度，所述用户第二输入为物理清晰程度与地质清晰程度都较高的断穿层位；
对选取的目的层顶面和底面开展层位追踪；
对走滑断裂两侧局部范围剥蚀量恢复。
6.根据权利要求1所述的方法，其特征在于，所述S4包括：按照同一方向获取已选取测线的残余地层厚度；其中，所述残余地层厚度为选取位置目的层的顶面深度与底面深度差的绝对值；
将残余地层厚度与选取位置的地层剥蚀量相加，按照同一方向获取已选取测线的原始地层厚度；
以获取测线地层厚度方向的起点为原点，编制“目的层原始厚度-距离测线一侧末端的距离”折线图；其中，同一组测线的“目的层原始厚度-距离测线一侧末端的距离”折线置于同一折线图中进行比较。
7.根据权利要求1所述的方法，其特征在于，所述S5包括：通过每一组中分别位于走滑断裂两侧测线的“目的层原始厚度-距离测线一侧末端的距离”折线之间的平移后是否重合性判定走滑断裂的滑动方向；其中，若位于走滑断裂某一侧测线的“目的层原始厚度-距离测线一侧末端的距离”折线经平移后可与另一侧测线的“目的层原始厚度-距离测线一侧末端的距离”折线重合，则相反方向为走滑断裂的滑动方向；
基于每一组中“目的层原始厚度-距离测线一侧末端的距离”折线的几何特征厘定特征点位置；根据每一组中“目的层原始厚度-距离测线一侧末端的距离”折线特征点的偏移量估算断裂的走滑距；其中，特征点的选取为“目的层原始厚度-距离测线一侧末端的距离”折线上的最大值点、最小值点、极大值点、极小值点中的一种或多种。
8.一种厘定地下走滑断裂滑动方向与走滑距离的装置，其特征在于，所述装置包括：
获取单元：用于获取走滑断裂的断穿层位以及走滑断裂的走向；根据走滑断裂的断穿层位以及走滑断裂的走向确定走滑断裂范围，根据走滑断裂范围确定走滑断裂的延伸长度；在走滑断裂两侧等间隔选取数组平行于走滑断裂走向的测线；确定目的层，并开展层位追踪与走滑断裂两侧局部范围剥蚀量恢复；
编制单元：用于按照同一方向获取已选取测线的目的层原始厚度，并编制“目的层原始厚度-距离测线一侧末端的距离”折线图；
估算单元：用于基于每组平行于走滑断裂走向的测线的“目的层原始厚度-距离测线一侧末端的距离”折线图判定断裂滑动方向，基于地层厚度曲线特征点位置偏移量估算走滑距离。
9.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述的方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至7中任一项所述的方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210396893.0A CN114996795B (zh) | 2022-04-15 | 2022-04-15 | 一种厘定地下走滑断裂滑动方向与走滑距离的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210396893.0A CN114996795B (zh) | 2022-04-15 | 2022-04-15 | 一种厘定地下走滑断裂滑动方向与走滑距离的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114996795A true CN114996795A (zh) | 2022-09-02 |
| CN114996795B CN114996795B (zh) | 2023-11-10 |
Family
ID=83023569
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210396893.0A Active CN114996795B (zh) | 2022-04-15 | 2022-04-15 | 一种厘定地下走滑断裂滑动方向与走滑距离的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114996795B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
| CN110275205A (zh) * | 2018-03-15 | 2019-09-24 | 中国石油化工股份有限公司 | 一种厘定盆地地下中小尺度走滑断裂活动期次的方法 |
| CN113050161A (zh) * | 2021-03-29 | 2021-06-29 | 中国矿业大学 | 一种基于断层交面线的断裂走向断距图解计算方法 |
| CN113589371A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | 碳酸盐岩内断层的封闭性评价方法 |
-
2022
- 2022-04-15 CN CN202210396893.0A patent/CN114996795B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104200039A (zh) * | 2014-09-17 | 2014-12-10 | 中国石油大学(华东) | 一种构造裂缝产状定量预测方法 |
| CN110275205A (zh) * | 2018-03-15 | 2019-09-24 | 中国石油化工股份有限公司 | 一种厘定盆地地下中小尺度走滑断裂活动期次的方法 |
| CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
| CN113589371A (zh) * | 2020-04-30 | 2021-11-02 | 中国石油化工股份有限公司 | 碳酸盐岩内断层的封闭性评价方法 |
| CN113050161A (zh) * | 2021-03-29 | 2021-06-29 | 中国矿业大学 | 一种基于断层交面线的断裂走向断距图解计算方法 |
Non-Patent Citations (2)
| Title |
|---|
| 童亨茂;宓荣三;于天才;刘宝鸿;孟令箭;杨景勇;: "渤海湾盆地辽河西部凹陷的走滑构造作用", 地质学报, no. 08 * |
| 黄超;余朝华;张桂林;傅良同;袁志云;范兴燕;: "郯庐断裂中段新生代右行走滑位移", 吉林大学学报(地球科学版), no. 03 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114996795B (zh) | 2023-11-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Agosta et al. | From fractures to flow: A field-based quantitative analysis of an outcropping carbonate reservoir | |
| US10459098B2 (en) | System and method for automatically correlating geologic tops | |
| US20180225778A1 (en) | System and method for automatically correlating geologic tops | |
| US8670944B2 (en) | Prospect stress prediction | |
| Anderson et al. | The Quaternary marine terraces of Santa Cruz, California: Evidence for coseismic uplift on two faults | |
| Castillo et al. | Trap integrity in the Lam In Aria High-Nancar Trough region, Timor Sea: prediction of fault seal failure using well-constrained stress tensors and fault surfaces interpreted from 3D seismic | |
| Asfahani | Geoelectrical investigation for characterizing the hydrogeological conditions in semi-arid region in Khanasser valley, Syria | |
| US6182013B1 (en) | Methods and apparatus for dynamically estimating the location of an oil-water interface in a petroleum reservoir | |
| Gartrell et al. | A new model for assessing trap integrity and oil preservation risks associated with postrift fault reactivation in the Timor Sea | |
| CN108106582B (zh) | 基于多期构造活动背景下走滑断裂的走滑位移计算方法 | |
| CN104597494A (zh) | 地震地层体分析方法及装置 | |
| Massiot et al. | Effects of regional and local stresses on fault slip tendency in the southern Taranaki Basin, New Zealand | |
| Wagreich et al. | Backstripping dip‐slip fault histories: apparent slip rates for the Miocene of the Vienna Basin | |
| CN114996795A (zh) | 一种厘定地下走滑断裂滑动方向与走滑距离的方法 | |
| CN110794477B (zh) | 一种井震结合识别超覆地层尖灭线的方法 | |
| CN108457646B (zh) | 确定储层流体性质的方法 | |
| CN114876454A (zh) | 一种大斜度井复杂岩性水淹层识别方法及系统 | |
| du Pisani et al. | Borehole radar delineation of the Ventersdorp Contact Reef in three dimensions | |
| CN113552623A (zh) | 曲流河沉积砂体边界的确定方法和确定装置 | |
| Karam et al. | Synthetic transfer zone characterization using seismic attributes: An example from the Parihaka fault system in the Taranaki Basin, New Zealand | |
| CN110439534A (zh) | 一种地震资料井漏预测方法 | |
| CA3000537A1 (en) | System and method for automatically correlating geologic tops | |
| Meyer et al. | Discovery of the Baldy toreva near urban areas along the southern Wasatch Range, Utah | |
| CN114647004B (zh) | 一种用于确认地下走滑断裂不同时期的滑动方向的方法 | |
| CN109709606B (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 |