一种地震属性优化方法及装置制造方法及图纸

本申请实施例公开了一种地震属性优化方法及装置。所述方法包括：基于微地震数据计算工区的第一裂缝特征；对于预设地震属性集合中的每个地震属性，基于该地震属性计算所述工区的第二裂缝特征；选取与所述第一裂缝特征相匹配的第二裂缝特征，将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。本申请实施例的方法和装置，可以提高优化结果的准确性。

【技术实现步骤摘要】

本申请涉及地球物理勘探
，特别涉及一种地震属性优化方法及装置。
技术介绍
裂缝已经称为人们认识和发现油气藏的重要手段。可靠的裂缝预测不仅能够有效的指导水平井位部署，还能够为钻井工程提供指导，优化后续压裂施工设计。地震数据携带有大量的储层地质信息。因此，目前通常基于地震数据来对裂缝进行预测。提高裂缝预测精度的关键在于提取并优选与地质特性有关的地震属性。优化选择地震属性而不是地震数据的原因在于，当工区的地质特性确定时，只有在储层或流体性质变化导致相应特征参数变化达到某一程度时，地震数据才会有明显变化，而地震属性却对这些特征参数敏感得多，并且许多地震属性是非线性的，可以增加裂缝预测的准确性。地震属性技术对裂缝预测的主要依据是储层物性和填充在储层中的流体性质的空间变化，造成地震波反射速度、振幅、频率等一系列基于几何的、运动学的、动力学的地震属性的变化。地震属性的种类繁多，与裂缝之间的关系复杂。不同的工区和/或不同的储层，对裂缝敏感的地震属性是不完全相同的。地震属性优化技术可以用来解决上述问题。地震属性优化方法可以利用人的经验或数学方法，优选出对裂缝最敏感(或最有效、最有代表性)的、个数最少的地震属性或地震属性组合，目的在于提高裂缝预测的精度，从而为认识和发现油气藏提供帮助。现有技术中，通常根据专家的知识挑选对裂缝预测最有影响的地震属性。油田专家对某个地区带有较多储层信息的特征是比较了解的，可以凭经验进行地震属性的选取。虽然专家有时可以提出几组较优的地震属性或者地震属性组合，但是很难进一步确定哪一组最优。现有技术中根据专家知识优化地震属性的方法多解性强、稳定性不够，并且选择过程主要根据个人经验，缺乏充分的理论依据，因此，亟需一种有充分理论证明的、稳定的优化地震属性的方法。
技术实现思路
本申请实施例的目的是提供一种地震属性优化方法及装置，以提高优化结果的准确性。为解决上述技术问题，本申请实施例提供一种地震属性优化方法及装置是这样实现的：一种地震属性优化方法，包括：基于微地震数据计算工区的第一裂缝特征；对于预设地震属性集合中的每个地震属性，基于该地震属性计算所述工区的第二裂缝特征；选取与所述第一裂缝特征相匹配的第二裂缝特征，将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。一种地震属性优化装置，包括：第一计算单元，用于基于微地震数据计算工区的第一裂缝特征；第二计算单元，用于对于预设地震属性集合中的每个地震属性，基于该地震属性计算所述工区的第二裂缝特征；选取单元，用于选取与所述第一裂缝特征相匹配的第二裂缝特征，将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。由以上本申请实施例提供的技术方案可见，本申请实施例可以基于微地震数据计算工区的第一裂缝特征；对于预设地震属性集合中的每个地震属性，可以基于该地震属性计算所述工区的第二裂缝特征；选取与所述第一裂缝特征相匹配的第二裂缝特征，可以将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。与现有技术相比，本申请实施例所采用的数学优化方法，比专家经验优化方法具有更充足的理论证明和严谨的推导过程，优化结果也更加准确。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为本申请地震属性优化方法一个实施例的流程图；图2为本申请利用蚂蚁体属性计算工区裂缝特征示意图；图3为本申请地震属性优化装置一个实施例的功能结构示意图。具体实施方式为了使本
的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。目标基于地震数据的裂缝预测技术主要包括横波分裂技术、多波多分量技术、多方位VSP(VerticalSeismicProfiling，垂直地震剖面)技术、各向异性AVO(Amplitudevariationwithoffset)及振幅方位各向异性等。随着地震属性技术的发展，基于地震数据的裂缝预测技术还包括了地球物理学家在叠后地震数据体上进行的一些特殊分析技术，例如，相干分析技术、谱分解技术、曲率分析技术、倾角/方位角分析技术等。上述基于地震数据的裂缝预测技术，关键在于提取并优选与地质特性有关的地震属性。地质储层可以包括常规储层和非常规储层。目前主要基于地震数据来对常规储层的裂缝进行预测。因此，对于常规储层，有必要挑选对裂缝预测最有影响的地震属性，以提高裂缝预测的精度。非常规储层中油气的赋存方式与常规储层存在很大差异。常规储层中的油气主要以游离方式赋存于储层中，而非常规储层中的油气主要以吸附方式赋存于储层中。因此，微裂缝或微裂隙会直接影响非常规储层中油气藏的富集程度。目前非常规储层的裂缝预测所采用的方法通常与常规储层相同。但非常规储层中裂缝的尺度比常规储层小的多，常规储层的裂缝预测技术，例如相干分析、边缘检测和曲率分析等，只能反映相对大尺度的裂缝发育带，对于中小尺度的裂缝带检测效果不佳。因此，在常规储层中能够取得较好预测效果的地震属性不一定适用于非常规储层。对于非常规储层，更有必要挑选对裂缝预测最有影响的地震属性，以对其裂缝进行准确地预测。在水力压裂井中，由压裂井中注水压力的变化所产生的地震数据通常称为微地震数据。水力压裂所产生的人工裂缝的展布特征与天然裂缝展布特征存在关系。水力压裂趋向于激活天然裂缝发育区，天然裂缝发育区往往出现较多的微地震事件。微地震事件的方向很大程度上也能表征天然裂缝方向。微地震监测能够很好的诊断水力压裂所产生的人工裂缝的形态。因此，在已知微地震展布的前提下，可以利用微地震事件的方向和强度作为约束条件，将微地震事件的方向和强度作为一种标定成果，从而达到求取最适合工区的地震属性。下面介绍本申请地震属性优化方法的一个实施例。请参阅图1，该实施例可以包括如下的步骤。步骤S10：基于微地震数据计算工区的第一裂缝特征。所述工区通常指施工的区域。所述裂缝特征通常可以包括裂缝的形态、方向、大小、倾角、以及强度等，其中，所述强度通常指单位面积内的裂缝数量。通常地，可以向工区的压裂井注水压裂，以形成微地震数据，所述微地震数据可以由监测井中的检波器接收到。所述监测井与所述压裂井的距离通常小于800m。所述基于微地震数据计算工区的第一裂缝特征，可以包括：基于所述微地震数据确定所述工区的震源位置；基于所述震源位置确定所述工区的天然裂缝区域；基于所述天然裂缝区域内的裂缝计算所述工区的第一裂缝特征。所述基于所述微地震数据确定所述工区的震源位置，可以包括：对所述微地震数据进行极化分析，以确定所述微地震数据中P波(P-wave，或primarywave)和S波(S-wave，或secondarywave)的传播方向；基于所述P波的传播方向，获取所述P波的初至时间；基于所述S波的传本文档来自技高网...

