大斜度井纵、横波测井速度各向异性同时校正方法技术

本发明专利技术公开了一种大斜度井纵、横波测井速度各向异性同时校正方法，该方法是，对含有泥质，裂缝发育，层状分布地层中的斜井中所测得的纵、横波速度同时进行校正，消除斜井纵、横波速度所受地层中泥质固有各向异性、裂缝各向异性以及地层层状各向异性的影响，将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。本发明专利技术可将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度，可将具有一定倾角的大斜度井中测得的纵、横波速度精确地校正为不受泥质固有各向异性，裂缝，地层层状各向异性影响的垂直于地层的纵、横波速度。

【技术实现步骤摘要】

本专利技术涉及油气勘探
，尤其涉及一种。
技术介绍
测井纵、横波速度是叠前反演、AVO分析等地震储层预测方法的基础资料。测井中所测得的速度是沿井筒方向的速度，而在地震中所得到的速度是垂直于地层的速度。当井筒垂直于地层时，两者一致；当井筒相对于地层有一定的倾角时，且地层为各向异性介质，即不同方向的速度不同时，在井中测得的速度不再是垂直于地层的速度，这样就与地震速度有一定的差异，从而对两者的联合研究，包括井震标定，子波估算，AVO分析，以及最终的叠前反演都有较大的影响。地层中的各向异性主要有泥质固有各向异性，以及地层层状分布的各向异性。Hornby (2003)提出了一种校正斜井纵波速度的办法，使用泥质的各向异性参数，基于实测纵波测井速度反演斜井测井速度对应的垂直纵波速度，但是他只考虑了泥质各向异性，没有考虑地层层状的影响，且仅仅是对纵波进行校正，因此，无法将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。所述AV0，英文为Amplitude Versus Offset,中文为振幅随偏移距变化,根据AVO分析,可以对油气直接进行检测。
技术实现思路
本专利技术要解决的技术问题是提供一种，将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。为了解决上述技术问题，本专利技术的是，对含有泥质，裂缝发育，层状分布地层中的斜井中所测得的纵、横波速度同时进行校正，消除斜井纵、横波速度所受地层中泥质固有各向异性、裂缝各向异性以及地层层状各向异性的影响，将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。具体步骤如下:步骤1:根据研究区域岩心测量数据，取泥质的固有各向异性参数，纵横波速度和密度；认为砂岩与灰岩为各向同性介质，各向异性参数为O,根据地区经验取砂岩,灰岩的纵横波速度和密度；根据复杂岩性测井解释程序，对测井曲线进行解释，得到页岩，砂岩，灰岩的相对体积含量；并根据电阻率解释得到裂缝密度；步骤2:在步骤I所得参数及岩性体积含量的基础之上，根据巴卡斯平均方程计算混合矿物的各向异性参数，根据哈德森各向异性岩石物理模型计算含裂缝岩石的各向异性参数；所述巴卡斯平均方程指Backus平均方程；所述哈德森各向异性岩石物理模型指Hudson各向异性岩石物理模型；步骤3:根据步骤2求得的各向异性参数，根据汤姆森相速度公式求取井筒角度下每一个采样点处的相速度，进而计算每一个采样点处的群速度；所述汤姆森相速度公式指Thomsen相速度公式；步骤4:使用步骤3计算的群速度与实测的速度构建反演目标函数，即速度校正目标函数，并迭代求解方程，得到垂直方向的纵波速度Vp (O)，横波速度Vs (O)。步骤2中，混合矿物的各向异性参数是在考虑泥质固有各向异性的基础之上，考虑了地层的层状分布所造成的各向异性以及裂缝造成的各向异性，基于每种矿物的各向异性参数和体积含量，由巴卡斯平均公式计算得到，并且经哈德森各向异性岩石物理模型计算了裂缝所造成的各向异性值。步骤4中的速度反演目标函数使用群速度公式与计算的群速度构成，而不是直接使用相速度公式，反演同时得到垂直方向的纵波速度和横波速度。以上三点使得该方法可以在已知井筒倾角、实测纵、横波的情况下，计算得到井筒垂直情况下的垂直速度。本专利技术的与现有技术相比具有以下有益效果。1、本技术方案由于采用了对含有泥质，裂缝发育，层状分布地层中的斜井中所测得的纵、横波速度同时进行校正，消除由于井筒倾角以及地层中泥质固有各向异性、裂缝各向异性以及地层层状各向异性对纵、横波速度造成的影响的技术手段，所以，可将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。2、本技术方案由于采用了考虑了混合矿物的各向异性参数是在考虑泥质固有各向异性的基础之上，考虑了地层的层状分布所造成的各向异性以及裂缝造成的各向异性，基于每种矿物的各向异性参数和体积含量，由巴卡斯平均公式计算得到，并且经哈德森各向异性岩石物理模型计算了裂缝所造成的各向异性值的技术手段，所以，可将具有一定倾角的大斜度井中测得的纵、横波速度精确地校正为不受泥质固有各向异性，裂缝，地层层状各向异性影响的垂直于地层的纵、横波速度。【附图说明】图1是大斜度井示意图。图2是在各向影响下的纵、横波速度受井斜角度影响示意图。图3是校正前后纵波速度、横波速度值随倾角变化图。图4为校正前后大斜度井信息及纵横波速度变化图。图5校正前后小斜度井信息及纵横波速度变化图。图6为校正前与校正后大斜度井与小斜度井纵横波慢度数据直方图。【具体实施方式】图1是大斜度井示意图。图2是在各向影响下的纵、横波速度受井斜角度影响示意图。各向异性相速度及由于井筒角度的变化，所造成的速度数据向理论各向异性相速度靠拢，而与垂直方向所能测得到的各向同性速度差距变大；图2中数据点代表具有一定井筒倾角Θ所对应的实测速度，其中纵波速度及对应倾角数据来自于Hornby等于2003年发表于Geophysics杂志第 68 卷第 2 期的 464-471 页的 Anisotropy correction for deviated-well soniclogs:Application to seismic well tie—文；横波数据是由Hornby的纵波速度数据，根据Castagna的Vp-Vs关系式计算得到，所述Castagna的Vp-Vs关系式指Castagna等在1985年发表于Geophy杂志第50卷第571-581页的Relationships betweencompressional wave and shear wave velocities in clastic silicate rocks 一文，公式为:Vp=1360+l.16*Vs，单位为 m/s。图3是校正前后纵波速度、横波速度值随倾角变化图，其中VpM，VsM为原始纵波速度和横波速度，VpC, VsC为校正后的纵波速度和横波速度，校正后，纵波速度和横波速度明显减小，与图2中所示的不受方向影响的各向同性速度所示规律一致。图4为校正前后大斜度井信息及纵横波速度变化图，图中自左往右分别为泥质含量(Vshale)，密度(P )，倾角(Dip)，井轨迹方位(Azimuth)，纵波速度(Vp)，校正前后纵波速度变化(Vp_Res)，横波速度(Vs).校正前后横波速度变化(Vs_Res)，其中，倾角显式，该井井筒倾角从5200米开始增大，在5500-井底之间，其角度基本在60°到80°之间，纵横波速度道中，黑色曲线为校正前数据，红色曲线为校正后数据，在井斜度小的区域，由于井的垂直，速度值为垂直速度，所以校正后速度变化不大，在井斜度较大的层段，校正后的速度比原始速度明显减小，与5100-5300米之间受角度影响较小的速度值靠拢，并且与图1所示各向同性与各向异性速度之间关系所示规律一致。图5为校正前后小斜度井信息及纵横波速度变化图，图中自左往右分别为泥质含量(Vshale)，密度(P )，倾角(Dip)，井轨迹方位(Azimuth)，纵波速度(Vp)，校正前后纵波速度变化(Vp_Res)，横波速度(Vs).校正前后横波速度变化(Vs_Res)，纵横波速度道中，黑本文档来自技高网...

