古构造图成图方法及装置制造方法
【专利摘要】本发明公开了一种古构造图成图方法及装置。该古构造图成图方法包括以下步骤：构造变速成图；井校得到深度误差；地震叠前三参数同步反演得到拉梅常数和密度的乘积；提取地震平均振幅；构建涉及深度误差、地震平均振幅、拉梅常数和密度的乘积的超定方程组：求出超定方程组的最小二乘解；根据公式计算未修正点的目的层的深度误差；生成古构造图。该古构造图成图方法及装置，通过结合井校数据和多个地震属性，估算了井校点的深度误差和地震属性之间的响应关系，进而利用地震属性对于无井区域的深度进行修正，从而能够在井控数量较少的情况下生成更准确可靠的古构造图。
【专利说明】古构造图成图方法及装置
【技术领域】
[0001]本发明涉及一种古构造图成图方法及装置。
【背景技术】
[0002]古构造图是指导勘探开发活动的重要基础性图件，是构造描述的重要内容，同时也是寻找构造圈闭、落实油气层含烃面积的主要依据。随着油气藏勘探开发的不断深入，对该类图件的精度及准确性提出了更高的要求，尤其是对于砂体微幅构造的描述的精确性，能够直接影响后期剩余油挖潜方向的决策。
[0003]目前较为成熟的古构造图生成技术中的一种为构造变速成图技术，其已广泛应用于Geoframe、Landmark等地震解释系统的相应构造成图模块之中。这一方法的进步性在于，相对于传统的构造常速成图技术，能够对由速度变化引起这部分深度估计误差得到一定校正。但是，精确的速度场建立主要依靠测井资料的点控，在实际运用中该方法往往受到井控数量的制约，导致应用效果达不到最佳效果，容易出现构造图在少井区深度失真的情况。然而这种少井区域往往是勘探或开发挖潜井的部署有利区域。因此，迫切需求一种能够生成更为准确的古构造图的方法。
[0004]通过进一步剖析构造图的深度误差产生原因可知，从根本上说导致岩石在地下产生差异压实与地震波传播速度的变化的根本原因应当是目的层的地质属性，即目的层岩石的岩性及物理性质，尤其是这些性质的横向非均质性。但在不均匀井网的客观条件限制下，构造图无法得到充分的校正，容易产生深度的失真，对于地质属性的获取也极为困难，通常情况下都只能依靠部分位置的井校数据进行校准，这里所说的井校数据即井点处的地质分层结果。现有技术中有多种处理方法均可利用少量井校数据来对更大区域的构造图进行修正，但均局限于从纯粹的数据处理的角度进行修正，其中的一个典型例子即采用插值处理。然而，现有的这些方法生成的古构造图，其准确性仍然不够高，并且在井校数据较少且地质条件的横向分布复杂的情况下其准确性不高的缺陷尤为明显。
【发明内容】
[0005]本发明要解决的技术问题是为了克服现有技术中的古构造图成图技术在井控数量较少的情况下生成的古构造图的不够准确的缺陷，提出一种古构造图成图方法及装置。
[0006]本发明是通过下述技术方案来解决上述技术问题的:
[0007]本发明提供了一种古构造图成图方法，其特点在于，包括以下步骤:
[0008]S1、采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
[0009]S2、获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
[0010]S3、根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式如=I； -2Z 2计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、P为密度、Zp为纵波阻抗、Zs为横波阻抗；
[0011]S4、根据目的层开取相应时窗，提取各点的地震平均振幅；
[0012]S5、构建以下超定方程组:
【权利要求】
1.一种古构造图成图方法，其特征在于，包括以下步骤: S1、采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度； S2、获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差； S3、根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式/lp = Zp2 -2Z计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、P为密度、Zp为纵波阻抗、Zs为横波阻抗； S4、根据目的层开取相应时窗，提取各点的地震平均振幅； S5、构建以下超定方程组:
2.一种古构造图成图装置，其特征在于，包括: 变速成图模块，用于采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度； 深度误差计算模块，用于获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差； 同步反演模块，用于根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式如=-2ZS2计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、P为密度、ΖΡ为纵波阻抗、Zs为横波阻抗； 振幅提取模块，用于根据目的层开取相应时窗，提取各点的地震平均振幅； 方程组构建模块，用于构建以下超定方程组:
【文档编号】G01V1/36GK103885083SQ201410043533
【公开日】2014年6月25日 申请日期:2014年1月29日 优先权日:2014年1月29日
【发明者】李 一 申请人:上海石油天然气有限公司