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

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

利用方位角道集处理方法预测煤层裂隙

煤层瓦斯含量受多种因素影响,而煤层中的节理裂隙密集带则容易导致瓦斯聚积。在煤田地质勘探中,依据振幅随炮点和检波点入射角的变化可以预测裂隙的发育方向和发育密度。在对淮南地区张集西三煤矿的三维地震数据资料处理时,按不同的方位角增量抽取方位角道集,并对其方位角道集进行速度分析、NMO校正、叠加和偏移,得到6个方位偏移数据体。根据Malick等人的理论模型对其进行计算,获得d、A、B值,并据此预测裂缝发育方向与发育密度。将其预测成果与常规地震解释成果对比,发现两者之间非常吻合,可见利用方位角道集处理方法预测煤层裂隙具有可行性。     

炮检分离浮动基准面处理方法在复杂山地地震成像中的应用——以准噶尔盆地南缘为例

准噶尔盆地南缘山地地表起伏剧烈、地下构造复杂且速度场纵横向变化大,地震数据原始单炮一致性差、信噪比低,这对山前带起伏地表地震成像造成了巨大挑战.为完善复杂条件下低信噪比山地资料的起伏地表处理技术,以提高山地资料地震成像精度,本文提出了保持炮检点一致性的浮动基准面处理方法.通过配套的时间域处理和深度域速度建模技术,实现了...     

应用折射层析方法解决塔里木盆地浅部地震速度的分布

首次提出了用于解决塔里木盆地地震资料数据处理中浅部地震速度分析的折射层析成象方法。该方法可以通过求解任意变速介质中的二维程函方程、射线追踪和ＳＩＲＴ迭代解法实现。通过对塔里木盆地地震资料的层析反演取得了较好的结果,揭示了浅部地震速度具有较大的横向变化。这一变化对建立地下速度模型及反射构造成象有较大影响,钻井资料证实了这一结果。这项研究的成功对研究塔里木盆地地震波速度的变化规律及进一步提高该区整体地震资料处理与解释的质量具有重要的实际意义。     

一种基于随钻VSP地震地质导向的钻井轨迹高效优化方法及其应用

塔里木盆地发现多个奥陶系碳酸盐岩缝洞型油气田,随着开发的不断推进,中小型储集规模的缝洞体已经成为产能建设的主力目标。勘探开发过程中,由于受到复杂地层速度的影响,缝洞体系偏移成像、归位存在一定误差,常导致钻井的失利。笔者提出一种基于随钻VSP地震地质导向的钻井轨迹高效优化方法,首先应用VSP速度与基于三角网格的新构造模型高效率地校正原始地震偏移速度场,再联合修正后的各向异性参数场对地震数据进行快速偏移成像。其次,对新数据的全层系地层的构造特征再分析,预测优质储层的发育位置,确定最终靶点。最后,根据工程、地质情况,采用最为经济的优化设计钻井轨迹,钻进入靶。应用表明,该技术在大沙漠腹地塔中地区应用取得了较好的效果,地震资料成像品质得到提升。调整轨迹后的储层钻遇率和钻完井直投率分别提升16%和13.2%,有效地促进了油田实施低成本效益开发方案。     

一种适用于欠压实成因的地层孔隙流体压力预测技术

为了提高基于欠压实成因机制的地层孔隙流体压力预测的准确性,这里从提高地震层速度场的计算精度的角度出发,提出了一种高精度地层孔隙流体压力预测技术。①该技术针对浅部地层(无测井曲线的地层段)采用三维约束DIX反演方法和叠加速度体、层位数据计算浅部地层的地震层速度体;②针对中深层地层(有测井曲线的地层段)采用叠后波阻抗反演方法和叠后三维地震数据体、测井数据、层位数据计算中深层地层的地震层速度;③沿着研究区的某一浅部地层分界面将浅部地层与中深层地层的地震层速度体进行拼接处理,得到一个完整的地震层速度体;④采用"Fillippone方法+Eaton方法"的组合方法进行三维地层孔隙流体压力预测。某区页岩气地层的实际计算结果验证了该技术的正确性与可靠性,其预测的地层孔隙流体压力与实钻井一致,且目的层的地层孔隙压力平面分布特征与地质认识一致,有助于指导研究区后续的勘探与开发。     