【技术保护点】
一种大斜度井纵、横波测井速度各向异性同时校正方法，其特征在于：对含有泥质，裂缝发育，层状分布地层中的斜井中所测得的纵、横波速度同时进行校正，消除斜井纵、横波速度所受地层中泥质固有各向异性、裂缝各向异性以及地层层状各向异性的影响，将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。

【技术特征摘要】
1.一种大斜度井纵、横波测井速度各向异性同时校正方法，其特征在于:对含有泥质，裂缝发育，层状分布地层中的斜井中所测得的纵、横波速度同时进行校正，消除斜井纵、横波速度所受地层中泥质固有各向异性、裂缝各向异性以及地层层状各向异性的影响，将具有一定倾角的大斜度井中测得的纵、横波速度校正为垂直于地层的纵、横波速度。2.根据权利要求1所述的方法，其特征在于:包括如下步骤: 步骤1:根据研究区域岩心测量数据，取泥质的固有各向异性参数，纵横波速度和密度；认为砂岩与灰岩为各向同性介质，各向异性参数为O，根据地区经验取砂岩，灰岩的纵横波速度和密度；根据复杂岩性测井解释程序，对测井曲线进行解释，得到页岩，砂岩，灰岩的相对体积含量；并根据电阻率解释得到裂缝密度； 步骤2:在步骤I所得参数及岩性体积含量的基础之上，根据巴卡斯平均方程计算混合矿物的各向异性参数，根据哈德森各向异性岩...

【专利技术属性】
技术研发人员：刘致水，孙赞东，田军，
申请(专利权)人：孙赞东，
类型：发明
国别省市：北京;11