CN110888163B - 微断裂发育位置的确定方法和装置 - Google Patents
微断裂发育位置的确定方法和装置 Download PDFInfo
- Publication number
- CN110888163B CN110888163B CN201811049433.0A CN201811049433A CN110888163B CN 110888163 B CN110888163 B CN 110888163B CN 201811049433 A CN201811049433 A CN 201811049433A CN 110888163 B CN110888163 B CN 110888163B
- Authority
- CN
- China
- Prior art keywords
- target
- relative height
- acquiring
- height value
- structure diagram
- 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
- 208000013201 Stress fracture Diseases 0.000 title claims abstract description 100
- 238000011161 development Methods 0.000 title claims abstract description 73
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000010586 diagram Methods 0.000 claims abstract description 146
- 238000010276 construction Methods 0.000 claims abstract description 42
- 238000009499 grossing Methods 0.000 claims abstract description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 31
- 208000010392 Bone Fractures Diseases 0.000 claims description 11
- 206010017076 Fracture Diseases 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 3
- 230000000875 corresponding effect Effects 0.000 description 54
- 238000012545 processing Methods 0.000 description 7
- 239000011435 rock Substances 0.000 description 6
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000005465 channeling Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 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
- 239000011148 porous material Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种微断裂发育位置的确定方法和装置，该方法包括：根据目标地层的多个地震剖面中每个地震剖面对应的地震成果层位解释曲线，获取目标地层的平面构造图，平面构造图用于表示目标地层中每个位置的相对高度值，相对高度值为相对于目标地层的零平面的高度值；对平面构造图进行平滑处理，获取平滑后的平面构造图；根据平面构造图和平滑后的平面构造图，获取目标地层对应的差值平面构造图；根据差值平面构造图，获取目标地层的微断裂发育位置。本发明根据目标地层的平面构造图和平滑后的平面构造图之间形成的差值平面构造图，能够有效地获取微断裂发育位置。
Description
技术领域
本发明涉及地球物理勘探开发技术领域，尤其涉及一种微断裂发育位置的确定方法和装置。
背景技术
地层断裂指的是岩层或岩体顺破裂面发生明显的位移，岩层呈现不连续、不完整状态。在地震资料中大的断裂具有明显的识别特征，如断裂在地震剖面上呈现同相轴明显错断，或者，呈现自下而上相同部位地震同相轴杂乱、连续性变差等异常特征。因为微断裂发育规模小，横向延展规律性不强，地震反射波组特征不清晰，地震剖面同相轴错断不明显，至多表现微挠曲特征，因此微断裂在地震剖面中的识别特征并不明显，特别是在地震资料品质较差的情况下，微断裂的识别就更加困难。
但由于在油气勘探开发过程中，微断裂对油气疏导和储层的改善至关重要，因此精细识别和刻画微断裂是油气勘探开发生产中的重要研究内容。而现有技术中鲜有对微断裂的精细识别的报道。
发明内容
本发明提供一种微断裂发育位置的确定方法和装置，本发明根据目标地层的平面构造图和平滑后的平面构造图之间形成的差值平面构造图，能够有效地获取微断裂发育位置。
本发明的方面提供一种微断裂发育位置的确定方法，包括：
根据目标地层的多个地震剖面中每个所述地震剖面对应的多个地震成果层位解释曲线，获取所述目标地层的平面构造图，所述平面构造图用于表示所述目标地层中每个位置的相对高度值，所述相对高度值为相对于所述目标地层的零平面的高度值；
对所述平面构造图进行平滑处理，获取平滑后的平面构造图；
根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图；
根据所述差值平面构造图，获取所述目标地层的微断裂发育位置。
可选的，所述根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图，包括：
获取每个所述地震成果层位解释曲线对应的每个剖面位置；
根据所述地震道的主测线、联络线、各所述剖面位置、以及每个所述剖面位置对应的地震成果层位解释曲线，获取所述目标地层的平面构造图。
可选的，所述对所述平面构造图进行平滑处理，获取平滑后的平面构造图，包括：
获取所述目标地层中所述每个位置的相对高度值；
将所述每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个所述相对高度值集合中包含多个位置的相对高度值；
获取每个所述相对高度值集合中的中间值；
使用所述每个位置对应的所述相对高度值集合中的中间值作为所述每个位置更新后的相对高度值，获取所述平滑后的平面构造图。
可选的，所述根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图，包括：
将所述每个位置的相对高度值、以及，所述每个位置更新后的相对高度值的差值作为所述每个位置的差值高度值，获取所述差值平面构造图。
可选的，所述根据差值平面构造图，获取所述目标地层的微断裂发育位置，包括：
根据所述差值高度值和预设高度阈值，获取所述目标地层的微断裂发育位置。
可选的，所述根据所述差值高度值和预设高度阈值，获取所述目标地层的微断裂发育位置，包括：
将大于所述预设高度阈值的所述差值高度值对应的位置确定为所述目标地层的微断裂发育位置；或者，
将大于所述预设高度阈值的所述差值高度值对应的位置确定为待选位置；
根据所述待选位置的断裂原因，获取所述目标地层的微断裂发育位置。
可选的，所述根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图之前，还包括：
获取所述目标地层。
可选的，所述获取所述目标地层，包括：
若目标区域中包含有已钻井，根据所述已钻井的测井资料，获取所述目标地层；或者，
若所述目标区域中不包含有已钻井，根据所述目标区域的相邻区域的已钻井的测井资料，获取所述目标地层；或者，
若所述目标区域中不包含有已钻井，且所述相邻区域不包含有已钻井，则根据野外露头数据，获取所述目标地层。
本发明的第二方面提供一种微断裂发育位置的确定装置，包括：
平面构造图获取模块，用于根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图，所述平面构造图用于表示所述目标地层中每个位置的相对高度值，所述相对高度值为相对于所述目标地层的零平面的高度值；
平滑后的平面构造图获取模块，用于对所述平面构造图进行平滑处理，获取平滑后的平面构造图；
差值平面构造图获取模块，用于根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图；
微断裂发育位置获取模块，用于根据差值平面构造图，获取所述目标地层的微断裂发育位置。
本发明的第三方面提供一种微断裂发育位置的确定装置，包括：至少一个处理设备和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理设备执行所述存储器存储的计算机执行指令，使得所述微断裂发育位置的确定装置执行上述微断裂发育位置的确定方法。
本发明的第四方面提供一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机执行指令，当所述计算机执行指令被处理设备执行时，实现上述微断裂发育位置的确定方法。
本发明提供一种微断裂发育位置的确定方法和装置，该方法包括：根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取目标地层的平面构造图，平面构造图用于表示目标地层中每个位置的相对高度值，相对高度值为相对于目标地层的零平面的高度值；对平面构造图进行平滑处理，获取平滑后的平面构造图；根据平面构造图和平滑后的平面构造图，获取目标地层对应的差值平面构造图；根据差值平面构造图，获取目标地层的微断裂发育位置。本发明根据目标地层的平面构造图和平滑后的平面构造图之间形成的差值平面构造图，能够有效地获取微断裂发育位置。
附图说明
图1为本发明提供的微断裂发育位置的确定方法流程示意图一；
图2为本发明提供的微断裂发育位置的确定方法流程示意图二；
图3为本发明提供的微断裂发育位置的确定装置的结构示意图一；
图4为本发明提供的微断裂发育位置的确定装置的结构示意图二；
图5为本发明提供的微断裂发育位置的确定装置的结构示意图三。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明的实施例，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1为本发明提供的微断裂发育位置的确定方法流程示意图一，图1所示方法流程的执行主体可以为微断裂发育位置的确定装置，该微断裂发育位置的确定装置可由任意的软件和/或硬件实现。如图1所示，本实施例提供的微断裂发育位置的确定方法可以包括：
S101，根据目标地层的多个地震剖面中每个地震剖面对应的地震成果层位解释曲线，获取目标地层的平面构造图，平面构造图用于表示目标地层中每个位置的相对高度值，相对高度值为相对于目标地层的零平面的高度值。
获取目标区域的地层微断裂情况对于研究目标区域的油气疏导和储层的改善至关重要，目标区域可以是勘探区域或者已钻井区域；本实施例中的目标地层为储层，储层为具有连通孔隙、允许油气在其中储存的岩层，为油气富集的岩层。
储层为处于某一深度的地层，该地层的深度可以为相对于地面具有不同深度的地层，根据目标区域的地震资料，获取目标地层时间域或者深度域的成果资料，本实施例以及下述实施例中采用的地震剖面地震成果层位解释曲线深度域的成果资料。
沿地震采集时主测线方向和联络线方向分别对目标地层进行逐道解释，通过解释，每个地震剖面上对应都有一条地震成果层位解释曲线，该地震成果层位解释曲线的起伏状态即表示该目标地层在该地震剖面上的起伏情况，多个地震剖面上的地震成果层位解释曲线组合一起就能表示目标地层的起伏状况。
将目标地层的多个地震成果层位解释曲线进行组合，绘制在一幅图中，即可得到目标地层的平面构造图，平面构造图用于表示目标地层中每个位置相对于目标地层的零相位的相对高度值，即平面构造图上的每个位置具有相对于目标地层的零平面的高低起伏状态；其中，目标地层的零平面可以预先设置，零平面即表示相对高度为0所在的“零平面”；具体的，平面构造图中有相对高度值为正值的位置，即表示该位置的高度为“零平面”之上；平面构造图中有相对高度值为负值的位置，即表示该位置的高度为“零平面”之下。
本实施例中，该平面构造图可采用不同的灰度对应不同的相对高度值，即采用灰度图的方式表示平面构造图；或者也可以采用不同等高线对应不同的相对高度值，即采用等高线的方式表示平面构造图；本实施例对平面构造图的具体形式不做限制，只要能够获取目标地层中每个位置对应的相对高度值即可。
S102，对平面构造图进行平滑处理，获取平滑后的平面构造图。
本实施例中，由于平面构造图是由多个地震成果层位解释曲线组合而成，而由于地震成果层位解释曲线呈现出的每个位置的相对高度值不足以直接获取微断裂发育位置，为了确定目标地层中的哪些位置发生微断裂，需要获取每个位置相对于零相位的挠曲程度，即相对于零相位的变形情况。
本实施例中，对平面构造图进行平滑处理，可以是采用现有技术中的如小波去噪等手段，去除地震成果层位解释曲线中的杂峰，使得地震成果层位解释曲线连续平滑，进而获取平滑后的平面构造图；此外，也可采用将相对高度值为正值的位置对应的高度设置为所有正值的相对高度值的中间值，将相对高度值为负值的位置对应的高度设置为所有负值的相对高度值的中间值，来获取平滑后的平面构造图。
S103，根据平面构造图和平滑后的平面构造图，获取目标地层对应的差值平面构造图。
由于平面构造图表示的是目标地层每个位置对应的相对高度值，而平滑后的平面构造图表示的为目标地层每个位置对应的平滑后的相对高度值，可以将二者中的各位置的相对高度值和平滑后的相对高度值做差，获取目标地层对应的差值平面构造图，该差值平面构造图用于表示目标地层的每个位置的新的相对高度值，该新的相对高度值可以反映目标地层同相轴挠曲变形的程度。
S104，根据差值平面构造图，获取目标地层的微断裂发育位置。
本实施例中，差值平面构造图用于表示目标地层的每个位置的新的相对高度值，具体的，可以用相对高度值和平滑后的相对高度值做差后获取的差值绝对值表示目标地层同相轴挠曲变形的程度，即该差值绝对值越大，表明该位置对应的同相轴挠曲变形的程度越重，微断裂发育的规模越大；该差值绝对值越小，表明该位置对应的同相轴挠曲变形的程度越轻，微断裂发育的规模越小。
本实施例中，还可以预先设置差值阈值，当差值绝对值大于该差值阈值时，确定该位置为微断裂发育位置。
本实施例提供的微断裂发育位置的确定方法包括：根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取目标地层的平面构造图，平面构造图用于表示目标地层中每个位置的相对高度值，相对高度值为相对于目标地层的零平面的高度值；对平面构造图进行平滑处理，获取平滑后的平面构造图；根据平面构造图和平滑后的平面构造图，获取目标地层对应的差值平面构造图；根据差值平面构造图，获取目标地层的微断裂发育位置。本实施例根据目标地层的平面构造图和平滑后的平面构造图，能够有效地获取微断裂发育位置。
下面结合图2对本发明提供的微断裂发育位置的确定方法进行详细说明，图2为本发明提供的微断裂发育位置的确定方法流程示意图二，如图2所示，本发明提供的微断裂发育位置的确定方法可以包括：
S201，获取目标地层。
在目标区域中获取目标地层，其中，若目标区域为已钻井区域，即目标区域中包含有已钻井，则根据已钻井的测井资料，获取目标地层；具体的，可根据测井资料中的岩性描述分析获取各地层的岩性，根据各岩性获取目标地层；或者可根据测井资料中的声波曲线或其他测井曲线，将符合储层特征的深度对应的地层确定为储层，本实施例对如何根据测井资料获取目标地层不做限制。或者，
若目标区域为勘探区域，目标区域中不包含有已钻井，则根据目标区域的相邻区域的已钻井的测井资料，获取目标地层。或者，
若目标区域中不包含有已钻井，且目标区域相邻区域不包含有已钻井，即目标区域以及与目标区域相邻的区域均为勘探区域，则根据野外露头数据，获取目标地层；其中，野外露头数据可以是与目标区域的各地层相同岩性岩石数据，根据该野外露头数据获取目标地层。
S202，获取每个地震成果层位解释曲线对应的每个剖面位置。
本实施例中，在获取地震剖面图时，便会获取剖面位置，剖面位置可以表示地震剖面所在的方向以及在目标区域中的位置坐标。本实施例中，地震剖面、地震成果层位解释曲线和剖面位置具有对应关系。获取目标地层的多个地震成果层位解释曲线对应的多个剖面位置。
S203，根据地震道的主测线、联络线、各剖面位置、以及每个剖面位置对应的地震成果层位解释曲线，获取目标地层的平面构造图。
地震道的主测线所在的方向为“北-东”方向，联络线是与主测线垂直，所在的方向为“北-西”方向。根据主测线、联络线的位置方向构建坐标系，根据各剖面位置，将各剖面位置出对应的地震成果层位解释曲线按照方向以及坐标位置进行对应，获取目标地层的平面构造图。
本实施例中可采用灰度图的方式作为实例对目标地层的平面构造图进行展示，其中，灰度相对高度值越大的地方代表相对高度值越小，灰度值越大的地方代表相对高度值越大。
S204，获取目标地层中每个位置的相对高度值；将每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个相对高度值集合中包含多个位置的相对高度值。
获取目标地层中每个位置对应的相对高度值，具体的，表现在地震成果层位解释曲线上即是获取地震成果层位解释曲线的波形上各点距离零相位的距离值；本实施例中，相对高度值有正值、有负值，将相对高度值按照大小顺序划分为多个相对高度值集合，每个相对高度值集合中包含多个相对高度值；具体的，大小顺序可按照数值的大小顺序进行排列。
示例性的，如相对高度值包含有“-10、-9.8、-9.5、-9.3、0、5.3、5.5、 6.8、7.1”，则可以将负值划分为一个集合{-10、-9.8、-9.5、-9.3}，将正值划分为一个集合{5.3、5.5、6.8、7.1}。
可以想到的是，目标地层中的相对高度值的总数量非常多，因此可预先设置相对高度值集合数量，可以将相对高度值的总数量与相对高度值集合数量的商值确定为每个相对高度值集合中包含的相对高度值集合的数量。
S205，获取每个相对高度值集合中的中间值；使用每个位置对应的相对高度值集合中的中间值作为每个位置更新后的相对高度值，获取平滑后的平面构造图。
本实施例中，在获取多个相对高度值集合后，每个相对高度值集合中包含多个相对高度值，具体的可以求取每个相对高度值集合中的中间值；当相对高度值集合中包含的相对高度值的数量为奇数个时，将大小处于中间位置的相对高度值确定为中间值；当相对高度值集合中包含的相对高度值的数量为偶数个时，将大小处于中间位置的两个相对高度值的平均值确定为中间值。
使用每个位置对应的相对高度值集合中的中间值作为每个位置更新后的相对高度值，获取平滑后的平面构造图。
本实施例中，为了提高获取的微断裂发育位置的准确性，可以对平面构造图进行多次平滑处理，具体的，平滑处理的过程可如上所示，在此不做赘述。
S206，将每个位置的相对高度值、以及，每个位置更新后的相对高度值的差值作为每个位置的差值高度值，获取差值平面构造图。
根据每个位置的相对高度值、以及，每个位置的更新后的相对高度值，即每个位置对应的中间值，获取每个位置的差值高度值，差值高度值为相对高度值和中间值的差值。
其中，平面构造图中每个位置对应的为相对高度值，平滑后的平面构造图中每个位置对应的为更新后的相对高度值，即中间值，将各位置的相对高度值和中间值做差获取各位置对应的差值高度值，其中，差值高度值为相对高度值和中间值的差值，具体的，该差值高度值可用于显示各位置的挠曲变形程度。
获取了目标层为每个位置对应的差值高度值后，将每个位置的相对高度值、以及，每个位置更新后的相对高度值的差值作为每个位置的差值高度值，获取差值平面构造图。本实施例中可以采用灰度图的方式对各位置不同的差值高度值，其中，灰度值越大，代表差值高度值越高，相应的，该位置的挠曲变形越重，微断裂发育规模越大；灰度值越小，代表差值高度值越低，相应的，该位置的挠曲变形越轻，微断裂发育规模越小。
S207，根据差值高度值和预设高度阈值，获取目标地层的微断裂发育位置。
本实施例中，可以将大于预设高度阈值的差值高度值对应的位置确定为目标地层的微断裂发育位置，即在上述实施例的基础上，可设置灰度阈值，将大于灰度阈值对应的位置确定为目标地层的微断裂发育位置。
另外，可以先将大于预设高度阈值的差值高度值对应的位置确定为待选位置；再根据待选位置的断裂原因，获取目标地层的微断裂发育位置，断裂原因可用于寻找造成待选位置处发生断裂的目标位置，具体的，可以将该目标位置作为目标地层的微断裂发育位置。
本实施例中，获取目标地层中每个位置的相对高度值；将每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个相对高度值集合中包含多个位置的相对高度值；获取每个相对高度值集合中的中间值；使用每个位置对应的相对高度值集合中的中间值作为每个位置更新后的相对高度值，获取平滑后的平面构造图，本实施例中获取平滑后的平面构造图可与原平面构造图产生差别，进而获取各位置的挠曲变形程度，可更为准确的获取目标地层的微断裂发育位置；进一步的，将每个位置的相对高度值、以及，每个位置更新后的相对高度值的差值作为每个位置的差值高度值，获取差值平面构造图，能更加直观的获取目标地层的微断裂发育位置；且本实施例中，将大于预设高度阈值的差值高度值对应的位置确定为目标地层的微断裂发育位置；或者，将大于预设高度阈值的差值高度值对应的位置确定为待选位置；根据待选位置的断裂原因，获取目标地层的微断裂发育位置，可进一步提高获取的微断裂发育位置的准确性。
图3为本发明提供的微断裂发育位置的确定装置的结构示意图一，如图 3所示，该微断裂发育位置的确定装置300包括：平面构造图获取模块301、平滑后的平面构造图获取模块302、差值平面构造图获取模块303和微断裂发育位置获取模块304。
平面构造图获取模块301，用于根据目标地层的多个地震剖面中每个地震剖面对应的果层位线解释曲线，获取目标地层的平面构造图，平面构造图用于表示目标地层中每个位置的相对高度值，相对高度值为相对于目标地层的零平面的高度值。
平滑后的平面构造图获取模块302，用于对平面构造图进行平滑处理，获取平滑后的平面构造图。
差值平面构造图获取模块303，用于根据平面构造图和平滑后的平面构造图，获取目标地层对应的差值平面构造图。
微断裂发育位置获取模块304，用于根据差值平面构造图，获取目标地层的微断裂发育位置。
本实施例提供的微断裂发育位置的确定装置与上述微断裂发育位置的确定方法实现的原理和技术效果类似，在此不作赘述。
可选的，图4为本发明提供的微断裂发育位置的确定装置的结构示意图二，如图4所示，该微断裂发育位置的确定装置300还包括：目标地层获取模块305。
目标地层获取模块305，用于获取目标地层。
可选的，目标地层获取模块305，具体用于若目标区域中包含有已钻井，根据已钻井的测井资料，获取目标地层；或者，
若目标区域中不包含有已钻井，根据目标区域的相邻区域的已钻井的测井资料，获取目标地层；或者，
若目标区域中不包含有已钻井，且相邻区域不包含有已钻井，则根据野外露头数据，获取目标地层。
可选的，平面构造图获取模块301，具体用于获取每个地震成果层位解释曲线对应的每个剖面位置；
根据地震道的主测线、联络线、各剖面位置、以及每个剖面位置对应的地震成果层位解释曲线，获取目标地层的平面构造图。
可选的，平滑后的平面构造图获取模块302，具体用于获取目标地层中每个位置的相对高度值；
将每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个相对高度值集合中包含多个位置的相对高度值；
获取每个相对高度值集合中的中间值；
使用每个位置对应的相对高度值集合中的中间值作为每个位置更新后的相对高度值，获取平滑后的平面构造图。
可选的，差值平面构造图获取模块303，具体用于将每个位置的相对高度值、以及，每个位置更新后的相对高度值的差值作为每个位置的差值高度值，获取差值平面构造图。
可选的，微断裂发育位置获取模块304，具体用于根据差值高度值和预设高度阈值，获取目标地层的微断裂发育位置。
可选的，微断裂发育位置获取模块304，具体用于将大于预设高度阈值的差值高度值对应的位置确定为目标地层的微断裂发育位置；或者，
将大于预设高度阈值的差值高度值对应的位置确定为待选位置；
根据待选位置的断裂原因，获取目标地层的微断裂发育位置。
图5为本发明提供的微断裂发育位置的确定装置的结构示意图三，该微断裂发育位置的确定装置例如可以是终端设备，比如智能手机、平板电脑、计算机等。如图5所示，该微断裂发育位置的确定装置400包括：存储器401 和至少一个处理器402。
存储器401，用于存储程序指令。
处理器402，用于在程序指令被执行时实现本实施例中的微断裂发育位置的确定方法，具体实现原理可参见上述实施例，本实施例此处不再赘述。
该微断裂发育位置的确定装置400还可以包括及输入/输出接口503。
输入/输出接口503可以包括独立的输出接口和输入接口，也可以为集成输入和输出的集成接口。其中，输出接口用于输出数据，输入接口用于获取输入的数据，上述输出的数据为上述方法实施例中输出的统称，输入的数据为上述方法实施例中输入的统称。
本发明还提供一种可读存储介质，可读存储介质中存储有执行指令，当微断裂发育位置的确定装置的至少一个处理器执行该执行指令时，当计算机执行指令被处理器执行时，实现上述实施例中的微断裂发育位置的确定方法。
本发明还提供一种程序产品，该程序产品包括执行指令，该执行指令存储在可读存储介质中。微断裂发育位置的确定装置的至少一个处理器可以从可读存储介质读取该执行指令，至少一个处理器执行该执行指令使得微断裂发育位置的确定装置实施上述的各种实施方式提供的微断裂发育位置的确定方法。
在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。
作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等) 或处理器(英文：processor)执行本发明各个实施例方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(英文：Read-Only Memory，简称：ROM)、随机存取存储器(英文：Random Access Memory，简称：RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
在上述网络设备或者终端设备的实施例中，应理解，处理器可以是中央处理单元(英文：Central Processing Unit，简称：CPU)，还可以是其他通用处理器、数字信号处理器(英文：Digital Signal Processor，简称：DSP)、专用集成电路(英文：ApplicationSpecific Integrated Circuit，简称：ASIC)等。通用处理器可以是微处理器或者该处理器也可以是任何常规的处理器等。结合本申请所公开的方法的步骤可以直接体现为硬件处理器执行完成，或者用处理器中的硬件及软件模块组合执行完成。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (8)
1.一种微断裂发育位置的确定方法，其特征在于，包括：
根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图，所述平面构造图用于表示所述目标地层中每个位置的相对高度值，所述相对高度值为相对于所述目标地层的零平面的高度值；
对所述平面构造图进行平滑处理，获取平滑后的平面构造图；
根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图；
根据所述差值平面构造图，获取所述目标地层的微断裂发育位置；
所述根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图，包括：
获取每个所述地震成果层位解释曲线对应的每个剖面位置；
根据地震道的主测线、联络线、各所述剖面位置、以及每个所述剖面位置对应的地震成果层位解释曲线，获取所述目标地层的平面构造图；
所述对所述平面构造图进行平滑处理，获取平滑后的平面构造图，包括：
获取所述目标地层中所述每个位置的相对高度值；
将所述每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个所述相对高度值集合中包含多个位置的相对高度值；
获取每个所述相对高度值集合中的中间值；
使用所述每个位置对应的所述相对高度值集合中的中间值作为所述每个位置更新后的相对高度值，获取所述平滑后的平面构造图。
2.根据权利要求1所述的方法，其特征在于，所述根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图，包括：
将所述每个位置的相对高度值、以及，所述每个位置更新后的相对高度值的差值作为所述每个位置的差值高度值，获取所述差值平面构造图。
3.根据权利要求2所述的方法，其特征在于，所述根据所述差值平面构造图，获取所述目标地层的微断裂发育位置，包括：
根据所述差值高度值和预设高度阈值，获取所述目标地层的微断裂发育位置。
4.根据权利要求3所述的方法，其特征在于，所述根据所述差值高度值和预设高度阈值，获取所述目标地层的微断裂发育位置，包括：
将大于所述预设高度阈值的所述差值高度值对应的位置确定为所述目标地层的微断裂发育位置；或者，
将大于所述预设高度阈值的所述差值高度值对应的位置确定为待选位置；
根据所述待选位置的断裂原因，获取所述目标地层的微断裂发育位置。
5.根据权利要求1所述的方法，其特征在于，所述根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图之前，还包括：
获取所述目标地层。
6.根据权利要求5所述的方法，其特征在于，所述获取所述目标地层，包括：
若目标区域中包含有已钻井，根据所述已钻井的测井资料，获取所述目标地层；或者，
若所述目标区域中不包含有已钻井，根据所述目标区域的相邻区域的已钻井的测井资料，获取所述目标地层；或者，
若所述目标区域中不包含有已钻井，且所述相邻区域不包含有已钻井，则根据野外露头数据，获取所述目标地层。
7.一种微断裂发育位置的确定装置，其特征在于，包括：
平面构造图获取模块，用于根据目标地层的多个地震剖面中每个所述地震剖面对应的地震成果层位解释曲线，获取所述目标地层的平面构造图，所述平面构造图用于表示所述目标地层中每个位置的相对高度值，所述相对高度值为相对于所述目标地层的零平面的高度值；
平滑后的平面构造图获取模块，用于对所述平面构造图进行平滑处理，获取平滑后的平面构造图；
差值平面构造图获取模块，用于根据所述平面构造图和所述平滑后的平面构造图，获取所述目标地层对应的差值平面构造图；
微断裂发育位置获取模块，用于根据差值平面构造图，获取所述目标地层的微断裂发育位置；
所述平面构造图获取模块，
具体用于获取每个所述地震成果层位解释曲线对应的每个剖面位置；
根据地震道的主测线、联络线、各所述剖面位置、以及每个所述剖面位置对应的地震成果层位解释曲线，获取所述目标地层的平面构造图；
所述平滑后的平面构造图获取模块，
具体用于获取所述目标地层中每个位置的相对高度值；
将所述每个位置的相对高度值按照大小顺序进行排序，并将排序后的相对高度值划分为多个相对高度值集合，每个所述相对高度值集合中包含多个位置的相对高度值；
获取每个所述相对高度值集合中的中间值；
使用所述每个位置对应的所述相对高度值集合中的中间值作为所述每个位置更新后的相对高度值，获取所述平滑后的平面构造图。
8.一种微断裂发育位置的确定装置，其特征在于，至少一个处理器和存储器；
所述存储器存储计算机执行指令；
所述至少一个处理器执行所述存储器存储的计算机执行指令，使得所述微断裂发育位置的确定装置执行权利要求1-6任一项所述的方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811049433.0A CN110888163B (zh) | 2018-09-10 | 2018-09-10 | 微断裂发育位置的确定方法和装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811049433.0A CN110888163B (zh) | 2018-09-10 | 2018-09-10 | 微断裂发育位置的确定方法和装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110888163A CN110888163A (zh) | 2020-03-17 |
| CN110888163B true CN110888163B (zh) | 2022-01-04 |
Family
ID=69744946
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811049433.0A Active CN110888163B (zh) | 2018-09-10 | 2018-09-10 | 微断裂发育位置的确定方法和装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110888163B (zh) |
Citations (4)
| 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 |
| CN103635831A (zh) * | 2011-02-24 | 2014-03-12 | 界标制图有限公司 | 在3d各向异性介质中的基于敏感度内核的偏移速度分析 |
| CN105717540A (zh) * | 2016-03-14 | 2016-06-29 | 中国海洋石油总公司 | 一种微幅度构造精细预测方法 |
| CN108121009A (zh) * | 2016-11-30 | 2018-06-05 | 中国石油化工股份有限公司 | 开发后期复杂断块变速构造成图方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6941228B2 (en) * | 2003-12-05 | 2005-09-06 | Schlumberger Technology Corporation | Method and system and program storage device for analyzing compressional 2D seismic data to identify zones of open natural fractures within rock formations |
-
2018
- 2018-09-10 CN CN201811049433.0A patent/CN110888163B/zh active Active
Patent Citations (4)
| 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 |
| CN103635831A (zh) * | 2011-02-24 | 2014-03-12 | 界标制图有限公司 | 在3d各向异性介质中的基于敏感度内核的偏移速度分析 |
| CN105717540A (zh) * | 2016-03-14 | 2016-06-29 | 中国海洋石油总公司 | 一种微幅度构造精细预测方法 |
| CN108121009A (zh) * | 2016-11-30 | 2018-06-05 | 中国石油化工股份有限公司 | 开发后期复杂断块变速构造成图方法 |
Non-Patent Citations (5)
| Title |
|---|
| AUTOMATIC DETECTION AND EXTRACTION OF FAULTS FROM THREE DIMENSIONAL SEISMIC DATA;Trygve Randen et al.;《SEG Technical Program Expanded Abstracts 2001》;20010101;第1-5页 * |
| Seismic attributes and integrated prediction of fractured and caved carbonate reservoirs in the Tarim Basin, China;Liu Lifeng et al.;《Pet.Sci》;20111231;第455-461页 * |
| 兴城 —丰乐地区深层火山岩储层识别;宋宗平等;《天然气工业》;20070831;第27卷(第8期);第29-30页 * |
| 塔里木盆地寒武系—奥陶系白云岩潜山勘探新发现与勘探意义;谢会文等;《中国石油勘探》;20170531;第22卷(第3期);第1-11页 * |
| 松辽盆地莺山地区营城组火山岩喷发特征研究;张鹏等;《地球物理学进展》;20161231;第31卷(第5期);第2264-2266页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110888163A (zh) | 2020-03-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111164462B (zh) | 一种人工源面波勘探方法、面波勘探装置及终端设备 | |
| CA2921340C (en) | Automated horizon auto-picking on multiple volumes | |
| EP3881105B1 (en) | Passive seismic imaging | |
| CN104950334A (zh) | 一种预测储层分布的方法及装置 | |
| CN111025384B (zh) | 基于波形分类交汇融合的储层预测方法和装置 | |
| CN103792576A (zh) | 一种基于梯度结构张量的储层非均质检测方法及设备 | |
| CN109188520B (zh) | 薄储层厚度预测方法及装置 | |
| WO2017132294A1 (en) | Feature index-based feature detection | |
| CA2876632A1 (en) | System and method for migration velocity modeling | |
| CN110675498B (zh) | 长大带状三维地质环境建模方法、装置及终端设备 | |
| CN109598068B (zh) | 古构造约束建模方法、装置和设备 | |
| CN103628866A (zh) | 获取地层的地层压力系数的方法 | |
| CN111859709B (zh) | 一种含水层结构变异转移概率的地质统计模拟方法及装置 | |
| CN110888163B (zh) | 微断裂发育位置的确定方法和装置 | |
| CN111208558A (zh) | 超深低幅度三维地质构造的建立方法及装置 | |
| CN114578445A (zh) | 基于重力资料确定断裂位置的方法及装置 | |
| CN113204855B (zh) | 地层缩短速率确定方法、系统、设备及可读介质 | |
| Parchkoohi et al. | Automatic detection of karstic sinkholes in seismic 3D images using circular Hough transform | |
| CN112782780A (zh) | 一种基于岩石物理相的储层评价方法、装置及设备 | |
| CN112987100B (zh) | 预测受层序界面控制的弱振幅岩溶层中储层的方法和装置 | |
| CN113050193B (zh) | 基底火成岩的分布识别方法、装置及系统 | |
| CN118444380A (zh) | 生物礁速度确定方法、装置、电子设备及可读介质 | |
| CN113075747A (zh) | 储层裂缝发育区域的预测方法及装置 | |
| CN112666607A (zh) | 重力反演黄土层厚度分布的方法及装置 | |
| CN114355453A (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 |