速度模型对地震波场偏移成像的影响

由于地震偏移速度的不确定性,偏移成像的质量很难得到保证。长期以来,人们对速度误差的分析大多是定性的,而且针对的是较简单的地质构造模型。为了揭示速度模型对地震波场偏移成像的影响特征,提高偏移成像的精度,对一些复杂地质模型进行了波动方程正演和偏移,并从理论模型的数值模拟结果中,推导出了速度模型对地震波场偏移成像的一些影响特征,提出了相应的改进建议和一些有益的结论和认识。     

GRISYS系统速度谱格式转换工具软件开发

本文通过分析ＧＲＩＳＹＳ系统速度谱格式的特点,利用一个特殊自定义格式（ＣＭＰ点号＋ＴＶ对序号＋ＴＶ值）和匹配搜索算法实现了ＧＲＩＳＹＳ系统速度谱数据的自动转换,压缩了现场处理周期,及时提供了现场处理剖面,保证了现场处理的及时质量监督     

反射槽波探测技术中速度分析方法研究

在巷道掘进过程中,利用反射槽波探测技术探查巷道侧帮构造发育情况,并准确解释槽波波阻抗反射界面的位置对矿井生产意义重大,但现阶段反射槽波地震解释断层的位置常存在偏差,主要原因是数据处理阶段选取速度存在较大误差。为降低反射槽波速度选取误差,系统分析国内13个煤矿15个煤层的槽波速度特征,总结出反射槽波速度与直达槽波、透射槽波速度的关系。研究发现反射槽波速度高于直达槽波速度3%~12%;低于透射槽波速度1%~6%。将其研究成果引入山西某矿反射槽波探测项目,槽波解释断层位置与探巷揭露位置偏差仅为3 m,取得了较为理想的应用效果。该研究为进一步提高反射槽波探测精度提供了新思路。      

江陵凹陷北部高速玄武岩下二维广角地震初探

江汉油田江陵凹陷北部高速玄武岩发育,利用常规地震资料很难得到高速玄武岩下的地震反射。分析研究表明,广角地震有利于玄武岩下的地震成像。通过综合正演和实际地震数据完成了广角地震资料的处理,获得了江陵凹陷北部高速玄武岩下的反射地震信息。在开展广角地震波场的识别和利用研究过程中,从正演模型出发,探讨了在不同地质条件下高速层对地震波场的影响,并针对靶区资料讨论了如何识别广角地震信息。在上述基础上,通过对实际地震资料做精细的分析和处理,最终完成了高速玄武岩下地震资料的成像。     

井-震标定参数分析及薄层标定——以塔里木盆地哈得四油田为例

新疆塔里木盆地哈得四（HD4）油田的主要储层东河砂岩属于薄储层,其厚度小于1/4地震波长,根据地震数据难以分辨储层的顶界和底界。为了利用声波速度对薄储层进行有效的井-震标定,从井-震标定的原理出发,以哈得四油田的井-震标定为例,分析了标定过程中子波的主频、类型和相位及声波速度对合成地震道的影响,同时分析了声波速度与叠加速度转换的层速度之间的差别。提出以参考层为标准,并对声波速度进行微调的薄储层的标定方法。参考层在研究区是一套沉积环境相对稳定、厚度较大的地层,它对应于地震剖面上的一个振幅较强、反射轴连续并易于追踪的反射界面。利用上述方法,以哈得四油田二叠系火成岩及石炭系双峰灰岩为参考层,选取与研究区地震数据吻合性最好的子波,合理地调整声波速度,对哈得四油田的主要储层东河砂岩进行井-震标定,获得了较为理想的效果。     

柴达木盆地东部石炭系地震层速度求取方法的研究

