基于矢量波场分离弹性波逆时偏移成像-李振春.pdf

《基于矢量波场分离弹性波逆时偏移成像-李振春.pdf》由会员分享，可在线阅读，更多相关《基于矢量波场分离弹性波逆时偏移成像-李振春.pdf（7页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、2016年第40卷第1期中国石油大学学报(自然科学版)Journal of China University of PetroleumV0140 No1Feb2016文章编-号：1673-5005(2016)01-0042-07 doi：103969jissn1673-5005201601006基于矢量波场分离弹性波逆时偏移成像李振春，雍 鹏，黄建平，李庆洋，谷丙洛(中国石油大学地球科学与技术学院。山东青岛266580)摘要：在传统的弹性波逆时偏移中。通过Helmholtz分解获取的横波的极性随入射方向改变，导致转换波的成像值也随入射方向发生极性变化(极性反转)，多炮叠加后转换波成像同相轴将受

2、到破坏。在传统波场分离基础上，对标量势(纵波)与矢量势(横波)分别进行梯度和旋度处理，得到矢量纵波与矢量横波。针对矢量纵横波波场，引入内积成像条件进行偏移成像，通过这种成像方法无需对转换波成像进行极性校正。波场分析结果表明，该方法既适用于同类波高精度成像，也适用于转换波成像。模型实验结果表明了本文方法的正确性和适应性。关键词：弹性波逆时偏移；极性反转；波场分离；内积成像条件中图分类号：P 6314 文献标志码：A引用格式：李振春，雍鹏，黄建平，等基于矢量波场分离弹性波逆时偏移成像J中国石油大学学报(自然科学版)，2016，40(1)：4248LI Zhenchun，YONG Peng，HUAN

3、G Jianping，et a1Elastic wave reverse time migration based on vector wavefield sepera-tionJJournal of China University of Petroleum(Edition of Natural Science)，2016，40(1)：4248Elastic wave reverse time migration based on vector wavefield seperationLI Zhenchun，YONG Peng，HUANG Jianping，LI Qingyang，GU Bi

4、ngluo(School of Geosciences in China University of Petroleum，Qingdao 266580，China)Abstract：With the advances of seismic acquisition，seismic ima季ng has migrated from using single-component data alone tomulticomponent elastic wave ima西ngBecause of its high precision imaging ability，reversetime migrati

5、on(RTM)has beenwidely appfied to complex structure imagingHowever，S-mode polarization changes via incident direction，which will causepolarity reversal，and this polarity reversal will seriously damage the continuity of reflection events in conventional PS imagingas well as SP imaging after stacking o

6、ver multiple shotsIn the conventional elastic RTMHelmholtz decomposition is used toobtain the scalar potential(Pwave)and the vector potential(S-wave)In this paper，we intmduce gradient and curl opera-tors tO the scalar potential and the vector potential，respectivelyConsequently，we can obtain vector P

7、-waves and S-wavesAs for the vector wavefield，we introduce all inner product imaging condition to the elastic RTM，which Can avoid using extrapolarity correctionWavefield analyses show that this method is applicable not only to primary imaging(PP and SS)but alsoto conveaed wave imaging(PS and SP)The

8、effectiveness and feasibility of our new imaging method in simple and complexmodels are illustrated in numerical examples using synthetic dataKeywords：elastic wave reverse time migration；polarity reversal；wavefield decomposition；inner product ima西ng condition逆时偏移1。1采用全波算子，能够处理强速度变化与陡倾角成像，已经逐步成为地下复杂构

9、造成像的重要手段。多分量地震数据更能代表地震波传播的物理规律，因而相比声波逆时偏移，弹性波逆时偏移能够提供更多的地质信息。逆时偏移主要包括两大步骤：构建震源波场与检波波场，利用成像条件提取成像值。波场的构建通常是通过数值求解波动方程，有限差分法是一种广泛使用的偏微分方程数值解法m。由于弹性波逆时偏移使用多分量数据进行成像，因此成像条件也比较多。其中一种简单的成像条收稿日期：20150721基金项目：大地测量与地球动力学国家重点实验室重点基金项目(SKLGED201552一Ez)；国家”973”项目(2014CB239006，201ICB202402)；国家自然科学基金项目(41104069，4

10、1274124)作者简介：李振春(1963一)，男，教授，博士，博士生导师，研究方向为地震波正演及偏移成像。E-mail：leonliupceduca。通讯作者：黄建平(1982一)，男，副教授，博士，研究方向为地震波正演及偏移成像。E-mall：jphuangchinagmailcorn。万方数据第40卷第1期 李振春，等：基于矢量波场分离弹性波逆时偏移成像 43件是直接将震源波场与检波波场的各分量分别进行互相关成像6，由于纵横波耦合在一起，导致串扰噪声7|。另外一种广泛使用的成像方式是先进行纵横波波场分离7，然后再使用成像条件进行成像。纵横波分离通常采用Helmholtz分解8，对原始波场

