发明名称 采用点谱技术计算吸收系数的方法
摘要 本发明涉及一种采用点谱技术计算吸收系数的方法,选取复数小波,即:Morlet小波,为高斯包络下的单频复正弦函数,利用连续小波变换对地震信号进行分析,求取每一点各尺度的瞬时振幅,并对着每一点的每一尺度都进行插值,得到每一点上每一个尺度的值,即:点谱,用高低频瞬时振幅比求取吸收系数。本发明所提供的采用点谱技术计算吸收系数的方法,其计算结果准确,不存在人为的异常;计算精度高;弥补传统基于傅立叶分析计算吸收系数的不足,能更好地预测油气,提高石油勘探开发的精度。
申请公布号 CN1873442A 申请公布日期 2006.12.06
申请号 CN200510074810.2 申请日期 2005.06.03
申请人 中国石油大学(北京) 发明人 曹思远;周鹏;李国福;张凤君;韩瑞冬
分类号 G01V1/28(2006.01) 主分类号 G01V1/28(2006.01)
代理机构 北京同立钧成知识产权代理有限公司 代理人 刘芳
主权项 1、一种采用点谱技术计算吸收系数的方法,其特征在于,它包括如下步骤:步骤1:选取复数小波,利用连续小波变换对地震信号进行分析,具体地为,假设x(t)为一道地震信号,对x(t)做小波变换:<math> <mrow> <msub> <mi>W</mi> <mi>x</mi> </msub> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>a</mi> </mfrac> <munderover> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mrow> <mo>+</mo> <mo>&infin;</mo> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>&phi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>dt</mi> </mrow> </math> <math> <mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mi>a</mi> </mfrac> <munderover> <mi>&Sigma;</mi> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mrow> <mo>+</mo> <mo>&infin;</mo> </mrow> </munderover> <mo>[</mo> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>i&phi;</mi> <mi>I</mi> </msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mo>]</mo> <mi>dt</mi> </mrow> </math> <math> <mrow> <mo>=</mo> <msub> <mi>W</mi> <mi>xR</mi> </msub> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>iW</mi> <mi>xI</mi> </msub> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> </mrow> </math> 式中,b为时间因子,且t,b∈R;a为尺度因子,且a∈R;为小波函数的共轭;WxR(b,a)为小波变换的实部;WxI(b,a)为小波变换的虚部;并且小波函数φ(t)满足φ∈L1(R)∩L2(R)<math> <mrow> <msub> <mi>c</mi> <mi>&phi;</mi> </msub> <mo>=</mo> <munderover> <mo>&Integral;</mo> <mn>0</mn> <mo>&infin;</mo> </munderover> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <mi>&omega;</mi> <mo>)</mo> </mrow> </mrow> <mi>&omega;</mi> </mfrac> <mi>d&omega;</mi> <mo>&lt;</mo> <mo>&infin;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <msub> <mi>c</mi> <mi>&phi;</mi> </msub> <mo>&NotEqual;</mo> <mn>0</mn> </mrow> </math> 这样便得到小波变换后的实部WxR(b,a)和虚部WxI(b,a);步骤2:利用步骤1对信号x(t)进行小波变换所得到的小波系数Wx(b,a)的实部和虚部相差90度的相位角的变换性质,直接求取地震信号x(t)在不同尺度上的瞬时振幅;定义对应不同尺度因子a下的瞬时振幅为:<math> <mrow> <mi>A</mi> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <msubsup> <mi>W</mi> <mi>xR</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> <mo>+</mo> <msubsup> <mi>W</mi> <mi>xI</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <mi>b</mi> <mo>,</mo> <mi>a</mi> <mo>)</mo> </mrow> </msqrt> </mrow> </math> 得到每一点的瞬时振幅;步骤3:根据步骤2得到的瞬时振幅,得到每一点的吸收系数:<math> <mrow> <msub> <mi>&xi;</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>&kappa;</mi> <mi>i</mi> </msub> <mfrac> <msub> <mi>A</mi> <mi>Hi</mi> </msub> <msub> <mi>A</mi> <mi>L</mi> </msub> </mfrac> </mrow> </math> 式中AHi,AL分别表示高、低频的瞬时振幅;κi为对应不同AHi校正系数;ξi为吸收系数;i为高频段的频数序号,a∈[ai-1,ai],AHi为此区间的振幅积分;AL具体低频频率由试验确定;然后根据井旁道测井显示油气层标定井旁道地震道吸收系数的κi校正系数,根据井信息,能够确定假定κi=1的时候,井旁道计算吸收系数的误差,根据已知的井旁的岩性,确定它的衰减系数,然后外推到井旁道的远道;通过上式得出对应每一个地层及气层的敏感的吸收频段a∈[ai-1,ai]。
地址 102249北京市昌平区府学路20号