一种虚拟井的自然伽马曲线重构方法技术

本发明专利技术提供了一种虚拟井的自然伽马曲线重构方法，用于在具有多个实钻井的目标区域获取虚拟井的自然伽玛曲线，包括：获取多个实钻井的自然伽玛数据和波阻抗数据以及虚拟井处的波阻抗数据；拟合出每一实钻井对应的波阻抗数据与自然伽玛数据之间的线性关系y＝k*x+b，其中，x为波阻抗数据，y为自然伽玛数据)，从而获得每一实钻井位置的k参数值和b参数值；通过拟合算法获得虚拟井处的k参数值和b参数值；根据虚拟井处的k参数值和b参数值以及波阻抗数据，获得虚拟井处的自然伽玛数据；和获得虚拟井的自然伽玛曲线。本发明专利技术的虚拟井构建方法提高了虚拟井的分辨能力和刻画细节的能力，同时，还具有更高的精度，更加符合真实的地质情况。

A natural gamma curve reconstruction method for virtual wells

The invention provides a natural gamma curve reconstruction method for a virtual well, which is used to obtain the natural gamma curve of a virtual well in a target area with multiple real drilling, including obtaining the natural gamma data and wave impedance data of multiple real drilling and the wave impedance data at the virtual well, and fitting the corresponding waves for each drilling well. The linear relationship between the impedance data and the natural gamma data is y = k*x+b, in which the X is the wave impedance data and the Y is the natural gamma data. The k parameter and the b parameter of each real drilling position are obtained, and the k parameter and the b parameter value of the virtual well are obtained by the fitting algorithm; the k parameter value and the b parameter value of the virtual well are based on the virtual well. And the wave impedance data, the natural gamma data in the virtual well are obtained, and the natural gamma curve of the virtual well is obtained. The virtual well construction method of the invention improves the ability of resolving and depicting the details of the virtual well. At the same time, it also has higher precision and is more in line with the real geological condition.

