地震资料约束下的品质因子计算方法
【技术领域】
[0001] 本发明涉及油气地球物理勘探领域中的数据处理技术领域，具体地，涉及一种地 震资料约束下的品质因子计算方法。
【背景技术】
[0002] 在油气田和矿产资源地震勘探中反射地震资料处理技术中，品质因子Q作为表征 介质吸收衰减特性的重要参量，在预测岩性、油气藏位置和范围等方面具有重要的意义，因 此如何有效的计算出Q值对于后续的地球物理勘探工作是相当重要的。现有的Q值计算方 法主要有李庆忠经验公式法、谱比法等。
[0003] 谱比法：频率为f的地震波的振幅A随时间t的变化关系可以用公式 J(/,Z2) =火进行描述，tl、t2为不同的时刻。对该公式两边求取自然对数可 得：
【主权项】
1. 地震资料约束下的品质因子计算方法，采用公式Q = a ·ν pb计算品质因子Q，V p为纵 波层速度，其特征在于，斜率a和指数b动态取值，且斜率a和指数b在不同取值下计算出 的品质因子值还进行反Q滤波计算，获得各自的滤波道集；所有滤波道集还与正演道集进 行相似性比较，将与正演道集最相似的滤波道集对应的品质因子值作为品质因子Q的最终 确定值。
2. 根据权利要求1所述的地震资料约束下的品质因子计算方法，其特征在于，具体包 括以下步骤： 51、 确定斜率a和指数b的取值区间和取值步长； 52、 初始化斜率a和指数b的值； 53、 将a、b、Vp代入Q = a · V pb计算Q值； 54、 根据求得的Q值对叠前道集或者叠后道集进行反Q滤波计算，得到滤去Q值的滤波 道集； 55、 将步骤S4中最新得到的滤波道集与正演道集进行相似性比较，记录比较结果； 56、 判断斜率a和指数b的取值是否完成所有取值组合的遍历，是则跳转到步骤S7 ;否 则改变a和/或b的取值，使a和b被赋予一组未使用过的值组合，跳转到步骤S3重新进 行循环； 57、 将与正演道集最相似的滤波道集对应的品质因子值作为品质因子Q的最终确定 值。
3. 根据权利要求2所述的地震资料约束下的品质因子计算方法，其特征在于，步骤S5 中基于拟合优度对正演道集与步骤S4中最新得到的滤波道集进行相似比较，计算出拟合 优度值并记录。
4. 根据权利要求3所述的地震资料约束下的品质因子计算方法，其特征在于，所述拟 合优度值为步骤S4中最新得到的滤波道集平方和与正演道集平方和的比值。
5. 根据权利要求1至4任一所述的地震资料约束下的品质因子计算方法，其特征在于， 所述正演道集由测井数据生成。
6. 根据权利要求2至4任一所述的地震资料约束下的品质因子计算方法，其特征在于， 步骤S2中，初始化斜率a和指数b的值时将a的取值区间的下限值赋给斜率a，将b的取值 区间的下限值赋给b。
7. 根据权利要求2至4任一所述的地震资料约束下的品质因子计算方法，其特征在于， 步骤Sl中，在确定了斜率a和指数b的取值区间和取值步长后，还包括构建a和b的 取值空间的步骤； 步骤S6中，判断斜率a和指数b的取值是否遍历完取值空间的所有取值组合，是则跳 转到步骤S7,否则在取值空间选取一个新的取值组合赋给a和b，并跳转到步骤S3重新进 行循环。
【专利摘要】地震资料约束下的品质因子计算方法，包括以下步骤：S1、确定a和b的取值区间和取值步长；S2、初始化a和b；S3、将a、b、纵波层速度Vp代入Q＝a·Vpb计算Q值；S4、根据求得的Q值进行反Q滤波计算，得到滤波道集；S5、将步骤S4中得到的滤波道集与正演道集进行相似行分析，记录比较结果；S6、判断a和b的取值是否完成所有取值组合的遍历，是则将与正演道集最相似的滤波道集对应的品质因子值作为品质因子Q的最终确定值，否则改变a和/或b的取值，跳转到步骤S3重新进行循环。本发明在李庆忠经验公式法的基础上动态选取斜率a和指数b的值，a、b取值灵活，在正演道集的控制下,计算的Q值更加符合实际的地质情况。
【IPC分类】G01V1-28
【公开号】CN104635268
【申请号】CN201510101505
【发明人】赵奇, 刘力辉
【申请人】成都晶石石油科技有限公司
【公开日】2015年5月20日
【申请日】2015年3月9日