专利名称：流体因子提取方法
技术领域：
本发明涉及地震勘探领域。更具体地讲，涉及一种流体因子提取方法。
背景技术：
在地震勘探中，流体因子的提取是一项非常重要的技术。精确的流体因子可以更准确地对地下的矿藏进行检测和识别。当前，可通过时频分析方法来对地震数据进行分析， 以获取流体因子。时频分析方法可提供时间域和频率域的联合分布信息，能够更清楚地描述信号的频率随时间变化的关系。下面，示出一种现有技术的流体因子提取方法。对于给定的某道地震数据χ (t)，其Wigner (维格纳)分布Wx (t，f)的定义如下
权利要求
1.一种流体因子提取方法，包括(a)检测工区内的至少一道地震数据；(b)设置用于分数阶傅里叶变换的分数阶数α的初始值；(C)使用分数阶数α，计算所述至少一道地震数据中的一道地震数据的分数阶傅里叶变换；(d)计算分数阶傅里叶变换的解析信号，并计算解析信号的带交叉项抑制的 Wigner-Ville 分布；(e)计算带交叉项抑制的Wigner-Ville分布的带窗的Radon变换；(f)利用带窗的Radon变换重构时频分布图；(g)从时频分布图提取流体因子。
2.如权利要求1所述的方法，其中通过将带窗的Radon变换所在的域顺时针旋转分数阶数α来重构时频分布图。
3.如权利要求1所述的方法，还包括(h)使用提取的流体因子构建地质结构；(i)对构建的地质结构与根据测井数据获取的地质结构进行比较，确定构建的地质结构与根据测井数据获取的地质结构是否在预定误差范围内吻合；(j)如果构建的地质结构与根据测井数据获取的地质结构在预定误差范围内吻合，则确定得到最优的分数阶数α。
4.如权利要求3所述的方法，其中，如果构建的地质结构与根据测井数据获取的地质结构没有在预定误差范围内吻合，则改变分数阶数α的值，返回步骤(C)。
5.如权利要求3所述的方法，还包括当所述至少一道地震数据为至少两道地震数据时，使用最优的分数阶数α对所述至少一道地震数据中的其他道地震数据执行步骤 (c)-(g)。
6.如权利要求1所述的方法，其中，从时频分布图提取流体因子的步骤包括从时频分布图提取频率属性因子。
7.如权利要求6所述的方法，其中，频率属性因子包括峰值瞬时值频率、均方根频率、 频率衰减梯度中的至少一个。
8.如权利要求1所述的方法，其中，根据当前已知的地震资料来设置分数阶数α的初始值。
9.如权利要求8所述的方法，其中，当前已知的地震资料是测井数据或者是邻近工区的地质结构资料。
10.如权利要求3所述的方法，其中，从所述工区内的测井获取所述测井数据。
11.如权利要求10所述的方法，其中，在检测所述至少一道地震数据之前或之后形成所述测井。
12.一种地质结构检测方法，其特征在于，使用根据权利要求1-11中的任意一项提取的流体因子来检测地质结构。
全文摘要
本发明提供一种流体因子提取方法。所述方法包括检测工区内的至少一道地震数据；设置用于分数阶傅里叶变换的分数阶数的初始值；使用分数阶数，计算所述至少一道地震数据中的一道地震数据的分数阶傅里叶变换；计算分数阶傅里叶变换的解析信号，并计算解析信号的带交叉项抑制的Wigner-Ville分布；计算带交叉项抑制的Wigner-Ville分布的带窗的Radon变换；利用带窗的Radon变换重构时频分布图；从时频分布图提取流体因子。根据本发明的流体因子提取方法，即使在地震波中存在色散情况，也能准确地从地震波中提取流体因子。
文档编号G01V1/28GK102305940SQ20111013496
公开日2012年1月4日 申请日期2011年5月24日 优先权日2011年5月24日
发明者吴秋波, 张洞君, 张靖, 彭真明, 李忠, 洪余刚, 邹文 申请人:中国石油集团川庆钻探工程有限公司