一种利用松弛因子叠前地震反演进行储层预测的方法

摘要

本发明提供了一种利用松弛因子叠前地震反演进行储层预测的方法，该方法包括：A、收集地震工区内的叠前地震角道集数据以及地震工区内相应的测井数据；B、根据先验地质层位以及所述的测井数据生成纵波速度、横波速度和密度的初始模型；C、将所述的叠前地震角道集数据叠加成不同角度的n个叠加地震道集；D、根据所述的n个叠加地震道集提取对应的n个子波；E、逐次从n个叠加地震道集及其对应的n个子波中提取一个叠加地震道集及其对应的1个子波等。通过将由纵波速度、横波速度和密度生成的松弛因子加入到反演过程中，进行储层的预测，提高了储层预测的精确性，对油气有利位置的指示更为明显。

主权项

1.一种利用松弛因子叠前地震反演进行储层预测的方法，其特征是，所述的方法包括：A、收集地震工区内的叠前地震角道集数据以及地震工区内相应的测井数据；B、根据先验地质层位以及所述的测井数据生成纵波速度、横波速度和密度的初始模型；C、将所述的叠前地震角道集数据叠加成不同角度的n个叠加地震道集；D、根据所述的n个叠加地震道集提取对应的n个子波；E、逐次从n个叠加地震道集及其对应的n个子波中提取一个叠加地震道集及其对应的1个子波；F、根据纵波速度、横波速度和密度的初始模型确定纵波反射系数；G、将所述的纵波反射系数与提取的1个叠加地震道集对应的子波褶积得到合成地震道集；H、确定所述的合成地震道集与叠加地震道集的误差；I、判断所述的误差是否小于预定的阈值，判断为否时，执行步骤J，否则执行步骤K；J、根据所述的误差修改纵波速度、横波速度和密度的初始模型，返回执行步骤F；K、根据所述的误差确定优化的纵波速度、横波速度和密度；L、判断当前的叠加地震道集是否为n个中的最后一个叠加地震道集，判断为否时，返回执行步骤E，否则执行步骤M；M、对n个优化的纵波速度、横波速度和密度求均值得到最优化的纵波速度、横波速度和密度；N、对所述最优化的纵波速度、横波速度和密度进行综合解释，生成当前区域的储层预测结果；其中，步骤B包括：根据先验地质层位以及所述的测井数据进行井震标定，将深度域测井的声波时差曲线标定为时间域以确定储层的目的层段；在所述的目的层段进行地震数据层位解释，拾取地质层位对应的地震层位；根据所述的地震层位对所述的测井数据进行插值以生成纵波速度、横波速度和密度的初始模型；步骤F包括：根据纵波速度、横波速度的初始模型生成松弛因子；根据纵波速度、横波速度和密度的初始模型以及松弛因子确定出纵波反射系数；其中，根据纵波速度、横波速度的初始模型生成松弛因子通过如下公式进行：$S_k=\left(\right|(\frac{{\sec }^2{i}_k}{2}-\overline{V_{\mathrm{ak}}})|/\overline{V_{\mathrm{ak}}}|(4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k-\overline{V_{\mathrm{ak}}})|/\overline{V_{\mathrm{ak}}}|(\frac{1}{2}(4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k-{\tan }^2{i}_k)-\overline{V_{\mathrm{ak}}})|/\stackrel{}{V_{\mathrm{ak}}})$其中，v_P,v_S分别为纵波速度、横波速度，v_S/v_P为横纵波速度比，k为叠加地震道集数，i为叠加地震道集对应的入射角度，通过如下公式进行：$\overline{V_{\mathrm{ak}}}=\left(\right|\frac{{\sec }^2{i}_k}{2}|+|4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k|+|\frac{1}{2}(4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k-{\tan }^2{i}_k)\left|\right)/3$根据纵波速度、横波速度和密度的初始模型以及松弛因子确定出纵波反射系数通过如下公式进行：$R_{\mathrm{PP}}\left(i_k\right)\approx (\frac{{\sec }^2{i}_k}{2}(\frac{\Delta v_P}{v_P}+\frac{\mathrm{\Delta \rho }}{\rho })-4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k(\frac{\Delta v_S}{v_S}+\frac{\mathrm{\Delta \rho }}{\rho })\frac{1}{2}(4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k-{\tan }^2{i}_k)\frac{\mathrm{\Delta \rho }}{\rho })·S_k^T·\overline{V_k},---\left(1\right)$其中，Δv_P/v_P,Δv_S/v_S,Δρ/ρ分别为纵波速度、横波速度和密度的相对变化量，通过如下公式进行：$\overline{V_{\mathrm{gk}}}={|\frac{{\sec }^2{i}_k}{2}·4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k·\frac{1}{2}(4{\left(\frac{v_S}{v_P}\right)}^2{\sin }^2{i}_k-{\tan }^2{i}_k)|}^{\frac{1}{3}}.$