基于地质统计先验信息的储层物性参数同步反演.pdf

2014 年09 月 APPLIED GEOPHYSICS （应用地球物理） 第11 卷，第3 期 基于地质统计先验信息的储层物性参数同步反演//Simultaneous inversion of petrophysical parameters based on geostatistical a priori information, APPLIED GEOPHYSICS, 2014, 11(3), P. 311-320. DOI: 10.1007/s11770-014-0445-1 基于地质统计先验信息的储层物性参数 同步反演 ♦ 印兴耀，孙瑞莹，王保丽，张广智 （中国石油大学（华东）地球科学与技术学院，山东 青岛 266580) 摘要：本文提出的储层物性参数同步反演是一种高分辨率的非线性反演方法，该 方法综合利用岩石物理和地质统计先验信息，在贝叶斯理论框架下，首先通过变 差结构分析得到合理的变差函数，进而利用快速傅里叶滑动平均模拟算法（Fast Fourier Transform-Moving Average, FFT-MA ）和逐渐变形算法（Gradual Deformation Method, GDM ）得到基于地质统计学的储层物性参数先验信息，然后根据统计岩 石物理模型建立弹性参数与储层物性参数之间的关系，构建似然函数，最终利用 Metropolis 算法实现后验概率密度的抽样，得到物性参数反演结果。并将此方法处 理了中国陆上探区的一块实际资料，本方法的反演结果具有较高的分辨率，与测 井数据吻合度较高；由于可以直接反演储层物性参数，避免了误差的累积，大大 减少了不确定性的传递，且计算效率较高。 关键词：地质统计先验信息，统计岩石物理，贝叶斯理论，物性参数同步反演 引言 等基于FFT 的离散谱模拟方法，通过快速傅里叶 变换（Fast Fourier Transform，简称FFT ）巧妙 地简化了Oliver （1995）提出的滑动平均模拟， 随着油气勘探开发难度加大及反演精度的 提出FFT-MA 模拟方法，并结合Hu （2000 ）提 提高，如何从已知资料获取储层物性参数成为储 出的逐步变形更新算法（GDM ），在石油工程领 层地球物理的研究热点。基于非线性理论的随机 域得到了成功应用，本文将FFT-MA 和GDM 引 反演方法能够在空间相关性和测井数据的约束 入地球物理反演中构建地质统计先验信息。由于 下，模拟出地震频带以外的信息，分辨率常高于 FFT-MA 是频率域的模拟方法，计算效率得到了 常规的确定性反演，因此更适用于岩性和物性高 较大的提高。 度非均质的复杂隐蔽性储层的研究（e.g.,Francis， 孔隙度、饱和度与泥质含量等储层物性参数 2005 ）。1992 年，Bortoli 等首次提出地质统计反 是进行储层描述的关键，而岩石物理是建立地震 演方法（Bortol