在柴达木盆地东部石炭系地震勘探过程中,地震层速度是一个重要的参数。层速度不仅仅可以应用于时深转换,还可以应用于反映地层岩性、构造和地层压力等方面的信息,地震剖面与地震速度也密切相关。柴达木盆地东部石炭系由于资料所限,速度问题一直未能很好地解决,时深转换过程中速度误差较大。因此,为了解决研究过程中速度误差较严重的问题,文章在前人工作基础上,探索性地提出了一种新的层速度求取方法,以提高柴达木盆地东部石炭系地震勘探过程中层速度的精度,从而更好的认识本区的构造和地质特征。通过对单井速度分析和二维地震叠加速度多次拾取,保证了速度谱拾取的精度,然后将地震叠加速度换算成均方根速度,进而利用基于单井速度趋势模型约束Dix反演层速度技术将均方根速度转换成层速度,拟合出全区的平均速度场。同时在拟合后,以井上测井速度为控制点,对平均速度进行了合理校正,从而建立比较合理的研究区层速度模型。利用新建立的层速度模型,对原来的构造重新进行时深转换,取得了较好的效果。新的改进Dix公式层速度求取方法求取的速度模型与前人成果相比较,利用了单井速度约束Dix公式进行速度反演,因而比前人成果更精确,更符合单井速度趋势,因而更符合地质认识规律。     

Analysis of Ground Vibrations Caused by Open Pit Production Blasts - A Case Study

In order to study the characteristics of ground vibrations caused by production blasts in an open pit gold mine and to evaluate the impact of the vibrations to buildings, a seismic survey was conducted at an open pit mine. Two monitoring lines with multiple seismic stations were surveyed. The first line had a length of 4,492 m and the second line runs approximately perpendicular to the first line with a length of 823 m. The seismographs recorded the particle acceleration, velocity, and displacement in longitudinal, transverse and vertical directions. The magnitudes of peak particle motions were calculated and compared with several established damage criteria used in mining and geotechnical/structure engineering. Empirical equations were established based on the field measurements. Analyses were also performed to exam the effects of geological structures on the attenuation of individual vibration components as well as the peak vector sum of the particle velocity.     

胜利探区复杂地质构造的地震波照明与观测系统优化

随着胜利探区勘探开发的逐步深入,主要勘探目标开始转向小断块、小砂体等微幅构造以及岩性隐蔽油气藏,传统的地震采集设计技术难以满足当前生产需求。为此,选取胜利探区复杂构造地质模型,开展基于波动方程理论的地震波照明和面向目标的观测系统优化设计技术的研究,形成了一套适合胜利探区二次采集的观测系统优化设计的技术系列。研究结果表明,该技术系列能够保证该区地质构造的形态完整,小断块、小砂体等微幅构造的成像质量也有较大幅度的提高。     

Statistical Characterization of a Random Velocity Field Using Stacking Velocity Profiles

In petroleum exploration and production, it is important to have good estimations of the uncertainties on the reserves. Uncertainties in the velocity model used during the data processing are of major importance in this estimation. The use of geo-statistical tools can help in dealing with these uncertainties. Up to now, a strong limitation has been the inability to properly merge velocity functions measured in the wells with seismic velocity data. This was due to the different “supports” among the two, i.e. the well velocity may be regarded as a direct measurement of the instantaneous velocity field, while the seismic velocities correspond to an “average along the travelled paths” of this field. The problem is that, apart from the well positions, the instantaneous velocity field is out of reach. Luckily, for many practical applications, it is enough to know just its covariance model. However, no algorithmic method is available in the literature to yield the covariance model, and geologists are forced to use arbitrary distributions. The present paper proposes an original method to obtain a good estimate of this covariance model, using widely available information, mainly seismic stacking velocities. This method was first developed in a simple one-layer case with constant velocity, and then extended to more realistic situations. Finally, a real data application is performed, which highlights the robustness of the resulting estimation.     

二维介质中的空间变量分离速度成像技术

二维介质速度成像的新方法是基于地震波的旅行时间，由射线路径沿线的速度分布特征所决定的，将二维介质中的速度分布用两个分别以二维空间坐标为自变量的勒让德多项式的线性组合来描述。根据走时积分方程，即可建立地震波走时与其射线路径沿线的速度分布之间的对应关系，进而重建出地震波射线覆盖区域的速度结构。该成像技术是一个迭代过程，且适用于折射波及反射波。     