【技术实现步骤摘要】
一种虚拟井的自然伽马曲线重构方法本申请是申请日为2015年12月01日、申请号为201510868288.9、题为“一种基于地质信息的虚拟井构建方法”的专利技术专利申请的分案申请。
本专利技术涉及石油地球物理探测领域，特别是涉及一种虚拟井的自然伽马曲线重构方法。
技术介绍
钻井(探井)可以直接获取最可靠的地质资料，但是钻井成本高昂，在石油探区勘探程度较低的情况下探井不适合大量钻探。地下岩性信息一般是通过井约束条件下对地震数据进行波阻抗反演得到的，某些石油探区内由于钻井稀少且分布不均等情况下的地震数据反演过程中缺少井信息约束，反演结果精度较低，导致对石油探区的地质认识受到限制。建立虚拟井约束地震反演是必要的，传统的虚拟井构建技术是基于地震反演波阻抗数据体得到的波阻抗信息，利用波阻抗信息来构建虚拟井有分辨率较低和忽略细节的不足，故传统的虚拟井构建技术存在较大的误差。因此寻求一种构建高精度的虚拟井的技术对沉积盆地勘探以及提高反演精度都具有重要意义。
技术实现思路
本专利技术的目的是要提供一种恢复精度高、分辨率高、符合地质规律的虚拟井构建方法。特别地，本专利技术提供了一种虚拟井的自然伽马曲线重构方法，用于在具有多个实钻井的目标区域获取所述目标区域中预定位置处的虚拟井的自然伽玛曲线，包括：以沿深度方向的第一采样密度获取所述多个实钻井各自的自然伽玛数据和波阻抗数据、以及所述虚拟井处的波阻抗数据；对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y＝k*x+b，其中，x为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据；和基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。可选地，所述拟合算法为插值算法；可选地，所述插值算法为克里金插值算法。可选地，还包括：以沿所述深度方向的第二采样密度获取所述虚拟井处的预定时间长度范围内的辐射量；其中，所述第二采样密度大于所述第一采样密度；对所述多个实钻井各自的自然伽玛数据进行去趋势化操作并建立其与所述辐射量之间的缩放关系，再基于所述虚拟井处的辐射量以及所述缩放关系获得对应的去趋势化的第二自然伽玛数据；和基于所述第二自然伽玛数据获得所述虚拟井的第二自然伽玛曲线。可选地，其中，将所述第一自然伽玛曲线与所述第二自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。可选地，还包括：根据所述虚拟井处的地质层结构信息以沿深度方向的第三采样密度获取所述虚拟井的自然伽玛数据的变化趋势，并获得反映所述变化趋势对应的第三自然伽玛曲线；其中，所述第三采样密度大于所述第一采样密度但小于第二采样密度。可选地，其中，用与所述虚拟井最近的所述实钻井处的地质层结构信息作为所述虚拟井处的地质层结构信息。可选地，所述地质层结构信息包括层序划分特征信息、时间地层格架信息及地震层位信息。可选地，其中，将所述第一自然伽玛曲线与所述第三自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。可选地，其中，将所述第一自然伽玛曲线、所述第二自然伽玛曲线与所述第三自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。本专利技术的虚拟井的自然伽马曲线重构方法，通过在由波阻抗数据获得的第一自然伽玛曲线的基础上叠加由天文周期控制的第二自然伽玛曲线和由地质模型信息确定的第三自然伽玛曲线，从而获得虚拟井的自然伽玛曲线，这充分考虑了地震数据在横向的高连续性、实钻井数据纵向的高分辨率性和地质模型的约束性。本专利技术的自然伽马曲线重构方法提高了虚拟井的分辨能力和刻画细节的能力，同时，还具有更高的精度，更加符合真实的地质情况。根据下文结合附图对本专利技术具体实施例的详细描述，本领域技术人员将会更加明了本专利技术的上述以及其他目的、优点和特征。附图说明后文将参照附图以示例性而非限制性的方式详细描述本专利技术的一些具体实施例。附图中相同的附图标记标示了相同或类似的部件或部分。本领域技术人员应该理解，这些附图未必是按比例绘制的。附图中：图1是根据本专利技术一个实施例的虚拟井的自然伽马曲线重构方法的流程图；图2是本专利技术实施例的石油探区的示意图；图3是本专利技术实施例的过实钻井w2和w4的地震数据的振幅剖面图；图4是本专利技术实施例的过实钻井w2和w4的波阻抗数据的剖面图；图5是本专利技术实施例的实钻井w2的波阻抗数据与自然伽玛数据的线性关系图；图6是本专利技术实施例的实钻井w4的波阻抗数据与自然伽玛数据的线性关系图；图7是本专利技术实施例的目的层段自然伽玛数据的频谱特征图；图8是本专利技术实施例的目的层段理论辐射量曲线图；图9是本专利技术实施例的自然伽玛曲线图。具体实施方式图1是根据本专利技术一个实施例的虚拟井的自然伽马曲线重构方法的流程图。其一般性地包括：S1、以沿深度方向的第一采样密度获取多个已有实钻井的各自的自然伽玛数据和波阻抗数据、以及虚拟井处的波阻抗数据，所述多个实钻井各自的自然伽马数据通过沿实钻井的井身测量岩层的天然伽马射线强度获得；地震数据是反映地层波阻抗变化界面的振幅数据，可以根据已有的数据记录或者对目标区域进行勘探测量获得目标区域内的多个实钻井处及虚拟井处的地震数据，对该地震数据进行波阻抗反演，可以获得它们对应的波阻抗数据。同时，还可以沿实钻井井身测量岩层的天然伽马射线强度，从而获得各个实钻井的自然伽玛数据。S2、对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y＝k*x+b，其中，x为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；提取每一实钻井的波阻抗数据的采样深度点a1、a2、a3、……、an，将由步骤1所获得的各实钻井的自然伽玛数据光滑后获得每一实钻井的光滑曲线，在每一实钻井的光滑曲线上以采样深度点a1、a2、a3、……、an获得对应的自然伽玛数据，通过采样深度点a1、a2、a3、……、an所对应的多组自然伽玛数据和波阻抗数据回归出每一实钻井的波阻抗数据与自然伽玛数据之间的线性关系y＝k*x+b(其中，x为波阻抗数据，y为自然伽玛数据)，即求出每一实钻井位置处的k参数值和b参数值。S3、根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；通过克里金插值算法求出各个实钻井位置与其对应的k参数值和b参数值的分布情况，拟合出其对应的曲线，从该曲线中提取虚拟井的位置，从而获得虚拟井处的k参数值和b参数值。S4、根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据；可以将虚拟井处的k参数值和b参数值以及波阻抗数据代入到线性关系y＝k*x+b(其中，x为波阻抗数据，y为自然伽玛数据)中，从而可以求出虚拟井处的自然伽玛数据，将该自然伽玛数据设定为第一自然伽玛数据。S5、基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。可以通过数据拟合的方法将第一自然伽玛数据拟合回归出第一自然伽玛曲线。第一自然伽玛本文档来自技高网...

【技术保护点】
1.一种虚拟井的自然伽马曲线重构方法，用于在具有多个实钻井的目标区域获取所述目标区域中预定位置处的虚拟井的自然伽玛曲线，包括：以沿深度方向的第一采样密度获取所述多个实钻井各自的自然伽玛数据和波阻抗数据、以及所述虚拟井处的波阻抗数据，所述多个实钻井各自的自然伽马数据通过沿实钻井的井身测量岩层的天然伽马射线强度获得；对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y＝k*x+b，其中，x为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据；和基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。

【技术特征摘要】
1.一种虚拟井的自然伽马曲线重构方法，用于在具有多个实钻井的目标区域获取所述目标区域中预定位置处的虚拟井的自然伽玛曲线，包括：以沿深度方向的第一采样密度获取所述多个实钻井各自的自然伽玛数据和波阻抗数据、以及所述虚拟井处的波阻抗数据，所述多个实钻井各自的自然伽马数据通过沿实钻井的井身测量岩层的天然伽马射线强度获得；对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y＝k*x+b，其中，x为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据；和基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。2.根据权利要求1所述的自然伽马曲线重构方法，其特征在于，所述拟合算法为插值算法。3.根据权利要求2所述的自然伽马曲线重构方法，其特征在于，所述插值算法为克里金插值算法。4.根据权利要求1-3中任一项所述的自然伽马曲线重构方法，其特征在于，还包括：以沿所述深度方向的第二采样密度获取所述虚拟井处的预定时间长度范围内的辐射量；其中，所述第二采样密度大于所述第一采样密度；对所述多个实钻井各自的自然伽玛数据进行去趋势化操作...

【专利技术属性】
技术研发人员：张元福，张海波，张高阳，姜在兴，王敏，
申请(专利权)人：中国地质大学北京，
类型：发明
国别省市：北京,11