主权项 |
本发明涉及一种基于优化系数的混合域傅立叶有限差分叠前深度偏移方法,其特征在于包括如下步骤:步骤1:将二维声波方程变换到频率‑波数域,并用波场的背景速度v<sub>0</sub>(x,z)代替介质的实际速度v(x,z)造成的频散方程的单平方根的误差,其误差如(1)式所示:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>ξ</mi><mo>=</mo><msup><mrow><mo>[</mo><mfrac><msup><mi>ω</mi><mn>2</mn></msup><mrow><msup><mi>v</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><msubsup><mi>k</mi><mi>x</mi><mn>2</mn></msubsup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><msup><mrow><mo>[</mo><mfrac><msup><mi>ω</mi><mn>2</mn></msup><mrow><msubsup><mi>v</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><msubsup><mi>k</mi><mi>x</mi><mn>2</mn></msubsup><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000817493160000011.GIF" wi="1316" he="166" /></maths>(1)式中:x、z为方向坐标,k<sub>x</sub>,k<sub>z</sub>分别为x、z方向上的波数,单位为:/m,v<sub>0</sub>(x,z)为波场的背景速度,单位为:m/s,v(x,z)为波场的实际速度,单位为:m/s,ω为角频率,单位为:rad/s;设d=v(x,z)·k<sub>x</sub>/ω,d∈(‑1,1),则(1)式中根式可以写为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>d</mi><mo>)</mo></mrow><mo>=</mo><msup><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>d</mi><mn>2</mn></msup><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000817493160000012.GIF" wi="1196" he="119" /></maths>步骤2:利用padé有理函数逼近(2)式,表达式如(3)式所示:<img file="FDA0000817493160000013.GIF" wi="1352" he="163" />(3)式中:d=v(x,z)·k<sub>x</sub>/ω,d∈(‑1,1),a<sub>0</sub>,a<sub>1</sub>,b<sub>1</sub>为逼近式系数,T<sub>0</sub>(d),T<sub>1</sub>(d)为第一类切比雪夫多项式;步骤3利用切比雪夫多项式函数对(3)式中的系数a<sub>0</sub>,a<sub>1</sub>,b<sub>1</sub>进行优化,取切比雪夫多项式函数的前三项与(3)式相等,系数a<sub>0</sub>=1.358,1a<sub>1</sub>=1.2866,b<sub>1</sub>=0.3684,则(3)式可表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>R</mi><mrow><mo>(</mo><mn>1</mn><mo>,</mo><mn>1</mn><mo>)</mo></mrow></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1.3581</mn><mo>-</mo><mn>1.2866</mn><msup><mrow><mo>(</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>k</mi><mi>x</mi></msub></mrow><mi>ω</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mn>1.3684</mn><mo>-</mo><mn>0.7368</mn><msup><mrow><mo>(</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><msub><mi>k</mi><mi>x</mi></msub></mrow><mi>ω</mi></mfrac><mo>)</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000817493160000014.GIF" wi="1411" he="327" /></maths>步骤4将(4)式代入到(1)式并变换到频率‑空间域,代入到优化系数的混合域叠前深度偏移下行波外推方程中,如(5)式所示:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>u</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>+</mo><mi>Δ</mi><mi>z</mi><mo>;</mo><mi>ω</mi><mo>)</mo></mrow><mo>=</mo><mi>u</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>ω</mi><mo>)</mo></mrow><mo>·</mo><msup><mi>e</mi><mrow><mo>-</mo><msub><mi>ik</mi><mi>z</mi></msub><mo>·</mo><mi>Δ</mi><mi>z</mi></mrow></msup></mrow></mtd></mtr><mtr><mtd><mrow><mo>=</mo><mi>u</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>;</mo><mi>ω</mi><mo>)</mo></mrow><mo>·</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><msup><mrow><mo>[</mo><mfrac><msup><mi>ω</mi><mn>2</mn></msup><mrow><msubsup><mi>v</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>+</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>]</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup><mo>·</mo><mi>Δ</mi><mi>z</mi></mrow></msup><mo>·</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mrow><mo>(</mo><mfrac><mi>ω</mi><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mi>ω</mi><mrow><msub><mi>v</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>)</mo></mrow><mo>·</mo><mi>Δ</mi><mi>z</mi></mrow></msup><mo>·</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>i</mi><mfrac><mrow><mn>0.5601</mn><mo>×</mo><mo>[</mo><mfrac><mi>ω</mi><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mi>ω</mi><mrow><msub><mi>v</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow></mfrac><mo>]</mo><mo>+</mo><mn>1.8364</mn><mo>×</mo><mfrac><mrow><mi>v</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>v</mi><mn>0</mn></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><mi>ω</mi></mfrac><mo>·</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac></mrow><mrow><mn>0.6316</mn><mo>+</mo><mn>0.7368</mn><mo>×</mo><mfrac><mrow><msup><mi>v</mi><mn>2</mn></msup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow><mo>+</mo><msubsup><mi>v</mi><mn>0</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>z</mi><mo>)</mo></mrow></mrow><msup><mi>ω</mi><mn>2</mn></msup></mfrac><mo>·</mo><mfrac><msup><mo>∂</mo><mn>2</mn></msup><mrow><mo>∂</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>)</mo></mrow></mfrac><mo>·</mo><mi>Δ</mi><mi>z</mi></mrow></msup></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000817493160000015.GIF" wi="1903" he="305" /></maths>步骤5按照(5)式并以相关型成像原理对地震数据进行偏移延拓成像,获得最终的偏移成像叠加剖面。 |