本发明属于油气勘探，具体涉及一种缝洞型碳酸盐岩储集体孔隙体积定量计算方法。
背景技术：
1、缝洞型碳酸盐岩油藏储层中，储集体是由沉积碳酸盐岩压实作用残余的孔隙、后期断裂及溶蚀产生的洞、缝组成，储集体发育具强非均质性，由于钻井成本高，井数量有限，且许多井需酸压才能沟通储层见产，仅靠单井资料去估算储量几乎不可能。地震具有很好的空间采样密度，通过地震振幅来刻画出洞穴的平面分布范围，但无法准确描述储层的顶面和底面，也无法用地震准确计算地下储集体的体积；常规方式是利用地震反演利用井资料插值形成低频模型开展工作，由于井数据的缺乏及井点数据对地层数据的代表性不足，造成误差很大，无法实现有效的储层预测，碳酸盐岩油藏的开发方案设计和井位设计遇到较大困难。
技术实现思路
1、为了解决现有技术中存在的上述技术问题，本发明提供一种缝洞型碳酸盐岩储集体孔隙体积定量计算方法，通过地震正演和反演结合，相互标定，可以准确刻画储集体的分布范围、顶底面位置和孔隙度大小，并计算储层孔隙空间，为储量计算奠定基础。
2、为实现上述目的，本发明的技术方案如下：
3、一种缝洞型碳酸盐岩储集体孔隙体积定量计算方法，包括：
4、s1、利用已钻井资料开展地震子波提取；
5、s2、利用地震子波和实际地震资料开展第一轮反演；
6、s3、建立地震模型，并开展正演，并对得到的正演资料进行处理；
7、s4、对于步骤s3中的正演资料进行反演，得到校正系数；
8、s5、结合步骤s2的反演结果以及校正系数，得到迭代后的低频模型；
9、s6、迭代反演得到孔隙度体解释储集体体积。
10、进一步地，步骤s1中，选择测井井径良好的井段的声波和密度资料，计算纵波阻抗及地层的反射系数，开展井震标定，结合井旁地震道提取地震平均子波以及地震资料的频带范围。
11、更进一步地，所述纵波阻抗zi为：
12、zi＝ρi×vi
13、所述反射系数ri为：
14、
15、其中，zi为第i层的纵波阻抗，ρi为第i层的岩石密度，vi为第i层的岩石纵波速度，ri为第i层和i+1层之间界面的反射系数。
16、进一步地，步骤s2中，利用实际地震资料开展稀疏脉冲反演，消除子波调谐效应，得到带限波阻抗数据体；然后开展缝洞体解释，得到主要储集体的空间包络顶面、底面位置以及平均储集体厚度。
17、进一步地，步骤s3中，利用实际测井资料中经过岩心标定的孔隙度和阻抗曲线，进行交会图分析，拟合建立孔隙度和纵波阻抗变换函数，基于平均储集体厚度建立从小到大的不同孔隙度变化地质模型，利用变换函数计算得到各孔隙度模型对应的纵波阻抗值，填入地质模型，利用实际提取子波开展波动方程正演。
18、更进一步地，步骤s3中对正演资料进行处理包括：对采集得到的正演地震单炮进行去噪和偏移处理，得到偏移后的地震全叠加数据。
19、进一步地，步骤s4中，利用地震子波和正演得到的地震资料，开展第一轮正演数据的约束稀疏脉冲反演，消除子波调谐效应，得到正演模型的带限波阻抗数据体，对比正演模型的反演带限波阻抗体及正演模型的实际波阻抗值，得到储层波阻抗校正系数及储集体的纵波阻抗截止值。
20、更进一步地，步骤s5中，利用步骤s4中得到的纵波阻抗截止值，将步骤s2中的带限反演体的纵波阻抗提取出来，并应用步骤s4中得到的校正系数进行校正，之后加上纯碳酸盐岩岩石纵波阻抗，得到反演的纵波阻抗低频模型。
21、进一步地，步骤s6中，利用步骤s1中的地震子波和步骤s5中的低频模型，开展第二轮稀疏脉冲反演，得到全频的纵波阻抗数据体；然后利用步骤s3中的岩石物理模型得到的孔隙度和纵波阻抗变换表，将纵波阻抗体转换为孔隙度体；对孔隙度体进行解释、雕刻，得到地下精确的孔洞顶底面位置及孔隙度大小；然后利用积分方式得到精确的目的层段孔隙度体积。
22、相对于现有技术，本发明具有以下有益效果：
23、本发明提供的缝洞型碳酸盐岩储集体孔隙体积定量计算方法，通过地震正演和反演结合，相互标定，可以准确刻画储集体的分布范围、顶底面位置和孔隙度大小，并计算储层孔隙空间，得到的孔隙度体和正演模型剖面具有很好的一致性，在不同厚度、不同孔隙度下预测精度显著提高。
1.一种缝洞型碳酸盐岩储集体孔隙体积定量计算方法，其特征在于，包括：
2.根据权利要求1所述的方法，其特征在于，步骤s1中，选择测井井径良好的井段的声波和密度资料，计算纵波阻抗及地层的反射系数，开展井震标定，结合井旁地震道提取地震平均子波以及地震资料的频带范围。
3.根据权利要求2所述的方法，其特征在于，所述纵波阻抗zi为：
4.根据权利要求1所述的方法，其特征在于，步骤s2中，利用实际地震资料开展稀疏脉冲反演，消除子波调谐效应，得到带限波阻抗数据体；然后开展缝洞体解释，得到主要储集体的空间包络顶面、底面位置以及平均储集体厚度。
5.根据权利要求1所述的方法，其特征在于，步骤s3中，利用实际测井资料中经过岩心标定的孔隙度和阻抗曲线，进行交会图分析，拟合建立孔隙度和纵波阻抗变换函数，基于平均储集体厚度建立从小到大的不同孔隙度变化地质模型，利用变换函数计算得到各孔隙度模型对应的纵波阻抗值，填入地质模型，利用实际提取子波开展波动方程正演。
6.根据权利要求5所述的方法，其特征在于，步骤s3中对正演资料进行处理包括：对采集得到的正演地震单炮进行去噪和偏移处理，得到偏移后的地震全叠加数据。
7.根据权利要求1所述的方法，其特征在于，步骤s4中，利用地震子波和正演得到的地震资料，开展第一轮正演数据的约束稀疏脉冲反演，消除子波调谐效应，得到正演模型的带限波阻抗数据体，对比正演模型的反演带限波阻抗体及正演模型的实际波阻抗值，得到储层波阻抗校正系数及储集体的纵波阻抗截止值。
8.根据权利要求7所述的方法，其特征在于，步骤s5中，利用步骤s4中得到的纵波阻抗截止值，将步骤s2中的带限反演体的纵波阻抗提取出来，并应用步骤s4中得到的校正系数进行校正，之后加上纯碳酸盐岩岩石纵波阻抗，得到反演的纵波阻抗低频模型。
9.根据权利要求5所述的方法，其特征在于，步骤s6中，利用步骤s1中的地震子波和步骤s5中的低频模型，开展第二轮稀疏脉冲反演，得到全频的纵波阻抗数据体；然后利用步骤s3中的岩石物理模型得到的孔隙度和纵波阻抗变换表，将纵波阻抗体转换为孔隙度体；对孔隙度体进行解释、雕刻，得到地下精确的孔洞顶底面位置及孔隙度大小；然后利用积分方式得到精确的目的层段孔隙度体积。
