CN113568046A - 裂缝方位约束的多维信息裂缝表征方法 - Google Patents

裂缝方位约束的多维信息裂缝表征方法 Download PDF

Info

Publication number
CN113568046A
CN113568046A CN202110854228.7A CN202110854228A CN113568046A CN 113568046 A CN113568046 A CN 113568046A CN 202110854228 A CN202110854228 A CN 202110854228A CN 113568046 A CN113568046 A CN 113568046A
Authority
CN
China
Prior art keywords
crack
fracture
azimuth
dimensional information
difference
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
CN202110854228.7A
Other languages
English (en)
Other versions
CN113568046B (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 Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
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 Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN202110854228.7A priority Critical patent/CN113568046B/zh
Publication of CN113568046A publication Critical patent/CN113568046A/zh
Application granted granted Critical
Publication of CN113568046B publication Critical patent/CN113568046B/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/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
    • 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
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • 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/63—Seismic attributes, e.g. amplitude, polarity, instant phase

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

一种裂缝方位约束的多维信息裂缝表征方法，包括以下步骤：一，选择目标工区钻至目的地层的参考井；二，对裂缝储层特征进行分析；三，对窄方位地震资料进行方位划分，并进行裂缝方向预测；四，基于方位属性差的裂缝密度预测；校正后的裂缝发育强度采用如下公式：
Figure DDA0003183502510000011
Figure DDA0003183502510000012
五，基于远近偏移距属性差的裂缝密度预测；六，叠前裂缝孔隙度反演；七，裂缝方向约束下的多维信息融合；融合公式为：
Figure DDA0003183502510000013
本发明不仅能够补充窄方位地震资料进行裂缝储层预测时的方位信息的不足的问题，全面准确的刻画潜山裂缝储层发育特征；而且，还能够提高窄方位地震资料进行潜山裂缝型储层预测的精度。

Description

裂缝方位约束的多维信息裂缝表征方法
技术领域
本发明属于油气田勘探领域，尤其涉及裂缝方位约束的多维信息裂缝表征方法。
背景技术
目前，裂缝型储层分布广泛，油气储量规模巨大，随着勘探的不断深入，中深层裂缝型储层逐渐成为勘探的主要目标，对裂缝型储层的研究也成为行业内热点。现有的对裂缝型储层预测的方法主要有：多分量转换横波预测方法、基于纵波方位各向异性的预测方法、叠后地震多属性分析方法以及应力场分析方法等。其中，由于多分量地震资料的采集和处理耗费成本较大，未能得到普遍应用。而基于纵波方位各向异性的预测方法、叠后多属性分析方法和应力场分析分析方法等技术却得到广泛的应用，并成为裂缝型储层研究的主流方法。
随着陆上全方位和宽方位资料采集和应用的增加，基于等效各向异性介质理论的纵波方位各向异性裂缝的预测方法也得到推广，并取得了良好的效果。但是，在海上勘探中，由于受到勘探风险、勘探成熟度和精度以及施工条件等限制，其宽方位资料的采集和应用相对较少。因此，导致现有的海上拖缆地震资料的方位信息不足，而单纯的利用纵波方位各向异性方法不能取得满意的效果。
目前，针对海上窄方位资料进行裂缝型储层预测的研究相对较少，如何充分挖掘和利用窄方位地震资料的有效信息，已成为海上深埋潜山裂缝型油气藏勘探评价的关键。
发明内容
本发明目的在于提供一种裂缝方位约束的多维信息裂缝表征方法,以解决利用海上拖缆采集的窄方位地震资料进行潜山裂缝型储层预测时的方位信息不充足的技术问题。
为实现上述目的，本发明的裂缝方位约束的多维信息裂缝表征方法的具体技术方案如下：
一种裂缝方位约束的多维信息裂缝表征方法，包括以下步骤：
第一步，选择目标工区钻至目的地层的参考井；
第二步，对裂缝储层特征进行分析：设有以下两种方式：
1)对成像测井裂缝储层特征进行统计分析：
2)对介质进行岩石物理正演；
第三步，对窄方位地震资料进行方位划分，并进行裂缝方向预测；
⑴对参考井旁的地震资料进行椭圆拟合分析；
⑵确定适合用于椭圆拟合的偏移距范围和地震属性；包括：振幅属性、阻抗属性、频率属性以及衰减属性；
第四步，基于方位属性差的裂缝密度预测；
1)反射波振幅变化是随着入射角的增加而增大；
2)当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感，而在平行和垂直裂缝方向敏感程度较低；
3)裂缝发育强度和观测方向与裂缝方向的夹角会同时对预测结果产生影响，从而，给裂缝预测结果带来不确定性；
4)为了消除由于采集方向与裂缝方向的夹角
Figure BDA0003183502490000021
给预测裂缝强度带来的影响，进一步对各向异性预测裂缝密度进行校正；通过引入校正因子κ，令
Figure BDA0003183502490000022
为了确保校正因子的稳定性，同时，引入γ作为稳定系数，于是，
Figure BDA0003183502490000023
最终，校正后的裂缝发育强度为：
Figure BDA0003183502490000031
第五步，基于远近偏移距属性差的裂缝密度预测；
第六步，叠前裂缝孔隙度反演；
第七步，裂缝方向约束下的多维信息融合；融合公式为：
Figure BDA0003183502490000032
进一步，所述第一步中，选择目标工区后，对参考井进行合成地震记录的制作和时深标定，且输出时深关系。
进一步，所述第二步中，1)对成像测井裂缝储层特征进行统计分析：利用成像测井的解释结果进行裂缝密度或裂缝孔隙度、裂缝产状，包括：走向、倾向、倾角的裂缝特征的统计分析；2)对介质进行岩石物理正演，包括：根据上面分析的裂缝特征，结合纵、横波速度和密度测井曲线，进行各向异性岩石物理正演，并分析裂缝等效介质的AVAZ特征。
所述第三步中，对窄方位地震资料进行方位划分，并进行裂缝方向预测的指导原则为：
进一步，所述第三步中，对窄方位地震资料进行方位划分，并进行裂缝方向预测的指导原则为：
①使得拟合的椭圆方向，能够基本与实际裂缝方向一致；
②指导确定方位划分方案；
③将地震资料进行分方位叠加，并计算分析得到的敏感属性，从而进行椭圆拟合得到最终的预测裂缝方向。
进一步，所述第四步中，公式(1)中的DA：为方位属性差异预测裂缝发育强度，Df：为校正后裂缝发育强度，
Figure BDA0003183502490000033
为裂缝方与采集方向的夹角，γ：为稳定系数，取值范围为0＜γ＝1。
进一步，所述第四步中，当稳定因子为0.01时，通过校正能够有效的减弱采集方向与裂缝方向夹角，对预测结果带来的影响，使得结果能够更大程度的反映裂缝发育强度。
进一步，所述第五步中，对裂缝密度预测：①先对地震资料进行针对性处理，并进行远近偏移距部分叠加；②结合测井资料，对井点处地震属性进行分析，并优选出敏感属性；求取远近偏移距的敏感属性差异，从而得到裂缝密度。
进一步，所述第六步中，对叠前裂缝孔隙度反演过程：⑴首先，以裂缝总孔隙度曲线作为目标曲线，横波阻抗曲线与密度曲线作为计算曲线，通过参与计算曲线的线性匹配关系模拟目标曲线，并与先验目标曲线进行误差分析，优选误差最小的组合关系作为最终的双参数拟合关系输出；⑵将双参数拟合关系应用于叠前反演的横波阻抗体和密度体从而实现叠前裂缝孔隙度反演。
进一步，所述第七步中，对多维信息融合如下：①当裂缝方向与观测方向垂直(90°)时，入射角30°与入射角10°的反射振幅差异达到最大；而当裂缝方向与观测方向平行(0°)时，入射角30°与入射角10°的反射振幅差异达到最小；这说明了，远近偏移距属性差异在垂直裂缝方向最为敏感，且随着观测方向与裂缝方向夹角的变小；因此，其能够利用远近偏移距属性差异进行垂直采集方向裂缝的检测；
②当裂缝与采集方向平行时，其AVO与常规各向同性介质是一致的，因此，其能够利用各向同性的反演方法；同时，在窄方位资料情况下，当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感；以实现多维信息的互补，从而实现裂缝的综合表征。
进一步，所述第七步中，公式(2)中的DA：为椭圆拟合结果，DO：为远近偏移距属性差结果，DI：为叠前裂缝孔隙度反演结果。
本发明的裂缝方位约束的多维信息裂缝表征方法具有以下优点：
1、本发明从岩石物理正演分析入手，明确了窄方位地震资料进行裂缝方向预测的可行性。
2、本发明通过岩石物理交汇分析，建立了弹性参数到裂缝孔隙度的关系，完成了弹性参数到物性参数的叠前反演，实现了叠前双参数裂缝孔隙度反演。
3、本发明从裂缝等效各项异性介质的机理研究出发，结合海上拖缆采集的窄方位地震资料的特征，根据方位各向异性、远近偏移距属性差异以及叠前裂缝孔隙度反演的敏感范围和互补特征，实现了在裂缝方向约束下的多维信息裂缝综合表征技术。
4、本发明能够有效的补充窄方位地震资料进行裂缝储层预测时的方位信息的不足的问题，全面准确的刻画潜山裂缝储层发育特征，实现了对裂缝型储层的精细描述。
5、本发明能够有效地提高窄方位地震资料进行潜山裂缝型储层预测的精度。
附图说明
图1为本发明的流程示意图；
图2A为本发明第一实施例中真实裂缝方向(成像测井)示意图；(其为屏幕上的实际图形)
图2B为本发明第一实施例中预测裂缝方向(道集)示意图；(其为屏幕上的实际图形)
图2C为本发明第一实施例中预测裂缝方向(分方位叠加)示意图；(其为屏幕上的实际图形)
图2D为本发明第二实施例中真实裂缝方向(成像测井)示意图；(其为屏幕上的实际图形)
图2E为本发明第二实施例中预测裂缝方向(道集)示意图；(其为屏幕上的实际图形)
图2F为本发明第二实施例中预测裂缝方向(分方位叠加)示意图；(其为屏幕上的实际图形)
图3为本发明的预测裂缝方向的平面玫瑰示意图；(其为屏幕上的实际图形)
图4A为本发明不同裂缝密度方位振幅曲线示意图；(其为屏幕上的实际图形)
图4B为本发明不同方位振幅差异度曲线示意图；(其为屏幕上的实际图形)
图4C为本发明校正结果曲线示意图；(其为屏幕上的实际图形)
图5A为本发明校正前方位各向异性预测裂缝密度平面示意图；(其为屏幕上的实际图形)
图5B为本发明校正后方位各向异性预测裂缝密度平面示意图；(其为屏幕上的实际图形)
图6为本发明利用远近偏移距属性差预测裂缝孔隙度平面示意图；(其为屏幕上的实际图形)
图7为本发明利用横波阻抗和密度曲线拟合裂缝孔隙度曲线示意图；(其为屏幕上的实际图形)
图8为本发明叠前裂缝孔隙度反演平面示意图；(其为屏幕上的实际图形)
图9A为本发明不同入射角的方位振幅曲线示意图；(其为屏幕上的实际图形)
图9B为本发明不同方位角AVO曲线示意图；(其为屏幕上的实际图形)
图10为本发明各属性敏感范围以互补融合示意图；(其为屏幕上的实际图形)
图11为本发明裂缝方位约束下多维信息融合预测裂缝密度平面示意图。
(其为屏幕上的实际图形)
具体实施方式
为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明裂缝方位约束的多维信息裂缝表征方法做进一步详细的描述。
如图1-图11所示，本发明包括以下步骤：
第一步，如图1所示，选择目标工区钻至目的地层，不仅要求井径质量较好，测井曲线完整，同时，具有电成像测井的已钻井，作为指导和验证的参考井；并对参考井进行合成地震记录的制作和时深标定，且输出时深关系；
第二步，对裂缝储层特征进行分析：
1)对成像测井裂缝储层特征进行统计分析：
利用成像测井的解释结果进行裂缝密度或裂缝孔隙度、裂缝产状(包括：走向、倾向、倾角)等裂缝特征的统计分析；
2)对介质进行岩石物理正演；
根据上面分析的裂缝特征，结合纵、横波速度和密度测井曲线，进行各向异性岩石物理正演，并分析裂缝等效介质的AVAZ特征；
第三步，对窄方位地震资料进行方位划分，并进行裂缝方向预测；
如图2A-图2F所示，⑴对参考井旁的地震资料进行椭圆拟合分析；⑵利用第二步中得到的裂缝特征和AVAZ特征作为指导，确定适合用于椭圆拟合的偏移距范围和地震属性；包括：但不限于振幅属性、阻抗属性、频率属性以及衰减属性等，其指导原则为：①使得拟合的椭圆方向，能够基本与实际裂缝方向一致；在此基础上，②指导确定方位划分方案；③根据分析得到的方位划分方案，将地震资料进行分方位叠加，并计算分析得到的敏感属性，从而进行椭圆拟合得到最终的预测裂缝方向(如图3所示)。
第四步，基于方位属性差的裂缝密度预测；
如图4A所示，1)由于各向异性导致的反射波振幅变化，随着入射角的增加而增大；如图4B所示，为方位差为11°时振幅最大差异度，2)当观测方向与裂缝方向夹角为45°时，其由于各向异性导致的振幅差异最大，而观测方向与裂缝方向平行和垂直情况下，振幅无差异；这说明，在窄方位资料情况下，当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感，而在平行和垂直裂缝方向敏感程度较低。3)同时，裂缝发育强度和观测方向与裂缝方向的夹角会同时对预测结果产生影响，从而，给裂缝预测结果带来不确定性；4)为了消除由于采集方向与裂缝方向的夹角
Figure BDA0003183502490000081
给预测裂缝强度带来的影响，进一步对各向异性预测裂缝密度进行校正；通过引入校正因子κ，令
Figure BDA0003183502490000082
为了确保校正因子的稳定性，同时，引入γ作为稳定系数，于是，
Figure BDA0003183502490000083
最终，校正后的裂缝发育强度为：
Figure BDA0003183502490000084
其中，DA：为方位属性差异预测裂缝发育强度，Df：为校正后裂缝发育强度，
Figure BDA0003183502490000085
为裂缝方与采集方向的夹角，γ：为稳定系数，取值范围为0＜γ＝1。
如图4C所示，为不同稳定系数下，对方位属性差异曲线进行校正的结果，从中可以看出，当稳定因子为0.01时，通过校正可以有效的减弱采集方向与裂缝方向夹角，对预测结果带来的影响，使得结果能够更大程度的反映裂缝发育强度。
如图5B所示，其为校正后方位各向异性预测裂缝密度平面图。
第五步，基于远近偏移距属性差的裂缝密度预测；
①对地震资料进行针对性处理，并进行远近偏移距部分叠加；
②结合测井资料，对井点处地震属性进行分析，并优选出敏感属性；求取远近偏移距的敏感属性差异，从而得到裂缝密度(如图6所示)。
第六步，叠前裂缝孔隙度反演；
如图7所示，⑴首先，以裂缝总孔隙度曲线作为目标曲线，横波阻抗曲线与密度曲线作为计算曲线，通过参与计算曲线的线性匹配关系模拟目标曲线，并与先验目标曲线进行误差分析，优选误差最小的组合关系作为最终的双参数拟合关系输出；
如图8所示，⑵将双参数拟合关系应用于叠前反演的横波阻抗体和密度体从而实现叠前裂缝孔隙度反演；
第七步，裂缝方向约束下的多维信息融合；
如图9A所示，①从图9A中，不同入射角的振幅曲线可以看出，当裂缝方向与观测方向垂直(90°)时，入射角30°与入射角10°的反射振幅差异达到最大；而当裂缝方向与观测方向平行(0°)时，入射角30°与入射角10°的反射振幅差异达到最小；这说明了，远近偏移距属性差异在垂直裂缝方向最为敏感，且随着观测方向与裂缝方向夹角的变小，该方法对裂缝的检测能力降低；因此，其可以利用远近偏移距属性差异进行垂直采集方向裂缝的检测；
如图9B所示，②从图9B不同方位角的AVO曲线可以看出，裂缝与采集方向平行时，其AVO与常规各向同性介质是一致的，因此，其可以利用各向同性的反演方法；同时，在上述步骤3中分析可知，在窄方位资料情况下，当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感；在这些认识的基础上，就可以在步骤3中预测的裂缝方位约束下，实现多维信息的互补，从而实现裂缝的综合表征；
如图10所示，根据各向异性、远近偏移距属性差以及叠前裂缝孔隙度反演的控制范围，在针对不同裂缝方位的裂缝检测结果的基础上，从各个方法和属性的可检测敏感范围和互补性出发，充分利用每个属性的优势，在裂缝方向的约束下进行融合，实现了多维信息综合裂缝表征的结果；其融合方法如下式：
Figure BDA0003183502490000091
其中，DA：为椭圆拟合结果，DO：为远近偏移距属性差结果，DI：为叠前裂缝孔隙度反演结果。
基于裂缝方位约束的多维信息裂缝表征方法，实现了窄方位资料裂缝带储层非均质性综合表征，精细描述了潜山裂缝储层发育规律，在海上拖缆采集窄方位地震资料进行潜山裂缝储层预测研究当中具有较好的应用前景。
上述未作说明的技术为现有技术，故不再赘述。
可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。

