地表沉陷空间预报及其有理分式解法的应用

通过阐述地面沉陷机理,分析地下岩体内应力状态和地下开挖空间的形成过程.引入数学积分模型建立地面沉陷理论模型,重点说明单元下沉盆地和特征点下沉,推导出地面沉陷预报公式.并利用有理分式解法求解地面沉陷预报数值.此法的计算精度高,适用面广,有效地解决了常规解法工作量大、计算精度低的缺点.为一种有效、实用的计算地面沉陷预报值的方法.     

地球物理学中的电磁场正演与反演

本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.     

Identifying effects of local and nonlocal factors of soil water storage using cyclical correlation analysis

There are various factors governing the spatial and temporal variability of soil water storage including soil properties, topography and vegetation. Some factors act locally, whereas others act nonlocally, which means that a factor measured at one location has effect on soil water storage at another location. The objective of this study was to examine the effects of local and nonlocal controls of soil water storage in a hummocky landscape using cyclical correlation analysis. Soil water storage, soil properties and terrain indices were measured along a 128‐point transect of 576 m long from the semiarid, hummocky, prairie pothole region of North America. There are large coefficients of determination (r²) between soil water storage and sand content (r² = 0.32–0.53), organic carbon content (r² = 0.22–0.56), depth to carbonate layer (r² = 0.13–0.63), wetness index (r² = 0.25–0.45) and other variables at the measurement scale at different times, indicating strong local effects from these variables. The correlation coefficients were also calculated by physically shifting the spatial series of soil water storage with respect to that of controlling factors. The shifting improves the correlation between the spatial series, and the length of shifting indicated the difference in the response of soil water to its controlling factors. For example, the value of r² increased more than eightfold (r² = 0.47–0.64) after shifting the spatial series of soil water storage by 54 m, almost equal to the average length of existing slope, compared with the very weak correlation (r² = 0.02–0.08) at the measurement scale. This indicated the nonlocal effect from the relative elevation. The identification of nonlocal effects from factors improves the prediction of soil water storage. Copyright © 2011 John Wiley & Sons, Ltd.     

Comparison of ordinary and disjunctive kriging in an Indian iron ore deposit

This paper presents the results of disjunctive kriging applied to a supergene iron ore deposit of Bailadila Range of India. Disjunctive kriging is applied firstly to compare estimates of the blocks by ordinary kriging and secondly to estimate benchwise local recoverable reserves of the orebody. Good agreement exists between block estimates by ordinary kriging and disjunctive kriging except for peripheral blocks with less borehole information. Estimation of benchwise reserves shows that the behavior of the distribution of grades is different in various benches. The study shows that disjunctive kriging can be applied successfully for estimation of local recoverable reserves in the case of a good grade hematite iron ore deposit.     

Evaluation of rainfall-runoff performance using the stochastic integral equation method

The stochastic integral equation method (S.I.E.M.) is used to evaluate the relative performance of a set of both calibrated and uncalibrated rainfall-runoff models with respect to prediction errors. The S.I.E.M. is also used to estimate confidence (prediction) interval values of a runoff criterion variable, given a prescribed rainfall-runoff model, and a similarity measure used to condition the storms that are utilized for model calibration purposes.Because of the increasing attention given to the issue of uncertainty in rainfall-runoff modeling estimates, the S.I.E.M. provides a promising tool for the hydrologist to consider in both research and design.     

由计算机产生的等距离区域分界线

本文介绍一种最近由作者在美国研制的,由计算机产生等距离区域分界线的程序及其数据结构,同时给出初步的试验结果。     

磁尾宽频带静电噪声的一种新的激发机制

本文提出了磁尾宽频带静电噪声强谱段的一种新的激发机制。在等离子体片边界层中存在着瞬时的局域性的晨昏电场。电子和离子对该电场的响应不同:在电场存在的区域内,电子作回旋漂移运动,而离子轨道在的时间尺度内可视为直线。由于等离子体片边界层中等离子体是非均匀的,这就导致了电荷分离的产生,从而可激发静电不稳定性。本文通过求解含电场的Vlasov方程,计算了由此产生的不稳定波的频率和增长率,考察了波矢方向。上述计算结果均与强谱段静电噪声的观测特征一致。     

Seismic wave attenuation in fluid-saturated porous media

Contrary to the traditional view, seismic attenuation in Biot's theory of fluid-saturated porous media is due to viscous damping of local (not global) pore-fluid motion. Since substantial inhomogeneities in fluid permeability of porous geological materials are to be expected, the regions of highest local permeability contribute most to the wave energy dissipation while those of lowest permeability dominate the fluid flow rate if they are uniformly distributed. This dichotomy can explain some of the observed discrepancies between computed and measured attenuation of compressional and shear waves in porous earth. One unfortunate consequence of this result is the fact that measured seismic wave attenuation in fluid-filled geological materials cannot be used directly as a diagnostic of the global fluid-flow permeability.     

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.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.