11、(位移场或速度场)分别取散度与旋度，得到原始波场的标量势与矢量势，分别代表纵波与横波9。对震源波场中的P波与检波波场中的s波进行互相关成像得到PS成像剖面。相似地，可以得到PP、SP和SS成像剖面。由于矢量势(S波)存在极性反转，在转换波(PS、SP)成像剖面中，同相轴也存在极性反转。多炮叠加后同相轴的连续性将遭到严重损坏7。很多学者对此提出了一些解决手段：如角度域成像71m11和计算反射波传播方向进行极性校正121 4|。前者计算较为准确，但计算量较大；后者计算量较小，但由于在复杂地区地震波传播方向难以精确计算，因而精度较差。Duan和Sava15提出了一种基于地层倾角与波传播方向进行极陛校

12、正的方法，改进了转换波成像效果，但提前计算地层倾角，增加了计算量。笔者利用一阶弹性波方程实现波场构建，在传统的Helmholtz分解得到的标量势与矢量势基础上，对标量势与矢量势分别进行梯度运算和旋度运算得到矢量纵波与矢量横波。传统分离方式得到的矢量势横波随入射方向发生极性变化，在矢量的纵横波波场中，纵波与横波的极性都随着入射角方向发生改变，且矢量纵横波极性变化是一致的。据此针对矢量的纵横波波场，提出一种内积的成像条件从矢量纵横波中提取成像值。1方法原理弹性逆时偏移主要包含两步：首先构建震源波场与检波波场；然后，使用成像条件提取成像值。11二维波动方程二维的弹性波一阶速度一应力方程为Ov。1 f

13、a盯。Oor。1百2石li+百l警=古陪+爿，百Oorzx=(A+驰)警+A誓， (1)等=A瓦Ovxm+驰)警，等=肛陪+钧式中，口，、叱为质点速度；oror。、or。为应力分量；p为密度；A以为拉梅常数。本文中采用时间2阶、空间10阶精度的交错网格有限差分法对式(1)进行差分离散，获取震源波场与检波波场。在进行波场模拟时，由于人工截断边界会导致边界反射，采用完全匹配层(PML)吸收边界16J 7I，以压制边界反射提高模拟区域的精度与信噪比。12弹性波波场分离利用式(1)得到震源波场与检波波场的矢量速度波场。矢量速度场的各分量均包含了P波与S波。若直接将其各分量分别进行互相关成像，会因为多种

