坐骨股骨撞击综合征的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.     

三维CT成像的数学模型

三维CT成像广泛应用于放疗模拟定位等临床实践。当扫描运动物体时三维图像存在的位移和形变,与扫描参数设置和被扫描物体的运动状态等关系密切。本文通过理论分析CT扫描过程,建立三维成像的数学模型,计算三维图像的时间分辨率、位移和形变等指标,并分析其变化规律,为将来进一步的实验研究奠定基础。     

医用PET系统技术概述

PET(正电子发射成像)是当今高层次的核医学技术。它综合了核技术、核电子学、计算机、数据处理等领域中的尖端技术,是当前医学界最先进的大型医疗诊断成像设备之一,在肿瘤、心血管疾病和神经系统疾病的临床诊断中发挥着越来越重要的作用。本文以国产32环高档临床PET为例,较详细地介绍了PET设备的工作原理、总体基本结构及各子系统基本功能,对PET整机性能及其临床应用也做了简述。     

应用核磁共振成像(NMRI)技术研究复合剂驱油特征

本文应用核磁共振成像技术,在均值、非均值模型中研究了复合剂驱油动态实验过程。观察及分析了复合剂驱油特征和效果,对复合剂驱油段塞前缘进行了实时跟踪,得到了复合剂驱油段塞前缘油墙的形成及剩余油分布与运移的二维图像,给出了复合剂驱油后,剩余油的饱和度分布曲线。实验表明:在层状非均值地层中,若渗透率差异大,复合剂段塞主要进入高渗层,而波及不到低渗层,即复合剂波及效率低。若渗透率差异小,复合剂段塞大部分进入高渗层,有少部分进入低渗层,即复合剂波及效率较高。     

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.