Three dimensional simulation of fluid flow in X‐ray CT images of porous media

A numerical scheme is developed in order to simulate fluid flow in three dimensional (3‐D) microstructures. The governing equations for steady incompressible flow are solved using the semi‐implicit method for pressure‐linked equations (SIMPLE) finite difference scheme within a non‐staggered grid system that represents the 3‐D microstructure. This system allows solving the governing equations using only one computational cell. The numerical scheme is verified through simulating fluid flow in idealized 3‐D microstructures with known closed form solutions for permeability. The numerical factors affecting the solution in terms of convergence and accuracy are also discussed. These factors include the resolution of the analysed microstructure and the truncation criterion. Fluid flow in 2‐D X‐ray computed tomography (CT) images of real porous media microstructure is also simulated using this numerical model. These real microstructures include field cores of asphalt mixes, laboratory linear kneading compactor (LKC) specimens, and laboratory Superpave gyratory compactor (SGC) specimens. The numerical results for the permeability of the real microstructures are compared with the results from closed form solutions. Copyright © 2004 John Wiley & Sons, Ltd.     

Implications of low-pass filtering on power spectra and autocorrelation functions of turbulent velocity signals

Filtering either through the electronics of an instrument or through digital procedure is performed routinely on geophysical data. When velocity fluctuations are measured in turbulent flows using electromagnetic current meters (ECMs), a builtin lowpass Butterworth filter of order n usually attenuates fluctuations at high frequencies. However, the effects of this filter may not be acknowledged in turbulence studies, thus impeding comparisons between data collected with different ECMs. This paper explores the implications of the filters on the characteristics of velocity signals, mainly on variance, power spectra, and correlation analyses. Variance losses resulting from filtering can be important but will vary with the order n of the Butterworth filter, decreasing as n increases. Knowing the filter response, it is possible to reconstruct the original signal spectrum to evaluate the effect of filtering on variance and to allow comparisons between data collected with different instruments. The autocorrelation function also is affected by filtering which increases the value of the coefficients in the first lags, resulting in an overestimation of the integral length scale of coherent structures. These important effects add to those related to size and shape differences in ECM sensors and must be taken into account in comparative studies.     

Electromagnetic signals associated with stick–slip of quartz-free rocks

We recorded clear transients in the electric and magnetic fields upon sudden slip in stick–slip experiments on dry, cylindrically shaped, quartz-free rock specimens of basalt and peridotite with a 30° saw-cut (representing a fault) at confining pressures of up to 120 MPa. The amplitudes of the measured electric field signals were always higher at the electrode pair oriented parallel to the strike of the fault than at the pair oriented perpendicular. This anisotropy suggests a preferred electric polarization normal to the slip surface. The transients in the electric and magnetic fields were observed only when the fault slip occurred by stick–slip mode, not by a stable mode of the sliding, and the amplitudes of the electric field signals increased with increasing stress drop. It is suggested that the generation process of the electromagnetic signals is closely related to the characteristic behavior of the fault at the time of the initiation of slip during stick–slip events, probably with respect to the intensity of the signals. We propose that one or both of the following two processes characteristic of the fault at the time of the initiation of slip during stick–slip events are essential for the generation of detectable electromagnetic signals: rapid slip along the simulated fault and separation of the rock masses across the fault.     

SOMATOM AR.C型CT机复位故障检修

本文介绍了由电源故障引起的不可逆转故障,并讨论了维修和故障的排除。     

Radiomagnetotelluric mapping, groundwater numerical modelling and 18-Oxygen isotopic data as combined tools to determine the hydrogeological system of a landslide prone area

Three methods were combined to determine the groundwater recharge and transfer processes of a landslide prone area. First, the radiomagnetotelluric method was used to investigate the distribution of electrical resistivity (ρ) of the subsurface and build a three-dimensional model of permeability (k), through an experimental relation between ρ and k. Second, this structural model of permeability and additional climatologic data were used to fix boundary and recharge conditions to perform a three-dimensional and transient numerical simulation of the groundwater flow. Finally 18-Oxygen time series observed at the main springs were used to validate the model. This association of methods led to an improved characterization of the groundwater flow system at local scale and a better understanding of the role of this system on the landslide phenomenon. This structured approach is thought to be useful to design specific remediation strategies to drain the unstable mass.     

Fast computation of Hankel Transform using orthonormal exponential approximation of complex kernel function

