主权项 |
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>∫</mo> <mrow> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>φ</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>Σ</mi> <mrow> <mo>-</mo> <mo>∞</mo> </mrow> <mrow> <mo>+</mo> <mo>∞</mo> </mrow> </munderover> <mo>[</mo> <msub> <mi>φ</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φ</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>φ</mi> </msub> <mo>=</mo> <munderover> <mo>∫</mo> <mn>0</mn> <mo>∞</mo> </munderover> <mfrac> <mrow> <msub> <mi>φ</mi> <mi>R</mi> </msub> <mrow> <mo>(</mo> <mi>ω</mi> <mo>)</mo> </mrow> </mrow> <mi>ω</mi> </mfrac> <mi>dω</mi> <mo><</mo> <mo>∞</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <msub> <mi>c</mi> <mi>φ</mi> </msub> <mo>≠</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>ξ</mi> <mi>i</mi> </msub> <mo>=</mo> <msub> <mi>κ</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]。 |