Claims (10)

1.一种裂缝方位约束的多维信息裂缝表征方法，其特征在于，包括以下步骤：
第一步，选择目标工区钻至目的地层的参考井；
第二步，对裂缝储层特征进行分析：设有以下两种方式：
1)对成像测井裂缝储层特征进行统计分析：
2)对介质进行岩石物理正演；
第三步，对窄方位地震资料进行方位划分，并进行裂缝方向预测；
⑴对参考井旁的地震资料进行椭圆拟合分析；
⑵确定适合用于椭圆拟合的偏移距范围和地震属性；包括：振幅属性、阻抗属性、频率属性以及衰减属性；
第四步，基于方位属性差的裂缝密度预测；
1)反射波振幅变化是随着入射角的增加而增大；
2)当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感，而在平行和垂直裂缝方向敏感程度较低；
3)裂缝发育强度和观测方向与裂缝方向的夹角会同时对预测结果产生影响，从而，给裂缝预测结果带来不确定性；
4)为了消除由于采集方向与裂缝方向的夹角
Figure FDA0003183502480000011
给预测裂缝强度带来的影响，进一步对各向异性预测裂缝密度进行校正；通过引入校正因子κ，令
Figure FDA0003183502480000012
为了确保校正因子的稳定性，同时，引入γ作为稳定系数，于是，
Figure FDA0003183502480000013
最终，校正后的裂缝发育强度为：
Figure FDA0003183502480000014
第五步，基于远近偏移距属性差的裂缝密度预测；
第六步，叠前裂缝孔隙度反演；
第七步，裂缝方向约束下的多维信息融合；融合公式为：
Figure FDA0003183502480000021
2.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第一步中，选择目标工区后，对参考井进行合成地震记录的制作和时深标定，且输出时深关系。
3.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第二步中，1)对成像测井裂缝储层特征进行统计分析：利用成像测井的解释结果进行裂缝密度或裂缝孔隙度、裂缝产状，包括：走向、倾向、倾角的裂缝特征的统计分析；2)对介质进行岩石物理正演，包括：根据上面分析的裂缝特征，结合纵、横波速度和密度测井曲线，进行各向异性岩石物理正演，并分析裂缝等效介质的AVAZ特征。
4.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第三步中，对窄方位地震资料进行方位划分，并进行裂缝方向预测的指导原则为：
①使得拟合的椭圆方向，能够基本与实际裂缝方向一致；
②指导确定方位划分方案；
③将地震资料进行分方位叠加，并计算分析得到的敏感属性，从而进行椭圆拟合得到最终的预测裂缝方向。
5.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第四步中，公式(1)中的DA：为方位属性差异预测裂缝发育强度，Df：为校正后裂缝发育强度，
Figure FDA0003183502480000022
为裂缝方与采集方向的夹角，γ：为稳定系数，取值范围为0＜γ＝1。
6.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第四步中，当稳定因子为0.01时，通过校正能够有效的减弱采集方向与裂缝方向夹角，对预测结果带来的影响，使得结果能够更大程度的反映裂缝发育强度。
7.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第五步中，对裂缝密度预测：①先对地震资料进行针对性处理，并进行远近偏移距部分叠加；②结合测井资料，对井点处地震属性进行分析，并优选出敏感属性；求取远近偏移距的敏感属性差异，从而得到裂缝密度。
8.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第六步中，对叠前裂缝孔隙度反演过程：⑴首先，以裂缝总孔隙度曲线作为目标曲线，横波阻抗曲线与密度曲线作为计算曲线，通过参与计算曲线的线性匹配关系模拟目标曲线，并与先验目标曲线进行误差分析，优选误差最小的组合关系作为最终的双参数拟合关系输出；⑵将双参数拟合关系应用于叠前反演的横波阻抗体和密度体从而实现叠前裂缝孔隙度反演。
9.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第七步中，对多维信息融合如下：①当裂缝方向与观测方向垂直(90°)时，入射角30°与入射角10°的反射振幅差异达到最大；而当裂缝方向与观测方向平行(0°)时，入射角30°与入射角10°的反射振幅差异达到最小；这说明了，远近偏移距属性差异在垂直裂缝方向最为敏感，且随着观测方向与裂缝方向夹角的变小；因此，其能够利用远近偏移距属性差异进行垂直采集方向裂缝的检测；
②当裂缝与采集方向平行时，其AVO与常规各向同性介质是一致的，因此，其能够利用各向同性的反演方法；同时，在窄方位资料情况下，当观测方向与裂缝方向存在一定夹角时，方位属性差异裂缝检测方法较敏感；以实现多维信息的互补，从而实现裂缝的综合表征。
10.根据权利要求1所述的裂缝方位约束的多维信息裂缝表征方法，其特征在于，所述第七步中，公式(2)中的DA：为椭圆拟合结果，DO：为远近偏移距属性差结果，DI：为叠前裂缝孔隙度反演结果。
CN202110854228.7A 2021-07-28 2021-07-28 裂缝方位约束的多维信息裂缝表征方法 Active CN113568046B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110854228.7A CN113568046B (zh) 2021-07-28 2021-07-28 裂缝方位约束的多维信息裂缝表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110854228.7A CN113568046B (zh) 2021-07-28 2021-07-28 裂缝方位约束的多维信息裂缝表征方法