14、混合的波型产生串扰假象。此外，这种成像方式难以给出明确的物理含义。为此，Yan和Sava提出在互相关成像之前。利用Helmholtz分解从矢量波提取纯纵波分量与横波分量。9=V1， (2)砂=V， (3)其中1，=(叱，吃)注意，采用散度算子得到的纯纵波妒实际上是标量势，采用旋度算子得到的纯横波砂为矢量势。为了方便讨论，在波数域对式(2)与(3)进行表达：(p=ik哥=ikP哥， (4)步=i七哥=iks驴 (5)式中，哥为原始波场l，在波数域的形式；与分别为9与砂在波数域的形式；k为地震波的传播单位矢量；k，与k。分别为纵、横波的传播单位矢量。图l所示的几何关系中，p，巾。分别代表纵横波中质

15、点的振动方向。在震源两侧波场和传播方向的戈分量的极性出现反转，不妨令左右两侧的反射波场分别为l，l=(。，口：)，l，=(一秽；，秽：)传播方向分别为kl=(k，k：)，k，=(-k，k。)二维情况下纵横波只含有一个分量，将波场变量带入式(4)、(5)，震源左右两侧的纵波数值相同，而两侧的横波数值相差一个负号。图1 平层模型波场极性反转解释Fig1 Polarity reversal analysis of layer-model若直接使用互相关成像，由于横波存在极性反万方数据44 中国石油大学学报(自然科学版) 2016年2月转，在多炮叠加时转换波的成像剖面(Ps、sP)中的反射同相轴将会因

16、极性反转受到破坏。需要注意的是，在三维情况横波有3个分量，纵波只有一个分量，因此通过这种波场分离方式难以给出转换波成像条件以及相应的物理含义。为了避免极性反转本文中在传统Helmholtz分解得到的标量势(纵波)与矢量势(横波)的基础上，提出使用梯度算子与旋度算子分别作用标量势妒与矢量势砂，得到矢量的纯P波波场和矢量的纯S波波场9|。P=7 9， (6)s=7砂 (7)式中，尸与S分别表示矢量纵波与矢量横波。通过式(6)与式(7)将得到矢量的P波与S波波场，后面将利用该矢量波场进行逆时偏移成像，以解决转换波成像中的极性反转。本文中波场分离方式在波数域可表示为P=-k(k帝)=一kP(kP哥P)

17、， (8)S=kx(七哥)=ksX(ksxOs) (9)式中。户与雪分别为矢量P波与矢量S波在波数域的表示形式。本文的波场分离方式可以看做Zhang和McMechan方法18在空间域的实现。直接在空间域实现避免了傅里叶变化带来的巨大计算量。为了方便讨论，将二维情况下的本文波场分离方式在波数域展开可得P，=一(k，2-，+矗，k：矛：)， (10)P；=一(k k：云。+磋云：)， (11)S，=一(移，-k。k：矛：)， (12)s：=一(彰舀：-k，k。云，) (13)针对图1所示的观测系统，带入波场变量与波矢量，左右两侧的矢量纵波与横波的各个分量的极性是一致的。事实上。在图l所示的几何系统

18、中，根据梯度、散度、旋度的性质可知，矢量势横波砂的方向垂直于剖面，而通过式(6)、(7)的分解，矢量纵波与横波的方向与原始的波场中的质点的振动方向p，妒。是一致的。13内积成像方式针对得到的矢量P波与矢量S波，同时为了补偿深部的照明不足与压制噪声，笔者给出一种震源归一化的内积成像条件=黯等篙端，(14，州=鼍等黄端， ，。=甓等黜，(16)，r1、一量：(兰!)：苎!(兰!1 2堕! r 1，7、155、4 7一fS3(x，t)S5(z，t)dt 、17式中，l(x)为成像值，其下标第一个字母代表震源波场，下标第二个字母代表检波波场；等式右边的上标s、r分别代表震源波场与检波波场。时间积分在数

19、值计算中表现为从激发时刻开始，对同一成像点的不同时刻的内积结果进行累加。对于复杂模型，深部通常存在照明不足的问题，使用归一化处理有利于增强深部的能量。例如，k(x)代表矢量P波震源波场与矢量S波在成像点戈处的归一化的成像值(Ps成像结果)。由于通过矢量的波场分解方式得到的纵横波的极性变化与质点偏振方向有关，这样成像值的极性也与质点偏振方向有关。由于地震勘探的主要能量的人射角都小于40。，因此采用内积成像方式不会出现极性反转。针对矢量纵横波场的内积成像方式将PP、Ps、sP、ss的成像统一于相同的成像框架中。传统的互相关成像方式中使用的纵横波是标量势与矢量势，通过Helmholtz分解得到S波的

20、极性与地震波入射到反射界面的几何关系方向有关。在矢量的纵横波波场中。纵波与横波的极性随着入射角的大小而改变。此外，内积成像条件不仅适用于二维的弹性波逆时偏移，同样适用于三维的弹性波逆时偏移。2正演波场分析为了说明极性反转现象，解释本文波场分离方法，本文中采用一个三层模型，该模型假定密度为常数，横波速度分别为2000、3000、2000 ms。纵横波速比为173。模型纵、横向网格点数分别为301、401，网格间隔为10 m。三层的层厚分别为600、1 800、600 m。采用炸药震源激发，震源位于模型中心，震源函数为雷克子波。图2(a)、(b)为通过有限差分正演模拟得到的钞。与钞：在035 s时

21、刻的波场快照。从图中可以看到，纵横波耦合在一起。采用Helmhohz分解后，可以得到标量势(纵波)、矢量势(横波)，如图2(c)、(d)所示。由于矢量势与地震波的入射到反射界面的几何关系有关，在图2(d)中，可以观察到反射横波在震源的左右两侧以及上下两侧出现极性反转。为了说明波场极性，设定震源左上方的反射波场为负，则其右上方的反射波的极性为正：同样可以看到震源下方左、右反射波场分别为正、负。注意到入射纵波(图2(c)在震源的上、下、左、右区域震源波场具有相同的极性。在传统方法中，直接使用互相关成像，在转换波万方数据第40卷第1期 李振春，等：基于矢量波场分离弹性波逆时偏移成像 。45成像剖面(

22、人射P波与反射s波互相关)中因为s波的极性变化将会出现极性反转。在实际的地震勘苦瑙翳苦恻翳0l230230230l23距离dkin0 1 2 3 4O(a)口。分量距离dkml 2 30(c)标量势9距离dkinl 2探中使用多炮叠加增加信噪比。由于极性反转，多炮叠加时将会损坏同相轴的连续性。距离dkin0 1 2 3 4苦魁翳O230230(b)口，分量距离dkl 2 3 44 00J苦憾聪023(d)矢量势距离dkinl 2 4(e)矢量P波t1分量距离dlkm1 2 3 4 0(g)矢量S波u分量_苦划残023(f)矢量P波坤分量距离dkinl 2 3 4(h)矢量S波埘分量图2三层模型

23、波场035 s时刻的波场快照Fig2 Snapshots of three-layer model at 035 s上盯趟聪上苦燃翳上趟聪万方数据46- 中国石油大学学报(自然科学版) 2016年2月对标量势、矢量势(图2(C)、(d)所示)分别取梯度算子与旋度算子，将得到矢量P波与矢量S波。在二维情况下，矢量P波具有M、W两个分量(图2(e)、(f)，矢量S波也包含、W两个分量(图2(g)、(h)，这里u分量与彬分量分别表示在石与z方向的投影。从图中可以看出，通过梯度算子、旋度算子的处理，纵、横波具有相似的极性分布。具体来讲，矢量P、S波的“分量均关于震源左右对称，矢量P、S波的W分量均关于

24、震源上下对称。在矢量纵横波波场中，波场的大小、极性与入射角的大小有关，与入射方向无关。期望通过内积算子从矢量纵、横波波场中获取反射信息，同时避开转换波成像中的极性反转。3数值模拟为了验证本文成像方法的有效性以及对复杂模距离dkl0 05 10 15 20上墓 -o嫩(a)传统方法型的适应性，本文中分别对平层模型以及Marmousi2模型进行了数值试算。31层状模型为了清楚地展示弹性波成像中的极性反转现象，首先使用简单平层模型进行单炮成像。平层模型纵、横向网格点数均为201，空间采样间隔为10 m。上、下层纵波速度分别为2 500、4000 ms，纵横波速比为173，假定密度为常数。震源纵横向网

25、格位置为(3，101)。由于正演模拟时采用炸药震源激发，只给出PS的成像结果。图3(a)与图3(b)分别为传统方法、本文方法的成像结果。从图中可以看到，传统方法成像结果在震源两侧极性发生了反转：而本文方法在震源两侧具有相同的极性。距离dk=0 05 10 15 20(b)本文方法图3平层模型单炮Ps成像结果Fig3 PS images of layer model with single shot图4为同一平层模型11炮PS成像结果叠加后 234。采样间隔为10 in。整个成像过程中采用全接的剖面。模拟过程中，横向排列网格点数为101，每 收方式，炮点所在网格点横向分布在277877范围一炮排

26、列挪动10个网格点，均采用中心激发的方 内，每炮的网格间隔数为10，共6l炮，采用炸药震式。从多炮的叠加剖面中可以看出，传统方法由于 源激发。成像过程中采用交错网格有限差分法实现极性反转叠加剖面中的同相轴连续性遭到损坏。 波场延拓。采用本文方法避开了极性反转，多炮叠加剖面中同 逆时偏移成像中使用全波动方程模拟波场，可相轴具有较好的连续性，有利于解释人员进行地质 以实现对陡倾角构造的成像，但这也带来低频噪声解释。 的问题。本文中采用Laplacian滤波算子1 9对成像32 Marmousi2模型 结果进行处理，图6(a)一(d)分别为滤去低频噪声图5显示了SEGEAGE提供的复杂的Marmou

27、一 后的PP、PS、sP、ss成像剖面，成像结果表明本文方si2模型的纵波速度。该模型较为复杂，在模型中部 法对复杂构造具有较好的适应性。从图6(a)、(d)有3条断层，中深部有两个背斜，在左上部有一气 中可以看出，采用本文提出的成像方法实现了同类藏。模拟中设定纵横波速度比为173，密度为常 波(PP、ss)对复杂的3条断层、两个背斜和气藏较数。模型横向网格点数为1 134，纵向网格点数为 为清楚的成像。由于在转换波成像过程中避开了极万方数据第40卷第1期 李振春，等：基于矢量波场分离弹性波逆时偏移成像 。47性反转。因此在多炮叠加成像结果中也能清楚地观察到复杂地质构造的形态(图6(b)、(C

28、)所示)。注意到在成像过程中由于使用炸药震源激发，在SP与距离dkm0 05 10 15 20O0510152O005苦 10魁聪15200sS剖面中深部构造的成像没有PP、PS的成像结果清楚。0 05距离dkin10 15 20(a)传统方法 (b)本文方法图4平层模型11炮叠加PS成像结果Fig4 PS images of layer model with 1 1 shots距离dk,2 3 4 5 6 7 8 9 10 11 速度，(-s一1)0苦1瑙酝2苦魁翳00图5 Marmousi2纵波速度模型Fig5 P-wave velocity of Marmousi2 model距离dk, 距离dk=4 6 8 10 0 2 4 6 8O 2(a)PP成像剖面距离dku4 6 8004 0003 0002 0000 2(b)PS成像削面距离dk,4 6 8(c)SP成像剖面 (d)SS成像剖面图6 Marmousi2模型多炮成像结果Fig6 Images of Marmousi2 modelIntd魁翳12InI避联l2上越聪万方数据万方数据