古构造图成图方法及装置
阅读：936发布：2020-05-15
专利汇可以提供古构造图成图方法及装置专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种古构造图成图方法及装置。该古构造图成图方法包括以下步骤：构造变速成图；井校得到深度误差； 地震 叠前三参数同步反演得到 拉梅常数 和 密度 的乘积；提取地震平均振幅；构建涉及深度误差、地震平均振幅、拉梅常数和密度的乘积的超定方程组：求出超定方程组的最小二乘解；根据公式计算未修正点的目的层的深度误差；生成古构造图。该古构造图成图方法及装置，通过结合井校数据和多个地震属性，估算了井校点的深度误差和地震属性之间的响应关系，进而利用地震属性对于无井区域的深度进行修正，从而能够在井控数量较少的情况下生成更准确可靠的古构造图。，下面是古构造图成图方法及装置专利的具体信息内容。
1.一种古构造图成图方法，其特征在于，包括以下步骤：
S1、采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
S2、获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
S3、根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗；
S4、根据目的层开取相应时窗，提取各点的地震平均振幅；
S5、构建以下超定方程组：
上述超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
S6、求出该超定方程组的最小二乘解c0、c1、c2；
S7、根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
2.一种古构造图成图装置，其特征在于，包括：
变速成图模块，用于采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
深度误差计算模块，用于获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
同步反演模块，用于根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗；
振幅提取模块，用于根据目的层开取相应时窗，提取各点的地震平均振幅；
方程组构建模块，用于构建以下超定方程组：
其中，该超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
方程组求解模块，用于求取该超定方程组的最小二乘解c0、c1、c2；
修正计算模块，用于根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
古构造图成图方法及装置
技术领域
[0001] 本发明涉及一种古构造图成图方法及装置。
背景技术
[0002] 古构造图是指导勘探开发活动的重要基础性图件，是构造描述的重要内容，同时也是寻找构造圈闭、落实油气层含烃面积的主要依据。随着油气藏勘探开发的不断深入，对该类图件的精度及准确性提出了更高的要求，尤其是对于砂体微幅构造的描述的精确性，能够直接影响后期剩余油挖潜方向的决策。
[0003] 目前较为成熟的古构造图生成技术中的一种为构造变速成图技术，其已广泛应用于Geoframe、Landmark等地震解释系统的相应构造成图模块之中。这一方法的进步性在于，相对于传统的构造常速成图技术，能够对由速度变化引起这部分深度估计误差得到一定校正。但是，精确的速度场建立主要依靠测井资料的点控，在实际运用中该方法往往受到井控数量的制约，导致应用效果达不到最佳效果，容易出现构造图在少井区深度失真的情况。然而这种少井区域往往是勘探或开发挖潜井的部署有利区域。因此，迫切需求一种能够生成更为准确的古构造图的方法。
[0004] 通过进一步剖析构造图的深度误差产生原因可知，从根本上说导致岩石在地下产生差异压实与地震波传播速度的变化的根本原因应当是目的层的地质属性，即目的层岩石的岩性及物理性质，尤其是这些性质的横向非均质性。但在不均匀井网的客观条件限制下，构造图无法得到充分的校正，容易产生深度的失真，对于地质属性的获取也极为困难，通常情况下都只能依靠部分位置的井校数据进行校准，这里所说的井校数据即井点处的地质分层结果。现有技术中有多种处理方法均可利用少量井校数据来对更大区域的构造图进行修正，但均局限于从纯粹的数据处理的角度进行修正，其中的一个典型例子即采用插值处理。然而，现有的这些方法生成的古构造图，其准确性仍然不够高，并且在井校数据较少且地质条件的横向分布复杂的情况下其准确性不高的缺陷尤为明显。
发明内容
[0005] 本发明要解决的技术问题是为了克服现有技术中的古构造图成图技术在井控数量较少的情况下生成的古构造图的不够准确的缺陷，提出一种古构造图成图方法及装置。
[0006] 本发明是通过下述技术方案来解决上述技术问题的：
[0007] 本发明提供了一种古构造图成图方法，其特点在于，包括以下步骤：
[0008] S1、采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
[0009] S2、获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
[0010] S3、根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗；
[0011] S4、根据目的层开取相应时窗，提取各点的地震平均振幅；
[0012] S5、构建以下超定方程组：
[0013]
[0014] 上述超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
[0015] S6、求出该超定方程组的最小二乘解c0、c1、c2；
[0016] S7、根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
[0017] S8、生成该区域的一古构造图，该古构造图中该多个点的目的层的深度为该多个点的目的层的井校深度，各个未修正点的目的层的深度为该原始构造图中各个未修正点的目的层的深度和S7计算得到的各个未修正点的目的层的深度误差之和。
[0018] 上述步骤S1中采用的构造变速成图技术为一种现有的构造成图方法，其原理大致为，应用沿层相干层速度反演方法计算层速度，然后使用图偏移的方法对时间构造图进行时深转换后形成深度构造图。
[0019] 更具体地，上述步骤S1包括了以下步骤：首先，对该区域内的各口实际钻井的测井资料，依据其声波时差与密度测井曲线求得纵波反射系数曲线，并与统计地震子波进行卷积生成合成地震记录，之后与实际地震资料的波形进行对比，多次进行时深关系调整，直至两者的相关系数达到最大值即实际地震与成地震记录的相似性达到最大；接着，利用上述标定过程获得的井控速度场对地震叠加速度场进行校正得到该区域的全区速度场；再根据井震标定结果，针对目标层位所属的地震反射进行地震解释，得到时间平面，利用地震解释成果进一步优化、微调井点处的地质分层结果，确定该目的层在钻井位置的真实深度；最后根据上述得到的时间平面和全区速度场的数据得到深度构造图，即该原始构造图。
[0020] 本领域技术人员应当理解，在步骤S2中的该多个点均为可以直接利用钻井实际测得的该目的层的深度进行校准的点，井校深度即该多个点经过钻井实测数据进行校准后得到的深度值。该多个点也可以理解为是井控区域或者井校区域内的点。由于该多个点可以利用钻井实际测得的深度值进行校准，可以认为这部分数据的准确性极高，并显著高于该原始构造图中的无法根据钻井实测数据进行直接校准的点的准确性。因此，在步骤S2中计算得到的深度误差近似于该原始构造图中的该多个点的深度数值和真实数值之间的误差。
[0021] 本领域技术人员也应当理解，地震反演是利用经过一系列处理过的地表（或海上）地震波观测资料，以地质规律和钻井、测井资料为约束，对地下岩层的物理性质进行成像的过程。而上述步骤S3中所采用的地震叠前三参数同步反演可以采用现有的叠前反演算法来实现。叠前反演的基本原理是当平面纵波以较大的入射角倾斜射入两种介质分界面的情况下，会产生反射纵波、反射横波、透射纵波和透射横波，应用弹性波动方程理论可以求得上述波的反射系数与界面上下介质的速度、密度与反射角的关系。因此，在纵波采集条件下，也可以得到横波的数据。其中的三维地震资料即为地震采集资料经过包括偏移成像与叠加等一系列处理后的地震振幅数据体，本发明中对三维地震资料的采集及处理采用常规技术手段。
[0022] 本发明中采用的地震叠前三参数同步反演的理论基础是Zoeppritz方程组，Zoeppritz方程组给出了反射系数、纵、横波速度、密度与入射角之间的关系。上述步骤S3中所采用的地震叠前三参数同步反演可以采用包括Hampson-Russel在内的地震解释软件依据Zoeppritz近似公式建立方程，利用叠前地震资料与测井纵、横波速度与密度，同步反演三个基本岩石弹性参数纵波阻抗、横波阻抗与密度。这里的叠前地震资料指经过限角度叠加的叠加前地震角度道集。在获取了纵波阻抗、横波阻抗与密度后，根据公式计算得到拉梅常数。本发明中的拉梅常数即拉梅第一参数，其可以表征岩石的压缩性，其对于流体性质尤为敏感。
[0023] 在步骤S4中根据目的层开取相应时窗，提取各点的地震平均振幅。应当理解的是，包括上述的三维地震资料在内的地震数据体均是时间域的，所以提取地震属性时需要针对相应的目的层来开取时窗，从而提取出目的层相对应的数据。也就是说，步骤S4中是根据开取的时窗内各地震采样的数据提取出平均振幅。举例来说，地震数据体在平面上某一点可能有n个采样对应于目的层，求取这些采样的地震振幅值的平均值即可得到这一点的地震平均振幅。
[0024] 步骤S5中构建的超定方程组以及S6的求解过程是基于对于深度误差对于拉梅常数和密度的乘积以及地震平均振幅的响应在一定程度上接近于线性响应，而对于线性响应的相关关系的近似求解可以采用求解超定方程组的最小二乘解的方式来进行。
[0025] 由步骤S5及S6可以看出，本发明和现有的古构造图的成图方法或者修正方法相比，最为显著的区别在于，本发明中并不是将井校数据简单地作为数值修正的基础，而是根据井校数据相应的区域以及这些区域的部分地震属性来估算变速成图得到的原始构造图和目的层的真实深度之间的误差如何响应地震属性。本领域技术人员应当理解，本发明中的井校数据可以理解为井点处的地质分层结果。这一区别使得本发明能够更好的修正由于地质属性的横向非均质性所带来的误差，即尽可能消除由于目的层岩石的岩性及物理性质而产生的误差。而在步骤S6求得的最小二乘解，即S7中采用的公式中的系数，因此利用这一公式来修正深度误差，进而生成的该古构造图的准确性不仅远远高于该原始构造图，也大大准确于现有技术中采用的插值法等数学分析方法对于构造图的修正。
[0026] 容易理解地，在该区域内的所有点中，该多个点可以利用井校数据进行修正，因此无需在步骤S7及S8中采用这一公式进行修正。除该多个点外的其他点均为未修正点，因而需要采用这一公式进行修正。
[0027] 本发明还提供了一种古构造图成图装置，其特点在于，包括：
[0028] 变速成图模块，用于采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
[0029] 深度误差计算模块，用于获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
[0030] 同步反演模块，用于根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗；
[0031] 振幅提取模块，用于根据目的层开取相应时窗，提取各点的地震平均振幅；
[0032] 方程组构建模块，用于构建以下超定方程组：
[0033]
[0034] 其中，该超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
[0035] 方程组求解模块，用于求取该超定方程组的最小二乘解c0、c1、c2；
[0036] 修正计算模块，用于根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
[0037] 成图模块，用于生成该区域的一古构造图，该古构造图中该多个点的目的层的深度为该多个点的目的层的井校深度，各个未修正点的目的层的深度为该原始构造图中各个未修正点的目的层的深度和修正计算模块计算得到的各个未修正点的目的层的深度误差之和。
[0038] 在符合本领域常识的基础上，上述各优选条件，可任意组合，即得本发明各较佳实例。
[0039] 本发明的积极进步效果在于：
[0040] 本发明的古构造图成图方法及装置，通过结合井点处的地质分层结果和多个地震属性，估算了井校点的深度误差和地震属性之间的响应关系，进而利用地震属性对于无井区域的深度进行修正，从而能够在井控数量较少的情况下生成更准确可靠的古构造图。附图说明
[0042] 图2a为本发明的应用实例中采用常规变速成图技术得到的构造图。
[0043] 图2b为本发明的应用实例中采用常规变速成图技术结合井校数据及常规数值处理得到的构造图。
[0044] 图2c为本发明的应用实例中采用实施例1的古构造图成图方法得到的构造图。
[0045] 图3a为本发明的应用实例中井校区域的深度误差分布图。
[0046] 图3b为本发明的应用实例中拉梅常数和密度的乘积的平面分布图。
[0047] 图3c为本发明的应用实例中地震平均振幅的平面分布图。
[0048] 图4为本发明实施例2的古构造图成图装置的示意图。
具体实施方式
[0049] 下面结合附图给出本发明较佳实施例，以详细说明本发明的技术方案，但并不因此将本发明限制在所述的实施例范围之中。
[0050] 实施例1
[0051] 如图1所示，本实施例的古构造图成图方法包括以下步骤：
[0052] S1、采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度；
[0053] S2、获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差；
[0054] S3、根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗；
[0055] S4、根据目的层开取相应时窗，提取各点的地震平均振幅；
[0056] S5、构建以下超定方程组：
[0057]
[0058] 上述超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
[0059] S6、求出该超定方程组的最小二乘解c0、c1、c2；
[0060] S7、根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
[0061] S8、生成该区域的一古构造图，该古构造图中该多个点的目的层的深度为该多个点的目的层的井校深度，各个未修正点的目的层的深度为该原始构造图中各个未修正点的目的层的深度和S7计算得到的各个未修正点的目的层的深度误差之和。
[0062] 具体来说，上述步骤S1包括以下步骤：首先，对该区域内的各口实际钻井的测井资料，依据其声波时差与密度测井曲线求得纵波反射系数曲线，并与统计地震子波进行卷积生成合成地震记录，之后与实际地震资料的波形进行对比，多次进行时深关系调整，直至两者的相关系数达到最大值即实际地震与成地震记录的相似性达到最大；接着，利用上述标定过程获得的井控速度场对地震叠加速度场进行校正得到该区域的全区速度场；再根据井震标定结果，针对目标层位所属的地震反射进行地震解释，得到时间平面，利用地震解释成果进一步优化、微调井点处的地质分层结果，确定该目的层在钻井位置的真实深度；最后根据上述得到的时间平面和全区速度场的数据得到深度构造图，即该原始构造图。
[0063] 其中，三维地震资料即为地震采集资料经过包括偏移成像与叠加等一系列处理后的地震振幅数据体。本实施例中采用的地震叠前三参数同步反演的理论基础是Zoeppritz方程组，Zoeppritz方程组给出了反射系数、纵、横波速度、密度与入射角之间的关系。上述步骤S3中所采用的地震叠前三参数同步反演采用Hampson-Russel地震解释软件依据Zoeppritz近似公式建立方程，利用叠前地震资料与测井纵、横波速度与密度，同步反演三个基本岩石弹性参数纵波阻抗、横波阻抗与密度。这里的叠前地震资料指经过限角度叠加的叠加前地震角度道集。在获取了纵波阻抗、横波阻抗与密度后，根据公式计算得到拉梅常数。本发明中的拉梅常数即拉梅第一参数，其可以表征岩石的压缩性，其对于流体性质尤为敏感。
[0064] 在步骤S4中根据目的层开取相应时窗，提取各点的地震平均振幅。应当理解的是，包括上述的三维地震资料在内的地震数据体均是时间域的，所以提取地震属性时需要针对相应的目的层来开取时窗，从而提取出目的层相对应的数据。也就是说，步骤S4中是根据开取的时窗内各地震采样的数据提取出平均振幅。举例来说，地震数据体在平面上某一点可能有n个采样对应于目的层，求取这些采样的地震振幅值的平均值即可得到这一点的地震平均振幅。
[0065] 在本实施例中，对于变速成图合成地震记录的生成、速度场的生成、通过地震解释得到时间平面、以及原始构造图的生成都采取了网格化处理的手段，也就是说，其中采用的网格均是相同的。并且，地震反演与属性提取结果的所属平面位置仍以三维地震资料地震道所在位置所决定，其间距为采集处理间距，与上述成图所用的网格可能是不一致的，但本实施例中仍然将其通过数学方法处理为相同的网格来表示。如此，步骤S3中通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，同样也采用相同的网格来表示，S4中提取的各点地震平均振幅也采用相同的网格来表示。
[0066] 本领域技术人员应当理解，本实施例中的上述网格化处理方法，只是为了实现数据处理中的一致性，提高处理效率。
[0067] 步骤S5中构建的超定方程组以及S6的求解过程是基于对于深度误差对于拉梅常数和密度的乘积以及地震平均振幅的响应在一定程度上接近于线性响应，而对于线性响应的相关关系的近似求解，本实施例采用求解超定方程组的最小二乘解的方式进行。然后根据求得的最小二乘解对深度误差进行修正。
[0068] 以下将对本实施例的古构造图成图方法的一个应用实例进行说明，这一应用实例涉及东海某油气产层砂体顶板构造图的制作成图。
[0069] 如图2a所示为本应用实例中采用常规的变速成图技术生成的构造图，即执行步骤S1后得到的原始构造图。在利用井校区域的井校深度，采用常规的数值处理方法校准后得到的构造图如图2b所示。图2a、图2b和图2c中的矩形所包围的区域即具有多个测井从而能够获取多个井校深度的井校区域，通常也称为井控区域。本领域技术人员应当理解，图2a、图2b和图2c的构造图中以连线表示平面中的等深度线（类似于等高线），在等深度线处标记的数字2670、2700、2730等均表示相应等深度线所对应的深度，因而这样的构造图能够全面形象地反应整个成图区域的平面深度分布。
[0070] 本应用实例中，在执行了步骤S1和S2后得到的井校区域的深度误差分布如图3a所示。进一步，在执行了步骤S3后得到的各点的目的层的拉梅常数和密度的乘积，其平面分布如图3c所示。执行了步骤S4后得到的各点的地震平均振幅，其平面分布如图3b所示。图3b和图3c中的曲线分别用于表示地震平均振幅以及拉梅常数和密度的乘积这两个参数在平面分布的等值线，图3b中的数值-1000、-2000等、图3c中的数值20、28、36等分别用于表示数值对应位置的等值线的取值。其中，地震平均振幅是由地震资料经过一系列处理后得到的，本应用实例中的地震平均振幅为无因次量，用于表示地震反射的相对强弱，是没有单位的物理量。本应用实例中，拉梅常数和密度的乘积的单位为GPa*g/cc。其中9
1GPa=1000MPa=1000000KPa=10Pa，GPa为压强单位，1GPa也就是1000兆帕。g/cc为密度单
3
位，1g/cc=1g/ml=1g/cm，即克每立方厘米。应当理解的是，本应用实例中的单位选取是和附图3b、3c相对应的，并不因此限制本发明的古构造成图方法必须采用上述单位。
1GPa=1000MPa=1000000KPa=10Pa，GPa为压强单位，1GPa也就是1000兆帕。g/cc为密度单
3
位，1g/cc=1g/ml=1g/cm，即克每立方厘米。应当理解的是，本应用实例中的单位选取是和附图3b、3c相对应的，并不因此限制本发明的古构造成图方法必须采用上述单位。
[0071] 根据井校区域的井校深度，以及图3b、3c所示的分布图得到相应区域（即井校深度相对应的平面点）的地震平均振幅以及拉梅常数和密度的乘积，并据此执行S5构建得到以下超定方程组：
[0072]
[0073] 其中，地震平均振幅在测量中的每个地震道为一个1ms或者2ms采样的波列图，其中含有有多个波峰波谷，上述超定方程组中地震平均振幅为负值即对应于波谷。
[0074] 执行S6求出上述超定方程组的最小二乘解c0=7.08；c1=0.0019；c2=0.122。然后执行S7根据公式ΔDm＝7.08+0.0019×Am+0.122×λmρm计算除井校深度相对应的平面点以外的未修正点的目的层的深度误差ΔDm，其中下标m表示遍历所有的未修正点。
[0075] 得出上述结果后执行S8生成该区域的一古构造图，其中各个未修正点的目的层的深度为该原始构造图中各个未修正点的目的层的深度和S7计算得到的各个未修正点的目的层的深度误差之和，本应用实例中生成的该古构造图如图2c所示。
[0076] 本应用实例中，以常规变速成图技术、常规变速成图技术结合井校数据和数值修正方法、以及本实施例的古构造图成图方法分别得到的构造图依次如图2a、图2b、图2c所示。并且，在图2a、图2b、图2c中的西北侧实际有一口测井，在三张图中均以PH7标志出，这一测井处于矩形的井校区域外，由图2a、图2b、图2c中的等深线分布就能够从图中获取图中任意位置处的深度，当然也就能够获取测井PH7处的深度。
[0077] 由图2a可以得到测井PH7所在点的深度为2730.06米，由图2b可以得到测井PH7所在点的深度为2728.38米（应当注意图2a和图2b中均没有深度为2740米或以上的等深线，深度最深的等深线为2730米），而由图2c得到的测井PH7所在点的深度为2733.44米（图2c中测井PH7的位置位于2730米等深线和2740米等深线之间）。
[0078] 本应用实例中测井PH7是实际存在的，由测井PH7实测得到的深度为2734.40米。也就是说，现有的变速成图方法以及利用数值方法校准得到的深度值的误差分别为4.34米和6.02米，而本发明的古构造图成图方法则将这一误差减小到了0.96米，相比于现有的成图方法减小误差达到70%以上，成图的准确性大幅提高。
[0079] 实施例2
[0080] 如图4所示，本实施例的古构造图成图装置包括变速成图模块1、深度误差计算模块2、同步反演模块3、振幅提取模块4、方程组构建模块5、方程组求解模块6、修正计算模块7、成图模块8。
[0081] 其中，变速成图模块1用于采用构造变速成图技术生成一区域的原始构造图，该原始构造图记录有该区域内各点的目的层的深度。深度误差计算模块2，用于获取该区域内的多个点的目的层的井校深度，并计算该多个点的目的层的深度误差，其中深度误差为井校深度和深度之差。
[0082] 同步反演模块3，用于根据三维地震资料，通过地震叠前三参数同步反演求得该区域内各点的目的层的纵波阻抗、横波阻抗与密度，然后根据公式 计算得到各点的目的层的拉梅常数和密度的乘积，公式中λ为拉梅常数、ρ为密度、Zp为纵波阻抗、Zs为横波阻抗。振幅提取模块4，用于根据目的层开取相应时窗，提取各点的地震平均振幅。
[0083] 方程组构建模块5，用于构建以下超定方程组：
[0084]
[0085] 其中，该超定方程组中方程的个数等于该多个点的个数n，c0、c1、c2为三个未知量，ΔDk、Ak、λkρk分别为该多个点中各点的目的层的深度误差、地震平均振幅、拉梅常数和密度的乘积，其中1≤k≤n；
[0086] 方程组求解模块6，用于求取该超定方程组的最小二乘解c0、c1、c2。修正计算模块7，用于根据公式ΔDm＝c0+c1Am+c2λmρm计算该区域内除该多个点外的未修正点的目的层的深度误差，其中m遍历满足n
[0087] 成图模块8，用于生成该区域的一古构造图，该古构造图中该多个点的目的层的深度为该多个点的目的层的井校深度，各个未修正点的目的层的深度为该原始构造图中各个未修正点的目的层的深度和修正计算模块计算得到的各个未修正点的目的层的深度误差之和。
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种预测煤层气甜点的方法 | 2020-05-13 | 794 |
| 一种裂缝型储层流体的识别方法及装置 | 2020-05-14 | 794 |
| 一种地震地层裂隙裂缝密度反演方法及系统 | 2020-05-20 | 201 |
| 建立铸造奥氏体不锈钢等轴晶声学特性计算模型的方法 | 2020-05-17 | 203 |
| 磁共振弹性成像重建方法及成像系统 | 2020-05-15 | 179 |
| 古构造图成图方法及装置 | 2020-05-15 | 936 |
| 确定、展示和使用地下岩石弹性模量和密度相对变化的方法 | 2020-05-12 | 163 |
| 一种阻尼材料黏弹性动态力学参数的测算方法 | 2020-05-18 | 309 |
| 基于非规则扇形网格的起伏隧道空间正演模拟系统及方法 | 2020-05-25 | 1005 |
| 一种直接探测地层中游离态气体的方法 | 2020-05-11 | 490 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
分析报告
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
QQ群二维码
意见反馈
意见反馈