CN110441818A - 单河道砂体的识别方法、装置和计算机存储介质 - Google Patents

单河道砂体的识别方法、装置和计算机存储介质 Download PDF

Info

Publication number
CN110441818A
CN110441818A CN201910703175.1A CN201910703175A CN110441818A CN 110441818 A CN110441818 A CN 110441818A CN 201910703175 A CN201910703175 A CN 201910703175A CN 110441818 A CN110441818 A CN 110441818A
Authority
CN
China
Prior art keywords
channel sand
formation
river channel
target locations
sampled 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.)
Granted
Application number
CN201910703175.1A
Other languages
English (en)
Other versions
CN110441818B (zh
Inventor
周连敏
宋舜尧
吴佳朋
王少波
黄国梁
王晶晶
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Natural Gas Co Ltd
Original Assignee
China Petroleum and Natural Gas Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Natural Gas Co Ltd filed Critical China Petroleum and Natural Gas Co Ltd
Priority to CN201910703175.1A priority Critical patent/CN110441818B/zh
Publication of CN110441818A publication Critical patent/CN110441818A/zh
Application granted granted Critical
Publication of CN110441818B publication Critical patent/CN110441818B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/62Physical property of subsurface
    • G01V2210/624Reservoir 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是本发明实施例提供的一种目标地层的结构示意图;
图6是本发明实施例提供的另一种目标地层的结构示意图;
图7是本发明实施例提供的一种采样区域的示意图;
图8是本发明实施例提供的一种复合河道砂体的结构示意图;
图9是本发明实施例提供的一种单河道砂体的识别装置的示意图;
图10是本发明实施例提供的另一种单河道砂体的识别装置的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种复合河道砂体内单河道砂体的接触位置图。如图1所示,该复合河道砂体中包括两条单河道砂体(C1、C2),两条单河道砂体拼合叠置,且两条单河道砂体接触位置为单河道砂体的边界。图2是本发明实施例提供的一种复合河道砂体内地震波波形变化示意图。如图2所示,两条单河道砂体中除接触位置地震波波形的变化是连续的,而在两条单河道砂体接触的边界位置,地震波波形出现了不连续的反射特征,因而在单河道砂体的边界处的地震属性参数的相似程度会与单河道砂体其他位置的地震属性参数存在差异。即通过判断地震属性参数的相似程度就可以确定出单河道砂体的边界。
图3是本发明实施例提供的一种单河道砂体的识别方法的流程图。如图3所示,该识别方法采用上位机执行,包括:
步骤101:获取目标地层区域的地震数据。
其中,目标地层区域为复合河道砂体所在区域。
步骤102:根据目标地层区域中采样点的地震数据,计算出采样点的至少两种地震属性参数的相干属性值。
其中,相干属性值反映采样点的地震属性参数与周围的采样点的地震属性参数的相似程度。
步骤103:将目标地层区域划分成多层切片区域。
步骤104:根据相干属性值确定每个切片区域中单河道砂体的边界线。
步骤105:将各层切片区域中单河道砂体的边界线相连,得到目标地层区域中单河道砂体的边界线。
本发明实施例中获取目标地层区域的地震数据,并根据目标地层区域中采样点的地震数据,计算出采样点的至少两种地震属性参数的相干属性值,同时将目标地层区域划分成多层切片区域,根据相干属性值确定每个切片区域中单河道砂体的边界线,将各层切片区域中单河道砂体的边界线相连,得到目标地层区域中单河道砂体的边界线。由于复合河道砂体内两单河道砂体接触的边界位置的地震波波形存在不连续反射特征,因而在单河道砂体的边界处的地震属性参数的相似程度会与单河道砂体其他位置的地震属性参数存在差异。因此,本发明通过采样点的地震属性参数的相干属性值就可以确定出单河道砂体的边界,同时采用切片区域的形式逐层确定目标地层区域中各个位置的边界线,也可以有效提高确定边界的准确性,并且采样点的相干属性值采用至少两个地震属性参数计算确定,因此能避免地震数据的准确性较差时仅依靠单一因素识别处的边界结果准确性较差的问题。
图4是本发明实施例提供的另一种单河道砂体的识别方法的流程图。如图4所示,该识别方法采用上位机执行,包括:
步骤201:获取目标地层区域的地震数据。
其中,目标地层区域为复合河道砂体所在区域。
获取地震数据的区域范围通常采用深度时窗范围作为计量单位。本实施例中,获取目标地层区域的地震数据为选取合适的深度时窗范围,该深度时窗范围能包含整个复合河道砂体所在的区域。以使得可以从深度时窗范围内包含的地震数据中提取出复合河道砂体的地震属性参数信息。
步骤201中,目标地层区域可以采用至少两种方式确定:
第一种:以同一深度时窗为起点,固定深度时窗长度的等时扫描确定目标地层区域。
示例性地,如图5所示,以某一深度时窗起点,且在以该深度时窗为起点向下扫描时,能扫描到目的地层G的层顶A和层底B。按照一定的深度时窗长度扫描,直至深度时窗范围满足包括整个目的地层G即可。该方法适用于目的地层G构造起伏较小的情况。
第二种:以目的地层的层顶或层底为起点,目的地层的深度时窗范围为长度确定目标地层区域。
示例性地,如图6所示,以目的地层G的层顶A或层底B为时间界面起点,按照目的地层的深度时窗范围为长度从层顶A向层底B或从层底B向层顶A扫描,直至深度时窗范围满足包括整个目的地层G即可。该方法适用于目的地层G构造起伏较大的情况。
可选地,在确定目的地层区域的深度时窗范围后,并获取该目的地层区域的地震数据后,还需要对该目的地层区域的地震数据进行过滤,保留复合河道砂体区域的地震数据,以提高识别的准确性。根据振幅高值区为河道砂岩发育区域的原理,判断砂泥岩的振幅值分界线,通过井上的河道砂体范围与地震振幅属性体相匹配,确定河道砂体的振幅值范围,过滤掉地震振幅数据体的低振幅。例如,河道砂体的振幅属性值为大于8000,则对低于8000的振幅值进行过滤处理。
可选地,根据复合河道砂体的深度时窗范围对深度视窗范围微调。例如,复合河道砂体的起始时窗深度为1800ms,终止时窗为1805ms,则可稍微放大一下深度时窗的范围,分别向上和向下放大1ms深度时窗,确保该复合河道砂体能显示完全。
步骤202:在目标地层区域的横切平面内,确定采样点对应的采样区域。
其中,横切平面是垂直于河道砂体延伸方向的平面(参见图7)。
如图7所示,采样区域(参见图7中虚线所示)包括:位于采样区域中心的基准采样点P和位于基准采样点周围的区域采样点p。
可选地,在确定采样区域时,根据当前选取的采样点,以该采样点为基准在垂直方向上向上和向下各增加一个深度时窗长度(例如1ms),在横切方向确保该基准采样点的两侧各包括2至4个区域采样点。设置采样范围可以精确分析地震数据的变化,克服数据太多,对计算结果形成干扰的问题。
步骤203:根据采样区域内的采样点的地震属性参数,计算基准采样点的不同地震属性参数对应的相干属性值。
其中,相干属性值反映采样点的地震属性参数与周围的采样点的地震属性参数的相似程度。
在确定各采样点的采样区域后,由于复合河道砂体内部两单河道砂体接触部位地震波波形的不连续反射特征,可以根据地震属性参数的相似程度确定单河道砂体的边界。其中各采样点的地震属性参数可以直接从地震数据中提取。
可选地,计算基准采样点的不同地震属性参数对应的相干属性值可以采用公式(1)确定。
其中,Xj(n)为地震属性参数,j=1、2…M,M为地震数据横切方向采样区域内的总道数,n=1、2…N,N为采样区域内的采样点。R为相干属性值,R的大小反映了地震属性参数的相似程度,R值大,相似性好。反之,则相似性差。
计算各个地震属性参数的相干属性值时,通过依次计算采样区域内各采样点的地震属性参数,最终完成相干属性值的计算。
步骤204:将目标地层区域划分成多层切片区域。
划分目标地层区域时可以采用以下两种划分方式:
第一种,将目标地层区域划分为多层与地层水平面平行的切片区域。也即是,各切片区域均与横切方向垂直。
可选地,切片区域的厚度可以跟目标地层区域的深度时窗长度确定。
示例性地,目标地层区域的深度时窗长度为10ms时,可以以1ms为间隔提取切片区域。
第二种,将目标地层区域划分为多层与地层水平面垂直的切片区域。也即是,各切片区域均与横切方向平行。
步骤205:将不同地震属性参数对应的相干属性值相加,得到复合相干属性值。
本实施例中,地震属性参数可以包括:地震道倾角和地震道振幅。从地震数据中提取出地震道倾角和地震道振幅后,可以根据公式(1)分别计算地震道倾角的相干属性值和地震道振幅的相干属性值,将地震道倾角的相干属性值和地震道振幅的相干属性值相加,计算得到复合相干属性值。
步骤206:在切片区域的横切平面内,将复合相干属性值最小的采样点确定为单河道砂体的边界点。
由于复合河道砂体内两单河道砂体接触的边界位置的地震波波形存在不连续反射特征,因而在单河道砂体的边界处的地震属性参数的相干属性值大于单河道砂体的其他位置的地震属性参数的相干属性值。而相干属性值R值越大,表明采样范围内地震属性参数的相似性好,相干属性值R值越小,表明采样范围内地震属性参数的相似性差。因此,在切片区域的横切平面内,复合相干属性值最小的采样点即可确定为单河道砂体的边界点。
步骤207:将切片区域的不同横切平面内的边界点相连,得到切片区域中单河道砂体的边界线。
如图8所示,将切片区域中不同横切平面内的边界点相连得到的连线即为切片区域中单河道砂体的边界线L1。
步骤208:将各层切片区域中单河道砂体的边界线相连,得到目标地层区域中单河道砂体的边界线。
如图8所示,在横切平面内,将各层切片区域中单河道砂体的边界线相连得到的连线即为横切方向上单河道砂体的边界线L2。边界线L1和L2共同构成目标地层区域中单河道砂体的边界线。
图9是本发明实施例提供的一种单河道砂体的识别装置的示意图。如图9所示,该识别装置包括:获取模块100、计算模块200、划分模块300、确定模块400。获取模块100用于获取目标地层区域的地震数据,目标地层区域为复合河道砂体所在区域。计算模块200用于根据目标地层区域中采样点的地震数据,计算出采样点的至少两种地震属性参数的相干属性值,相干属性值反映采样点的地震属性参数与周围的采样点的地震属性参数的相似程度。划分模块300用于将目标地层区域划分成多层切片区域。确定模块400用于根据相干属性值确定每个切片区域中单河道砂体的边界线,确定模块400用于将各层切片区域中单河道砂体的边界线相连,得到目标地层区域中单河道砂体的边界线。
可选地,确定模块400还用于在目标地层区域的横切平面内,确定采样点对应的采样区域。其中,采样区域包括:位于采样区域中心的基准采样点和位于基准采样点周围的区域采样点。计算模块200还用于根据采样区域内的采样点的地震属性参数,计算基准采样点的不同地震属性参数对应的相干属性值。
可选地,确定模块400用于将不同地震属性参数对应的相干属性值相加,得到复合相干属性值。确定模块400用于在切片区域的横切平面内,将复合相干属性值最小的采样点确定为单河道砂体的边界点。确定模块400用于将切片区域的不同横切平面内的边界点相连,得到切片区域中单河道砂体的边界线。
可选地,划分模块300用于将目标地层区域划分为多层与地层水平面平行的切片区域。或者,划分模块300用于将目标地层区域划分为多层与地层水平面垂直的切片区域。
可选地,地震属性参数包括:地震道倾角和地震道振幅。
图10是本发明实施例提供的另一种单河道砂体的识别装置的示意图。如图10所示,该单河道砂体的识别装置700可以是电脑等。
通常,单河道砂体的识别装置700包括有:处理器701和存储器702。
处理器701可以包括一个或多个处理核心,比如4核心处理器、8核心处理器等。处理器701可以采用DSP(Digital Signal Processing,数字信号处理)、FPGA(Field-Programmable Gate Array,现场可编程门阵列)、PLA(Programmable Logic Array,可编程逻辑阵列)中的至少一种硬件形式来实现。处理器701也可以包括主处理器和协处理器,主处理器是用于对在唤醒状态下的数据进行处理的处理器,也称CPU(Central ProcessingUnit,中央处理器);协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。在一些实施例中,处理器701可以在集成有GPU(Graphics Processing Unit,图像处理器),GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中,处理器701还可以包括AI(Artificial Intelligence,人工智能)处理器,该AI处理器用于处理有关机器学习的计算操作。
存储器702可以包括一个或多个计算机可读存储介质,该计算机可读存储介质可以是非暂态的。存储器702还可包括高速随机存取存储器,以及非易失性存储器,比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中,存储器702中的非暂态的计算机可读存储介质用于存储至少一个指令,该至少一个指令用于被处理器701所执行以实现本申请中方法实施例提供的单河道砂体的识别方法。
在一些实施例中,单河道砂体的识别装置700还可选包括有:外围设备接口703和至少一个外围设备。处理器701、存储器702和外围设备接口703之间可以通过总线或信号线相连。各个外围设备可以通过总线、信号线或电路板与外围设备接口703相连。具体地,外围设备包括:射频电路704、触摸显示屏705、摄像头706、音频电路707、定位组件708和电源709中的至少一种。
外围设备接口703可被用于将I/O(Input/Output,输入/输出)相关的至少一个外围设备连接到处理器701和存储器702。在一些实施例中,处理器701、存储器702和外围设备接口703被集成在同一芯片或电路板上;在一些其他实施例中,处理器701、存储器702和外围设备接口703中的任意一个或两个可以在单独的芯片或电路板上实现,本实施例对此不加以限定。
显示屏705用于显示UI(User Interface,用户界面)。该UI可以包括图形、文本、图标、视频及其它们的任意组合。当显示屏705是触摸显示屏时,显示屏705还具有采集在显示屏705的表面或表面上方的触摸信号的能力。该触摸信号可以作为控制信号输入至处理器701进行处理。此时,显示屏705还可以用于提供虚拟按钮和/或虚拟键盘,也称软按钮和/或软键盘。在一些实施例中,显示屏705可以为一个,设置单河道砂体的识别装置700的前面板;在另一些实施例中,显示屏705可以为至少两个,分别设置在单河道砂体的识别装置700的不同表面或呈折叠设计;在再一些实施例中,显示屏705可以是柔性显示屏,设置在单河道砂体的识别装置700的弯曲表面上或折叠面上。甚至,显示屏705还可以设置成非矩形的不规则图形,也即异形屏。显示屏705可以采用LCD(Liquid Crystal Display,液晶显示屏)、OLED(Organic Light-Emitting Diode,有机发光二极管)等材质制备。
电源709用于为单河道砂体的识别装置700中的各个组件进行供电。电源709可以是交流电、直流电、一次性电池或可充电电池。当电源709包括可充电电池时,该可充电电池可以支持有线充电或无线充电。该可充电电池还可以用于支持快充技术。
本领域技术人员可以理解,图10中示出的结构并不构成对单河道砂体的识别装置700的限定,可以包括比图示更多或更少的组件,或者组合某些组件,或者采用不同的组件布置。
本发明实施例还提供了一种计算机存储介质,当所述计算机存储介质中的指令由单河道砂体的识别装置的处理器执行时,使得单河道砂体的识别装置能够执行图3或图4所示实施例提供的单河道砂体的识别方法。
一种包含指令的计算机程序产品,当其在计算机上运行时,使得计算机执行上述图3或图4所示实施例提供的单河道砂体的识别方法。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (11)

1.一种单河道砂体的识别方法,其特征在于,所述识别方法包括:
获取目标地层区域的地震数据,所述目标地层区域为复合河道砂体所在区域;
根据所述目标地层区域中采样点的地震数据,计算出所述采样点的至少两种地震属性参数的相干属性值,所述相干属性值用于反映所述采样点的地震属性参数与周围的所述采样点的地震属性参数的相似程度;
将所述目标地层区域划分成多层切片区域;
根据所述相干属性值确定每个所述切片区域中单河道砂体的边界线;
将各层所述切片区域中单河道砂体的边界线相连,得到所述目标地层区域中单河道砂体的边界线。
2.根据权利要求1所述的单河道砂体的识别方法,其特征在于,所述根据所述目标地层区域中采样点的地震数据,计算出所述采样点的至少两种地震属性参数的相干属性值,包括:
在所述目标地层区域的横切平面内,确定所述采样点对应的采样区域,所述采样区域包括:位于所述采样区域中心的基准采样点和位于所述基准采样点周围的区域采样点;
根据所述采样区域内的采样点的地震属性参数,计算所述基准采样点的不同地震属性参数对应的所述相干属性值。
3.根据权利要求1所述的单河道砂体的识别方法,其特征在于,所述根据所述相干属性值确定每个所述切片区域中单河道砂体的边界线,包括:
将所述采样点的至少两种地震属性参数分别对应的所述相干属性值相加,得到所述采样点的复合相干属性值;
在所述切片区域的横切平面内,将所述复合相干属性值最小的采样点确定为单河道砂体的边界点;
将所述切片区域的不同横切平面内的所述边界点相连,得到所述切片区域中单河道砂体的边界线。
4.根据权利要求1至3任一项所述的单河道砂体的识别方法,其特征在于,所述将所述目标地层区域划分成多层切片区域,包括:
将所述目标地层区域划分为多层与地层水平面平行的切片区域;或者,
将所述目标地层区域划分为多层与地层水平面垂直的切片区域。
5.根据权利要求1至3任一项所述的单河道砂体的识别方法,其特征在于,所述地震属性参数包括:地震道倾角和地震道振幅。
6.一种单河道砂体的识别装置,其特征在于,所述识别装置包括:
获取模块,用于获取目标地层区域的地震数据,所述目标地层区域为复合河道砂体所在区域;
计算模块,用于根据所述目标地层区域中采样点的地震数据,计算出所述采样点的至少两种地震属性参数的相干属性值,所述相干属性值用于反映所述采样点的地震属性参数与周围的所述采样点的地震属性参数的相似程度;
划分模块,用于将所述目标地层区域划分成多层切片区域;
确定模块,用于根据所述相干属性值确定每个所述切片区域中单河道砂体的边界线;以及用于将各层所述切片区域中单河道砂体的边界线相连,得到所述目标地层区域中单河道砂体的边界线。
7.根据权利要求6所述的单河道砂体的识别装置,其特征在于,所述确定模块:用于在所述目标地层区域的横切平面内,确定所述采样点对应的采样区域,所述采样区域包括:位于所述采样区域中心的基准采样点和位于所述基准采样点周围的区域采样点;
所述计算模块:用于根据所述采样区域内的采样点的地震属性参数,计算所述基准采样点的不同地震属性参数对应的所述相干属性值。
8.根据权利要求6所述的单河道砂体的识别装置,其特征在于,所述确定模块:用于将不同地震属性参数对应的所述相干属性值相加,得到复合相干属性值;用于在所述切片区域的横切平面内,将所述复合相干属性值最小的采样点确定为单河道砂体的边界点;用于将所述切片区域的不同横切平面内的所述边界点相连,得到所述切片区域中单河道砂体的边界线。
9.根据权利要求6至8任一项所述的单河道砂体的识别装置,其特征在于,所述划分模块:用于将所述目标地层区域划分为多层与地层水平面平行的切片区域;或者,将所述目标地层区域划分为多层与地层水平面垂直的切片区域。
10.一种单河道砂体的识别装置,其特征在于,所述识别装置包括:处理器;被配置为存储处理器可执行指令的存储器;其中,所述处理器被配置为执行权利要求1至5任一项所述的识别方法。
11.一种计算机存储介质,其上存储有计算机指令,其特征在于,所述计算机指令被处理器执行时实现权利要求1至5任一项所述的识别方法。
CN201910703175.1A 2019-07-31 2019-07-31 单河道砂体的识别方法、装置和计算机存储介质 Active CN110441818B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910703175.1A CN110441818B (zh) 2019-07-31 2019-07-31 单河道砂体的识别方法、装置和计算机存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910703175.1A CN110441818B (zh) 2019-07-31 2019-07-31 单河道砂体的识别方法、装置和计算机存储介质

Publications (2)

Publication Number Publication Date
CN110441818A true CN110441818A (zh) 2019-11-12
CN110441818B CN110441818B (zh) 2021-03-30

Family

ID=68432646

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910703175.1A Active CN110441818B (zh) 2019-07-31 2019-07-31 单河道砂体的识别方法、装置和计算机存储介质

Country Status (1)

Country Link
CN (1) CN110441818B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103941288A (zh) * 2014-04-16 2014-07-23 中国海洋石油总公司 一种基于平面导航的三维目标地质体解释方法
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
CN105093303A (zh) * 2015-08-11 2015-11-25 中国海洋石油总公司 一种河流相砂体不连续界线精细表征的地震解释方法
CN105372703A (zh) * 2015-08-11 2016-03-02 中国石油化工股份有限公司 一种河道砂体精细识别方法
US9354338B1 (en) * 2012-02-22 2016-05-31 Westerngeco L.L.C. Generating synthetic seismic traces
CN109116414A (zh) * 2018-08-29 2019-01-01 中国海洋石油集团有限公司 一种基于谱分解和相干属性融合流程的水道复合体沉积要素精细识别方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9354338B1 (en) * 2012-02-22 2016-05-31 Westerngeco L.L.C. Generating synthetic seismic traces
CN103941288A (zh) * 2014-04-16 2014-07-23 中国海洋石油总公司 一种基于平面导航的三维目标地质体解释方法
CN104502969A (zh) * 2014-12-30 2015-04-08 中国石油化工股份有限公司 河道砂岩性油藏的识别方法
CN105093303A (zh) * 2015-08-11 2015-11-25 中国海洋石油总公司 一种河流相砂体不连续界线精细表征的地震解释方法
CN105372703A (zh) * 2015-08-11 2016-03-02 中国石油化工股份有限公司 一种河道砂体精细识别方法
CN109116414A (zh) * 2018-08-29 2019-01-01 中国海洋石油集团有限公司 一种基于谱分解和相干属性融合流程的水道复合体沉积要素精细识别方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
CUI XINKAI ET AL: "Study on the Identification of Composite Sand Body Single Channel by Well Seismic Combination", 《IOP CONFERENCE SERIES: EARTH AND ENVIRONMENTAL SCIENCE》 *
刘金连等: "济阳探区单一河道砂体边界地质建模及其地震正演响应特征分析", 《石油物探》 *
周连敏等: "复合曲流河道内的单河道识别", 《石油地球物理勘探》 *
徐立恒等: "井震反演识别复合砂体内单一河道", 《西南石油大学学报(自然科学版)》 *
王永刚: "《地震资料综合解释方法》", 31 May 2007, 中国石油大学出版社 *
石巨业等: "利用地层切片研究河道砂体发育特征—以南图尔盖盆地 Doshan地区为例", 《海相油气地质》 *
肖佃师等: "井震联合识别复合砂体中单一河道—以朝44区块扶余油层为例", 《石油地球物理勘探》 *

Also Published As

Publication number Publication date
CN110441818B (zh) 2021-03-30

Similar Documents

Publication Publication Date Title
CN109976930A (zh) 异常数据的检测方法、系统及存储介质
CN106537283A (zh) 智能装置唤醒
CN110400304A (zh) 基于深度学习的物体检测方法、装置、设备及存储介质
EP2756335B1 (en) Method and system of subsurface horizon assignment
CN106353818B (zh) 一种地震数据时差校正的方法及系统
CN106338763A (zh) 一种确定地震剖面显示数据的方法及装置
US20210173117A1 (en) Characterization of subsurface regions using moving-window based analysis of unsegmented continuous data
CN108375785A (zh) 裂缝带位置校正方法及装置
CN109696706A (zh) 盐构造检测方法、装置及存储介质
CN108897036A (zh) 一种地震数据处理方法及装置
CN110441818A (zh) 单河道砂体的识别方法、装置和计算机存储介质
WO2016115437A1 (en) Approaches to tying well logs to seismic data
CN108693557A (zh) 一种地震数据处理方法及装置
CN108052677A (zh) 页面处理方法、装置及可读存储介质
FR3061236A1 (fr) Segmentation des signatures de frequence temporelle pour la determination automatisee de defaut de tuyau
US9869799B2 (en) Object-based well correlation
CN115929292A (zh) 基于成像测井和常规测井融合的测井方法、装置及设备
CN106094026B (zh) 一种获取垂直地震数据走廊叠加剖面的方法及装置
CN106504107B (zh) 一种确定油层射孔层位的方法及装置
CN111426813B (zh) 砂体连通质量的确定方法、装置及存储介质
CN106054249B (zh) 储层评价模型建立方法及装置、储层评价方法及装置
CN108181656B (zh) 一种近偏移距转换波动校正方法及装置
CN107976712B (zh) 地质界面的确定方法和装置
CN111736215A (zh) 断层断距确定方法及装置
CN106353810B (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