黏滞声波高斯束叠前深度偏移

高斯束偏移不仅具有接近于波动方程偏移的成像精度,而且保留了Kirchhoff 积分法高效、灵活的优点,可以对复杂介质准确成像。由于实际地下介质具有黏滞性,因此研究黏滞声波叠前深度偏移具有一定的现实意义。笔者采用高斯束偏移方法对地震数据进行吸收衰减补偿。首先给出共炮域高斯束叠前深度偏移原理;然后在此基础上推导补偿吸收衰减的表达式,校正品质因子Q引起的振幅衰减和相位畸变,实现基于吸收衰减补偿的高斯束叠前深度偏移;最后用两层模型和气云模型对偏移方法进行了测试。结果表明,在考虑地下介质的黏滞性时,黏滞声波高斯束叠前深度偏移比声波高斯束叠前深度偏移具有更高的成像分辨率。     

基于方位—反射角度道集的高斯束层析速度建模方法及实现

层析反演是速度建模中最重要的方法之一,结合偏移成像在成像域进行波动方程线性化走时层析速度建模是当前比较实用有效且精度较高的技术组合。文中首先给出了高斯束偏移提取方位—反射角度道集的方法,之后从高斯束偏移角度道集出发,在波动方程的一阶Born近似和Rytov近似下,推导了成像域波动方程线性化走时层析方程及其显式表达的层析核函数,并利用高斯束传播算子计算该核函数。基于高斯束传播算子的偏移成像与层析成像相结合进行深度域速度建模迭代及偏移成像,体现了速度建模与成像一体化的思想。数值计算及实际数据应用证明了基于高斯束传播算子的层析成像与偏移成像方法的有效性。     

曙光4000A上三维叠前深度偏移并行计算的应用设计

波动方程叠前深度偏移在地震勘探成像处理方面起着不可替代的作用。随着高性能大规模并行计算机技术的发展,波动方程叠前深度偏移计算在地震勘探中的应用有了很大进步。在波动方程叠前深度偏移处理中,庞大的数据规模与海量计算对计算性能提出了很高的要求。曙光4000A超级计算机系统是我国目前峰值速度最快的商用超级计算机系统,无论是硬件平台建设还是应用软件的配置方面,都具有良好的应用性能。基于该系统设计的三维波动方程叠前深度偏移(炮域)PSDM软件,采用动态负载平衡并行计算模式,具有较高的计算效率,高度的可扩展性和可靠性。     

泌阳凹陷南部陡坡带地震资料的炮域波动方程叠前深度偏移

泌阳凹陷南部陡坡带是河南油田的重点勘探老区,其构造产状不仅陡,断层也多,叠后时间偏移难以得到较好的地震成像效果,基于Kirchhoff积分偏移方法具有运算速度快、效率高的特点,以及炮域波动方程方法能保持地震波的动力学特征,应用Kirchhoff积分偏移方法进行了偏移速度分析,建立了精细的速度模型.然后基于该模型在深度域进行了剩余速度分析以及层析成像,进一步提高了速度模型的精度. 通过分析偏移孔径、去假频参数以及延拓步长等成像处理参数对成像效果的影响,确定出最佳偏移处理参数,最后利用炮域波动方程方法对工区的三维地震数据进行了叠前偏移成像.成像结果表明,该方法能够使该区陡坡带反射层位得到较好的偏移,信噪比和横向分辨率都得到提高.     

波动方程角度域共成像道集

针对基于波场深度递推的波动方程叠前深度偏移成像方法难以准确地给出如Kirchhoff积分叠前深度偏移成像方法的偏移距域共成像道集的不足,在分析研究现有的各种波动方程偏移成像共成像道集方法的基础上,利用波场外推的单平方根算子和波场的窗口Fourier框架展开与重构方法,提出一种局部角度域共成像道集方法——成像点处反射角共成像道集方法.把这种角度域共成像道集方法应用于国际标准的Marmousi模型数据和一条实际二维地震数据都取得了理想的结果.提出的波动方程反射角度域共成像道集方法可为进一步的叠前偏移数据振幅随角度变化分析和偏移速度分析工作提供基础.     

