专利名称：频散avo反演方法
技术领域：
本发明涉及地震勘探技术领域，是一种新的地震属性，用于描述地震信号的频散程度，具体讲是一种频散AVO (Amplitude Versus Offset)反演方法。
背景技术：
目前的纵波AVO反演技术仅仅考虑了地震波的振幅信息随偏移距的变化情况，认为储层中流体性质和流体饱和度的变化导致地层中波阻抗的差异，从而引起地震波反射振幅的变化，而这种变化是与频率无关的。然而无论砂岩储层还是碳酸盐岩储层，都是由固体和液体两种部分组成的双相介质，流体的存在导致地震波发生不同程度的频散和衰减。目前的AVO分析忽略了振幅随频率的变化特征，而这一特征受储层中所含流体性质的影响， 是描述储层中流体性质的依据之一。Zo印pritz方程是AVO技术的理论基础。当弹性波传播到弹性分界面时，要发生波型转换和能量重新分配，在不同介质分界面上，根据位移和应力的边界条件，通过解波动方程，可以得到弹性分界面能量分配关系。Knott和heppritz利用反射界面两侧位移和应力的连续性作为边界条件，得到了反射系数和透射系数，它们是入射角和介质弹性参数(密度、体积模量和剪切模量)的函数，即Zc^ppritz方程
权利要求
1.一种频散AVO反演方法，其特征在于，所述频散AVO反演方法包括 获取叠前道集，并对所述叠前道集进行球面扩散补偿和动校正；对所述进行球面扩散补偿和动校正的叠前道集进行时频分析，生成不同频率下的振幅谱，所述振幅谱为采样时间、接收道号和频率的函数； 对所述不同频率下的振幅谱进行振幅谱均衡处理；获取所述叠前道集中的采样点的叠加速度，并根据所述叠加速度生成所述采样点的层速度；根据所述谱均衡处理后的振幅谱和所述采样点的层速度，采用最小二乘反演，生成设定频率下的纵横波速度变化率和所述采样点的频散程度，所述频散程度为纵横波速度变化率随频率变化的快慢；将所述生成的频散程度按照道集序号进行排列，生成频散程度属性剖面。
2.如权利要求1所述的频散AVO反演方法，其特征在于，所述对进行球面扩散补偿和动校正的叠前道集进行时频分析，生成不同频率下的振幅谱，包括所述叠前道集有η个接收道，表示为矩阵形式为s (t，η)，采用短时傅里叶变换或者小波变换的时频分析方法，进行频谱分解得到不同频率f下的振幅谱S (t，η, f)。
3.如权利要求1所述的频散AVO反演方法，其特征在于，所述对不同频率下的振幅谱进行振幅谱均衡处理，包括对不同频率的振幅谱通过加权函数ω进行谱均衡Sb(t,n,f) = S(t，n，f) ω (f)，其中，S (t，n，f)为不同频率下的振幅谱，&(t，n，f)为进行谱均衡处理后的振幅谱，ω (f)为进行谱均衡处理的加权函数。
4.如权利要求1所述的频散AVO反演方法，其特征在于，所述根据所述叠加速度生成所述采样点的层速度，包括获取所述叠前道集中的采样点的叠加速度，根据Dix公式
5.如权利要求3所述的频散AVO反演方法，其特征在于，所述根据所述谱均衡处理后的振幅谱和所述采样点的层速度，根据最小二乘反演，生成设定频率下的纵横波速度变化率和所述采样点的频散程度，包括根据公式下述两个公式获取所述纵横波速度变化率和频散程度
6.如权利要求1所述的频散AVO反演方法，其特征在于，所述叠前道集包括共反射点道集CRP、共中心点道集CMP和共深度点道集⑶P。
全文摘要
本发明公开了一种频散AVO反演方法，包括获取叠前道集，并对叠前道集进行球面扩散补偿和动校正；对叠前道集进行瞬时频谱分析，生成不同频率下的振幅谱；对振幅谱进行振幅谱均衡处理；获取叠前道集中的采样点的叠加速度，并根据叠加速度生成采样点的层速度；根据谱均衡处理后的振幅谱和采样点的层速度，根据最小二乘反演，生成设定频率下的纵横波速度变化率和所述采样点的频散程度；将生成的频散程度按照道集序号进行排列，生成频散程度属性剖面。本发明实施例的频散AVO反演方法，消除了由于弹性界面波阻抗差异导致的强震幅能量团的影响，也可避免叠后地震资料频谱分析时叠加处理可能导致的虚假“频率异常”，使得到的结果更加可靠。
文档编号G01V1/28GK102478667SQ201010567538
公开日2012年5月30日 申请日期2010年11月30日 优先权日2010年11月30日
发明者吴小羊, 张峰, 李向阳, 陈双全 申请人:中国石油大学(北京), 中国石油天然气集团公司