根据渐进谱的统计规律生成地震加速度时程

基于美国西部80条基岩上的近场强震记录, 采用Nakayama方法生成记录的渐进谱, 并参照Kameda方式,用统计方法建立了根据震级、距离等地震参数预测渐进谱的统计模型. 提出一种以渐进谱为目标谱的生成幅值和频率非平稳地震加速度时程的迭代方法. 由于考虑了渐进谱幅值和相位的相互影响,所生成的时程的相位也是时频非平稳的,并在相位调整中识别了相位谱增量符号以加速迭代收敛进程. 最后根据统计回归的目标渐进谱模型和本文提出的拟合目标渐进谱的方法,可生成不同震级、距离条件下的幅值和频率均非平稳的地震加速度时程.      

利用克里金技术标定新疆及附近地区P/S震相幅值比及其在地震事件识别中的应用

针对中国的乌鲁木齐台（WMQ）和哈萨克斯坦的马坎奇台（MAK）,采用贝叶斯克里金技术建立了新疆及附近地区地震事件的震相幅值比的校正曲面,并分析了校正曲面与地震波传播路径的关系,以及校正曲面对于改善震相幅值比判据识别效果的有效性.分析结果表明,建立的校正曲面基本合理,结合传播路径差异和地震波传播规律,能够较好地解释校正曲面的起伏特征. 而且,在传播路径差异较大,使得震相幅值比经震中距校正后离散度仍然较高的情况下,利用校正曲面修正震相幅值比,可以进一步降低其离散度,从而改善判据的识别效果.经校正曲面修正后,对天然地震的误识率由16.3%下降为5.2％.      

Spatial and multivariate analysis of geochemical data from metavolcanic rocks in the Ben Nevis area, Ontario

A study of the lithogeochemistry of metavolcanics in the Ben Nevis area of Ontario, Canada has shown that factor analysis methods can distinguish lithogeochemical trends related to different geological processes, most notably, the principal compositional variation related to the volcanic stratigraphy and zones of carbonate alteration associated with the presence of sulphides and gold. Auto- and cross-correlation functions have been estimated for the two-dimensional distribution of various elements in the area. These functions allow computation of spatial factors in which patterns of multivariate relationships are dependent upon the spatial auto- and cross-correlation of the components. Because of the anisotropy of primary compositions of the volcanics, some spatial factor patterns are difficult to interpret. Isotropically distributed variables such as CO ₂ are delineated clearly in spatial factor maps. For anisotropically distributed variables (SiO ₂ ), as the neighborhood becomes smaller, the spacial factor maps becomes better. Interpretation of spatial factors requires computation of the corresponding amplitude vectors from the eigenvalue solution. This vector reflects relative amplitudes by which the variables follow the spatial factors. Instability of some eigenvalue solutions requires that caution be used in interpreting the resulting factor patterns. A measure of the predictive power of the spatial factors can be determined from autocorrelation coefficients and squared multiple correlation coefficients that indicate which variables are significant in any given factor. The spatial factor approach utilizes spatial relationships of variables in conjunction with systematic variation of variables representing geological processes. This approach can yield potential exploration targets based on the spatial continuity of alteration haloes that reflect mineralization.List of symbols c _i Scalar factor that minimizes the discrepancy between andU _i - D Radius of circular neighborhood used for estimating auto- and cross-correlation coefficients - d Distance for which transition matrixU is estimated - d _ij Distance between observed valuesi andj - E Expected value - E _i Row vector of residuals in the standardized model - F(d _ij) Quadratic function of distanced _ij F(d _ij)=a+bd _ij+cd _ij ² - L Diagonal matrix of the eigenvalues ofU - _i Eigenvalue of the matrixU;ith diagonal element ofL - N Number of observations - p Number of variables - Q Total predictive power ofU - R Correlation matrix of the variables - R _0j Variance-covariance signal matrix of the standardized variables at origin;j is the index related tod andD (e.g.,j=1 ford=500 m,D=1000 m) - R _1j Matrix of auto- and cross-correlation coefficients evaluated at a given distance within the neighborhood - R _m ² Multiple correlation coefficient squared for themth variable - S _i Column vectori of the signal values - s _k ² Residual variance for variablek - T _i Amplitude vector corresponding toV _i;ith row ofT=V ^–1 - T Total variation in the system - U Nonsymmetric transition matrix formed by post-multiplyingR ₀₁ ^–1 byR _ij - U _i Componenti of the matrixU, corresponding to theith eigenvectorV _i;U _i= _iV_iT_i - U* _i ComponentU _i multiplied byc _i - U _ij Sum of componentsU _i+U _j - V _i Eigenvector of the matrixU;ith column ofV withUV=VL - w Weighting factor; equal to the ratio of two eigenvalues - X _i Random variable at pointi - x _i Value of random variable at pointi - y _i Residual ofx _i - Z _i Row vectori for the standardized variables - z _i Standardized value of variable     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Development of the ultrasonic measurement technique and its application in the observation of rock mass stability

Summary A method of measuring ultrasonic wave properties in rocks during the complete stress-strain process is described. The relevant experimental laboratory study, to reveal the change of the amplitude spectrum parameters with strain (or stress) has been carried out. A preliminary study was conducted on the application of the ultrasonic measurement technique at a belt conveyor roadway of the north wing in Baodian coal mine, Shandong province. A borehole ultrasonic device with dry coupling was developed to provide better coupling and more accurate measurement data in comparison with those of a water coupled situation. Based on the interrelations between the complete stress-strain properties of specimens and their wave properties, the structural properties of surrounding rocks, the range of yield zones, and the change of stresses within the rocks surrounding a roadway which was subjected to mining influence of upper longwall face were analysed. Amplitude spectrum parameters were used to give a better prediction of the physico-mechanical state of the surrounding rocks.     

The Fourier Integral Method: An efficient spectral method for simulation of random fields

The Fourier Integral Method (FIM) of spectral simulation, adapted to generate realizations of a random function in one, two, or three dimensions, is shown to be an efficient technique of non-conditional geostatistical simulation. The main contribution is the use of the fast Fourier transform for both numerical calculus of the density spectral function and as generator of random finite multidimensional sequences with imposed covariance. Results obtained with the FIM are compared with those obtained by other classic methods: Shinozuka and Jan Method in 1D and Turning Bands Method in 2D and 3D, the points for and against different methodologies are discussed. Moreover, with the FIM the simulation of nested structures, one of which can be a nugget effect and the simulation of both zonal and geometric anisotropy is straightforward. All steps taken to implement the FIM methodology are discussed.     

国产WSP—2半米平面光栅测定超痕量金方法介绍

本文叙述了国产半米平面光栅光程短、光强好、谱线锐度佳的特点,通过最佳条件选择,改进分析方法,使其测金灵敏度为3×10~(-9)g,取样10g可测定含金量范围为0.3—1000ppb,取得了很好的效果。     

油气资源遥感勘查新方法

地下深处的油气藏中的烃类物质,在压力梯度下会向地表产生微渗漏,导致地表物质产生一系列理化异常,在地表形成吸附烃异常、红层退色、粘土化、碳酸盐化、放射性异常和地植物异常等标志。根据这些地表标志的地物波谱特征,选择出油气资源遥感勘查的最佳工作波段,在此基础上,以机载遥感信息为主,以星载遥感信息为辅,以地表地球化学、地植物学、地面遥感勘查作为验证手段,最终施行多元信息复合分析,可达到寻找油气藏的目的。