Laplace-transform analytic element solution of transient flow in porous media

A Laplace-transform analytic element method (LT-AEM) is described for the solution of transient flow problems in porous media. Following Laplace transformation of the original flow problem, the analytic element method (AEM) is used to solve the resultant time-independent modified Helmholtz equation, and the solution is inverted numerically back into the time domain. The solution is entirely general, retaining the mathematical elegance and computational efficiency of the AEM while being amenable to parallel computation. It is especially well suited for problems in which a solution is required at a limited number of points in space–time, and for problems involving materials with sharply contrasting hydraulic properties. We illustrate the LT-AEM on transient flow through a uniform confined aquifer with a circular inclusion of contrasting hydraulic conductivity and specific storage. Our results compare well with published analytical solutions in the special case of radial flow.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 2: Complex modes

A method, based on the Hilbert–Huang spectral analysis, has been proposed by the authors to identify linear structures in which normal modes exist (i.e., real eigenvalues and eigenvectors). Frequently, all the eigenvalues and eigenvectors of linear structures are complex. In this paper, the method is extended further to identify general linear structures with complex modes using the free vibration response data polluted by noise. Measured response signals are first decomposed into modal responses using the method of Empirical Mode Decomposition with intermittency criteria. Each modal response contains the contribution of a complex conjugate pair of modes with a unique frequency and a damping ratio. Then, each modal response is decomposed in the frequency–time domain to yield instantaneous phase angle and amplitude using the Hilbert transform. Based on a single measurement of the impulse response time history at one appropriate location, the complex eigenvalues of the linear structure can be identified using a simple analysis procedure. When the response time histories are measured at all locations, the proposed methodology is capable of identifying the complex mode shapes as well as the mass, damping and stiffness matrices of the structure. The effectiveness and accuracy of the method presented are illustrated through numerical simulations. It is demonstrated that dynamic characteristics of linear structures with complex modes can be identified effectively using the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

利用MATLAB实现图像处理与分析

本文通过对图像的均衡操作、二值图像的处理、小波去噪等实例，介绍了MATLAB进行图像处理和分析的实用方法与技巧，它对各种软钉包、进行图像处理和分析等，提供了高效便捷的方法。     

Laser Ablation ICP-MS Analysis of Geological Materials Prepared as Lithium Borate Glasses

A simple, rapid and precise method is described for determining trace elements by laser ablation (LA)-ICP-MS analysis in bulk geological materials that have been prepared as lithium borate glasses following standard procedures for XRF analysis. This approach reliably achieves complete sample digestion and provides for complementary XRF and LA-ICP-MS analysis of a full suite of major and trace elements from a single sample preparation. Highly precise analysis is enabled by rastering an ArF excimer laser (λ= 193nm) across fused samples to deliver a constant sample yield to the mass spectrometer without inter-element fractionation effects during each analysis. Capabilities of the method are demonstrated by determination of twenty five trace elements (Sc, Ti, V, Ga, Rb, Sr, Y, Zr, Nb, Cs, Ba, REE, Hf, Ta, Pb, Th and U) in a diverse range of geological reference materials that includes peridotites, basalts, granites, metamorphic rocks and sediments. More than 90% of determinations are indistinguishable from published reference values at the 95% confidence level. Systematic bias greater than 5% is observed for only a handful of elements (Zr, Nb and U) and may be attributed in part to inaccurate calibration values used for the NIST SRM 612 glass in the case of Zr and Nb. Detection limits for several elements, most notably La, are compromised at ultra-trace levels by impurities in the lithium borate flux but can be corrected for by subtracting appropriate procedural blanks. Reliable Pb analysis has proved problematic due to variable degrees of contamination introduced during sample polishing prior to analysis and from Pt-crucibles previously used to fuse Pb-rich samples. Scope exists for extending the method to include internal standard element/isotope spiking, particularly where integrated XRF analysis is not available to characterise major and trace elements in the fused lithium borate glasses prior to LA-ICP-MS analysis.     

A Reappraisal of Rb, Y, Zr, Pb and Th Values in Geochemical Reference Material BHVO-1

The geochemical reference material BHVO-1 was analysed by a variety of techniques over a six year period. These techniques included inductively coupled plasma-mass spectrometry and atomic emission spectroscopy (ICP-MS and ICP-AES, respectively), laser ablation ICP-MS and spark source mass spectroscopy. Inconsistencies between the published consensus values reported by Gladney and Roelandts (1988, Geostandards Newsletter) and the results of our study are noted for Rb, Y, Zr, Pb and Th. The values reported here for Rb, Y, Zr and Pb are generally lower, while Th is higher than the consensus value. This is not an analytical artefact unique to the University of Notre Dame ICP-MS facility, as most of the BHVO-1 analyses reported over the last ten to twenty years are in agreement with our results. We propose new consensus values for each of these elements as follows: Rb = 9.3 ± 0.2 μg g-1 (compared to 11 ± 2 μg g-1), Y = 24.4 ± 1.3 μg g-1 (compared to 27.6 ± 1.7 μg g-1), Zr = 172 ± 10 μg g-1 (compared to 179 ± 21 μg g-1), Pb = 2.2 ± 0.2 μg g-1 (compared to 2.6 ± 0.9 μg g-1) and Th = 1.22 ± 0.02 μg g-1 (compared to 1.08 ± 0.15 μg g-1).     

Hilbert-Huang变换在密频结构阻尼识别中的应用

Hilbert—Huang变换是一种新的数据处理方法，由经验模分解(Empirical Mode Decomposition)技术及Hilbert变换两部分组成。本文研究此方法对于密频结构阻尼识别的应用。首先对于两自由度系统模型，说明该方法用于阻尼识别的步骤。进而研究存在频率密集现象的高层建筑的阻尼识别问题。上述结果与理论值及由半功率带宽法的识别值进行了比较，对比显示Hilbert．Huang方法较传统方法具有良好的识别密频结构阻尼的性能，适用于大型结构的系统识别。     

小波变换在识别储层流体性质中的应用

根据小波变换的特点提出一种识别储层流体的新方法. 本文假设测井所得的信号为地层微观孔隙结构、地层流体性质及岩性和岩相等各部分的综合贡献之和，不同的储层中小波分析的能量谱特征不同，因而通过对已知储层模型进行小波能量谱分析，可得到各种储层的特征能量谱，根据这些储层特征能量谱就可以对待划分的储层段进行储层流体性质判别.     

河北省黄壁庄水库坝基渗漏单井示踪试验研究

以河北省黄壁庄水库副坝坝基的渗漏塌坑为例，利用“九五”国家重点科技成果推广项目“智能化地下水动态参数测量仪”在天然流场下，测量出水库自然渗漏量地下水的渗透流速流向、垂直向流速流向，每米含水层的渗透系数等工程必须的水文地质参数。确定性地诊断和检测出黄壁水库坝基渗漏隐患出现的具体位置和高程，从而为大坝的渗漏处理和安全处理提供详细的现场测量数据。     

地震资料小波剖面制作系统

对不同的小波变换算法和小波基函数进行了研究，得到了适用于地震资料处理的最佳小波基函数，并用Visual C 语言开发了基于Windows操作平台下的地震资料小波剖面制作系统。利用小波剖面制作系统和最佳小波基函数对煤田实际地震资料进行了处理，取得了令人满意的地质效果。