Imaging mantle transition zone thickness with SdS-SS finite-frequency sensitivity kernels

We invert differential SdS-SS traveltime residuals measured from stacked waveforms and finite-frequency sensitivity kernels for topography on the 410- and 660-km discontinuities. This approach yields higher resolution images of transition zone thickness than previous stacking methods, which simply average/smooth over topographic features. Apparent structure measured using simple stacking is highly dependent upon the bin size of each stack. By inverting for discontinuity topography with a variety of bin sizes, we can more accurately calculate the true structure. The inverted transition zone model is similar to simple stack models with an average thickness of 242 km, but the lateral variations in thickness are larger in amplitude and smaller in scale. Fast seismic velocities in 3-D mantle models such as SB4L18 correlate with areas of thicker transition zone. The elongated curvilinear regions of thickened transition zone that occur near subduction zones are narrow and high amplitude, which suggests relatively little lateral spreading and warming of subducted lithosphere within the transition zone. The anomalously thin transition zone regions are laterally narrow, and not broadly continuous. If these variations in transition zone thickness are interpreted as thermal in nature, then this model suggests significant temperature variations on small lateral scales.     

Parametrizing surface wave tomographic models with harmonic spherical splines

We present a mathematical framework and a new methodology for the parametrization of surface wave phase-speed models, based on traveltime data. Our method is neither purely local, like block-based approaches, nor is it purely global, like those based on spherical harmonic basis functions. Rather, it combines the well-known theory and practical utility of the spherical harmonics with the spatial localization properties of spline basis functions. We derive the theoretical foundations for the application of harmonic spherical splines to surface wave tomography and summarize the results of numerous numerical tests illustrating the performance of a practical inversion scheme based upon them. Our presentation is based on the notion of reproducing-kernel Hilbert spaces, which lends itself to the parametrization of fully 3-D tomographic earth models that include body waves as well.     

Quantitative analysis of soil pores under natural vegetation successions on the Loess Plateau

Soil pore is a key attribute of the soil structure that affects soil reservoir under natural vegetation recovery on the Loess Plateau.This study is to quantitatively analyze soil pore parameters,measured with Computed Tomography(CT) at 15-57 mm depths under five different vegetation succession stages using a concept of substituting space for time in the Ziwuling Forest Region of the Loess Plateau.The results showed that the soil pore parameters,such as pore number,porosity,circularity,and fractal dimension,...     

应用地震层析成象技术研究深部速度结构

罗马尼亚Ｖａｎｃｅａ地震区是大陆上发生与板块磁撞和削减有关的中深部地震活动的地区之一。本文介绍了应用地震层析成象方法研究该地区深部速度结构成果。在研究中使用了地方和区域地震所记录的４３３个浅源和中部地震的到时资料反演求解深至２００ｋｍ的三维速度结构,在走时和射线路径的计算中利用了有效的三维射线跟踪技术,在反演中采用ＬＳＱＲ算法,高分辨率的地震层析图象揭示了速度结构的广泛不均匀性,结果表明,地震层析     

肺癌CT诊断与分析

本文目的是探讨ＣＴ在肺癌诊断中的价值,所用的材料与方法是搜集赤峰市第六医院的肺癌病例４６例,所有病例均手术病理,支气管镜检或细胞这检查证实,常规ＣＴ扫描后对选择病例在选择层面行ＨＲＣＴ扫描,结果得到了ＣＴ在肺癌诊断上,特别在肺癌分期上优于常规Ｘ线检查,此外,ＨＲＣＴ在显示细微病变或病变伪细微结构方面有重要价值,结论是ＣＴ不仅可用于肺癌诊断还可用于肺癌分期,因而更有益于临床医生制订治疗方案。     

Image‐ray Tomography

