Impact of user parameters in inversion velocity analysis

Migration velocity analysis is a method devoted to the evaluation of both reflectivity and background velocity models, associated with the high and low wavenumber components of the model, respectively. Inversion velocity analysis is one of its improved versions, leading to more stable background velocity updates. Still, the impact of the user parameters should be understood for an optimal update of the background velocity. We show that a sign reversal of the background velocity gradient could occur when the selected surface offset range or the space lag range is too small. We derive the theoretical limits and check their consistency through simulations in a simple model with a single interface. These guidelines determine the necessary ranges of surface offsets and space lags for a proper update of the background velocity model. We discuss their applicability on the Marmousi model. Artefacts in the retrieved background velocity model are observed when the guidelines are not satisfied.     

全干涉成像的微地震定位方法研究

基于偏移成像的微地震定位方法由于可避免走时拾取误差以及可实现自动化定位等优点被广泛应用.绕射叠加方法将能量沿走时曲线聚焦到空间网格点上,进行成像时需要搜索发震时刻.干涉成像方法利用互相关提取的走时差信息可避免搜索发震时刻,但定位结果受数据信噪比的影响较大.为了进一步提高干涉成像法的定位精度,本文提出一种同时使用互相关和自相关道集的全干涉成像方法,增加的自相关道集提取的S-P走时差可降低震源-检波器方向的定位误差,提高定位精度.单井监测的理论测试和实际数据的研究结果表明,全干涉成像的定位精度高于仅使用互相关道集的干涉成像方法,同时计算效率高于绕射叠加方法.     

Effects of aperture on Marchenko focussing functions and their radiation behaviour at depth

A focussing function is a specially constructed field that focusses on to a purely downgoing pulse at a specified subsurface position upon injection into the medium. Such focussing functions are key ingredients in the Marchenko method and in its applications such as retrieving Green's functions, redatuming, imaging with multiples and synthesizing the response of virtual sources/receiver arrays at depth. In this study, we show how the focussing function and its corresponding focussed response at a specified subsurface position are heavily influenced by the aperture of the source/receiver array at the surface. We describe such effects by considering focussing functions in the context of time-domain imaging, offering explicit connections between time processing and Marchenko focussing. In particular, we show that the focussed response radiates in the direction perpendicular to the line drawn from the centre of the surface data array aperture to the focussed position in the time-imaging domain, that is, in time-migration coordinates. The corresponding direction in the Cartesian domain follows from the sum (superposition) of the time-domain direction and the directional change due to time-to-depth conversion. Therefore, the result from this study provides a better understanding of focussing functions and has implications in applications such as the construction of amplitude-preserving redatuming and imaging, where the directional dependence of the focussed response plays a key role in controlling amplitude distortions.     

坐骨股骨撞击综合征的MRI表现

目的:总结坐骨股骨撞击综合征(IFIS)的MRI表现,提高对该病的认识。资料与方法:回顾性分析经我院确诊的IFIS 84个髋关节(病例组)及正常志愿者髋关节116个(对照组)的MRI图像,通过轴位T1WI序列测量两组人群坐骨股骨间隙(IFS)、股方肌间隙(QFS)、坐骨结节间距,并通过质子加权脂肪抑制序列观察股方肌形态及信号。结果:独立样本t检验分析表明对照组IFS及QFS((2.36±0.47)cm、(1.58±0.39)cm)明显大于病例组((1.56±0.42)cm、(0.89±0.24)cm),对照组坐骨结节间距((10.51±1.53)cm)小于病例组((11.39±1.01)cm),两组之间差异均有统计学意义(P <0.05)。Pearson相关分析表明髋关节IFS与QFS之间呈显著正相关,r=0.813,P <0.05;坐骨结节间距与IFS之间存在负相关关系,r=-0.263,P <0.05。在IFIS患者中,股方肌均有不同程度的变形及水肿。所测全部数据中,中青年组IFS、QFS宽度((2.09±0.55)cm,(1.32±0.47)cm)均大于老年组((1.90±0.58)cm,(1.19±0.39)cm),两组之间差异有统计学意义(P <0.05);而两组之间坐骨结节间距无显著性差异。结论:IFS和QFS变窄,同时股方肌出现变形、水肿,是坐骨股骨撞击综合征的典型MRI表现。     

基于两步弥散算子的共反射面元射线束成像方法研究

基于克希霍夫统一成像理论,对共反射面元射线束成像方法进行了深入分析.通过这种分析可以得到一种新的共反射面元射线束类成像实现方式:基于两步弥散算子的共反射面元射线束类成像方法.当仅实施第一步弥散时,即可获得信噪比大幅提高同时更为规则的叠前道集.二维、三维的理论和实际数据算例证实了上述观点.     

青藏高原岩石圈三维电性结构

