一种基于地质信息的虚拟井构建方法
【技术领域】
[0001]本发明涉及石油地球物理探测领域，特别是涉及一种基于地质信息的虚拟井构建方法。
【背景技术】
[0002]钻井(探井)可以直接获取最可靠的地质资料，但是钻井成本尚昂，在石油探区勘探程度较低的情况下探井不适合大量钻探。地下岩性信息一般是通过井约束条件下对地震数据进行波阻抗反演得到的，某些石油探区内由于钻井稀少且分布不均等情况下的地震数据反演过程中缺少井信息约束，反演结果精度较低，导致对石油探区的地质认识受到限制。建立虚拟井约束地震反演是必要的，传统的虚拟井构建技术是基于地震反演波阻抗数据体得到的波阻抗信息，利用波阻抗信息来构建虚拟井有分辨率较低和忽略细节的不足，故传统的虚拟井构建技术存在较大的误差。因此寻求一种构建高精度的虚拟井的技术对沉积盆地勘探以及提高反演精度都具有重要意义。
【发明内容】
[0003]本发明的目的是要提供一种恢复精度高、分辨率高、符合地质规律的虚拟井构建方法。
[0004]特别地，本发明提供了一种基于地质信息的虚拟井构建方法，用于在具有多个实钻井的目标区域获取所述目标区域中预定位置处的虚拟井的自然伽玛曲线，包括:
[0005]以沿深度方向的第一采样密度获取所述多个实钻井各自的自然伽玛数据和波阻抗数据、以及所述虚拟井处的波阻抗数据；
[0006]对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y = k*x+b，其中，X为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；
[0007]根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；
[0008]根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据;和
[0009]基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。
[0010]可选地，所述拟合算法为插值算法；
[0011 ]可选地，所述插值算法为克里金插值算法。
[0012]可选地，还包括:以沿所述深度方向的第二采样密度获取所述虚拟井处的预定时间长度范围内的辐射量;其中，所述第二采样密度大于所述第一采样密度；
[0013]对所述各实钻井的自然伽玛数据进行去趋势化操作并建立其与所述辐射量之间的缩放关系，再基于所述虚拟井处的辐射量以及所述缩放关系获得对应的去趋势化的第二自然伽玛数据;和
[0014]基于所述第二自然伽玛数据获得所述虚拟井的第二自然伽玛曲线。
[0015]可选地，其中，将所述第一自然伽玛曲线与所述第二自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。
[0016]可选地，还包括:根据所述虚拟井处的地质层结构信息以沿深度方向的第三采样密度获取所述虚拟井的自然伽玛数据的变化趋势，并获得反映所述变化趋势对应的第三自然伽玛曲线;其中，所述第三采样密度大于所述第一采样密度但小于第二采样密度。
[0017]可选地，其中，用与所述虚拟井最近的所述实钻井处的地质层结构信息作为所述虚拟井处的地质层结构信息。
[0018]可选地，所述地质层结构信息包括层序划分特征信息、时间地层格架信息及地震层位信息。
[0019]可选地，其中，将所述第一自然伽玛曲线与所述第三自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。
[0020]可选地，其中，将所述第一自然伽玛曲线、所述第二自然伽玛曲线与所述第三自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。
[0021]本发明的基于地质信息的虚拟井构建方法，通过在由波阻抗数据获得的第一自然伽玛曲线的基础上叠加由天文周期控制的第二自然伽玛曲线和由地质模型信息确定的第三自然伽玛曲线，从而获得虚拟井的自然伽玛曲线，这充分考虑了地震数据在横向的高连续性、实钻井数据纵向的高分辨率性和地质模型的约束性。本发明的虚拟井构建方法提高了虚拟井的分辨能力和刻画细节的能力，同时，还具有更高的精度，更加符合真实的地质情况。
[0022]根据下文结合附图对本发明具体实施例的详细描述，本领域技术人员将会更加明了本发明的上述以及其他目的、优点和特征。
【附图说明】
[0023]后文将参照附图以示例性而非限制性的方式详细描述本发明的一些具体实施例。附图中相同的附图标记标示了相同或类似的部件或部分。本领域技术人员应该理解，这些附图未必是按比例绘制的。附图中:
[0024]图1是根据本发明一个实施例的基于地质信息的虚拟井构建方法的流程图；
[0025]图2是本发明实施例的石油探区的示意图；
[0026]图3是本发明实施例的过实钻井w2和w4的地震数据的振幅剖面图；
[0027]图4是本发明实施例的过实钻井w2和w4的波阻抗数据的剖面图；
[0028]图5是本发明实施例的实钻井w2的波阻抗数据与自然伽玛数据的线性关系图；
[0029]图6是本发明实施例的实钻井w4的波阻抗数据与自然伽玛数据的线性关系图；
[0030]图7是本发明实施例的目的层段自然伽玛数据的频谱特征图；
[0031 ]图8是本发明实施例的目的层段理论辐射量曲线图；
[0032]图9是本发明实施例的自然伽玛曲线图。
【具体实施方式】
[0033]图1是根据本发明一个实施例的基于地质信息的虚拟井构建方法的流程图。其一般性地包括:
[0034]S1、以沿深度方向的第一采样密度获取多个已有实钻井的各自的自然伽玛数据和波阻抗数据、以及虚拟井处的波阻抗数据；
[0035]地震数据是反映地层波阻抗变化界面的振幅数据，可以根据已有的数据记录或者对目标区域进行勘探测量获得目标区域内的多个实钻井处及虚拟井处的地震数据，对该地震数据进行波阻抗反演，可以获得它们对应的波阻抗数据。同时，还可以沿实钻井井身测量岩层的天然伽马射线强度，从而获得各个实钻井的自然伽玛数据。
[0036]S2、对于每一实钻井，拟合出其对应的所述波阻抗数据与所述自然伽玛数据之间的线性关系y = k*x+b，其中，X为波阻抗数据，y为自然伽玛数据，从而获得每一实钻井位置处的k参数值和b参数值；
[0037]提取每一实钻井的波阻抗数据的采样深度点a1、a2、a3、……、an，将由步骤I所获得的各实钻井的自然伽玛数据光滑后获得每一实钻井的光滑曲线，在每一实钻井的光滑曲线上以采样深度点a1、a2、a3、......、an获得对应的自然伽玛数据，通过采样深度点a1、a2、
a3、……、an所对应的多组自然伽玛数据和波阻抗数据回归出每一实钻井的波阻抗数据与自然伽玛数据之间的线性关系y = k*x+b(其中，X为波阻抗数据，y为自然伽玛数据)，即求出每一实钻井位置处的k参数值和b参数值。
[0038]S3、根据各个实钻井位置处的k参数值和b参数值，通过拟合算法获得所述虚拟井处的k参数值和b参数值；
[0039]通过克里金插值算法求出各个实钻井位置与其对应的k参数值和b参数值的分布情况，拟合出其对应的曲线，从该曲线中提取虚拟井的位置，从而获得虚拟井处的k参数值和b参数值。
[0040]S4、根据所述虚拟井处的k参数值和b参数值以及波阻抗数据，基于所述线性关系获得所述虚拟井处的第一自然伽玛数据；
[0041]可以将虚拟井处的k参数值和b参数值以及波阻抗数据代入到线性关系y= k*x+b(其中，X为波阻抗数据，y为自然伽玛数据)中，从而可以求出虚拟井处的自然伽玛数据，将该自然伽玛数据设定为第一自然伽玛数据。
[0042]S5、基于所述第一自然伽玛数据获得所述虚拟井的第一自然伽玛曲线，所述第一自然伽玛曲线即为所述虚拟井的所述自然伽玛曲线。
[0043]可以通过数据拟合的方法将第一自然伽玛数据拟合回归出第一自然伽玛曲线。第一自然伽玛曲线是基于第一采样密度由虚拟井处的波阻抗数据确定的低频信息。需要说明的是，本实施例中所述第一自然伽玛数据即为虚拟井处的自然伽玛数据，所述第一自然伽玛曲线即为虚拟井的自然伽玛曲线。
[0044]在本发明的另一个实施例中，还包括:
[0045]以沿所述深度方向的第二采样密度获取所述虚拟井处的预定时间长度范围内的辐射量;其中，所述第二采样密度大于所述第一采样密度；
[0046]本实施例的前提条件是确定由天文周期引起的辐射量是否对该目标区域的自然伽玛数据产生影响。可以根据该目标区域的玮度和研究目的层段持续的时间范围得出其理论辐射量曲线。
[0047]对所述各实钻井的自然伽玛数据进行去趋势化操作并建立其与所述辐射量之间的缩放关系，再基于所述虚拟井处的辐射量以及所述缩放关系获得对应的去趋势化的第二自然伽玛数据;和
[0048]基于所述第二自然伽玛数据获得所述虚拟井的第二自然伽玛曲线。
[0049]可以通过数据拟合的方法将第二自然伽玛数据拟合回归出第二自然伽玛曲线。第二自然伽玛曲线是基于第二采样密度由天文周期辐射量变化得出的高频信息。
[0050]需要说明的是，将所述第一自然伽玛曲线与所述第二自然伽玛曲线相叠加以获得所述虚拟井的所述自然伽玛曲线。
[0051]在本发明的另一个实施例中，还包括:
[0052]根据所述虚拟井处的地质层结构信息以沿深度方向的第三采样密度获取所述虚拟井的自然伽玛数据的变化趋势；
[0053]地质层结构信息反映了在一定地质条件下发育的沉积岩分布、演化及特征的综合，因此，在相同的地质层结构中岩性特征及演化规律具有一定的相似性。所以，本实施例采用与所述虚拟井最邻近的所述实钻井处的地质层结构信息作为所述虚拟井处的地质层结构信息。具体地，所述地质层结构信息包括层序划分特征信息、时