弹性模量、波速与应力的关系及其应用

为充分了解不同尺度地质系统的岩石动力学特征,更好利用不同规模勘探资料有效解释地质系统的特征,在岩石物理实验对岩石弹性模量、波速与应力的关系进行研究总结的基础上,探讨野外实际地震勘探中地震波速与地应力的关系。研究表明,岩石、岩体或地质系统是多相(固液气)、多种矿物、各种构造的综合体,不论尺度如何均可以用多重介质孔隙裂隙地质模型描述,而波速是一定尺度不均匀地质体的等效波速,在一定地质条件下,各种尺度介质波速、岩石弹性模量与有效应力存在依存关系。坚硬的岩石通常很脆,如果岩石处于很大的应力之下,应力的释放将在岩石中产生微裂隙,微裂隙会降低地震速度。只要把岩石恢复到原来的应力状态就可能消除所引发的微裂隙。在岩体未达到破坏时,应力越大,波速愈高。压力释放引发的微裂隙有助于现场确定最大/最小水平应力方向,波速的椭圆长轴即为主应力方向。而上覆地层压力不变(或变化很小),净储层压力的影响正好与储层压力对地震特性的影响相反。随着净储层压力的增加,所有岩石的地震特性都增高。这种增高的量级取决于若干其它的因素(孔隙形状、孔隙度、孔隙流体、岩性等等)。不论野外地震勘探还是实验研究,地震波速与地应力的一致性为地质系统的正反演解释提供了宏观控制依据:通过测量弄清     

A Simple Algorithm for Sequentially Incorporating Gravity Observations in Seismic Traveltime Tomography

The geologic structure of the Earth's upper crust can be revealed by modeling variation in seismic arrival times and in potential field measurements. We demonstrate a simple method for sequentially satisfying seismic traveltime and observed gravity residuals in an iterative 3-D inversion. The algorithm is portable to any seismic analysis method that uses a gridded representation of velocity structure. Our technique calculates the gravity anomaly resulting from a velocity model by converting to density with Gardner's rule. The residual between calculated and observed gravity is minimized by weighted adjustments to the model velocity-depth gradient where the gradient is steepest and where seismic coverage is least. The adjustments are scaled by the sign and magnitude of the gravity residuals, and a smoothing step is performed to minimize vertical streaking. The adjusted model is then used as a starting model in the next seismic traveltime iteration. The process is repeated until one velocity model can simultaneously satisfy both the gravity anomaly and seismic traveltime observations within acceptable misfits. We test our algorithm with data gathered in the Puget Lowland of Washington state, USA (Seismic Hazards Investigation in Puget Sound [SHIPS] experiment). We perform resolution tests with synthetic traveltime and gravity observations calculated with a checkerboard velocity model using the SHIPS experiment geometry, and show that the addition of gravity significantly enhances resolution. We calculate a new velocity model for the region using SHIPS traveltimes and observed gravity, and show examples where correlation between surface geology and modeled subsurface velocity structure is enhanced.     

低信噪比地震资料处理中的精细速度分析方法

为了提高低信噪比资料速度分析的精度,得到中深层可靠的速度信息。针对关中地区深部隐伏断裂构造调查的地质目标,在二维地震勘探的基础上,对低信噪比地震资料采用"优选道集速度分析、动校正、F-K滤波、反动校正、高密度精细速度分析"的技术流程,使得速度谱中能量团聚集,不仅提高了低信噪比地震资料速度分析的精度,而且达到了提高信噪比的效果,进而获得了较高品质的中深层速度场信息与成像剖面资料。      

双谱速度分析下的地层压力预测

为提高依靠地震速度谱数据进行压力预测的准确度,提出将主成分分析( PCA) 判断准则应用到双参数速度扫描中来拾取地震速度,从而减少扫描数据时随着反射层数的增加及偏移距的增大而造成的偏差,并将所得速度参数及时间参数应用到改进的Fillippone 公式中,改善了压力预测的精度。