CN114063154B - 走滑断层位移量计算方法、装置、电子设备及介质 - Google Patents
走滑断层位移量计算方法、装置、电子设备及介质 Download PDFInfo
- Publication number
- CN114063154B CN114063154B CN202010746666.7A CN202010746666A CN114063154B CN 114063154 B CN114063154 B CN 114063154B CN 202010746666 A CN202010746666 A CN 202010746666A CN 114063154 B CN114063154 B CN 114063154B
- Authority
- CN
- China
- Prior art keywords
- seismic
- calculating
- point
- walk
- reference point
- 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
- 238000006073 displacement reaction Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- 230000001427 coherent effect Effects 0.000 claims abstract description 22
- 230000000149 penetrating effect Effects 0.000 claims abstract description 9
- 230000007480 spreading Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 238000000605 extraction Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 10
- 238000005755 formation reaction Methods 0.000 description 10
- 238000004062 sedimentation Methods 0.000 description 6
- 239000013049 sediment Substances 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 150000003839 salts Chemical class 0.000 description 3
- 230000009471 action Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- 239000011707 mineral 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
- 230000008569 process Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012360 testing method 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
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/64—Geostructures, e.g. in 3D data cubes
- G01V2210/642—Faults
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
技术领域
本发明涉及地球科学与地震勘探领域，更具体地，涉及一种走滑断层位移量计算方法、装置、电子设备及介质。
背景技术
走滑断层是在剪切应力作用下的一种普遍地质构造现象，所谓“走滑断层(Strike-Slip Fault)”是指两盘断块体基本上沿走向发生相对位移运动的断层，也称走向滑动断层或平移断层。走滑断层的分类方法很多，通常，根据两盘的相对运动方向可以分为左旋走滑断层和右旋走滑断层；根据倾向应力分量的特征可以分为正走滑断层和负走滑断层；根据与同一应力场中其他地层变形之间的相互关系可以分为转换断层、调解断层等。Sylvester(1988年)提议将走滑断层分为板块间和板块内两类，其中板块间的断层称为平移断层，而板块内部只限于地壳的称为横推断层。夏义平等人(2007年)提出根据构造运动规模大小，依次把走滑断层分为板块级走滑断层、盆地级走滑断层、区带级走滑断层、圈闭级走滑断层和显微级走滑断层。其中后三类中小型的走滑断层都属于盆内低序级的走滑断层。
图1示出了走滑断层地层错断的示意图。
在剪切应力作用下，走滑断层两盘原来连接在一体的地层被错断并移动，而两盘体地层和地质体以及所包含的各种地质现象却仍然留存于断层两侧，如图1所示。这些地层和地质体的所有地质现象与特征，包括构造带(构造单元、地层起伏、断层、褶皱等)、沉积相带(扇体、河道、礁体、煤层、盐底辟、泥底辟等)、地层变化(地层厚度、地层超覆、地层不整合、地层尖灭等)、成岩带(岩溶、压实带、白云岩体等)、古地貌、岩脉(火山岩体、侵入体和岩浆底辟等)、变质体和矿带体等，凡此种种所有被错开的特征都是走滑位移量估算的重要依据。
走滑断层相对位移量(滑移距)可以通过断层两盘的地质特征与识别标志的错开距离来确定，但由于难以直接准确找到两侧的地质识别标志，使得走滑断层相对位移量的估算变得十分困难。因此，长期以来走滑位移量估算一直是困扰走滑构造研究相关领域的难点问题。文献检索表明，走滑位移量的计算方法目前主要有四种方法：两盘地质参考点对比法、古地磁学研究方法、地壳变形速度估算法和沉降速率与走滑速率关系法。
两盘地质参考点对比法是沿走向对比断层两盘的断层带、化石带、岩相、岩性带或航磁异常带等特征，从而发现相似识别标志，据此推算出两盘相对位移量。这种方法是目前使用最广泛、最直接的一种估算方法。但在运用这种方法时，必须保证两盘地质体的可观测性，包括直接裸露在地表如断层带、化石带等，或者能通过地球物理方法间接反映出来，如航磁异常带等。然而，对于深埋地下的地层来说，钻井岩芯资料极其有限，也只能借助于地球物理的各类观测资料来实现。但航磁异常探测精度很低，尤其是对于深部地层更是如此。
古地磁学方法是通过测定两盘不同时期地层的古地磁数值，据此换算恢复两盘的古纬度，再由古纬度差值对比推算出该时期内断层两盘的相对移动方向和距离。这种方法比较适用于近南北向、走滑断距较大的断层。应用这种方法时，必须保证取样、测试的精度较高而且各取样岩石有比较准确的年代数据。
地壳变形速度估算法认为走滑断层两盘地壳变形的速度差与时间的累积就是断层两盘的相对位移量的大小。显然，这种方法不仅须要全面弄清走滑断层的构造演化过程，而且还须要准确确定各构造期次中两盘的地壳变形速度的大小。
沉降速率与走滑速率关系法是通过理论模型或实验模拟研究，推导出走滑拉分盆地中盆地的走滑速率与沉降速率之间的关系，即走滑速率同盆地的几何形状参数、最大沉降深度和沉降速率存在一定的数值关系。由此可在恢复走滑拉分盆地的沉降史的前提下估算出走滑盆地的累计走滑位移量。显然，这种方法对模型和实验模拟都具有很高的要求。
上述四种方法中，后两种方法主要是针对盆地级以上的巨型走滑断裂来说的，而前两种方法适用于不同规模的走滑断层，其中“两盘地质参考点对比法”应用最为普遍，但其参考特征标志往往难以获取。总体上，现有的技术方法都具有其局限性和繁苛的前提条件，实用性较差。
因此，有必要开发一种盆内低序级走滑断层位移量计算方法、装置、电子设备及介质。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种走滑断层位移量计算方法、装置、电子设备及介质，其能够通过地震反射波组特征与多种地震属性特征，采用地震道的相关性计算与分析，定量计算出盆内低序级走滑断层两盘地层的相对位移量，为研究走滑断层提供了科学依据。
第一方面，本公开实施例提供了一种走滑断层位移量计算方法，包括：
通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；
沿所述地震反射层位提取地震特征；
根据所述地震特征，确定走滑位移量计算的参考点；
根据所述参考点确定对应点，进而计算位移量。
优选地，通过三维叠偏数据体计算地震相干数据体包括：
分别计算所述三维叠偏数据体在纵线方向与横线方向的互相关系数；
计算所述三维叠偏数据体在每个时间点的地震相干值，获得所述地震相干数据体。
优选地，通过公式(1)计算所述三维叠偏数据体在每个时间点的地震相干值：
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值。
优选地，沿所述地震反射层位提取地震特征包括：
沿所述地震反射层位从所述地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；
沿所述地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
优选地，根据所述参考点，确定对应点包括：
根据所述参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；
若无法直接确定所述对应点，则通过地震道相关系数曲线图确定所述对应点。
优选地，通过地震道相关系数曲线图确定所述对应点包括：
根据所述参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；
根据所述地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为所述参考点的对应点。
优选地，通过所述参考点与所述对应点，通过坐标距离计算公式计算所述位移量。
作为本公开实施例的一种具体实现方式，
第二方面，本公开实施例还提供了一种走滑断层位移量计算装置，包括：
数据输入模块，通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；
地震特征提取模块，沿所述地震反射层位提取地震特征；
参考点确定模块，根据所述地震特征，确定走滑位移量计算的参考点；
计算模块，根据所述参考点确定对应点，进而计算位移量。
优选地，通过三维叠偏数据体计算地震相干数据体包括：
分别计算所述三维叠偏数据体在纵线方向与横线方向的互相关系数；
计算所述三维叠偏数据体在每个时间点的地震相干值，获得所述地震相干数据体。
优选地，通过公式(1)计算所述三维叠偏数据体在每个时间点的地震相干值：
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值。
优选地，沿所述地震反射层位提取地震特征包括：
沿所述地震反射层位从所述地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；
沿所述地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
优选地，根据所述参考点，确定对应点包括：
根据所述参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；
若无法直接确定所述对应点，则通过地震道相关系数曲线图确定所述对应点。
优选地，通过地震道相关系数曲线图确定所述对应点包括：
根据所述参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；
根据所述地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为所述参考点的对应点。
优选地，通过所述参考点与所述对应点，通过坐标距离计算公式计算所述位移量。
第三方面，本公开实施例还提供了一种电子设备，该电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的走滑断层位移量计算方法。
第四方面，本公开实施例还提供了一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的走滑断层位移量计算方法。
本发明的方法和装置具有其它的特性和优点，这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的，或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述，这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述，本发明的上述以及其它目的、特征和优势将变得更加明显，其中，在本发明示例性实施例中，相同的参考标号通常代表相同部件。
图1示出了走滑断层地层错断的示意图。
图2示出了根据本发明的一个实施例的利用三维地震数据实现走滑断层位移量估算的示意图。
图3示出了根据本发明的一个实施例的走滑断层位移量计算方法的步骤的流程图。
图4示出了根据本发明的一个实施例的三维地震数据沿层相干切片的示意图。
图5a与图5b分别示出了根据本发明的一个实施例的采用2个不同的参考点计算的地震道最大相关系数结果的示意图。
图6示出了根据本发明的一个实施例的计算地震数据中实际走滑断层的位移量结果的示意图。
图7示出了根据本发明的一个实施例的一种走滑断层位移量计算装置的框图。
附图标记说明：
201、数据输入模块；202、地震特征提取模块；203、参考点确定模块；204、计算模块。
具体实施方式
下面将更详细地描述本发明的优选实施方式。虽然以下描述了本发明的优选实施方式，然而应该理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。
走滑断层两盘相对位移量可以通过断层两盘地层的地质特征与识别标志的错开距离来估算。但其实现的关键是在两盘都能够准确找到可供参考的地质特征标识，而现有的技术方法都很难满足这个需求，尤其对于深埋地下的地层来说。
被走滑断层错开之前相邻的地下地层都是在同一个环境下、同样的地质条件下所形成的地质体，因此具有相同的地质特征；而被走滑断层错断位移之后分属两盘的地质体依然保留着原来“同根同源”的地质特征。而地震反射响应特征是对地下各种地质特征的间接反映，在地震可分辨的前提下，不同的地质特征往往都能找到相应的地震反射响应特征，因此，被走滑断层所错开的具有相同地质特征的地层，相应地，也会在对应地层的地震反射中搜寻到相似的特征标识。
图2示出了根据本发明的一个实施例的利用三维地震数据实现走滑断层位移量估算的示意图。
本发明就是采用技术方法，在三维地震数据中确定走滑断层位置，并根据走滑断层任一盘中选定的地震参考点及其特征从另一盘的地震数据中搜寻与之最相似点的位置，这个最相似点就是与参考点同处一套地层的相邻点被错开后的位置，依据参考点与原相邻点的位置座标，就能计算出两盘的相对走滑位移量。如图2所示，它截取了三维地震数据体的一个局部及其侧切面，顶面两个垂直方向分别代表三维地震数据的线(Inline)和道(Crossline)方向，顶面每个黑圆点都对应于一道地震数据，而其中加粗的横线表示地震相干识别出的走滑断层，其侧面就是地震反射道，侧面加粗的2个十分相似的地震反射波组的位置就对应于原来同处一套地层的相邻2个点被走滑断层两盘错开后各自所处的位置点(星形标注点为参考点，带圆环星形标注点为原来与参考点相邻的对应点)。
本发明提供一种走滑断层位移量计算方法，包括：
通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；在一个示例中，通过三维叠偏数据体计算地震相干数据体包括：分别计算三维叠偏数据体在纵线方向与横线方向的互相关系数；计算三维叠偏数据体在每个时间点的地震相干值，获得地震相干数据体。
在一个示例中，通过公式(1)计算三维叠偏数据体在每个时间点的地震相干值：
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值。
具体地，三维地震叠偏数据体纵线方向(对应X轴)上任意两点A(xi,yi)和B(xi+1,yi)在某时间t的地震道s，在时移量为n时的互相关系数Cx为：
同样，三维地震叠偏数据体横线方向(对应Y轴)上任意两点A(xi,yi)和C(xi,yi+1)在某时间t的地震道s，在时移量为m时的互相关系数Cy为：
上述两式中，相关时窗长度为2ω+1，时移量n与m分别近似为纵横线上地层时间视倾角。
再利用公式(1)对上述纵横向相邻地震道相关值取几何平均值Cxy，就得到三维地震叠偏数据体A点在t时间的地震相干值，如此依次逐点求取整个三维数据体各个时间点的相干值，就能得到三维地震相干数据体。
沿地震反射层位提取地震特征；在一个示例中，沿地震反射层位提取地震特征包括：沿地震反射层位从地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；沿地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
根据地震特征，确定走滑位移量计算的参考点。
具体地，根据地震相干切片得到的走滑断裂平面展布特征并综合地层沉积-地震相平面特征来确定参考点，参考点一般选定在走滑断裂附近一侧，越靠近断裂越好，但注意不要选在断裂破碎带上。而且，参考点最好选在地层沉积-地震相平面图上具有典型沉积现象(如河道、火山、盐丘等)的位置处，这样便于后续快速、准确地在断裂另外一侧找到对应点。
根据参考点确定对应点，进而计算位移量。在一个示例中，根据参考点，确定对应点包括：根据参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；若无法直接确定对应点，则通过地震道相关系数曲线图确定对应点。
在一个示例中，通过地震道相关系数曲线图确定对应点包括：根据参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；根据地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为参考点的对应点。
在一个示例中，通过参考点与对应点，通过坐标距离计算公式计算位移量。
具体地，根据走滑断裂和地层地震沉积分布特征，在参考点所在走滑断裂的另外一盘，寻找沉积地震相特征最相近的对应点；如果能够准确找到参考点的对应点，那么直接可以根据两点的位置，用两点坐标距离计算公式得出该走滑断层对所选地层上造成的位移量；如果无法直接找到参考点的对应点，那么就提取参考点位置处的地震波组特征，沿参考点的另一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，得到一条地震道相关系数曲线图，找到最大相关系数值对应的地震道位置点，记为参考点的对应点，相关系数一般应该大于60％以上。同样，由两点坐标距离公式计算得到该走滑断层对所选地层上造成的位移量。
本发明还提供一种走滑断层位移量计算装置，包括：
数据输入模块，通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位。在一个示例中，通过三维叠偏数据体计算地震相干数据体包括：分别计算三维叠偏数据体在纵线方向与横线方向的互相关系数；计算三维叠偏数据体在每个时间点的地震相干值，获得地震相干数据体。
在一个示例中，通过公式(1)计算三维叠偏数据体在每个时间点的地震相干值：
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值。
具体地，三维地震叠偏数据体纵线方向(对应X轴)上任意两点A(xi,yi)和B(xi+1,yi)在某时间t的地震道s，在时移量为n时的互相关系数Cx为公式(2)。
同样，三维地震叠偏数据体横线方向(对应Y轴)上任意两点A(xi,yi)和C(xi,yi+1)在某时间t的地震道s，在时移量为m时的互相关系数Cy为公式(3)。
上述两式中，相关时窗长度为2ω+1，时移量n与m分别近似为纵横线上地层时间视倾角。
再利用公式(1)对上述纵横向相邻地震道相关值取几何平均值Cxy，就得到三维地震叠偏数据体A点在t时间的地震相干值，如此依次逐点求取整个三维数据体各个时间点的相干值，就能得到三维地震相干数据体。
地震特征提取模块，沿地震反射层位提取地震特征；在一个示例中，沿地震反射层位提取地震特征包括：沿地震反射层位从地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；沿地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
参考点确定模块，根据地震特征，确定走滑位移量计算的参考点。
具体地，根据地震相干切片得到的走滑断裂平面展布特征并综合地层沉积-地震相平面特征来确定参考点，参考点一般选定在走滑断裂附近一侧，越靠近断裂越好，但注意不要选在断裂破碎带上。而且，参考点最好选在地层沉积-地震相平面图上具有典型沉积现象(如河道、火山、盐丘等)的位置处，这样便于后续快速、准确地在断裂另外一侧找到对应点。
计算模块，根据参考点确定对应点，进而计算位移量。在一个示例中，根据参考点，确定对应点包括：根据参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；若无法直接确定对应点，则通过地震道相关系数曲线图确定对应点。
在一个示例中，通过地震道相关系数曲线图确定对应点包括：根据参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；根据地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为参考点的对应点。
在一个示例中，通过参考点与对应点，通过坐标距离计算公式计算位移量。
具体地，根据走滑断裂和地层地震沉积分布特征，在参考点所在走滑断裂的另外一盘，寻找沉积地震相特征最相近的对应点；如果能够准确找到参考点的对应点，那么直接可以根据两点的位置，用两点坐标距离计算公式得出该走滑断层对所选地层上造成的位移量；如果无法直接找到参考点的对应点，那么就提取参考点位置处的地震波组特征，沿参考点的另一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，得到一条地震道相关系数曲线图，找到最大相关系数值对应的地震道位置点，记为参考点的对应点，相关系数一般应该大于60％以上。同样，由两点坐标距离公式计算得到该走滑断层对所选地层上造成的位移量。
本发明还提供一种电子设备，电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的走滑断层位移量计算方法。
本发明还提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现上述的走滑断层位移量计算方法。
为便于理解本发明实施例的方案及其效果，以下给出四个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
实施例1
图3示出了根据本发明的一个实施例的走滑断层位移量计算方法的步骤的流程图。
如图3所示，该走滑断层位移量计算方法包括：步骤101，通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；步骤102，沿地震反射层位提取地震特征；步骤103，根据地震特征，确定走滑位移量计算的参考点；步骤104，根据参考点确定对应点，进而计算位移量。
图4示出了根据本发明的一个实施例的三维地震数据沿层相干切片的示意图。
图5a与图5b分别示出了根据本发明的一个实施例的采用2个不同的参考点计算的地震道最大相关系数结果的示意图。
图4为三维地震数据体沿某层(如H1)地震反射层位的相干切片，图中右侧竖直的黑线为人为设定的一条走滑断层，沿该走滑断层把三维地震数据左右两侧的地震道认为错位移动了50道。图5a与图5b是应用本发明方法，选定2个不同的参考点分别计算出的地震道相关系数曲线，不难看出，根据2条相关系数曲线找到的最大相关系数值的对应点所计算出的位移量均为1250米，50个道，地震道间隔为25米。这与预先人为设定的结果完全吻合，从而，验证了本发明方法估算结果的可靠性与准确性。
图6示出了根据本发明的一个实施例的计算地震数据中实际走滑断层的位移量结果的示意图。
图4中左侧斜黑线为地震数据相干体切片得到的实际走滑断层，计算结果表明该走滑断层两盘相对位移了4个地震道，地震道间隔为25米，因此，估算该走滑断层位移量为100米，如图6所示。本发明所估算的位移量与该走滑断层的级别及该区地质构造活动强度相符。
实施例2
图7示出了根据本发明的一个实施例的一种走滑断层位移量计算装置的框图。
如图7所示，该走滑断层位移量计算装置，包括：
数据输入模块201，通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；
地震特征提取模块202，沿地震反射层位提取地震特征；
参考点确定模块203，根据地震特征，确定走滑位移量计算的参考点；
计算模块204，根据参考点确定对应点，进而计算位移量。
作为可选方案，通过三维叠偏数据体计算地震相干数据体包括：
分别计算三维叠偏数据体在纵线方向与横线方向的互相关系数；
计算三维叠偏数据体在每个时间点的地震相干值，获得地震相干数据体。
作为可选方案，通过公式(1)计算三维叠偏数据体在每个时间点的地震相干值：
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值。
作为可选方案，沿地震反射层位提取地震特征包括：
沿地震反射层位从地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；
沿地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
作为可选方案，根据参考点，确定对应点包括：
根据参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；
若无法直接确定对应点，则通过地震道相关系数曲线图确定对应点。
作为可选方案，通过地震道相关系数曲线图确定对应点包括：
根据参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；
根据地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为参考点的对应点。
作为可选方案，通过参考点与对应点，通过坐标距离计算公式计算位移量。
实施例3
本公开提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述走滑断层位移量计算方法。
根据本公开实施例的电子设备包括存储器和处理器。
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(RAM)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(ROM)、硬盘、闪存等。
该处理器可以是中央处理单元(CPU)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
实施例4
本公开实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的走滑断层位移量计算方法。
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
上述计算机可读存储介质包括但不限于：光存储介质(例如：CD－ROM和DVD)、磁光存储介质(例如：MO)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置ROM的媒体(例如：ROM盒)。
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (7)
1.一种走滑断层位移量计算方法，其特征在于，包括：
通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；
沿所述地震反射层位提取地震特征；
根据所述地震特征，确定走滑位移量计算的参考点；
根据所述参考点确定对应点，进而计算位移量；
其中，通过三维叠偏数据体计算地震相干数据体包括：
分别计算所述三维叠偏数据体在纵线方向与横线方向的互相关系数；
计算所述三维叠偏数据体在每个时间点的地震相干值，获得所述地震相干数据体；
其中，通过公式（1）计算所述三维叠偏数据体在每个时间点的地震相干值：
（1）
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值；
其中，根据所述参考点，确定对应点包括：
根据所述参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；
若无法直接确定所述对应点，则通过地震道相关系数曲线图确定所述对应点。
2.根据权利要求1所述的走滑断层位移量计算方法，其中，沿所述地震反射层位提取地震特征包括：
沿所述地震反射层位从所述地震相干数据体中提取沿层相干切片，获得走滑断裂的平面展布情况；
沿所述地震反射层位提取沿层地层沉积地震相切片，获得地层沉积-地震相特征的平面展布情况。
3.根据权利要求1所述的走滑断层位移量计算方法，其中，通过地震道相关系数曲线图确定所述对应点包括：
根据所述参考点所在走滑断裂的另外一盘，从地震叠偏数据体中逐点计算与参考点位置地震道反射波组的地震相关系数，获得地震道相关系数曲线图；
根据所述地震道相关系数曲线图，确定最大相关系数值对应的地震道位置为所述参考点的对应点。
4.根据权利要求1或3所述的走滑断层位移量计算方法，其中，通过所述参考点与所述对应点，通过坐标距离计算公式计算所述位移量。
5.一种走滑断层位移量计算装置，其特征在于，包括：
数据输入模块，通过三维叠偏数据体计算地震相干数据体，输入走滑断层断穿地层的地震反射解释层位；
地震特征提取模块，沿所述地震反射层位提取地震特征；
参考点确定模块，根据所述地震特征，确定走滑位移量计算的参考点；
计算模块，根据所述参考点确定对应点，进而计算位移量；
其中，通过三维叠偏数据体计算地震相干数据体包括：
分别计算所述三维叠偏数据体在纵线方向与横线方向的互相关系数；
计算所述三维叠偏数据体在每个时间点的地震相干值，获得所述地震相干数据体；
其中，通过公式（1）计算所述三维叠偏数据体在每个时间点的地震相干值：
（1）
其中，Cxy为三维叠偏数据体在时间点t的地震相干值，maxCx(t,n,xi,yi)为纵线方向搜索时窗内时移量n时的最大相关值，maxCy(t,m,xi,yi)为横线方向搜索时窗内时移量m时的最大相关值；
其中，根据所述参考点，确定对应点包括：
根据所述参考点所在走滑断裂的另外一盘，直接寻找沉积地震相特征最相近的对应点；
若无法直接确定所述对应点，则通过地震道相关系数曲线图确定所述对应点。
6.一种电子设备，其特征在于，所述电子设备包括：
存储器，存储有可执行指令；
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现权利要求1-4中任一项所述的走滑断层位移量计算方法。
7.一种计算机可读存储介质，其特征在于，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现权利要求1-4中任一项所述的走滑断层位移量计算方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010746666.7A CN114063154B (zh) | 2020-07-29 | 2020-07-29 | 走滑断层位移量计算方法、装置、电子设备及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010746666.7A CN114063154B (zh) | 2020-07-29 | 2020-07-29 | 走滑断层位移量计算方法、装置、电子设备及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114063154A CN114063154A (zh) | 2022-02-18 |
| CN114063154B true CN114063154B (zh) | 2024-05-07 |
Family
ID=80226879
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010746666.7A Active CN114063154B (zh) | 2020-07-29 | 2020-07-29 | 走滑断层位移量计算方法、装置、电子设备及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114063154B (zh) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5537320A (en) * | 1994-10-31 | 1996-07-16 | Landmark Graphics Corporation | Method and apparatus for identifying fault curves in seismic data |
| US5892732A (en) * | 1996-04-12 | 1999-04-06 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
| CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
| CN109856678A (zh) * | 2019-03-07 | 2019-06-07 | 中国石油化工股份有限公司 | 一种压扭性走滑断层走滑速率的计算方法 |
| RU2698551C1 (ru) * | 2018-04-12 | 2019-08-28 | Федеральное государственное бюджетное учреждение науки Институт земной коры Сибирского отделения Российской академии наук (ИЗК СО РАН) | Способ определения параметров структуры разломной трещиноватости литосферы |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7069149B2 (en) * | 2001-12-14 | 2006-06-27 | Chevron U.S.A. Inc. | Process for interpreting faults from a fault-enhanced 3-dimensional seismic attribute volume |
| RU2594377C1 (ru) * | 2013-05-15 | 2016-08-20 | Лэндмарк Графикс Корпорейшн | Стратиграфия сейсмических трасс для интерпретации геологических разломов |
| FR3043481A1 (zh) * | 2015-11-05 | 2017-05-12 | Services Petroliers Schlumberger |
-
2020
- 2020-07-29 CN CN202010746666.7A patent/CN114063154B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5537320A (en) * | 1994-10-31 | 1996-07-16 | Landmark Graphics Corporation | Method and apparatus for identifying fault curves in seismic data |
| US5892732A (en) * | 1996-04-12 | 1999-04-06 | Amoco Corporation | Method and apparatus for seismic signal processing and exploration |
| CN101452081A (zh) * | 2007-12-05 | 2009-06-10 | 中国科学院地质与地球物理研究所 | 消除地震多次波的方法 |
| CN102033242A (zh) * | 2010-10-22 | 2011-04-27 | 中国石油化工股份有限公司 | 一种深层倾斜裂缝储层地震振幅预测方法 |
| RU2698551C1 (ru) * | 2018-04-12 | 2019-08-28 | Федеральное государственное бюджетное учреждение науки Институт земной коры Сибирского отделения Российской академии наук (ИЗК СО РАН) | Способ определения параметров структуры разломной трещиноватости литосферы |
| CN109856678A (zh) * | 2019-03-07 | 2019-06-07 | 中国石油化工股份有限公司 | 一种压扭性走滑断层走滑速率的计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114063154A (zh) | 2022-02-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104977618B (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| US6498989B1 (en) | Method for predicting dynamic parameters of fluids in a subterranean reservoir | |
| US6721661B2 (en) | Method of distinguishing types of geologic sedimentation | |
| CN109425896A (zh) | 白云岩油气储层分布预测方法及装置 | |
| WO2015042103A1 (en) | Identifying geological formation depth structure using well log data | |
| CN104853822A (zh) | 一种评价页岩气储层及寻找甜点区的方法 | |
| Caers et al. | Stochastic integration of seismic data and geologic scenarios: A West Africa submarine channel saga | |
| CN104636980B (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
| CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
| Aarre et al. | Seismic detection of subtle faults and fractures | |
| Li et al. | Part II: A demonstration of integrating multiple-scale 3D modelling into GIS-based prospectivity analysis: A case study of the Huayuan-Malichang district, China | |
| Li | Geophysical exploration technology: Applications in lithological and stratigraphic reservoirs | |
| Harris | The role of geology in reservoir simulation studies | |
| Chen et al. | Interpretation of fault system in the Tana Sag, Kenya, using edge recognition techniques and Euler deconvolution | |
| Aboud | Determination of sedimentary cover and structural trends in the Central Sinai area using gravity and magnetic data analysis | |
| CN115327663A (zh) | 深部矿产资源勘查用空-地-井立体地球物理探测方法 | |
| Elebiju et al. | Investigation of links between Precambrian basement structure and Paleozoic strata in the Fort Worth Basin, Texas, USA, using high-resolution aeromagnetic (HRAM) data and seismic attributes | |
| Abdel-Fattah et al. | Diacritical seismic signatures for complex geological structures: case studies from Shushan Basin (Egypt) and Arkoma Basin (USA) | |
| CN105628576B (zh) | 灰质白云岩孔隙度测量方法及装置 | |
| Vasco | Seismic imaging of reservoir flow properties: Time-lapse pressure changes | |
| CN104820237B (zh) | 计算地层的方差体的方法 | |
| Ibanez et al. | An integrated remote sensing and biomarker maturity parameter-based approach to trace petroleum migration in a complex fault zone lineament of the Brazilian Equatorial Margin | |
| CN114063154B (zh) | 走滑断层位移量计算方法、装置、电子设备及介质 | |
| Elmasry et al. | Exploration of geothermal resources utilizing geophysical and borehole data in the Abu Gharadig basin of Egypt’s Northern Western desert | |
| Djunin et al. | Hydrogeodynamics of oil and gas basins |
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 |