本文报道通过综合大地电磁调查数据研究青藏高原岩石圈三维电阻率模型的初步成果.大地电磁法调查区域已经覆盖了高原大部分面积,为全区三维电阻率成像研究打下了可靠的基础.对多个测区大地电磁数据进行精细的同化处理和反演成像,取得了青藏高原可靠的岩石圈三维电阻率结构图像.成像的区域为28°N—35°N,80°E—104°E.三维反演计算时采用的网格尺寸为20km×20km,垂直方向不等间距剖分为26层.结果表明,青藏高原现今岩石圈电阻率扰动主要反映印度克拉通对亚欧大陆板块俯冲引起的热流体运动和大陆碰撞和拆离产生的构造.在岩石圈地幔,察隅地块、喜马拉雅地块和拉萨地块东部联成统一的高电阻率地块,它们反映了向北东俯冲的印度克拉通.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,表明岩石圈深处有热流体活动.雅鲁藏布江、班公—怒江和金沙江缝合带都有明显的低电阻率异常,也表明它们的岩石圈还有流体活动.青藏高原东部的低阻区沿100°E向地幔下方扩大,反映了金沙江断裂带有切穿岩石圈的趋势.地幔电阻率平面扰动的模式显示,青藏高原东西部的地体碰撞拼合形式和方向是不同的.在青藏高原西部,羌塘、拉萨和喜马拉雅等地体从北到南碰撞拼合.在青藏高原东部,羌塘—拉萨、察隅、印支、雅安和扬子等地体多方向拆离拼合,在地壳造成不正交的拆离带和压扭构造系.从高阻-低阻区的分布看,东部的地体拼合有地幔的根源,今后还会进一步发展.察隅地块岩石圈对青藏高原东部的楔入,使其北部和东部地块的岩石圈发生拆离撕裂,也造成热流体上涌的低电阻率异常.     

逆时偏移在声反射成像测井中的应用方法研究

声反射成像测井中常用的基于射线理论的Kirchhoff积分偏移算法和基于单程波理论的F-K偏移算法均可实现井旁缝洞反射体的快速偏移成像,但其仅适用于地层垂向变化较弱的速度场和高陡角裂缝的偏移成像,无法实现低角度反射体的准确偏移归位,产生偏移假象误导测井解释.逆时偏移基于全波动方程,可适应强垂向变化速度场,实现近似水平反射体的偏移成像.本文详细分析了将逆时偏移应用于声反射成像测井时存在的数据准备、时间采样间隔匹配和成像条件改进等若干问题,通过设置多组理论模型来说明算法对井旁不同反射体的识别能力.模拟资料和实际资料处理结果证实,较F-K偏移算法,逆时偏移算法成像精度更高、收敛性更好,可有效实现近似水平构造偏移归位.改进的归一化互相关成像条件可解决深部地层的远井壁成像衰减问题,降低测井解释的多解性.逆时偏移将成为声反射成像测井高精度偏移技术的发展方向.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

Deep massive sulphide exploration using 2D and 3D geoelectrical and induced polarization data in Skellefte mining district,northern Sweden

Geoelectrical and induced polarization data from measurements along three profiles and from one 3D survey are acquired and processed in the central Skellefte District, northern Sweden. The data were collected during two field campaigns in 2009 and 2010 in order to delineate the structures related to volcanogenic massive sulphide deposits and to model lithological contacts down to a maximum depth of 1.5 km. The 2009 data were inverted previously, and their joint interpretation with potential field data indicated several anomalous zones. The 2010 data not only provide additional information from greater depths compared with the 2009 data but also cover a larger surface area. Several high‐chargeability low‐resistivity zones, interpreted as possible massive sulphide mineralization and associated hydrothermal alteration, are revealed. The 3D survey data provide a detailed high‐resolution image of the top ～450 m of the upper crust around the Maurliden East, North, and Central deposits. Several anomalies are interpreted as new potential prospects in the Maurliden area, which are mainly concentrated in the central conductive zone. In addition, the contact relationship between the major geological units, e.g., the contact between the Skellefte Group and the Jörn Intrusive Complex, is better understood with the help of 2010 deep‐resistivity/chargeability data. The bottommost part of the Vargfors basin is imaged using the 2010 geoelectrical and induced polarization data down to ～1‐km depth.     

Point‐spread functions for interferometric imaging

Interferometric redatuming is a data‐driven method to transform seismic responses with sources at one level and receivers at a deeper level into virtual reflection data with both sources and receivers at the deeper level. Although this method has traditionally been applied by cross‐correlation, accurate redatuming through a heterogeneous overburden requires solving a multidimensional deconvolution problem. Input data can be obtained either by direct observation (for instance in a horizontal borehole), by modelling or by a novel iterative scheme that is currently being developed. The output of interferometric redatuming can be used for imaging below the redatuming level, resulting in a so‐called interferometric image. Internal multiples from above the redatuming level are eliminated during this process. In the past, we introduced point‐spread functions for interferometric redatuming by cross‐correlation. These point‐spread functions quantify distortions in the redatumed data, caused by internal multiple reflections in the overburden. In this paper, we define point‐spread functions for interferometric imaging to quantify these distortions in the image domain. These point‐spread functions are similar to conventional resolution functions for seismic migration but they contain additional information on the internal multiples in the overburden and they are partly data‐driven. We show how these point‐spread functions can be visualized to diagnose image defocusing and artefacts. Finally, we illustrate how point‐spread functions can also be defined for interferometric imaging with passive noise sources in the subsurface or with simultaneous‐source acquisition at the surface.