Publications (2)

Publication Number Publication Date
CN113568046A true CN113568046A (zh) 2021-10-29
CN113568046B CN113568046B (zh) 2023-02-03

Family

ID=78168228

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110854228.7A Active CN113568046B (zh) 2021-07-28 2021-07-28 裂缝方位约束的多维信息裂缝表征方法

Country Status (1)

Country Link
CN (1) CN113568046B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114019565A (zh) * 2021-11-02 2022-02-08 中海石油(中国)有限公司 快速叠前裂缝预测方法、装置、计算机存储介质及设备
CN114459911A (zh) * 2022-01-24 2022-05-10 湖南继善高科技有限公司 一种油气压裂裂缝的方位解释方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040117121A1 (en) * 2002-09-27 2004-06-17 Veritas Dgc Inc. Reservoir fracture characterization
WO2011005353A1 (en) * 2009-07-06 2011-01-13 Exxonmobil Upstream Research Company Method for seismic interpretation using seismic texture attributes
CN102478662A (zh) * 2010-11-30 2012-05-30 中国石油天然气集团公司 纵波裂缝预测与剥层方法
WO2013134427A2 (en) * 2012-03-06 2013-09-12 Ion Geophysical Corporation Model predicting fracturing of shale
CN105445814A (zh) * 2015-11-09 2016-03-30 恒泰艾普石油天然气技术服务股份有限公司 海上窄方位角资料的方位各向异性属性差裂缝预测方法
CN106772601A (zh) * 2016-12-21 2017-05-31 成都理工大学 一种基于叠前地震数据的裂缝储层快速预测方法
US20170160413A1 (en) * 2015-12-04 2017-06-08 Cgg Services Sas Method and apparatus for analyzing fractures using avoaz inversion
CN106873033A (zh) * 2017-03-13 2017-06-20 中国海洋石油总公司 一种基于神经网络的多属性地震信息融合裂缝预测的方法
CN107728204A (zh) * 2016-08-11 2018-02-23 中国石油化工股份有限公司 基于叠前纵波各向异性的裂缝预测方法及系统
US20190369279A1 (en) * 2018-05-31 2019-12-05 Cgg Services Sas METHODS AND SYSTEMS PERFORMING JOINT VVAz AND AVAz INVERSION

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040117121A1 (en) * 2002-09-27 2004-06-17 Veritas Dgc Inc. Reservoir fracture characterization
WO2011005353A1 (en) * 2009-07-06 2011-01-13 Exxonmobil Upstream Research Company Method for seismic interpretation using seismic texture attributes
CN102478662A (zh) * 2010-11-30 2012-05-30 中国石油天然气集团公司 纵波裂缝预测与剥层方法
WO2013134427A2 (en) * 2012-03-06 2013-09-12 Ion Geophysical Corporation Model predicting fracturing of shale
CN105445814A (zh) * 2015-11-09 2016-03-30 恒泰艾普石油天然气技术服务股份有限公司 海上窄方位角资料的方位各向异性属性差裂缝预测方法
US20170160413A1 (en) * 2015-12-04 2017-06-08 Cgg Services Sas Method and apparatus for analyzing fractures using avoaz inversion
CN107728204A (zh) * 2016-08-11 2018-02-23 中国石油化工股份有限公司 基于叠前纵波各向异性的裂缝预测方法及系统
CN106772601A (zh) * 2016-12-21 2017-05-31 成都理工大学 一种基于叠前地震数据的裂缝储层快速预测方法
CN106873033A (zh) * 2017-03-13 2017-06-20 中国海洋石油总公司 一种基于神经网络的多属性地震信息融合裂缝预测的方法
US20190369279A1 (en) * 2018-05-31 2019-12-05 Cgg Services Sas METHODS AND SYSTEMS PERFORMING JOINT VVAz AND AVAz INVERSION