【技术保护点】
一种地震属性优化方法，其特征在于，包括：基于微地震数据计算工区的第一裂缝特征；对于预设地震属性集合中的每个地震属性，基于该地震属性计算所述工区的第二裂缝特征；选取与所述第一裂缝特征相匹配的第二裂缝特征，将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。

【技术特征摘要】
1.一种地震属性优化方法，其特征在于，包括：基于微地震数据计算工区的第一裂缝特征；对于预设地震属性集合中的每个地震属性，基于该地震属性计算所述工区的第二裂缝特征；选取与所述第一裂缝特征相匹配的第二裂缝特征，将选取的第二裂缝特征所对应的地震属性作为所述工区的优化地震属性。2.如权利要求1所述的方法，其特征在于，所述基于微地震数据计算工区的第一裂缝特征，包括：基于所述微地震数据确定所述工区的震源位置；基于所述震源位置确定所述工区的天然裂缝区域；基于所述天然裂缝区域内的裂缝计算所述工区的第一裂缝特征。3.如权利要求2所述的方法，其特征在于，所述基于所述微地震数据确定所述工区的震源位置，包括：对所述微地震数据进行极化分析，以确定所述微地震数据中P波和S波的传播方向；基于所述P波的传播方向，获取所述P波的初至时间；基于所述S波的传播方向，获取所述S波的初至时间；基于所述P波的初至时间、S波的初至时间、以及预设的速度模型，计算所述工区的震源位置。4.如权利要求3所述的方法，其特征在于，所述计算所述工区的震源位置，包括：根据如下的公式计算所述工区的震源位置，[(xpi-xqk)2+(ypi-yqk)2+(zpi-zqk)2]12=ΔTkivPvSvP-vS,]]>其中，xpi、ypi、以及zpi分别为检测点pi在三维直角坐标系中的x轴坐标值、y轴坐标值、以及z轴坐标值；xqk、yqk、以及zqk分别为震源位置qk在三维直角坐标系中的x轴坐标值、y轴坐标值、以及z轴坐标值；ΔTki为P波和S波之间初至时间的差值；vp为P波在所述预设速度模型中的传播速度；vS为S波在所述预设速度模型中的传播速度。5.如权利要求3所述的方法，其特征在于，所述计算所述工区的震源位置，包括：根据如下的公式计算所述工区的震源位置，[(xpi-xqk)2+(ypi-yqk)2+(zpi-zqk)2]12-[(xpl-xqk)2+(ypl-yqk)2+(zpl-zqk)2]12=vP(Tki-Tkl),]]>其中，xpi、ypi、以及zpi分别为检测点pi在三维直角坐标系中的x轴坐标值、y轴坐标值、以及z轴坐标值；xpl、ypl、以及zpl分别为检测点pl在三维直角坐标系中的x轴坐标值、y轴坐标值、以及z轴坐标值；...

【专利技术属性】
技术研发人员：刘伟，王熙明，张宇生，徐刚，万小平，张晨，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京;11