基于深度加权的三维地震斜率层析成像

地震斜率层析成像是一种联合使用反射波局部相干同相轴的走时和斜率来建立宏观速度模型的有效方法。对于低信噪比资料,往往能够拾取的深层有效反射波数据远远少于浅层反射波数据,使深部地层速度的层析反演效果差。因此,提出了一种基于深度加权的三维地震斜率层析成像方法。在每次迭代的线性层析反演方程中,对观测数据关于离散模型节点速度的核函数进行深度加权,权系数根据离散模型节点深度和当前迭代的各个炮—检对的反射点深度来确定,实现了加大深部反射波数据对深层速度的约束作用、浅层速度主要由浅部反射波数据约束的目的,在保持浅层速度反演精度的同时,改善对深层速度的反演效果。对理论模型数据和实际资料应用测试的良好效果,表明了该方法的有效性。     

共偏移距可分表示法叠前深度偏移

基于双平方根方程的共偏移距可分表示法叠前深度偏移用于复杂介质成像, 该方法在中点-偏移距坐标中同时向下延拓炮点和检波点波场, 实现方式采用正反傅立叶变换.构造的双平方根方程波场延拓算子能够使波数域变量与空间(速度) 域变量分离, 波数域内进行相移计算, 在空间域对因介质横向变速引起的时移作修正.地震数据偏移不需要逐炮计算, 具有较高的效率.对Marmousi模型数据的偏移成像结果显示, 该方法较好地成像强横向变速介质中的复杂构造.      

基于保幅波动方程的广义高阶屏地震偏移成像方法

原始地震数据中饱含丰富的反映相位的走时信息和反映反射率的振幅信息。基于波场延拓理论的波动方程保幅地震偏移成像,是在给出正确构造成像位置的同时也给出真实反射振幅的有效完善。基于全声波方程,利用严格的解耦理论进行单程波动保幅分解,得到一个由波场传播项与振幅补偿项构成的,在走时与振幅上满足全声波方程对应的程函方程与输运方程的保幅单程波动方程;利用摄动理论进行单平方根算子渐进展开,推导出基于保幅波动方程的广义高阶屏地震偏移算子方程。模型测试和实际资料处理表明,该方法不但可以凭借更准确的相位归位和散射能量聚焦提高构造成像精度,而且输出了能更正确反映地下反射属性的能量信息,从而可以为更深层次的勘探开发,提供地球物理技术支撑。     

转换波波动方程叠前深度偏移

克希霍夫叠前时间偏移被广泛运用于转换波地震成像领域,该方法具有计算效率高、速度模型易于获取的优点。但是,纵横波速度差使纵波和转换波的成像在时间域不匹配,且克希霍夫偏移也无法较好地处理多次到达和保幅等问题。这里提出一种基于波动方程的转换波叠前深度偏移方法,以改善克希霍夫时间偏移的不足。     

波动方程辛几何算法三维叠前深度偏移在集群式并行机上的实现

波动方程三维叠前深度偏移研究,已成为近年国内外众多油气勘探科研单位的研究焦点,其关键问题之一是借助于软硬件及其匹配的进步,以求解决计算量庞大的障碍。作者在本文中从油田实际出发,提出了一套并行波动方程辛几何算法三维叠前深度偏移方法与处理流程,运用32节 点集群式并行机,完成了XJWZ地区150km的二维和WJT地区126km^2的三维地震资料波动方程叠前深度偏移处理。其成像的效果,在信噪比、分辨率、深层复杂构造成像等方面均获得了显著改善。     

基于CMP域的山地表层速度模型构建

