本发明涉及油田开发,特别是涉及到一种改进型的储层物性参数同步反演方法。
背景技术:
1、储层物性参数地震叠前反演通过研究地震叠前反射信号振幅随入射角的变化规律、以及弹性参数随物性参数的变化规律,可以直接从地震叠前反射信号振幅中反演获得多种反映地下岩石及其饱和流体的物性参数,为进行储层预测和油气藏描述提供了定量化的依据。由于以下三个方面的原因,地震叠前反演的适定性受到很大的挑战:(1)由于大地的滤波作用,我们能获得的地震信号是带限的;(2)地震信号中包含了多种类型的噪音,通常信噪比不高;(3)由于认识的差距,在地震反演过程中使用了大量的近似公式和模型。上述原因一方面带来了储层物性参数地震反演问题的解具有不唯一性。
2、针对这些不足,国内外的学者们开展了基于岩石物理的储层物性反演研究。斯坦福大学的mukerj i等(2001)和eidsvik等(2004)等人率先将岩石物理和叠前反演相结合,估计储层物性参数,bachrach(2006)进一步地将岩石物理模型用于储层孔隙度和饱和度反演,larsen等(2006)等人则提出了一种基于马尔科夫链的储层流体反演方法,随后grana等(2010)等人发展了这些前人的研究成果,将岩石物理和贝叶斯反演相结合,实现了物性的反演,取得了较好的效果。国内在基于岩石物理的储层物性反演方面的研究较少,郑继新等(2009)首次将统计岩石物理用于含气储层饱和度与孔隙度联合反演,取得了一定的成果,但是其对岩石物理模型和实际应用的研究不够深入,属于一种前瞻性的研究,胡华锋等(2012)进一步发展了这种方法,提出了基于贝叶斯分类的储层物性参数联合反演方法,也取得了较好的应用效果,但是方法研究仍需要进一步深入。杨培杰等(2017)在前人研究的基础上,针对砂泥岩储层,基于岩石物理理论、贝叶斯理论、蒙特卡罗和遗传算法,提出了一种孔隙度和含水饱和度同步反演方法。
3、如何从地震数据中有效地提取孔隙度、泥质含量、含水饱和度等储层物性参数是储层地球物理研究的热点和难点,目前的物性提取方法主要是通过两步法实现的,即首先反演弹性参数,然后通过弹性参数再次计算物性参数,该方法的不足是效率低,并且准确度不高。为提高物性提取的效率和准确度,以贝叶斯理论为基础,开展一种改进型的储层物性参数同步反演方法。
4、在申请号:cn201410513635.1的中国专利申请中,涉及到一种岩相约束储层物性参数反演方法及装置,其中方法包括:根据贝叶斯反演框架,结合岩相信息,建立基于岩相约束的目标反演函数;引入岩相信息,建立统计岩石物理岩相模型;利用随机模拟技术,基于统计岩石物理岩相模型产生包含岩相信息的弹性参数与储层物性参数训练样本集;输入叠前地震反演成果数据及岩相数据,利用训练样本集,对基于岩相约束的目标反演函数进行求解。该发明将岩相信息作为反演约束信息,通过加入岩相约束的手段使得反演原理更加严谨、反演过程更加合理,能大大减少反演结果的多解性;推导出的目标反演函数纳入了岩相信息,求解过程更具针对性,结果更加可靠;能够精准地识别出储层的有利位置。
5、在申请号:cn201710753641.8的中国专利申请中,涉及到一种储层物性参数直接提取新方法,该储层物性参数直接提取新方法包括:步骤1,输入地震数据、地震子波、测井数据以及地质层位;步骤2,以统计岩石物理分析为基础,从测井数据出发,建立起储层物性参数和弹性参数之间的定量关系;步骤3,推导物性提取的近似公式,里面包含有孔隙度、泥质含量、密度三个参数;步骤4,以叠前分析为主要手段,结合物性提取的近似公式和点对点方程,从叠前数据中直接提取物性参数。该储层物性参数直接提取新方法可以同时反演储层孔隙度、泥质含量和密度等参数,提高储层物性估计结果的客观性和准确性,降低投资风险,为地质人员估算油气储量、确定井位提供更为可靠的依据。
6、在申请号:cn201710501254.5的中国专利申请中,涉及到一种确定储层参数的方法及装置。所述方法包括:获取所述目的储层中钻井位置处岩石的属性信息;基于所述岩石的属性信息,建立与所述目的储层中钻井位置处的储层参数、纵横波速度比、纵波阻抗相关联的目标岩石物理模板;获取所述目的储层的叠前地震数据,基于所述叠前地震数据,确定所述目的储层的反演结果数据体;所述反演结果数据体包括多个数据点;基于所述目标岩石物理模板,从所述反演结果数据体中确定种子数据点;基于所述目标岩石物理模板和所述种子数据点,确定所述反演结果数据体中每一个数据点对应的储层参数值。该申请提供的技术方案,可以提高所确定的储层参数的准确度。
7、以上现有技术均与本发明有较大区别,未能解决我们想要解决的技术问题,为此我们发明了一种新的改进型的储层物性参数同步反演方法。
技术实现思路
1、本发明的目的是提供一种可以同步反演储层孔隙度、泥质含量和密度等参数,提高储层物性估计结果的客观性和准确性的叠前地震资料直接同步反演得到储层物性参数的改进型的储层物性参数同步反演方法。
2、本发明的目的可通过如下技术措施来实现:改进型的储层物性参数同步反演方法,该改进型的储层物性参数同步反演方法包括:
3、步骤1,进行叠前道集处理及优化;
4、步骤2,进行精细的合成记录标定;
5、步骤3,以统计岩石物理分析为基础,从井、地震数据出发,利用贝叶斯理论将地震反射信号与物性模型参数建立岩石物理模型;
6、步骤4,以叠前分析为主要手段,结合储层物性提取的近似公式和点对点方程,从叠前地震资料直接同步反演得到储层物性参数。
7、本发明的目的还可通过如下技术措施来实现:
8、在步骤1中,地震数据为叠前偏移道集或叠前角道集数据,叠前道集是储层物性参数反演的前提和基础,叠前道集质量的好坏关系到参数反演的可靠性。
9、在步骤1中,开展叠前道集数据的处理及优化;开展叠前共反射点道集或是共中心点道集的经过保幅、去多次、去噪、拉平这些优化处理,并分取不同的角道集部分叠加数据。
10、在步骤2中,通过步骤1得到可靠的叠前道集数据的基础上,完成精细的合成记录标定与子波提取,确定地震叠前道集数据充分满足叠前反演的需求,这是利用岩石物理模型建立地震反射信号与储层物性参数之间的桥梁关系,也是开展储层物性参数反演的重要手段。
11、在步骤3中,物性参数和地震反射信号建立的定量关系的表达式为
12、g.f(χ)+ε=d (1)
13、其中,g是地震子波矩阵,f是关于χ的函数,χ是储层物性参数构成的反射系数向量,d是地震反射信号,ε是随机噪音向量。
14、在步骤3中,为了同时反演得到储层物性参数χ,在贝叶斯理论框架下将地震反射信号d与模型参数(χ)同时作为随机变量,反演目标函数是以地震反射信号d为条件的关于模型参数(χ)的条件概率,反演结果是条件概率最大时对应的模型参数,可以表示为:
15、χ=argmax{p(χ|d)} (2)
16、根据贝叶斯原理,式(2)可以表示为地震似然函数p(χ)与模型先验函数的乘积:
17、χ=argmax{p(χ|d).p(χ)} (3)。
18、在步骤4中,将物性参数近似方程和叠前反演有效结合,开发了一种点对点的叠前反演策略,通过点对点方程实现约束。
19、在步骤4中,建立的点对点方程为:
20、
21、其中,ξ表示待提取的物性参数,n是约束点的数量,是变化的,m≥n>1,m是待反演的地震道的总点数,r表示地层反射系数。
22、本发明中的改进型的储层物性参数同步反演方法,提高物性提取的效率和准确度,以贝叶斯理论为基础,开展一种改进型的储层物性参数同步反演方法。本发明利用叠前地震资料直接同步反演得到储层物性参数,提高储层物性估计结果的客观性和准确性,降低了子波提取准确性较差、计算效率偏低、分步获取导致的累计误差,为地质人员估算油气储量、确定井位提供更为可靠的依据。本发明可以同步反演储层孔隙度、泥质含量和密度等参数,提高储层物性估计结果的客观性和准确性的叠前地震资料直接同步反演得到储层物性参数。