Seismic tomography is a well‐established approach to invert smooth macro‐velocity models from kinematic parameters, such as traveltimes and their derivatives, which can be directly estimated from data. Tomographic methods differ more with respect to data domains than in the specifications of inverse‐problem solving schemes. Typical examples are stereotomography, which is applied to prestack data and Normal‐Incidence‐Point‐wave tomography, which is applied to common midpoint stacked data. One of the main challenges within the tomographic approach is the reliable estimation of the kinematic attributes from the data that are used in the inversion process. Estimations in the prestack domain (weak and noisy signals), as well as in the post‐stack domain (occurrence of triplications and diffractions leading to numerous conflicting dip situations) may lead to parameter inaccuracies that will adversely impact the resulting velocity models. To overcome the above limitations, a new tomographic procedure applied in the time‐migrated domain is proposed. We call this method Image‐Incident‐Point‐wave tomography. The new scheme can be seen as an alternative to Normal‐Incidence‐Point‐wave tomography. The latter method is based on traveltime attributes associated with normal rays, whereas the Image‐Incidence‐Point‐wave technique is based on the corresponding quantities for the image rays. Compared to Normal‐Incidence‐Point‐wave tomography the proposed method eases the selection of the tomography attributes, which is shown by synthetic and field data examples. Moreover, the method provides a direct way to convert time‐migration velocities into depth‐migration velocities without the need of any Dix‐style inversion.     

Misfit functionals in Laplace‐Fourier domain waveform inversion,with application to wide‐angle ocean bottom seismograph data

In seismic waveform inversion, non‐linearity and non‐uniqueness require appropriate strategies. We formulate four types of L2 normed misfit functionals for Laplace‐Fourier domain waveform inversion: i) subtraction of complex‐valued observed data from complex‐valued predicted data (the ‘conventional phase‐amplitude’ residual), ii) a ‘conventional phase‐only’ residual in which amplitude variations are normalized, iii) a ‘logarithmic phase‐amplitude’ residual and finally iv) a ‘logarithmic phase‐only’ residual in which the only imaginary part of the logarithmic residual is used. We evaluate these misfit functionals by using a wide‐angle field Ocean Bottom Seismograph (OBS) data set with a maximum offset of 55 km. The conventional phase‐amplitude approach is restricted in illumination and delineates only shallow velocity structures. In contrast, the other three misfit functionals retrieve detailed velocity structures with clear lithological boundaries down to the deeper part of the model. We also test the performance of additional phase‐amplitude inversions starting from the logarithmic phase‐only inversion result. The resulting velocity updates are prominent only in the high‐wavenumber components, sharpening the lithological boundaries. We argue that the discrepancies in the behaviours of the misfit functionals are primarily caused by the sensitivities of the model gradient to strong amplitude variations in the data. As the observed data amplitudes are dominated by the near‐offset traces, the conventional phase‐amplitude inversion primarily updates the shallow structures as a result. In contrast, the other three misfit functionals eliminate the strong dependence on amplitude variation naturally and enhance the depth of illumination. We further suggest that the phase‐only inversions are sufficient to obtain robust and reliable velocity structures and the amplitude information is of secondary importance in constraining subsurface velocity models.     

Forced imbibition into a limestone: measuring P‐wave velocity and water saturation dependence on injection rate

Quantitative interpretation of time‐lapse seismic data requires knowledge of the relationship between elastic wave velocities and fluid saturation. This relationship is not unique but depends on the spatial distribution of the fluid in the pore‐space of the rock. In turn, the fluid distribution depends on the injection rate. To study this dependency, forced imbibition experiments with variable injection rates have been performed on an air‐dry limestone sample. Water was injected into a cylindrical sample and was monitored by X‐Ray Computed Tomography and ultrasonic time‐of‐flight measurements across the sample. The measurements show that the P‐wave velocity decreases well before the saturation front approaches the ultrasonic raypath. This decrease is followed by an increase as the saturation front crosses the raypath. The observed patterns of the acoustic response and water saturation as functions of the injection rate are consistent with previous observations on sandstone. The results confirm that the injection rate has significant influence on fluid distribution and the corresponding acoustic response. The complexity of the acoustic response —‐ that is not monotonic with changes in saturation, and which at the same saturation varies between hydrostatic conditions and states of dynamic fluid flow – may have implications for the interpretation of time‐lapse seismic responses.     

沿河开发建筑场地之隐伏沙洞的电阻率层析成像

采用电阻率层析成像二维测量的方法,对河南省三门峡市沿河开发的建筑场地进行了勘探,获取了地下隐伏沙洞的准确位置、埋藏深度、含水情况、尺寸大小等。说明电阻率层析成像二维测量的方法在地下隐伏沙洞的探测中是行之有效的;地下隐伏沙洞在没有水体存在时,电阻率层析成像二维测量成果图中反映的一般都是高阻异常封闭圈。如果有水体存在,其异常则表现为低阻异常封闭圈。