对于山地地震勘探而言,地形起伏,不能直接使用Herglotz-Wiechert公式反演近地表速度模型,这里采用不同于前人在炮域处理的思路,而是在CMP域建立了一组基于CMP浮动基准面对起伏地表上的炮检对走时和偏移距进行校正的公式,首先用本文提出的公式去除起伏地表对初至波的影响,然后应用Herglotz-Wiechert公式快速反演地形起伏地区近地表速度模型。该方法采用地震反射波的CMP浮动基准面校正思路,包含三个关键步骤:①产生CMP浮动基准面;②实测时距曲线的校正;③计算射线参数和速度-深度函数的计算。此方法无需射线追踪,无需初始模型,计算速度快,反演得到的模型具有和地震数据相同的水平分辨率,该模型可用于计算基准面静校正量,或用于叠前深度偏移速度场的表层建模,也可以作为回转初至波层析反演或全波形反演的初始模型。     

波动方程叠前深度偏移的GPU技术

偏移成像是地震资料处理流程中重要的一环,目前工业界所采用的方法是克希霍夫时间偏移,该方法的优点在于计算效率高,模型易于获取。但因其基于高频近似,难以适应多次到达,欠照明影响及横向变速剧烈的地质情况。波动方程叠前深度偏移则能较好地解决这些问题,但其计算成本远高于克希霍夫偏移方法,这也成为其在实际生产应用中的瓶颈之一。这里将利用GPU技术,研究波动方程叠前深度偏移的高性能计算方法,使之能够满足生产中海量数据处理的需要。     

基于PC机群波动方程叠前深度偏移的并行计算策略

基于波动方程的叠前深度偏移技术是解决复杂地质地貌地区地震成像的一种非常有效的手段,但是在实际应用中面临着数据量巨大和计算量巨大的双重困难.以PC 机集群为硬件环境,以MPI消息传递并行编程环境为并行程序设计平台,研究设计出请求分配作业的主从模式来实现波动方程叠前深度偏移的并行计算,从而解决了动态负载均衡难题,采用的作业登记与分析技术解决了容错处理问题.理论模型数据和实际地震资料测试结果表明:程序运行稳定,并行效率高.对我国东部地区某三维地震数据进行了处理,用请求分配作业的主从模式比平均分配作业的主从模式节省了23.68%的时间.     

煤矿采区高密度三维地震深度域资料解释方法

时间域地震资料解释比较成熟,深度域处理技术已经走向煤炭领域,但煤矿采区高密度三维地震深度域资料解释实际应用还存在很多问题。通过煤炭地震深度域层位标定、深度域断层解释、深度域底板成图的摸索应用,参考时间域解释的流程,初步建立在煤炭高密度三维地震深度域资料中直接解释煤田地质成果的方法。以淮北祁南矿三维叠前深度偏移地震资料应用为例,通过深度域、时间域地震数据的对比剖析,断层解释、回采面地震属性显示及底板成图,取得精度更高的结果,利于一线技术人员直接运用深度域地震资料来指导煤矿生产。      

基于角度域共成像点道集的微分相似优化法偏移速度分析

微分相似优化(Differential Semblance Optimization,DSO)法是一种可以有效避免非线性最小二乘反演过程中不聚焦问题的速度反演方法。基于DSO原理的波动方程偏移速度分析将成像道集的聚焦或者拉平程度作为速度判定的准则,其目标泛函具有良好的凸性,可以有效避免局部极值问题,目标泛函的梯度较为光滑,并且对地震数据中低频成分的要求较低,对速度模型中背景速度的估计较为准确。一般情况下,DSO方法利用偏移距域共成像点道集建立目标泛函,文中利用角度域共成像点道集建立目标泛函对速度的准确性进行判定。模型试算结果表明,通过角道集可以更加直观地对速度的准确性进行判定,并且准确反映速度与深度的对应关系。     

叠前地震全波形反演实际应用的关键影响因素