The computation of electromagnetic (EM) fields, for 1-D layered earth model, requires evaluation of Hankel Transform (HT) of the EM kernel function. The digital filtering is the most widely used technique to evaluate HT integrals. However, it has some obvious shortcomings. We present an alternative scheme, based on an orthonormal exponential approximation of the kernel function, for evaluating HT integrals. This approximation of the kernel function was chosen because the analytical solution of HT of an exponential function is readily available in literature. This expansion reduces the integral to a simple algebraic sum. The implementation of such a scheme requires that the weights and the exponents of the exponential function be estimated. The exponents were estimated through a guided search algorithm while the weights were obtained using Marquardt matrix inversion method. The algorithm was tested on analytical HT pairs available in literature. The results are compared with those obtained using the digital filtering technique with Anderson filters. The field curves for four types (A-, K-, H-and Q-type) of 3-layer earth models are generated using the present scheme and compared with the corresponding curves obtained using the Anderson sc heme. It is concluded that the present scheme is more accurate than the Anderson scheme     

X线及CT检查在外伤性骨质疏松中的应用及讨论

目的 讨论CT及X线检查在外伤性骨质疏松中的应用。方法 1998－2000年200例外伤复查病人行CT及X光检查。结果 2－3个月后出现骨疏松的病人113例，3－4个月后出现骨疏松的66例，未见骨疏松的21例。结论 骨折2－3个月后复查病人，应适量降低CT的窗宽，窗位和X线的照射剂量。     

X线侧位片与高分辨率CT对鼻骨骨折的对比研究

目的　通过86例鼻骨骨折的鼻骨X线侧位片与高分辨率CT(HRCT）对比研究，确定HRCT检查对鼻骨骨折的诊断价值。材料与方祛　86例鼻外伤或/和He、面外伤的患者均进行了鼻骨X线侧位摄影、HRCT扫描检查，经集体讨论，由有经验的放射医师对全部影像资料进行分析并进行统计学处理。结果HRCT显示所有的86例鼻骨骨折，但在此86例骨折中X线平片只显示了35例阳性，包括15例鼻骨汇合部骨折，7例双侧鼻骨骨折，3例复合骨折(鼻骨汇合部及双侧鼻骨骨折)，l0例鼻区骨折。X线平片显示5l例为阴性。结论　对鼻骨骨折的诊断，HRCT明显优于鼻骨X线侧位片，HRCT扫描检查可提高X线隐性鼻骨骨折的诊断，并可明确鼻骨骨折的类型、范围和程度，为临床诊断和治疗提供可靠的依据。     

螺旋CT三维重建在骨盆骨折中的应用

目的：探讨用螺旋CT三维重建评价骨盆骨折的临床应用价值。方法：收集骨盆骨折17例，采用螺旋CT三维重建检查。结果：17例检查均取得清晰立体图像，且可消除某一骨性结构，单独观察骨盆某一组成骨的骨折情况。结论：螺旋CT三维重建在骨盆骨折的诊断、手术方案的设计及教学等有较高价值。     

2002年CT扫描和三维成像学术年会评述

由CT理论与应用研究编辑部和北京信息工程学院主办、有中国体视学学会多位理事和会员参加的“2002年CT扫描和三维成像研讨会（CT科技学术年会）”于2002年10月18－20日，在北京大学逸夫苑理科楼学术报告厅举行。中国体视学学会副理事长刘国权教授出席致辞，对于这次CT科技年会的召开表示祝贺，并介绍了2003年在美国圣迭戈和北京将要召开第11届国际体视学会大会的有关情况，希望有关专家准备论文出席，来自全国各地48位专家，教授出席了会议，本次会议涉及CT理论、计算方法和软件及硬件技术；医学CT扫描成像和疾病诊断；工业CT和三维成像：地球层析和工程CT。美国爱荷华大学王革教授，纽约罗彻斯特大学宁若拉教授专程到会，分别做了精彩的特邀主题报告：显微CT的锥束重建新方法和近年来在这方面的进展和展望”，“锥光束三维CT技术在美国的过去、现在和未来”“国内专家也报告了前沿研究的新成果，在会前，CT理论与应用研究编辑部，根据收到的论文和摘要汇编和印制了专集，在会议最后一天，召开了“CT理论与应用研究”期刊编委会，传达了中国科协对于办好科技期刊的新要求，会议开得很成功。