CN108646292B - 裂缝密度预测方法、装置及计算机存储介质 - Google Patents

裂缝密度预测方法、装置及计算机存储介质 Download PDF

Info

Publication number
CN108646292B
CN108646292B CN201810170528.1A CN201810170528A CN108646292B CN 108646292 B CN108646292 B CN 108646292B CN 201810170528 A CN201810170528 A CN 201810170528A CN 108646292 B CN108646292 B CN 108646292B
Authority
CN
China
Prior art keywords
fracture
stratum
data volume
gather data
target area
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
Application number
CN201810170528.1A
Other languages
English (en)
Other versions
CN108646292A (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 National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201810170528.1A priority Critical patent/CN108646292B/zh
Publication of CN108646292A publication Critical patent/CN108646292A/zh
Application granted granted Critical
Publication of CN108646292B publication Critical patent/CN108646292B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • 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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/16—Survey configurations
    • G01V2210/167—Very long offset

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

裂缝密度预测方法、装置及计算机存储介质
技术领域
本说明书涉及石油及天然气勘探、开发技术领域，特别是裂缝密度预测方法、装置及计算机存储介质。
背景技术
在石油和天然气的勘探、开发行业中，把保存油气的圈闭称为油藏。塔里木盆地碳酸盐岩是塔里木油田油气勘探开发的重要领域，而塔里木碳酸盐岩储层非均质性强，储集层特征主要表现为三大类：裂缝-洞穴型、裂缝-孔洞型和裂缝型储层，进入油田开发阶段，裂缝-孔洞型储层是最主要的研究目标，优选高产稳产目标即是要找到连通性好且有一定储集规模的缝洞体系。实际生产中，只有发育达到一定强度的裂缝时才能有效沟通周边孔隙，且裂缝强度越强，其能沟通周边孔隙范围就越大，所形成的有效储集空间规模则越大，因此，井震联合有效裂缝密度预测在碳酸盐岩油气藏开发中尤为重要。
目前针对深埋缝洞型碳酸盐岩连通性分析，仅仅是进行层面分析或者简单地提取一个层段范围内属性进行分析，前人采用的常规叠后与常规叠前裂缝预测方法的研究精度远远不能满足生产的需要。
最常用的裂缝预测方法主要有叠前预测方法和叠后预测方法，叠后预测方法主要是相干技术和曲率技术等，该方法对断裂的解释有一定指导作用，但对微裂缝发育方向和密度的预测精度为足；叠前预测方法主要是各向异性裂缝预测方法，是基于椭圆拟合算法进行裂缝预测，其认为椭圆长短轴之间的关系反映裂缝强度，椭圆长短轴方向指示裂缝发育方向。该方法具有较大有局限性，即其只能对发育单组方向裂缝或发育方向相对一致的裂缝带进行有效预测，对以发育多组方向裂缝则不适用，比如，当发育有四组方向裂缝且各组裂缝密度相近，每组裂缝方向之间的方位角度数为45度，其椭圆拟合结果近似为圆，指示出各向同性特征，没有裂缝密度一说，而实际上其裂缝密度远远大于发育单一方向裂缝时的裂缝密度。另外，椭圆长短轴之间的直接关系反映的只是地震数据中的各向异性强弱程度，并不是裂缝密度，需要根据裂缝发育情况将裂缝各向异性强度的纯物理意义有效的转换为有地质意义裂缝密度。
综上所述，如何提出一种有效预测裂缝密度的有效性以及其分布特征，从而指导裂缝洞体系空间连通性分析是本领域亟需解决的问题。
发明内容
本说明书实施方式提供裂缝密度预测方法、装置及计算机存储介质，通过对已有数据进行分析处理，聚类地震反射特征，根据目标区域地震道记录中的偏移距信息和方位角信息，对特定地层的裂缝密度进行预测，确定所述特定地层裂缝密度的有效性，提高缝洞体系空间连通性分析的精度，为缝洞型油气勘探开发提供重要依据。
本说明书实施方式提供一种裂缝密度预测方法，包括，提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
本说明书实施方式提供一种裂缝密度预测装置，包括：道集数据体生成模块：用于提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；区域划分模块：用于根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；裂缝地质数据获取模块：用于根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；裂缝密度预测模块：用于根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
本说明书实施方式提供一种计算机存储介质，所述存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
由以上本说明书实施方式提供的技术方案可见，本说明书实施方式通过提取目标区域中偏移距和方位角相同或趋于相同的地震道记录生成道集数据体。根据所述道集数据体对所述目标区域进行分区，根据所述道集数据体和所述分区的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度，在处理过程中考虑到裂缝发育各向异性的特征，保留地震道记录的偏移距和方位角信息，提高了目标区域裂缝密度预测的精度。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本说明书实施方式提供的裂缝密度预测方法的执行流程图；
图2为本说明书实施方式提供的道集数据体示意图；
图3为本说明书实施方式提供的OVT道集规则化处理的前后对比图；
图4为本说明书实施方式提供的单组裂缝发育区裂缝各向异性特征示意图；
图5为本说明书实施方式提供的两组裂缝发育区裂缝各向异性特征示意图；
图6为本说明书实施方式提供的多组裂缝发育区裂缝各向异性特征示意图；
图7为本说明书实施方式提供的根据炮检距和方向角选择样点的示意图；
图8为本说明书实施方式提供的椭圆拟合后椭圆离心率与裂缝密度关系示意图；
图9为本说明书实施方式提供的两组不同方向裂缝发育区样点分布特征示意图；
图10为本说明书实施方式提供的蜘蛛网图与椭圆拟合效果对比图；
图11为本说明书实施方式提供的测井曲线井眼环境校正与VSP曲线校正前后对比图；
图12为本说明书实施方式提供的曲线校正前后井震标定结果对比示意图；
图13为本说明书实施方式提供的单组裂缝发育区井震标定后获取有效裂缝密度示意图；
图14为本说明书实施方式提供的两组裂缝发育区井震标定后获取有效裂缝密度示意图；
图15为本说明书实施方式提供的多组裂缝发育区井震标定后获取有效裂缝密度示意图；
图16为本说明书实施方式提供的有效裂缝密度分布示意图；
图17为本说明书实施方式提供的裂缝密度预测结果示意图；
图18为本说明书实施方式提供的地震剖面示意图；
图19为本说明书实施方式提供的生产曲线示意图
图20为本说明书实施方式提供的不同裂缝密度正演示意图；
图21为本说明书实施方式提供的裂缝密度预测装置的结构示意图。
具体实施方式
下面将结合本说明书实施例中的附图，对本说明书实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本说明书一部分实施例，而不是全部的实施例。基于本说明书中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书保护的范围。
在本说明书提供的一个场景示例中，使用本说明书提供的方法某山区进行裂缝密度预测。
在本场景示例中，根据所获取的目标区域所对应的观测系统参数，可以包括炮线距、接收线距、面元、最小覆盖次数以及偏移距等，获取炮线和检波线相互垂直的所有地震道集的集合，即十字排列子集；然后对所述十字排列子集中每一个十字排列按炮线距和检波线距进行等距离划分，得到许多小矩形，每个小矩形为OVT(偏移距向量片技术)向量片，从所有的OVT向量片中获取偏移距和方位角大致相同的OVT向量片，所获取的OVT向量片构成了OVT道集数据体。
在本场景示例中，请参与图2，在获得所述OVT向量片后，可以对所述OVT向量片进行OVT域数据规则化处理，首先对所述OVT向量片进行面元中心化和数据插值处理，然后对所述OVT向量片进行方位角规则化处理，最后进行OVT域偏移处理，为了保留其方位信息，进行全部OVT向量片偏移处理，直到完成所有的OVT向量片偏移处理，获得如图2所述的道集地震数据。
在本场景示例中，请参阅图3，所述道集数据体样点分布离散，同时对所述样点进行沿偏移距增大方向和掩方位角度增大方向进行分布位置的规则化处理。
在本场景示例中，请参阅图20，可以假设地层中裂缝为高角度或垂直裂缝，基本裂缝模型为硬币型模型，当发育单组方向裂缝时其裂缝密度为ρf，当发育正交两组方向裂缝时其裂缝密度为2ρf，当发育N组方向间隔相同的裂缝时其裂缝密度为Nρf(一般N≤4)。由于有效地震波为反射波，其他条件一定的情况下，可以认为裂缝密度越大则波阻抗差异越大，所以参与分析的道集数据体样点振幅越大则相应的裂缝密度越大。
在本场景示例中，请参阅图7，根据山区目标区域的所述道集数据体的方位角和偏移距选择在相同方位角和偏移距范围内单次覆盖的样点进行椭圆拟合，确定目标区域的单组裂缝发育区。
在本场景示例中，请参阅图9，针对两组裂缝，选择各方位覆盖次数均匀的样点时，样点并不是椭圆形态。请参阅图10，进行蜘蛛图分析，在给定提取800-5000m范围的偏移距范围内的数据时，按照10度方位角间隔(如10度，20度，30度等)，可以直接将各个相邻的样点直接连线而成一个闭合的不规则的图形，进而可以确定所述目标区域两组裂缝发育区。在数据点信息可信的基础上，数据点离原点距离最大的方位可能是裂缝发育方向。
在本场景示例中，可以将所述目标区域除了单组裂缝发育区和两组裂缝发育区以外的区域确定为多组裂缝发育区。
在本场景示例中，单组裂缝发育区，裂缝各向异性特征如图4所示；两组不同方向裂缝发育区，裂缝各向异性特征如图5所示；多组不同方向裂缝发育区，如图6所示，各裂缝之间方向夹角相等，表现各向同性特征，只是在反射强度上有所差异。
在本场景示例中，对所述已知地层的OVT道集数据体进行椭圆拟合，根据椭圆离心率获取所述已知地层的各向异性数据。建立获取的已知地层的裂缝密度和各向异性数据的关系，如图8所述所述裂缝密度和椭圆离心率(各向异性)之间满足的关系为y＝0.0362e2.6346x，其中，y标识裂缝密度，x表示椭圆离心率(各向异性)。根据待测地层的OVT道集数据体的方位角和偏移距，舍去偏移距较远和较近的样点，选择各方位覆盖次数均匀的样点进行椭圆拟合，计算得到椭圆的离心率得到所述待测区域的各向异性数据，根据y＝0.0362e2.6346x计算得到对应的裂缝密度。
在本场景示例中，在其他条件一定的情况下，裂缝密度越大则波阻抗差异越大，所以样点的能量值越大所对应的裂缝密度越大。根据测井资料和钻井资料，获取已知地层的裂缝密度，计算所述已知地层道集数据体的地震波反射能量，建立获取的已知地层的裂缝密度和所述地震波反射能量之间的关系。根据待测地层蜘蛛图分析结果得到待测地层道集数据体的地震波反射能量，根据所述待测地层地震波反射能量以及已经建立的已知地层的裂缝密度和所述振幅值之间的关系，计算得到待测地层的裂缝密度。
在本场景示例中，对多组不同方向裂缝发育区，根据待测地层的OVT道集数据体中的方位角数据和偏移距数据，舍去偏移距较远和较近的样点，选择各方位覆盖次数均匀的样点，去除道集样点中的奇异值点，根据相干、曲率等算法得到待测地层的裂缝密度。
在本场景示例中，请参阅图11，对目标区域的测井曲线进行环境校正处理等工作，对所有测井资料进行检查并对不同年份不同测井探管系列的数据进行标准化处理。应用VSP资料对声波曲线进行校正，从而完成对测井资料和地震信息的井震标定，如图12所示，为井震标定后的结果对比图。如图13至15所示，将所述单组方向裂缝发育区、两组方向裂缝发育区、多组方向裂缝发育区求得的裂缝密度与测井曲线进行匹配，得到各分区有效裂缝所对应的阈值，进而可以得到如图16所示的有效裂缝密度分布图。
在本场景示例中，请参阅图17至图19，应用本说明书实施方式提供的方法，裂缝连通性预测吻合率提高了83％。
请参阅图1，本说明书实施方式提供裂缝密度预测方法。所述裂缝密度预测方法具体可以包括以下步骤。
步骤S10：提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体。
在本实施方式中，所述符合指定条件的地震道记录可以包括，所述目标区域中，炮检距和方位角相同，或者地震道记录之间所述炮检距和方位角的差别在一定范围内。具体地，例如，设定方位角差别在5度以内，偏移距差别在10m以内，则可以将方位角45度，偏移距50m的地震道记录和方位角42度，偏移距52m的地震道记录提取出来生成道集数据体
在本实施方式中，所述道集数据体可以包括根据目标区域中偏移距和方位角符合指定条件的地震道记录而生成的道集记录中所包括的地震数据。具体地，可以包括炮检距、振幅、频率、方位角、时间、深度、东西南北坐标等地震数据。
在本实施方式中，所述目标区域可以是指整个勘探区域，也可以包括整个勘探区域中的部分区域或开发区域。所述目标区域可以包括已经勘探的已知地层，也可以包括待勘探的待测地层。待测地层也可以称为目的层。
在本实施方式中，所述道集数据体可以包括目标区域已知地层的道集数据体也可以包括待测地层的道集数据体。
在本实施方式中，所述偏移距可以包括在所述目标区域的观测系统中，炮点距离检波点之间的距离。
在本实施方式中，所述方位角可以包括在所述目标区域的观测系统中，在炮点与检波点位置关系与正北方向的夹角。可以根据不同方位角地震波的反射振幅差异获得较为精细的各向异性地震属性。在裂缝和断层的区域，不同方位角地震波的反射振幅可以是有差异的，而空间各点与周围点的反射振幅差异可以是裂缝引起的，也可能是岩性横向变化带来的。也就是说可以利用空间各点与周围点的反射振幅差异进行裂缝预测有更强的多解性，但是利用不同方位角地震波的反射振幅差异可以有利于消除这种多解性。
在本实施方式中，所述提取目标区域中偏移距和方位角相同或趋于相同的地震道记录可以包括根据所获取的目标区域所对应的观测系统参数，可以包括炮线距、接收线距、面元、最小覆盖次数以及偏移距等，获取炮线和检波线相互垂直的所有地震道集的集合，即十字排列子集；然后对所述十字排列子集中每一个十字排列按炮线距和检波线距进行等距离划分，得到许多小矩形，每个小矩形为OVT(偏移距向量片技术)向量片，从所有的OVT向量片中获取偏移距和方位角大致相同的OVT向量片，所获取的OVT向量片构成了OVT道集。根据偏移距和方位角，可以将目标区域地震道记录进行细分，编排成不同的OVT号，抽取不同CMP(共中心点道集)的具有相同OVT号的地震道就可以形成新的OVT道集。具体地，例如，将目标区域观测系统划为十字排子集，有OVT号为1的地震道记录，其偏移距范围大约为370m到410m，方位角大约为20度～40度。将目标区域内所有CMP位于1点位置的地震道记录抽取出来，重新排序，可以形成一个OVT道集，那么OVT片号为1的OVT道集内所有地震道偏移距都约为370m～410m，方位角约为20度～40度。每个OVT道集从每个CMP中抽取了一道，因此每个OVT道集又是一个单次覆盖的非零偏数据体，每道有相应的Inline和Xline线号，对应特定的物理位置。
在本实施方式中，可以对提取目标区域中偏移距和方位角相同或趋于相同的地震道记录进行规则化处理。所述规则化处理可以包括面元中心化处理、数值插值处理、方位角规则化处理等。规则化处理后可以对所述地震道记录进行偏移处理以获得所述道集数据体的道集数据。所述道集数据可以包括能量、振幅等反映所述道集数据体性质的数据。
步骤S12：根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区。
在本实施方式中，所述目标区域地层中的裂缝发育情况可以不同，可以包括发育单组裂缝，区域中仅包含发育方向一致的裂缝；可以包括发育两组裂缝，区域中包含发育方向不一致的两组裂缝；可以包括发育多组裂缝，区域中包含发育方向互不相同的多组裂缝。裂缝组数不同，裂缝密度预测方法也可以不同，具体地，例如，对于单组裂缝区域，可以进行椭圆拟合，根据椭圆离心率可以得到对应的裂缝密度，但对于多组裂缝区域，其椭圆拟合结果近似于圆，无法预测其裂缝密度。可以根据裂缝的发育情况对所述目标区域进行分区。
在本实施方式中，所述分区的类型可以不止一个，所述各个分区内裂缝的各向异性特征不同。例如，可以将1组裂缝的区域分成一组，2组裂缝的区域分成一组，3组裂缝的区域分成一组等。不同的裂缝组合，所述裂缝的各向异性特征也有所不同。另外，由于如果是沿方位角均匀分布的裂缝组合的地震响应可能就是一个各向同性的结果，所以裂缝组数小于等于4组。
在本实施方式中，可以根据裂缝的发育情况对所述目标区域进行分区。已知，对于单组发育裂缝可以根据所述道集数据体中的炮检距信息和方位角信息对其进行椭圆拟合，可以对目标区域内的道集数据体进行椭圆拟合，将拟合结果符合椭圆的道集数据体所在的区域可以确定为单组裂缝发育区。
在本实施方式中，可以将目标区域中确定的单组裂缝发育区去除，对剩余区域根据所述道集数据体中的炮检距信息和方位角信息对其进行进行蜘蛛图分析，在可以将相邻的样点直接连线而成一个闭合的不规则的图形道集数据体所在的区域可以确定为两组裂缝发育区。
在本实施方式中，可以将目标区域中除了单组裂缝发育区和两组裂缝发育区的区域确定为多组裂缝发育区。
步骤S14：根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据。
在本实施方式中，所述地质资料可以包括所述目标区域已有的测井资料和钻井资料。其中所述地质资料可以包括声波测井数据、成像测井数据、电阻率测井数据、自然伽马测井数据、中子测井数据等数据。
在本实施方式中，所述已知地层的裂缝地质数据可以包括反映所述已知地层裂缝性质的数据。可以包括裂缝密度、裂缝产状、裂缝组数、裂缝孔隙度等。具体地，例如，可以根据FMI(成像测井成果)统计得出裂缝方向的组数。
步骤S16：根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
在本实施方式中，所述根据所述道集数据体中的偏移距信息和方位角信息和所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度可以包括，针对不同的分区内已知地层的裂缝地质数据和对应的保留有方位角和偏移距信息的所述道集数据体，分别对各分区内的裂缝密度进行预测。
在本实施方式中，针对于单组裂缝发育区，可以根据待测地层道集数据体的方位角和偏移距选择合适的数据样点进行椭圆拟合，根据椭圆拟合结果得到椭圆的离心率进而可以得到待测地层所述裂缝的各向异性数据。可以根据已知地层的道集数据体进行椭圆拟合，可以得到已知地层的裂缝各向异性数据，从而可以建立已知地层裂缝各向异性与裂缝密度之间的关系，根据所述已知地层裂缝各向异性与裂缝密度之间的关系以及待测地层的裂缝的各向异性数据可以预测得到所述单组裂缝发育区待测地层的裂缝密度。
在本实施方式中，针对于两组裂缝发育区，可以根据所述待测地层道集数据体的方位角和偏移距以及蜘蛛图分析方法选择一定方位角范围内和偏移距范围内覆盖均匀的样点，计算待测地层的地震波反射能量，根据待测地层的地震波反射能量以及已知地层裂缝密度和地震波反射能量的关系可以预测得到所述两组裂缝发育区待测地层的裂缝密度。
在本实施方式中，针对具有多组方向的裂缝，对其进行椭圆拟合结果趋近于圆，指示出各向同性特征，没有表现出裂缝密度。可以根据所述道集数据体的方位角和偏移距选择各方位分布均匀的样点数据，去除数据奇异的样点，可以根据相干或者曲率等算法预测得到所述待测地层的裂缝密度。
在本实施方式中，通过提取目标区域中偏移距和方位角相同或者趋于相同的地震道记录生成地震道记录，保留了道集数据体的偏移距和方位角信息。根据目标区域已有地质资料得到目标区域已知地层的裂缝地质数据，根据所述裂缝地质数据对所述目标区域进行分区，根据所述道集数据体的偏移距和方位角信息，并结合目标区域已知地层的裂缝地质数据和道集地震数据综合预测所述待测地层的裂缝密度
在一个实施方式中，所述提取目标区域中偏移距和方位角相同或趋于相同的地震道记录的步骤包括，将所述目标区域的观测系统划分为多个十字排子集；其中，所述十字排子集包括至少一个向量片，其中，所述向量片内包括一定偏移距范围内以及一定方位角范围内工区的地震道记录；在所述多个十字排子集中，提取相同偏移距范围内以及相同方位角范围内的向量片。
在本实施方式中，在所述十字排子集可以包括目标区域的观测系统中炮线和检波线相互垂直的所有地震道集的集合。可以以一定的炮线距和接收线距将所述十字排子集进行等距离划分，可以得到多个小矩形，每个小矩形可以是一个向量片，所述向量片中可以包括至少一个具有一定偏移距和方位角的地震道记录。所述目标区域观测系统中可以包括至少一个所述十字排子集。可以将所述多个十字排子集中具有相同或趋于相同方位角和偏移距向量片，进而可以提取目标区域内偏移距和方位角相同或者趋于相同的地震道记录。
在本实施方式中，通过将所述目标区域的观测系统划分为多个十字排子集；所述十字排子集中可以包括至少一个向量片，其中，所述向量片内包括一定偏移距范围内以及一定方位角范围内工区的地震道记录；提取相同偏移距范围内以及相同方位角范围内的向量片，从而将目标区域内偏移距和方位角相同或者趋于相同的地震道记录提取出来，方案易于实现、使用灵活。
在一个实施方式中，所述生成所述道集数据体的步骤中包括，对所述相同或趋于相同偏移距范围以及方位角范围内的向量片进行规则化处理；在所述规则化处理后进行偏移处理，以生成所述道集数据体。
在本实施方式中，由于观测系统、观测设备等各种因素的影响，实际采集的资料可能会造成地下面元属性的非均匀性。例如，面元内各道偏离中心点位置、方位角和偏移距偏以及覆盖次数的不均等。
在本实施方式中，所述相同或趋于相同偏移距范围以及方位角范围内的向量片进行规则化处理的步骤可以包括，先对所述向量片进行面元中心化处理和数据插值处理，在指定的所述向量片中，实际地震道记录的信息可能不在指定向量片的中心，甚至可能没有数据信息。所述中心化处理可以是将信息点规则到中心点，插值可以是通过各个方向的向量片对没有数据信息的向量片进行插值并赋予插值结果。另外，由于受地震采集时观测系统的影响，检波装置接收到的地震信号并不是布满整个空间的，地震采集观测系统可以是有网格特征的，但又不是严格的正方形网格，在面元中心化和插值处理后一个十字排列子集可以时是四方形的，如果按照方位角取值，在相同的偏移距一定时，有的方位可能是没有值的，甚至是相邻的几个偏移距值处都取不到值，需要进行方位角规则化处理。
在本实施方式中，所述偏移处理可以包括，在规则化处理后将所述道集数据体中的反射波归位到产生它们的反射界面上并使绕射波收敛到产生它的绕射点上的。最终可以得到能够反映界面反射系数特点的并正确归位了的地震波形剖面。具体地，例如，由于所述目标区域中OVT向量片数一般等于覆盖次数，为了保留其方位信息，可以进行全部OVT向量片偏移处理，直到所有的OVT向量片偏移完成，从而实现OVT域的偏移成像。
在本实施方式中，通过对所述目标区域中偏移距和方位角相同或趋于相同的地震道记录进行规则化处理，在所述规则化处理后对所述地震道记录进行偏移处理，获得了品质较高的道集数据体。
在一个实施方式中，所述方法还包括，根据所述方位角和偏移距对所述道集数据体进行分布位置规则化处理，以使所述道集数据体由在直角坐标系内表示的形式转化为在极坐标内表示的形式。
在本实施方式中，由于受地震采集时观测系统的影响，接收到的地震信号可以不是布满整个偏移距和方位角范围内的，而是一系列规则分布的离散样点，但这些样点并不是轴对称分布的，请参阅图3，可以根据炮检距增大方向和方位角增大方向对所述道集数据体进行分布位置的规则化处理以使所述道集数据由在直角坐标系内表示的形式转化为在极坐标内表示的形式。保留了道集数据体中的偏移距和方位角信息。
在本实施方式中，通过对所述道集数据体进行分布位置的规则化处理，将所述道集数据体中的方位角信息和炮检距信息同时保留下来，方便直观的确定道集数据体的炮检距和方位角信息。
在一个实施方式中，所述根据所述道集数据体对所述目标区域进行区域划分的步骤包括，根据所述目标区域道集数据体的方位角数据和偏移距数据选择指定样点进行椭圆拟合，以确定所述目标区域的单组裂缝发育区；其中，所述样点用于表示处于一定方位角和偏移距范围内的道集数据体；根据所述目标区域除了单组裂缝发育区以外的道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析，以确定所述目标区域的两组裂缝发育区；将目标区域除了单组裂缝发育区和两组裂缝发育区以外的区域确定为多组裂缝发育区。
在本实施方式中，所述根据所述目标区域道集数据体的方位角数据和偏移距数据选择指定样点进行椭圆拟合，以确定所述目标区域的单组裂缝发育区可以包括，将所述目标区域内的道集数据体分别根据方位角和偏移距选择在一定方位角和偏移距范围内覆盖次数均匀的样点，对与选择出的样点可以对其进行椭圆拟合，其拟合结果如果满足椭圆要求，则可以将拟合后为椭圆的道集数据体所在区域确定为单组裂缝发育区。
在本实施方式中，所述根据所述目标区域除了单组裂缝发育区以外的道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析，以确定所述目标区域的两组裂缝发育区可以包括，排除目标区域中的单组裂缝发育区，将所述目标区域中除单组裂缝发育区以外的道集数据体分别根据方位角和偏移距选择在一定方位角和偏移距范围内覆盖次数均匀的样点，确定可以直接将各个相邻的样点直接连线而成一个闭合的不规则的图形的道集数据体所在区域确定为两组裂缝发育区。
在本实施方式中，可以将所述目标区域除单组裂缝发育区和两组裂缝发育区以外的区域确定为多组裂缝发育区。
在本实施方式中，由于不同裂缝发育情况对应的裂缝密度预测方法可以不同，因此对所述目标区域进行分区有利于提高裂缝密度预测的精度。
在一个实施方式中，对所述单组裂缝发育区内待测地层的裂缝密度进行预测，包括，根据目标区域已知地层的地质资料获取所述单组裂缝发育区内已知地层的裂缝密度和对应的各向异性数据；根据所述已知地层的裂缝密度和对应的各向异性数据得到所述裂缝密度和所述对应的各向异性数据的关系；根据待测地层道集数据体椭圆拟合结果得到待测地层裂缝的各向异性数据；根据所述待测地层的各向异性数据以及已知地层裂缝密度和对应的各向异性数据的关系预测得到单组裂缝发育区内待测地层的裂缝密度。
在本实施方式中，所述各向异性可以包括，当裂缝发育具有明显的方向性时，由于裂缝面与原地层物理性质的不同，地震波在裂缝介质中的传播速度、反射系数、频谱衰减特征、AVO特征等物理特性可以随传播方向与裂缝走向的夹角而变化。
在本实施方式中，单组裂缝发育区，由于裂缝主要表现为某一方向的裂缝发育，可以呈现明显的各向异性特征。地震波传播速度随方位角变化，沿裂缝主方向传播时速度最快，垂直于这个方向传播时速度最慢，整体表现为一个椭圆，所述椭圆的离心率可以表示地层的各向异性数据。
在本实施方式中，可以根据所述已知地层的地质资料获取已知地层的裂缝密度。可以根据已知地层的道集数据体，进行椭圆拟合，计算椭圆离心率可以得到所述已知地层的各向异性数据。可以根据所述已知地层各向异性数据和所述已知地层裂缝密度进行关系拟合，得到已知地层各向异性和裂缝密度的关系式，所述关系式可以是一次方程关系也可以是指数方程关系，可以根据拟合结果得到。
在本实施方式中，所述根据待测地层道集数据体椭圆拟合结果得到待测地层裂缝的各向异性数据包括，根据待测地层道集数据体中的方位角数据和偏移距数据可以去除偏移距，舍去偏移距较远和较近的样点，选择各方位覆盖次数均匀的样点进行椭圆拟合。例如在一定方位角和偏移距范围内选择只被覆盖一次的面元的中样点进行椭圆拟合。计算椭圆离心率可以得到待测地层的各向异性数据。
在本实施方式中，所述根据所述待测地层的各向异性数据和已知地层裂缝密度和对应的各向异性数据的关系预测得到所述待测地层的裂缝密度可以包括，根据得到的待测地层裂缝各向异性数据，可以将其带入已经建立的已知地层各向异性数据与裂缝密度的拟合关系中，可以计算得到对应的裂缝密度。所述裂缝密度可以为所述待测地层的裂缝密度。
在本实施方式中，通过以待测地层的方位角数据和炮检距数据选择样点进行椭圆拟合的方式求得各向异性数据，并根据已知地层的测井数据得到已知地层裂缝密度和对应各向异性数据的关系预测得出待测地层的裂缝密度。在原有椭圆拟合方法的基础上，充分利用了方位角信息，并且与已有的实际数据相结合，提高了单组裂缝发育区裂缝的预测精度。
在一个实施方式中，对所述两组裂缝发育区内待测地层的裂缝密度进行预测，包括，获取所述两组裂缝发育区内已知地层的裂缝密度和对应的道集数据体的地震波反射能量；根据所述已知地层裂缝密度和所述对应的道集数据体的叠加振幅得到所述裂缝密度和所述对应的道集数据体的地震波反射能量的关系；根据待测地层蜘蛛图分析结果得到待测地层道集数据体的地震波反射能量；根据所述待测地层道集数据体的地震波反射能量以及已知地层裂缝密度和对应的道集数据体的地震波反射能量的关系预测得到两组裂缝发育区内待测地层的裂缝密度。
在本实施方式中，所述得到所述地震波反射能量包括，对道集数据体进行叠加处理后可以得到所述道集数据体的叠加振幅，可以根据所述叠加振幅可以计算得到所述地震波反射能量。
在本实施方式中，实际地层中发育有两组方向不同的裂缝，可以根据待测地层道集数据体的方位角和偏移距选择各方位覆盖次数均匀的样点，样点拟合后并不是椭圆形态。所述样点用于表示处于一定方位角和偏移距范围内的道集数据体。
在本实施方式中，所述蜘蛛图分析可以包括，在给定提取一定的偏移距范围内的数据时，按照一定方位角间隔，在极坐标系下或是一个圆形态，或是一个椭圆形态，可以将直接将各个相邻的数据点直接连线而成一个闭合的不规则的图形，在数据点信息可信的基础上所述样点离原点距离最大的方位可能是裂缝发育方向。
在本实施方式中，所述根据所述蜘蛛图分析结果得到待测地层的道集数据体的地震波反射能量可以包括由于在其他条件一定的情况下，裂缝密度越大则波阻抗差异越大，所以道集数据体的能量值越大所对应的裂缝密度越大。可以根据蜘蛛图连接根据方位角数据和偏移距数据选择出的样点，计算样点的地震波反射能量。
在本实施方式中，可以根据已知地层的道集数据体获得已知地层的地震波反射能量，根据已知地层的地质数据可以得到已知地层的裂缝密度，根据已知地层道集数据体的地震波反射能量和裂缝密度可以进行关系拟合，得到已知地层道集数据体的地震波反射能量和裂缝密度的关系式，所述关系式可以是一次方程关系也可以是指数方程关系，可以根据拟合结果得到。
在本实施方式中，根据所述待测地层道集数据体的地震波反射能量以及已知地层裂缝密度和对应的道集数据体的地震波反射能量的关系预测得到两组裂缝发育区内待测地层的裂缝密度可以包括，根据得到的待测地层的地震波反射能量，可以将其带入已经建立的已知地层地震波反射能量与裂缝密度的拟合关系中，可以计算得到对应的裂缝密度。所述裂缝密度可以为所述待测地层的裂缝密度。
在本实施方式中，通过根据待测地层道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析的方式求得待测地层地震波反射能量，并根据已知地层的测井数据得到已知地层裂缝密度和对应地震波反射能量的关系预测得出待测地层的裂缝密度。充分利用了方位角信息，并且与已有的实际数据相结合，提高了两组裂缝发育区裂缝的预测精度。
在一个实施方式中，对所述多组裂缝发育区内待测地层的裂缝密度进行预测，包括，根据待测地层道集数据体的方位角数据和偏移距数据选择指定样点；其中，所述样点用于表示处于一定方位角和偏移距范围内的道集数据体；根据所述指定样点进行叠后处理预测所述待测地层的裂缝密度。
在本实施方式中，根据待测地层道集数据体的方位角数据和偏移距数据选择指定样点可以包括，根据待测地层道集数据体中的方位角数据和偏移距数据可以去除偏移距，舍去偏移距较远和较近的样点，以及奇异的数值点，选择各方位覆盖次数均匀的样点。其中所述奇异的数值点可以包括数值过大或者过小的点。
在本实施方式中，所述叠后处理可以包括相干处理和曲率处理，所述相干处理或曲率处理的处理过程可以参考现有技术中的方法，在此不再赘述。
在本实施方式中，由于裂缝组数达到一定条件时基本可以等效于各向同性介质，通过对多组裂缝发育区进行叠后处理得到对应的裂缝密度，相对于现有技术中无论裂缝组数均使用叠后处理的方法，提高了裂缝密度预测精度。
在一个实施方式中，在所述根据所述道集数据体和所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度的步骤中还包括，对所述目标区域内的已知地层的测井数据和所述道集数据体进行井震标定处理；根据井震标定结果确定所述分区对应的有效裂缝密度阈值；根据所述有效裂缝密度阈值获取所述目标区域有效裂缝的空间分布范围。
在本实施方式中，所述井震标定可以包括，将测井数据与道集数据进行匹配的过程。例如测井数据得到裂缝孔隙度，根据道集数据得到裂缝密度，将裂缝密度与所述裂缝孔隙度进行匹配的过程可以是井震标定的过程。在井震标定之前，由于测井资料中存在噪音以及不同年代的处理资料存在差异，可以对测井资料进行环境校正和标准化处理，以使不同井的测井曲线在数据上保持统一。在标准化处理后可以应用VSP资料对声波曲线进行校正，提高标定效果。
在本实施方式中，所述根据井震标定结果确定所述分区对应的有效裂缝密度阈值可以包括，根据井震标定结果可以在所述分区内的道集数据体中提取对应的裂缝密度值，可以将提取的裂缝密度值与测井曲线进行匹配，在测井曲线中可以选择裂缝密度值与测井曲线匹配度较高的区域对应的裂缝密度范围确定为所述有效裂缝密度的阈值区间。具体地，例如，请参阅图13，在单组裂缝发育区，将道集数据体得到的裂缝密度与测井数据得到的裂缝孔隙度进行井震标定处理，提取对应的裂缝密度，并与裂缝孔隙度进行交汇，确定有效裂缝密度的阈值为0.18。
在本实施方式中，通过井震标定，将道集数据与实际测井数据相融合，提高了裂缝预测的精度。
请参阅图21，本说明书实施方式还提供一种裂缝密度预测装置，包括，道集数据体生成模块，用于提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；区域划分模块，用于根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；裂缝地质数据获取模块，用于根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；裂缝密度预测模块，用于根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
本实施方式中提供的裂缝密度预测装置，其实现的功能和效果可以参见其它实施方式对照解释。
本说明书实施方式中所描述的各种说明性的模块都可以通过通用处理器，数字信号处理器，专用集成电路，现场可编辑门阵列或其他可编程逻辑装置，离散硬部件，或上述任何组合的设计来实现或操作所描述的功能。通用处理器可以为微处理器，可选地，该通用处理器也可以为任何传统的处理器、控制器、微控制器或状态机。处理器也可以通过计算装置的组合来实现，例如数字信号处理器，多个微处理器，一个或多个微处理器联合一个数字信号处理器核，或其他类似的配置来实现。
本说明书实施方式还提供一种计算机存储介质，所述存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现，提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度。
本实施方式中提供的存储介质，其程序指令被执行时实现的功能和效果可以参见其它实施方式对照解释。
在本实施方式中，所述存储介质包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。
本说明书实施方式提供的裂缝密度预测方法、装置及计算机存储介质，通过提取目标区域中偏移距和方位角相同或趋于相同的地震道记录生成道集数据体。根据已有地质资料对所述目标区域进行分区，根据所述道集数据体和所述分区的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度，在处理过程中考虑到裂缝发育各向异性的特征，保留地震道记录的偏移距和方位角信息，提高了目标区域裂缝密度预测的精度。
本说明书实施方式中所描述的方法或算法的步骤可以直接嵌入硬件、处理器执行的软硬模快、或者两者的结合。软件模块可以存储于RAM存储器、闪存、ROM存储器、EPROM存储器，寄存器、硬盘、可移动磁盘、CD-ROM或本领域中其他任意形式的存储媒介中。
本说明书实施方式所描述的上述功能可以在硬件、软件、固件或这三者的任意组合来实现。如果在软件中实现，这些功能可以存储与电脑可读的媒介上，或以一个或多个指令或代码的形式传输于电脑可读媒介上。电脑可读媒介包括电脑存储媒介和便于使得让电脑程序从一个地方转移到其他地方的通信媒介。存储媒介可以是任何通用或特殊电脑可以接入访问的可用媒体。例如，这样的电脑可读媒体可以包括但不限于RAM，ROM，EPROM，CD-ROM或其他光盘存储、磁盘存储或其他磁性存储装置，或其他任何可以用于承载或存储以指令或数据结构和其它可被通用或特殊电脑、或者通用或特殊处理器读取形式的程序代码的媒介。此外，任何连接都可以被适当地定义为电脑可读媒介，例如，如果软件是从一个网站站点、服务器或其他远程资源通过一个同轴电缆、光纤电缆、双绞线、数字用户线(DSL)或以例如红外、无线和微波等无线传输的也被包含在所定义的电脑可读媒介中。所述的碟片和磁盘，包括压缩磁盘、镭射盘、光盘、DVD、软盘和蓝光光盘，磁盘通常以磁性复制数据，而碟片通常以激光进行光学复制数据。上述组合也可以包含在电脑可读媒介中。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于装置及存储介质实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的说明书部分即可。
虽然通过实施例描绘了本说明书，本领域普通技术人员知道，本说明书有许多变形和变化而不脱离本说明书的精神，希望所附的权利要求包括这些变形和变化而不脱离本说明书的精神。

Claims (10)

1.一种裂缝密度预测方法，其特征在于，所述方法包括：
提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；
根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；其中，根据所述道集数据体对所述目标区域进行区域划分的步骤包括：根据所述目标区域除了单组裂缝发育区以外的道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析，以确定所述目标区域的两组裂缝发育区；
根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；
根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度；其中，对所述两组裂缝发育区内待测地层的裂缝密度进行预测，包括：获取所述两组裂缝发育区内已知地层的裂缝密度和对应的道集数据体的地震波反射能量；根据所述已知地层裂缝密度和所述对应的道集数据体的叠加振幅得到所述裂缝密度和所述对应的道集数据体的地震波反射能量的关系；根据待测地层蜘蛛图分析结果得到待测地层道集数据体的地震波反射能量；根据所述待测地层道集数据体的地震波反射能量以及已知地层裂缝密度和对应的道集数据体的地震波反射能量的关系预测得到两组裂缝发育区内待测地层的裂缝密度。
2.根据权利要求1所述的方法，其特征在于，提取目标区域中偏移距和方位角符合指定条件的地震道记录的步骤包括：
将所述目标区域的观测系统划分为多个十字排子集；其中，所述十字排子集包括至少一个向量片，其中，所述向量片内包括一定偏移距范围内以及一定方位角范围内工区的地震道记录；
在所述多个十字排子集中，提取相同或趋于相同偏移距范围以及方位角范围内的向量片。
3.根据权利要求2所述的方法，其特征在于，所述生成所述道集数据体的步骤中包括：
对所述相同或趋于相同偏移距范围以及方位角范围内的向量片进行规则化处理；
在所述规则化处理后进行偏移处理，以生成所述道集数据体。
4.根据权利要求3所述的方法，其特征在于，所述方法还包括：
根据所述方位角和偏移距对所述道集数据体进行分布位置规则化处理，以使所述道集数据体由在直角坐标系内表示的形式转化为在极坐标内表示的形式。
5.根据权利要求1所述的方法，其特征在于，所述根据所述道集数据体对所述目标区域进行区域划分的步骤还包括：
根据所述目标区域道集数据体的方位角数据和偏移距数据选择指定样点进行椭圆拟合，以确定所述目标区域的单组裂缝发育区；其中，所述样点用于表示处于一定方位角和偏移距范围内的道集数据体；
将所述目标区域除了单组裂缝发育区和两组裂缝发育区以外的区域确定为多组裂缝发育区。
6.根据权利要求5所述的方法，其特征在于，对所述单组裂缝发育区内待测地层的裂缝密度进行预测，包括：
根据目标区域已知地层的地质资料获取所述单组裂缝发育区内已知地层的裂缝密度和对应的各向异性数据；
根据所述已知地层的裂缝密度和对应的各向异性数据得到所述裂缝密度和所述对应的各向异性数据的关系；
根据待测地层道集数据体椭圆拟合结果得到待测地层裂缝的各向异性数据；
根据所述待测地层的各向异性数据以及已知地层裂缝密度和对应的各向异性数据的关系预测得到单组裂缝发育区内待测地层的裂缝密度。
7.根据权利要求5所述的方法，其特征在于，对所述多组裂缝发育区内待测地层的裂缝密度进行预测，包括：
根据待测地层道集数据体的方位角数据和偏移距数据选择指定样点；其中，所述样点用于表示处于一定方位角和偏移距范围内的道集数据体；
根据所述指定样点进行叠后处理预测所述待测地层的裂缝密度。
8.根据权利要求1所述的方法，其特征在于，在所述根据所述道集数据体和所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度的步骤中还包括：
对所述目标区域内的已知地层的测井数据和所述道集数据体进行井震标定处理；
根据井震标定结果确定所述分区对应的有效裂缝密度阈值；
根据所述有效裂缝密度阈值获取所述目标区域有效裂缝的空间分布范围。
9.一种裂缝密度预测装置，其特征在于，包括：
道集数据体生成模块：用于提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；
区域划分模块：用于根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；其中，根据所述道集数据体对所述目标区域进行区域划分的步骤包括：根据所述目标区域除了单组裂缝发育区以外的道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析，以确定所述目标区域的两组裂缝发育区；
裂缝地质数据获取模块：用于根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；
裂缝密度预测模块：用于根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度；其中，对所述两组裂缝发育区内待测地层的裂缝密度进行预测，包括：获取所述两组裂缝发育区内已知地层的裂缝密度和对应的道集数据体的地震波反射能量；根据所述已知地层裂缝密度和所述对应的道集数据体的叠加振幅得到所述裂缝密度和所述对应的道集数据体的地震波反射能量的关系；根据待测地层蜘蛛图分析结果得到待测地层道集数据体的地震波反射能量；根据所述待测地层道集数据体的地震波反射能量以及已知地层裂缝密度和对应的道集数据体的地震波反射能量的关系预测得到两组裂缝发育区内待测地层的裂缝密度。
10.一种计算机存储介质，所述存储介质存储有计算机程序指令，在所述计算机程序指令被执行时实现：
提取目标区域中偏移距和方位角符合指定条件的地震道记录生成至少一个道集数据体；
根据所述道集数据体中的偏移距信息和方位角信息对所述目标区域的裂缝发育情况进行区域划分；其中，所述目标区域被划分为至少一个分区；其中，根据所述道集数据体对所述目标区域进行区域划分的步骤包括：根据所述目标区域除了单组裂缝发育区以外的道集数据体的方位角数据和偏移距数据选择指定样点进行蜘蛛图分析，以确定所述目标区域的两组裂缝发育区；
根据目标区域已知地层的地质资料，获取所述已知地层的裂缝地质数据；
根据所述道集数据体中的偏移距信息和方位角信息以及所述分区内已知地层的所述裂缝地质数据预测所述目标区域的待测地层中裂缝的裂缝密度；其中，对所述两组裂缝发育区内待测地层的裂缝密度进行预测，包括：获取所述两组裂缝发育区内已知地层的裂缝密度和对应的道集数据体的地震波反射能量；根据所述已知地层裂缝密度和所述对应的道集数据体的叠加振幅得到所述裂缝密度和所述对应的道集数据体的地震波反射能量的关系；根据待测地层蜘蛛图分析结果得到待测地层道集数据体的地震波反射能量；根据所述待测地层道集数据体的地震波反射能量以及已知地层裂缝密度和对应的道集数据体的地震波反射能量的关系预测得到两组裂缝发育区内待测地层的裂缝密度。
CN201810170528.1A 2018-03-01 2018-03-01 裂缝密度预测方法、装置及计算机存储介质 Active CN108646292B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810170528.1A CN108646292B (zh) 2018-03-01 2018-03-01 裂缝密度预测方法、装置及计算机存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810170528.1A CN108646292B (zh) 2018-03-01 2018-03-01 裂缝密度预测方法、装置及计算机存储介质

Publications (2)

Publication Number Publication Date
CN108646292A CN108646292A (zh) 2018-10-12
CN108646292B true CN108646292B (zh) 2020-01-07

Family

ID=63744233

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810170528.1A Active CN108646292B (zh) 2018-03-01 2018-03-01 裂缝密度预测方法、装置及计算机存储介质

Country Status (1)

Country Link
CN (1) CN108646292B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112578432B (zh) * 2019-09-27 2023-05-26 中国石油化工股份有限公司 一种正交各向异性介质ovt域成像道集提取方法及系统
CN112780253B (zh) * 2020-01-20 2022-05-10 中国石油天然气集团有限公司 裂缝型储层的预测及评价方法
CN111506861B (zh) * 2020-04-24 2020-11-13 成都捷科思石油天然气技术发展有限公司 一种目的层有利区域裂缝强度计算方法
CN114814916B (zh) * 2022-04-01 2023-09-26 煤炭科学研究总院有限公司 地表裂缝定位方法、装置及电子设备

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103513277B (zh) * 2013-09-27 2016-11-16 中国石油天然气股份有限公司 一种地震地层裂隙裂缝密度反演方法及系统
CN105425292B (zh) * 2015-11-20 2019-04-12 中国石油天然气集团公司 一种油气预测方法及装置
CN107045143B (zh) * 2017-02-22 2018-12-28 中国石油天然气股份有限公司 一种预测裂缝发育的方法及装置

Also Published As

Publication number Publication date
CN108646292A (zh) 2018-10-12

Similar Documents

Publication Publication Date Title
CN108646292B (zh) 裂缝密度预测方法、装置及计算机存储介质
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
Verdon et al. Microseismic monitoring using a fiber-optic distributed acoustic sensor array
US8724429B2 (en) System and method for performing time-lapse monitor surverying using sparse monitor data
CA2964893C (en) Structure tensor constrained tomographic velocity analysis
CN106556861B (zh) 一种基于全方位地震资料的方位avo反演方法
US9268776B2 (en) Methods and apparatus for data collection
RU2339056C2 (ru) Обобщенное трехмерное прогнозирование кратных волн от поверхности
CA2961168A1 (en) Integrating vertical seismic profile data for microseismic anisotropy velocity analysis
CN109469477B (zh) 一种人工裂缝延伸方向的预测方法和装置
CN108594300B (zh) 地貌成像方法、装置及计算机存储介质
US10705241B2 (en) Determining sea water resistivity
Feng et al. The application of high-resolution 3D seismic acquisition techniques for carbonate reservoir characterization in China
CN106597549B (zh) 一种反射波提取方法及装置、反射声波成像测井仪器
CN110579799B (zh) 一种等旅行时间间隔的地震采集观测方法及系统
Li et al. Gaussian beam imaging of fractures near the wellbore using sonic logging tools after removing dispersive borehole waves
Jiang et al. Velocity model optimization for surface microseismic monitoring via amplitude stacking
CN104297786B (zh) 检测地层裂缝走向方位角度的方法和装置
CN106526668B (zh) 原始波形提取与成像方法
CN111474580B (zh) 一种基于炮检距矢量片的方位角道集提取方法和系统
Misaghi et al. Overburden complexity and repeatability of seismic data: Impacts of positioning errors at the Oseberg field, North Sea
CN110579798B (zh) 一种等反射角度间隔的地震采集观测方法及系统
CN113700466B (zh) 深层碳酸盐岩油气藏油气检测方法、装置、设备及介质
CN105989240A (zh) 砂层预测方法及装置
CN112630832A (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