一种地震振幅类属性勘探中补偿差异压实的方法
【专利摘要】本发明公开了一种地震振幅类属性勘探中补偿差异压实的方法，该方法的步骤包括：对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释，沿地震解释层提取地震属性平面图；统计样品点坐标与样品点之间的深度差异数据；按x，y大地坐标与所提取坐标寻求对应关系；得到在一个地震属性数值随深度增加而增加的关系式；运用得到的关系式对导出的地震属性数据体进行二次运算，得到补偿压实作用之后的平面地震属性矢量图。本发明通过这种补偿压实作用对地震属性工作影响的方法，最大限度的体现了地层由于岩性变化所引起阻抗值的变化，排除由于构造压实作用影响岩层密度从而引起的反演结果不准确，从而能真实反映地层发育特征与分布范围。
【专利说明】一种地震振幅类属性勘探中补偿差异压实的方法
【技术领域】
[0001]本发明属于地球物理【技术领域】，尤其涉及一种地震振幅类属性勘探中补偿差异压实的方法。
【背景技术】
[0002]根据目前国内外已知的补偿压实的方法，整体上说只是一个通过缩小研究工区面积降低构造压实对岩层密度的影响，而目前这类方法并不能真正起到真实反映波阻抗值由岩性变化而变化的趋势，存在相对较大的误差。如果只是传统的勘探方法，不仅投资大，而且延误工期，甚至有时因场地无法布置钻探而影响勘探精度。因此单一地质勘察方法已不能满足需求，需要应用新的勘察方法，以求的更高的经济效益。
[0003]地震属性是指那些由叠前或叠后地震数据经过数学变换而导出的有关地震波的几何形态、运动学特征和统计学特征的特殊测量值。根据不同的算法、不同的属性具有不同的地质意义，有些属性值可能擅长于揭示不易探测到的岩性变化，而有些属性则可以直接用于烃类检测。
[0004]地震属性的研究和使用始于20世纪70年代，最早主要是以振幅为基础的瞬时属性，用来直接指示油气与地下不同岩层的展布情况。70年代早期，Nigel Anstey发现了含气砂岩波阻抗的异常变化，使用了反射波振幅变化特征-“亮点”与“暗点”对含气砂岩储集体进行预测。80年代，地震属性的数量迅速增力卩，多元属性分析技术首次被使用。80年代晚期，多维属性如倾角和方位角的初步发展，导致了三维连续属性在90年代的广泛应用。究其原因，它的成功得益于清楚和明确的地质含义。90年代以来，基于精细储层描述以及三维地震数据体解释的需要，地震属性分析技术更是急剧发展。现今，地震属性技术已广泛应用于储层预测、油气藏动态监测、油气藏特征描述等领域，并取得了很好的效果。
[0005]综合地震属性的数学、物理及地质意义，也可以将地震属性分类两类:第一类是基本属性，包括反射时间、振幅及相位，他们各自具有独立的物理意义和地质意义；第二类是延伸属性，如对多种属性进行处理(叠合、求差，相关、求主成分等)得到的复合属性以及根据单一地震属性再提取出的新的级联属性。其中振幅类属性主要是用于体现地下不同岩性地层的分布特征与展布规律。然而传统的勘探技术不易探测到岩性的变化，且不能直接用于烃类检测，对于寻找剩余油、提高采收率成为油田开发工作的重点及难点，也是稳产的关键。
【发明内容】
[0006]本发明的目的在于提供一种地震振幅类属性勘探中补偿差异压实的方法，旨在解决传统的勘探方法，不仅投资大，而且延误工期，甚至有时因场地无法布置钻探而影响勘探精度的问题。
[0007]本发明是这样实现的，一种地震振幅类属性勘探中补偿差异压实的方法，该方法的步骤包括:[0008]步骤一、对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释，沿地震解释层提取地震属性平面图；
[0009]根据测井、录井数据完成井上层位标定；
[0010]根据井上层位标定结合地震反射特征完成地震解释；
[0011]根据地震数据体进行属性体提取的二次运算；
[0012]根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片，得到沿地震解释层位地震属性平面图。
[0013]步骤二、在对研究区进行岩心观察、录井解释的基础上进行小层对比，挑选相同岩石类型，相同地层沉积厚度但埋深不同的样品点，统计样品点坐标与样品点之间的深度差异数据；
[0014]根据对取心井观察与数据统计，挑选相同岩石类型与沉积厚度但不同埋深的样品
占.[0015]根据所述样品点选择，统计其大地坐标和样品点埋深数据；
[0016]根据所述埋深数值二次运算得到埋深数据之间的差值。
[0017]步骤三、将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的地震属性数值，按X，y大地坐标与步骤二中所提取样品点坐标寻求对应关系；
[0018]将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的包括大地坐标的地震属性数值；
[0019]根据所述地震属性数值中的大地坐标和井位坐标寻求相同坐标的样品点,将这部分样品点作为有效样品点挑选出来；
[0020]根据所挑选的有效样品点，计算其地震属性数值差异。
[0021]步骤四、根据工区研究对象精度要求，设定补偿压实作用差值最小计算单元，然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析，得到在一个地震属性数值随深度增加而增加的关系式；
[0022]根据研究要求，设定补偿差异压实作用最小计算单元；
[0023]根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析，得到一个关于两组差值数据之间的关系式:Y=aX+b，其中Y为地震属性差异数值，X为深度差异数值，a、b为根据交汇分析所得常量；
[0024]步骤五、运用得到的关系式对导出的地震属性数据体进行二次运算，对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面地震属性矢量图。
[0025]进一步，沿地震解释层提取地震属性平面图中的反演工作要求较为精确的地震层位以及断层解释，然后合理的合成记录制作以及时深转换；
[0026]进一步，统计样品点坐标与样品点之间的深度差异数据中数据提取工作要求合理的小层对比方案，尽可能多的提取数据点，保证样品点的代表性与非偶然性；
[0027]进一步，沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值中寻求数据对应的工作要求阻抗数值的大地坐标与井点坐标两者尽可能的接近甚至重合，要求平面上距离最短，需考虑方位角和偏移距等可能引起误差的各种因素。方位角和偏移距有误差使研究区沿地震解释层的波阻抗数值不精确，且导致阻抗数值的大地坐标与井点坐标两者出现偏差，不能达到平面上距离最短的要求；
[0028]进一步，为得到在一个地震属性数值随深度增加而增加的关系式需根据研究目标精度要求设定数据单元，且数据单元应该不小于研究目标要求精度；
[0029]进一步，为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰，尽可能的反映由于岩性变化而引起的属性变化空间展布和变化规律。
[0030]本发明提供的地震振幅类属性勘探中补偿差异压实的方法，该实验方法真实性高，通过这种补偿压实作用对地震属性工作影响的方法，最大限度的体现了地层由于岩性变化所引起阻抗值的变化，排除由于构造压实作用影响岩层密度从而引起的反演结果不准确，从而能真实反映地层发育特征与分布范围。
[0031]本发明提供的地震振幅类属性勘探中补偿差异压实的方法，该实验工作效率高，与目前国内外主流的研究区分块，分段地震解释后提取属性的方法相比(目前国内一般是20平方公里，美国一般是20平方英里)，现有方法在针对较大工区整体工作时具有无与伦比的效率优势和技术优势。首先，整个工区集中在一起研究能够从总体上把握整个研究区的构造与沉积背景，在地质上有一个整体的认识，从而避免了从局部无法把握全局的场面；其次，在地震属性提取的过程中整个工区一起进行能够最大限度的避免对地震体分块而导致的数据丢失，保证数据的完整性与精度；再次，集中属性提取能够在曲线归一化的基础上尽可能多的运用井资料和地震资料结合，相互影响，控制因素较多，最大限度的保证了研究工作的质量和准确性；最后，整体属性分析能够避免最后全区综合评价时期的一个拼接工作，提高工作效率。
【专利附图】
【附图说明】
[0032]图1是本发明提供的地震振幅类属性勘探中补偿差异压实的方法。
【具体实施方式】
[0033]本发明是这样实现的，一种地震振幅类属性勘探中补偿差异压实的方法，该方法的步骤包括:
[0034]在步骤SlOl中:对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释，沿地震解释层提取地震属性平面图；
[0035]根据测井、录井数据完成井上层位标定；
[0036]根据井上层位标定结合地震反射特征完成地震解释；
[0037]根据地震数据体进行属性体提取的二次运算；
[0038]根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片，得到沿地震解释层位地震属性平面图。
[0039]在步骤S102中:在对研究区进行岩心观察、录井解释的基础上进行小层对比，挑选相同岩石类型，相同地层沉积厚度但埋深不同的样品点，统计样品点坐标与样品点之间的深度差异数据；
[0040]根据对取心井观察与数据统计，挑选相同岩石类型与沉积厚度但不同埋深的样品占.[0041]根据所述样品点选择，统计其大地坐标和样品点埋深数据；
[0042]根据所述埋深数值二次运算得到埋深数据之间的差值。
[0043]在步骤S103中:将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的地震属性数值，按X，I大地坐标与步骤二中所提取样品点坐标寻求对应关系;
[0044]将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的包括大地坐标的地震属性数值；
[0045]根据所述地震属性数值中的大地坐标和井位坐标寻求相同坐标的样品点，将这部分样品点作为有效样品点挑选出来；
[0046]根据所挑选的有效样品点，计算其地震属性数值差异。
[0047]在步骤S104中:根据工区研究对象精度要求，设定补偿压实作用差值最小计算单元，然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析，得到在一个地震属性数值随深度增加而增加的关系式；
[0048]根据研究要求，设定补偿差异压实作用最小计算单元；
[0049]根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析，得到一个关于两组差值数据之间的关系式:Y=aX+b，其中Y为地震属性差异数值，X为深度差异数值，a、b为根据交汇分析所得常量；
[0050]在步骤S105中:运用得到的关系式对导出的地震属性数据体进行二次运算，对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面地震属性
矢量图。
[0051]进一步，沿地震解释层提取地震属性平面图中的反演工作要求较为精确的地震层位以及断层解释，然后合理的合成记录制作以及时深转换；
[0052]进一步，统计样品点坐标与样品点之间的深度差异数据中数据提取工作要求合理的小层对比方案，尽可能多的提取数据点，保证样品点的代表性与非偶然性；
[0053]进一步，沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值中寻求数据对应的工作要求阻抗数值的大地坐标与井点坐标两者尽可能的接近甚至重合，要求平面上距离最短，需考虑方位角和偏移距等可能引起误差的各种因素。方位角和偏移距有误差使研究区沿地震解释层的波阻抗数值不精确，且导致阻抗数值的大地坐标与井点坐标两者出现偏差，不能达到平面上距离最短的要求。
[0054]进一步，为得到在一个地震属性数值随深度增加而增加的关系式需根据研究目标精度要求设定数据单元，且数据单元应该不小于研究目标要求精度。
[0055]进一步，为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰，尽可能的反映由于岩性变化而引起的属性变化空间展布和变化规律。
[0056]本发明提供的地震振幅类属性勘探中补偿差异压实的方法，该方法包括:对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释，沿地震解释层提取地震属性平面图；在对研究区进行岩心观察、录井解释的基础上进行小层对比，挑选相同岩石类型，相同地层沉积厚度但埋深不同的样品点，统计样品点坐标与样品点之间的深度差异数据；将沿地震解释层位提取的属性图按数值格式导出得到一组研究区沿地震解释层的波阻抗数值，按X，y大地坐标与所提取坐标寻求对应关系；根据工区研究对象精度要求，设定补偿压实作用差值单元，然后将所对应好的数据进行交汇分析，得到在一个地震属性数值随深度增加而增加的关系式；运用得到的关系式对导出的地震属性数据体进行二次运算，对全区不同深度点的地震属性数据乘以不同系数之后得出新的属性值，然后重新平面成图，从而得到补偿压实作用之后的平面地震属性矢量图。
[0057]以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
【权利要求】
1.一种地震振幅类属性勘探中补偿差异压实的方法，其特征在于，该方法的步骤包括: 步骤一、对研究区利用现有地震体segy数据、测井数据与录井数据完成地震解释，沿地震解释层提取地震属性平面图； 根据测井、录井数据完成井上层位标定； 根据井上层位标定结合地震反射特征完成地震解释； 根据地震数据体进行属性体提取的二次运算； 根据地震解释层位对二次运算后的属性数据体进行沿地震解释层切片，得到沿地震解释层位地震属性平面图； 步骤二、在对研究区进行岩心观察、录井解释的基础上进行小层对比，挑选相同岩石类型，相同地层沉积厚度但埋深不同的样品点，统计样品点坐标与样品点之间的深度差异数据； 根据对取心井观察与数据统计，挑选相同岩石类型与沉积厚度但不同埋深的样品点； 根据所述样品点选择，统计其大地坐标和样品点埋深数据； 根据所述埋深数值二次运算得到埋深数据之间的差值； 步骤三、将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的地震属性数值，按X，y大地坐标与步骤二中所提取样品点坐标寻求对应关系； 将沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的包括大地坐标的地震属性数值； 根据所述地震属性数值中的大地 坐标和井位坐标寻求相同坐标的样品点,将这部分样品点作为有效样品点挑选出来； 根据所挑选的有效样品点，计算其地震属性数值差异； 步骤四、根据工区研究对象精度要求，设定补偿压实作用差值最小计算单元，然后将所对应好的数据将深度差异数据与地震属性差异数据进行交汇分析，得到在一个地震属性数值随深度增加而增加的关系式； 根据研究要求，设定补偿差异压实作用最小计算单元； 根据步骤三所挑选的有效样品点中的深度差异数据与地震属性差异数据交汇分析，得到一个关于两组差值数据之间的关系式:Y=aX+b，其中Y为地震属性差异数值，X为深度差异数值，a、b为根据交汇分析所得常量； 步骤五、运用得到的关系式对导出的地震属性数据体进行二次运算，对全区不同深度点的地震属性数据根据步骤四所得出的关系式乘以不同系数之后得出新的属性值，然后重新平面成图，从而得到补偿压实作用之后的沿地震解释层位的平面地震属性矢量图。
2.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法，其特征在于，沿地震解释层提取地震属性平面图，该步工作要求地震层位以及断层解释，然后合成记录制作以及时深转换公式。
3.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法，其特征在于，统计样品点坐标与样品点之间的深度差异数据，该步数据提取工作要求合理的小层对比方案；提取样品的时候需要在不同岩性区域，不同构造单元中提取样品点。
4.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法，其特征在于，沿地震解释层位提取的属性图按数值格式导出，得到一组研究区沿地震解释层的地震属性数值，然后将该组数值中的大地坐标寻求与权利要求1步骤2中的样品大地坐标进行对应，该步工作要求地震属性数值的大地坐标与样品坐标两者尽可能的接近甚至重合，要求平面上距离最短，需考虑样品坐标中方位角和偏移距能引起误差的各种因素。
5.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法，其特征在于，为得到在一个地震振幅类属性数值随深度增加而增加的关系式，需根据研究目标精度要求设定数据单元，且数据单元应该不小于研究目标要求精度。
6.如权利要求1所述的地震振幅类属性勘探中补偿差异压实的方法，其特征在于，为得到补偿压实作用之后的平面地震属性矢量图必须最大限度的排除由于压实作用对地震属性值的干扰，尽可能的 反映由于岩性变化而引起的振幅类属性变化的空间展布和变化规律。
【文档编号】G01V1/28GK103487830SQ201310421916
【公开日】2014年1月1日 申请日期:2013年9月17日 优先权日:2013年9月17日
【发明者】邢翔, 胡望水, 谢锐杰, 雷中英, 牛世忠, 王志坤, 孙立勋, 李相明, 李涛, 王炜, 胡志鹏 申请人:长江大学