Non-Patent Citations (8)

* Cited by examiner, † Cited by third party
Title
LIXIN TIAN ET AL.: "Application of Narrow Azimuth Seismic data in Metamorphic Fractured Reservoir Prediction", 《SEG DENVER 2010 ANNUAL MEETING》 *
YAO LI ET AL.: "Research on Seismic Fluid Detection Driven by Rock Physics: A case study in Bohai Bay", 《2018 SEG INTERNATIONAL EXPOSITION AND 88TH ANNUAL MEETING》 *
YONG XIAO ET AL.: "Fracture Prediction Method with Narrow Azimuth Seismic Data in Tazhong District of Tarim Basin", 《INTERNATIONAL JOURNAL OF ADVANCED MATERIALS CHEMISTRY AND PHYSICS》 *
刘军等: "远近偏移距属性差在碳酸盐岩储层裂缝检测中的应用—以流花１１－１油田为例", 《中国海上油气》 *
明君等: "锦州南油田潜山裂缝储层地震预测技术研究", 《中国石油勘探》 *
李尧等: "渤海油田渤中 A 构造太古宙潜山裂缝储层预测", 《物探与化探》 *
李梅等: "基于海上窄方位采集地震资料的裂缝预测技术及应用", 《成都理工大学学报（ 自然科学版）》 *
田立新等: "窄方位角地震资料在裂缝储层预测中的应用", 《成都理工大学学报(自然科学版)》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114019565A (zh) * 2021-11-02 2022-02-08 中海石油(中国)有限公司 快速叠前裂缝预测方法、装置、计算机存储介质及设备
CN114019565B (zh) * 2021-11-02 2023-10-13 中海石油(中国)有限公司 快速叠前裂缝预测方法、装置、计算机存储介质及设备
CN114459911A (zh) * 2022-01-24 2022-05-10 湖南继善高科技有限公司 一种油气压裂裂缝的方位解释方法
CN114459911B (zh) * 2022-01-24 2023-06-20 湖南继善高科技有限公司 一种油气压裂裂缝的方位解释方法

