Directional‐oriented wavefield imaging: a new wave‐based subsurface illumination imaging condition for reverse time migration

The key objective of an imaging algorithm is to produce accurate and high‐resolution images of the subsurface geology. However, significant wavefield distortions occur due to wave propagation through complex structures and irregular acquisition geometries causing uneven wavefield illumination at the target. Therefore, conventional imaging conditions are unable to correctly compensate for variable illumination effects. We propose a generalised wave‐based imaging condition, which incorporates a weighting function based on energy illumination at each subsurface reflection and azimuth angles. Our proposed imaging kernel, named as the directional‐oriented wavefield imaging, compensates for illumination effects produced by possible surface obstructions during acquisition, sparse geometries employed in the field, and complex velocity models. An integral part of the directional‐oriented wavefield imaging condition is a methodology for applying down‐going/up‐going wavefield decomposition to both source and receiver extrapolated wavefields. This type of wavefield decomposition eliminates low‐frequency artefacts and scattering noise caused by the two‐way wave equation and can facilitate the robust estimation for energy fluxes of wavefields required for the seismic illumination analysis. Then, based on the estimation of the respective wavefield propagation vectors and associated directions, we evaluate the illumination energy for each subsurface location as a function of image depth point and subsurface azimuth and reflection angles. Thus, the final directional‐oriented wavefield imaging kernel is a cross‐correlation of the decomposed source and receiver wavefields weighted by the illuminated energy estimated at each depth location. The application of the directional‐oriented wavefield imaging condition can be employed during the generation of both depth‐stacked images and azimuth–reflection angle‐domain common image gathers. Numerical examples using synthetic and real data demonstrate that the new imaging condition can properly image complex wave paths and produce high‐fidelity depth sections.     

Multipathing in three‐parameter common‐image gathers from reverse‐time migration

Reverse‐time migration gives high‐quality, complete images by using full‐wave extrapolations. It is thus not subject to important limitations of other migrations that are based on high‐frequency or one‐way approximations. The cross‐correlation imaging condition in two‐dimensional pre‐stack reverse‐time migration of common‐source data explicitly sums the product of the (forward‐propagating) source and (backward‐propagating) receiver wavefields over all image times. The primary contribution at any image point travels a minimum‐time path that has only one (specular) reflection, and it usually corresponds to a local maximum amplitude. All other contributions at the same image point are various types of multipaths, including prismatic multi‐arrivals, free‐surface and internal multiples, converted waves, and all crosstalk noise, which are imaged at later times, and potentially create migration artefacts. A solution that facilitates inclusion of correctly imaged, non‐primary arrivals and removal of the related artefacts, is to save the depth versus incident angle slice at each image time (rather than automatically summing them). This results in a three‐parameter (incident angle, depth, and image time) common‐image volume that integrates, into a single unified representation, attributes that were previously computed by separate processes. The volume can be post‐processed by selecting any desired combination of primary and/or multipath data before stacking over image time. Separate images (with or without artifacts) and various projections can then be produced without having to remigrate the data, providing an efficient tool for optimization of migration images. A numerical example for a simple model shows how primary and prismatic multipath contributions merge into a single incident angle versus image time trajectory. A second example, using synthetic data from the Sigsbee2 model, shows that the contributions to subsalt images of primary and multipath (in this case, turning wave) reflections are different. The primary reflections contain most of the information in regions away from the salt, but both primary and multipath data contribute in the subsalt region.     

人工源宽角地震实验的谱元法模拟和逆时偏移成像方法研究

人工源宽角反射(折射)地震资料具有偏移距较大、信噪比较低等特点,通常用于地震波走时反演重建地壳速度结构。逆时偏移成像方法作为勘探地震学领域获取地下构造形态的有效手段之一,可以有效弥补走时反演方法的不足。本文针对大偏移距宽角反射(折射)地震实验,利用四边形网格谱元法进行波场模拟,结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度获取模型合成地震记录,后采用逆时偏移成像方法将合成地震记录偏移归位,获取地壳几何结构,验证了逆时偏移成像方法在宽角地震资料处理及结果解释中的适用性,为后期实际地震资料的偏移成像提供了理论依据和支持。     

Fluid flow reconstruction in hanging and footwall carbonates: Compartmentalization by Cenozoic reverse faulting in the Northern Oman Mountains (UAE)