叠前地震全波形反演是利用波动方程正演模拟来描述地震波传播特征,通过非线性局部寻优方法来获取高精度反演结果的方法。与其他反演方法相比,全波形反演利用的是叠前炮域地震数据,没有经过叠加处理,更好地保留了振幅、相位、频率等波形参数随偏移距变化的特征。在反演结果中得到的不再是简单的界面和旅行时信息,而是反映地层速度和构造特征的深度域数据体,有助于对地质体的分析和认识。从叠前全波形反演方法的优势出发,通过理论模型的全波形反演,分析做好全波形反演的技术关键,然后将全波形反演方法应用于实际资料处理,并将反演结果应用于偏移处理和河流相储层识别中,取得了良好的应用效果。最后针对全波形反演实际资料应用中应该注意的问题进行了总结,为叠前全波形反演的进一步发展和应用提供借鉴。     

任意广角方程逆时偏移的脉冲响应及模型试算

采用高阶差分法求解任意广角声波方程,推导了各向同性介质中任意广角声波方程叠前逆时深度偏移的差分算子,给出了适用于任意广角声波方程逆时延拓的吸收边界条件和其计算方法,在此基础上实现了任意广角声波方程的高精度逆时延拓。分析了不同模型条件下任意广角声波方程叠前逆时深度偏移的脉冲响应,分析结果表明,任意广角方程逆时偏移算子可压制层间反射波,减少偏移噪声,理论模型的成像结果表明该方法能够对复杂构造模型精确成像。     

基于单程波波动方程的起伏地表叠前深度偏移技术

直接从起伏地表开始进行叠前深度偏移的波场初始化时,应用了同步边界处理策略和深度误差校正策略。前者可以消除起伏地表非波场传播区域的成像,后者可以消除起伏地表规则采样时深度误差对地下构造成像带来的影响。利用以上处理策略,优化了常规基于单程波波动方程延拓算子的起伏地表叠前深度偏移算法,从而更好实现了复杂条件下的构造成像。基于FFD延拓算子的波场下延法,用SEG起伏地表模型数据验证了该方法的有效性。该方法同样适用于各种基于频率域波动方程深度延拓算子的起伏地表深度偏移方法。      

地震数据炮域规则化方法研究及应用

叠前深度逆时偏移成像技术是目前处理复杂高陡构造成像的有效方法之一,作为一种高精度的炮域波动方程偏移方法,需要规则采集的高质量炮集数据作为输入。受复杂地质条件的限制,实际采集到的地震数据在空间方向往往是稀疏不规则采样的,存在严重的空间假频和噪声干扰,不满足炮域偏移算法对数据的要求。为了提高逆时偏移成像精度,笔者开展了叠前地震数据炮域规则化方法研究,依据偏移距和方位角将炮集数据分选为OVT(offset vector tile)道集,在OVT域进行地震数据全波数带的迭代加权反假频插值,再基于几何观测系统中炮检点坐标映射关系,从OVT域插值结果中提取规则化后的炮集数据并应用于逆时偏移成像。实际资料的应用结果表明对地震数据进行叠前炮域规则化处理可以有效提高逆时偏移的成像精度。     

叠前深度偏移技术在QMQ地区成像中的应用

位于东部QMQ地区的地震资料,其地表地质条件较好,但因勘探目的层是奥陶系复杂的古潜山,因此,适合应用叠前深度偏移技术进行潜山顶及其内幕的成像。通过Kirchhoff叠前深度偏移技术中精细的数据准备、速度-深度模型的建立与优化、偏移孔径的试验与选择和偏移成像4个关键环节的应用研究,完成了QMQ地区地震资料的Kirchhoff叠前深度偏移。通过与叠前时间偏移剖面对比,叠前深度成像改善了QMQ地区地震数据的信噪比和分辨率,同时提高了奥陶系潜山及其内幕的成像精度,为后续的深度域地质解释和构造成图提供了可靠的偏移数据,研究成果对深层煤勘探所面临的底板奥灰水的研究具有借鉴意义。