井下浅层地震探测技术的应用

钱 家 营矿 区内 8、12 煤 层顶 底 板为 细砂 岩 或粉 砂岩 ,煤 层与 顶 底板 之间 物 理特 性差 异 较大 ,有 利于 开 展 井下 -2浅 层地 震 勘探 。结合 矿 井开 采地 质 条件 及生 产 实际 ,建 立数 学物 理 模型 ,推 算时 距 曲线 、求 解地 层速 度 。在分 析选 用 仪器 的工 作 原理 及井 下 施工 方法 的 基础 上,对在 上 覆煤 层巷 道 中使 用该 技 术探 测 8、12-2煤层 开采 地 质条 件的 可 行 性及应 用前 景 进行 讨论 。     

Wavelets and the Generalization of the Variogram

The experimental variogram computed in the usual way by the method of moments and the Haar wavelet transform are similar in that they filter data and yield informative summaries that may be interpreted. The variogram filters out constant values; wavelets can filter variation at several spatial scales and thereby provide a richer repertoire for analysis and demand no assumptions other than that of finite variance. This paper compares the two functions, identifying that part of the Haar wavelet transform that gives it its advantages. It goes on to show that the generalized variogram of order k=1, 2, and 3 filters linear, quadratic, and cubic polynomials from the data, respectively, which correspond with more complex wavelets in Daubechies's family. The additional filter coefficients of the latter can reveal features of the data that are not evident in its usual form. Three examples in which data recorded at regular intervals on transects are analyzed illustrate the extended form of the variogram. The apparent periodicity of gilgais in Australia seems to be accentuated as filter coefficients are added, but otherwise the analysis provides no new insight. Analysis of hyerpsectral data with a strong linear trend showed that the wavelet-based variograms filtered it out. Adding filter coefficients in the analysis of the topsoil across the Jurassic scarplands of England changed the upper bound of the variogram; it then resembled the within-class variogram computed by the method of moments. To elucidate these results, we simulated several series of data to represent a random process with values fluctuating about a mean, data with long-range linear trend, data with local trend, and data with stepped transitions. The results suggest that the wavelet variogram can filter out the effects of long-range trend, but not local trend, and of transitions from one class to another, as across boundaries.     

On the efficient calculation of ordinary and generalized spherical harmonics

Algorithms for the stable computation of generalized and ordinary spherical harmonics are presented. The algorithms are fast and have the useful property that they can compute harmonics for isolated harmonic degrees. fortran and C programs implementing these algorithms are available from the authors.     

A Bootstrap Test Using Maximum Likelihood Ratio Statistics to Check the Similarity of Two 3-Dimensionally Oriented Data Samples

Comparing three-dimensionally oriented datasets is a problem encountered in various branches of earth science. A simple statistical tool for the comparison of two 3-dimensionally oriented datasets using the bootstrap method in line with the usual nonparametric permutation test is described here. This bootstrap test involves the estimation of maximum likelihood ratio statistic for properly constructed joint frequency tables of the datasets to be compared. This test does not use asymptotic result and will work well even for small sample sizes. Also this test does not make any specific distributional assumptions.     

The Second-Order Stationary Universal Kriging Model Revisited

Universal kriging originally was developed for problems of spatial interpolation if a drift seemed to be justified to model the experimental data. But its use has been questioned in relation to the bias of the estimated underlying variogram (variogram of the residuals), and furthermore universal kriging came to be considered an old-fashioned method after the theory of intrinsic random functions was developed. In this paper the model is reexamined together with methods for handling problems in the inference of parameters. The efficiency of the inference of covariance parameters is shown in terms of bias, variance, and mean square error of the sampling distribution obtained by Monte Carlo simulation for three different estimators (maximum likelihood, bias corrected maximum likelihood, and restricted maximum likelihood). It is shown that unbiased estimates for the covariance parameters may be obtained but if the number of samples is small there can be no guarantee of good estimates (estimates close to the true value) because the sampling variance usually is large. This problem is not specific to the universal kriging model but rather arises in any model where parameters are inferred from experimental data. The validity of the estimates may be evaluated statistically as a risk function as is shown in this paper.     

Regionalization of floods in New Brunswick (Canada)

This study uses the method of peaks over threshold (P.O.T.) to estimate the flood flow quantiles for a number of hydrometric stations in the province of New Brunswick, Canada. The peak values exceeding the base level (threshold), or `exceedances', are fitted by a generalized Pareto distribution. It is known that under the assumption of Poisson process arrival for flood exceedances, the P.O.T. model leads to a generalized extreme value distribution (GEV) for yearly maximum discharge values. The P.O.T. model can then be applied to calculate the quantiles X <sub>T</sub> corresponding to different return periods T, in years. A regionalization of floods in New Brunswick, which consists of dividing the province into `homogeneous regions', is performed using the method of the `region of influence'. The 100-year flood is subsequently estimated using a regionally estimated value of the shape parameter of the generalized Pareto distribution and a regression of the 100-year flood on the drainage area. The jackknife sampling method is then used to contrast the regional results with the values estimated at site. The variability of these results is presented in box-plot form. Received: June 1, 1997     

Regionalization of floods in New Brunswick (Canada)

This study uses the method of peaks over threshold (P.O.T.) to estimate the flood flow quantiles for a number of hydrometric stations in the province of New Brunswick, Canada. The peak values exceeding the base level (threshold), or `exceedances', are fitted by a generalized Pareto distribution. It is known that under the assumption of Poisson process arrival for flood exceedances, the P.O.T. model leads to a generalized extreme value distribution (GEV) for yearly maximum discharge values. The P.O.T. model can then be applied to calculate the quantiles X <sub>T</sub> corresponding to different return periods T, in years. A regionalization of floods in New Brunswick, which consists of dividing the province into `homogeneous regions', is performed using the method of the `region of influence'. The 100-year flood is subsequently estimated using a regionally estimated value of the shape parameter of the generalized Pareto distribution and a regression of the 100-year flood on the drainage area. The jackknife sampling method is then used to contrast the regional results with the values estimated at site. The variability of these results is presented in box-plot form. Received: June 1, 1997     

Generalized‐α methods for seismic structural testing

This paper presents novel predictor–corrector time‐integration algorithms based on the Generalized‐α method to perform pseudo‐dynamic tests with substructuring. The implicit Generalized‐α algorithm was implemented in a predictor–one corrector form giving rise to the implicit IPC–ρ∞ method, able to avoid expensive iterative corrections in view of high‐speed applications. Moreover, the scheme embodies a secant stiffness formula that can closely approximate the actual stiffness of a structure. Also an explicit algorithm endowed with user‐controlled dissipation properties, the EPC–ρb method, was implemented. The resulting schemes were tested experimentally both on a two‐ and on a six‐degrees‐of‐freedom system, using substructuring. The tests indicated that the numerical strategies enhance the fidelity of the pseudo‐dynamic test results even in an environment characterized by considerable experimental errors. Moreover, the schemes were tested numerically on severe non‐linear substructured multiple‐degrees‐of‐freedom systems reproduced with the Bouc–Wen model, showing the reliability of the seismic tests under these conditions. Copyright © 2004 John Wiley & Sons, Ltd.     

全球气温变化的多分形谱

运用计算多分形谱的多分形非趋势波动分析法，研究全球、北半球和中国月平均温度距平的多分形特征。全球气温、北半球气温和中国气温的变化是自相似的多分形；它们都表现出一定的正长程相关性，全球气温的长程相关性最强，北半球次之，中国气温相对最弱；中国气温的最强涨落较北半球气温大，而全球气温的最强涨落最小。