In this paper, the diagenesis from either side of a major Cenozoic reverse fault in the Northern Oman Mountains is documented. Detailed petrographical and geochemical analysis of calcite-filled fractures in carbonate strata of Late Triassic and Early Cretaceous age in the hanging wall and footwall in Wadi Ghalilah reflect a different diagenetic history. In both hanging wall and footwall most of the fractures are pre-burial, extensional in origin, formed by a crack-seal mechanism, and the calcite vein infill has a host-rock buffered signature. In the hanging wall, the fluid responsible for calcite precipitation of these extensional fractures was a marine fluid at 60 °C. These veins predate deep burial and contractional tectonic deformation and consequently do not provide any information about syntectonic fluid flow. Neither do the pre-burial extension fractures in the footwall which are also host-rock buffered. The fractures post-dating the tectonic stylolitization in the footwall, by contrast, show evidence of syntectonic migration of saline formation waters at temperatures between 80 and 160 °C during contractional deformation. These fluids probably were sourced from the subsurface via the reverse fault, which acted as a fluid conduit. At the same time, however, this fault functioned as a permeability barrier towards the hanging wall, since no evidence of syntectonic fluid flow is present here. In this way compartmentalization of the hanging wall and footwall block was realized.     

CT在航空深弹结构分析和拆卸中的应用研究

通过对逆向工程研究方法的分析,提出了利用工业CT断层扫描对航空深弹进行结构分析,精确测量其零部件位置、尺寸变化的方法。并据此绘制图纸、制定拆卸方案,从而保证了深弹拆卸和研究的安全性和可靠性。     

逆序断裂的发育特征与地震反射识别

在具有盆地基底和未固结盖层的双层结构的盆地里,由于松散盖层对基底断裂应变的吸收,晚期发育的基底断裂在进入盖层后,应变会逐渐减小甚至消失,形成下部错动大、上部错动小的特点,本文把这样的后期断裂称之为逆序断裂。当基底断裂为走滑断裂时,逆序断裂具有双层空间构造样式,即在盆地基底发育主干走滑断裂、而在盖层优先发育另外走向的次级雁行断裂。逆序断裂在渤海湾盆地普遍发育,本文以歧口凹陷的张北断裂带和白水头断裂带为例(这两处断裂带都在下构造层发育北东向基底走滑断裂、而在上构造层发育近东西向盖层次级正断裂),运用三维地震资料,结合沿层属性和垂直剖面,阐述了逆序断裂的发育特征和对地震反射剖面的识别。由于之前晚期的逆序断裂一般被当作早期同沉积生长断裂对待,本文希望能帮助重新认识渤海湾地区的区域断裂期次和构造发育演化过程。     

海拉尔盆地乌尔逊凹陷中的三角形断块：演化及对油气的控制作用

在乌尔逊凹陷中部发育一个三角形断块,因边界断裂的同沉积活动而使断块沉陷深,沉积厚度大.晚期的构造反转使整个三角形断块向上冲断,尤以北界断裂的冲断活动表现强烈.三角形断块为各时期的沉积中心,控制着深湖相和湖底扇沉积的发育;烃源岩的类型也与三角形断块有关,块内为Ⅰ型和Ⅱ型干酪根,边部和外部为Ⅱ型和Ⅲ型十酪根.南屯组的有效烃源岩在乌尔逊断陷分布范围较大,主力烃源岩只局限于三角形断块内.晚期反转时北界断裂有利于构造圈闭的形成,反转断层也为油气的运移提供了通道.因此,三角形断块在生储盖层的发育和油气生成、运移与聚集的空间配置关系上可谓得天独厚,是最有利的油气聚集区.     

近百年来南京降水变化的趋势和特征

本文对近百年来（１９０５－１９９４年）南京实测降水量资料序列和降水的其他统计特征作了详尽的诊断分析，得到如下十分有益的气候变化信息。     

逆断层对致密岩石构造裂缝发育的约束控制

致密岩石(如火山岩和碳酸盐岩)可广泛发育构造裂缝,这类致密岩石的构造裂缝数值建模是构造裂缝地震解释、地震反演和油藏数值模拟的基础.我国西北各盆地普遍发育逆断层,通过在新疆巴楚地区微波塔逆断层和一间房逆断层分别布置实测剖面对野外测量的构造裂缝面密度进行定量分析,发现受逆断层控制的裂缝面密度与逆断层面的曲率呈线性关系;裂缝...     

易门三家厂铜矿深部采区涌水量预测

利用上部开采坑道涌水量记录资料,对深部新采区水文地质模型进行合理概化,用上部涌水量反算矿床充水层的渗透系数,进而求算坑道涌水量。通过实践证实,预测结果正确。