Also Published As

Publication number Publication date
CN113568046B (zh) 2023-02-03

Similar Documents

Publication Publication Date Title
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
CN101907725B (zh) 裂缝预测方法和装置
WO2016041189A1 (zh) 一种评价页岩气储层及寻找甜点区的方法
CN108680951A (zh) 一种基于地震信息判断煤层气富集沉积控制作用的方法
CN104502971A (zh) 页岩储层总有机碳含量和脆性空间展布的检测方法
CN105277982B (zh) 一种泥页岩总有机碳含量地震预测方法
CN113568046B (zh) 裂缝方位约束的多维信息裂缝表征方法
CN104375182B (zh) 一种裂缝型储层流体的识别方法及装置
CN108663713B (zh) 一种建立深度域构造模型的方法
CN105445814B (zh) 海上窄方位角资料的方位各向异性属性差裂缝预测方法
CN104005760A (zh) 基于方位各向异性弹性阻抗的裂缝检测方法
CN105719433A (zh) 一种基于孔内地震波的超前预报方法
CN106324682B (zh) 一种应用于冻土层地区的表层结构调查方法
CN110703354A (zh) 一种基于井震结合的岩性边界识别方法
CN111208564A (zh) 一种深度域层位标定方法及装置
CN111308558B (zh) 页岩气水平井纵波时差校正方法
CN112946782B (zh) 一种致密油气储渗体地震精细刻画方法
CN112305602A (zh) 基于叠前多属性与古地貌融合技术的碳酸岩储层预测方法
Wang et al. Determination of the microstructure of a lithologic interface using the delayed response characteristics of horizontal well gamma ray logging curves: a case study of the Daqingzijing Oilfield, Songliao Basin, Northeast China
CN115857047B (zh) 一种地震储层综合预测方法
CN114876454A (zh) 一种大斜度井复杂岩性水淹层识别方法及系统
CN113806674B (zh) 古河道纵向尺度的量化方法、装置、电子设备及存储介质
CN112305594B (zh) 非均质储层的油气分布确定方法及系统
CN116500679A (zh) 一种井震标定方法、装置和相关设备
CN110780345A (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