利用地层空间的相对稳定性、井资料、地震资料模拟高分辨率波阻抗资料

由于地震资料是带限的,而井资料约束反演过程中的地震反演频带的展宽也是有一定限度的,使得反演地震波阻抗也是带限的。如何进一步拓宽反演地震波阻抗的频带宽度,提高地震资料的分辨率是文中研究的内容。利用井约束下的地震反演波阻抗资料、地震资料及井资料,借鉴地质模拟技术,研究如何利用地层的空间相对稳定性特征,通过相关函数的分析,估计反演更高分辨率的地震波阻抗资料。考虑到地质体空间的相关性及分辨率,相关半径取相关函数从最大点到第1个极值点的范围。在目标函数的计算过程中,既考虑当前计算点的权系数信息,又考虑井外推时的权系数信息,这样使计算结果更加稳定可靠。通过综合模拟计算,能够识别较小砂体     

阿尔金北缘大平沟金矿床成因

大平沟金矿是阿尔金北缘地区近年来新发现的一个中型金矿床。它受韧性剪切带构造的控制,矿化类型为蚀变糜棱岩型夹少量钾长石石英脉型。文章通过地质地球化学的研究,探讨了该矿床的成因。稀土元素、硫同位素、铅同位素和锶同位素地球化学特征反映成矿物质主要来源于变质岩;流体包裹体和氢氧同位素反映成矿流体主要来源于变质水,有少量重熔岩浆水和大气降水的混入,成矿温度198～2 90℃,成矿流体密度0 .791～0 .971g/cm3 ;盐度w(NaCleq) =2 .0 3%～5 .5 % ,成矿压力(42 0～72 0 )×10 5Pa ,成矿深度为1.6 1～2 .6 8km。石英流体包裹体Rb_Sr等时线年龄表明成矿时代为加里东期。大平沟金矿床属于中温动力变质热液成因矿床。     

Reliability analysis of drilled shaft behavior using finite difference method and <Emphasis Type="Italic">Monte Carlo</Emphasis> simulation

Load displacement analysis of drilled shafts can be accomplished by utilizing the “t-z” method, which models soil resistance along the length and tip of the drilled shaft as a series of springs. For non-linear soil springs, the governing differential equation that describes the soil-structure interaction may be discretized into a set of algebraic equations based upon finite difference methods. This system of algebraic equations may be solved to determine the load–displacement behavior of the drilled shaft when subjected to compression or pullout. By combining the finite difference method with Monte Carlo simulation techniques, a probabilistic load–displacement analysis can be conducted. The probabilistic analysis is advantageous compared to standard factor of safety design because uncertainties with the shaft–soil interface and tip properties can be independently quantified. This paper presents a reliability analysis of drilled shaft behavior by combining the finite difference technique for analyzing non-linear load–displacement behavior with Monte Carlo simulation method. As a result we develop probabilistic relationships for drilled shaft design for both total stress (undrained) and effective stress (drained) parameters. The results are presented in the form of factor of safety or resistance factors suitable for serviceability design of drilled shafts.     

Digital Image Based Approach for Three-Dimensional Mechanical Analysis of Heterogeneous Rocks

Summary This paper presents a digital image based approach for three-dimensional (3-D) numerical simulation and failure analysis of rocks by taking into account the actual 3-D heterogeneity. Digital image techniques are adopted to extract two-dimensional (2-D) material heterogeneity from material surface images. The 2-D image mesostructures are further extrapolated to 3-D cuboid mesostructures by assuming the material surface as a representation of the inner material heterogeneity within a very small depth. The iterative milling and scanning system is set up to generate the 3-D rock mesostructures. A Hong Kong granite specimen is used as an example to demonstrate the procedure of 3-D mesostructure establishment. The mechanical responses and failure process under the conventional Brazilian tensile test condition are examined through numerical analyses. The stress distribution, crack propagation process and failure model of heterogeneous material cases are simulated with a finite difference software. The numerical results indicate that material heterogeneity plays an important role in determining the failure behavior of rocks under external loading.     

Far-field simulation of the 1946 Aleutian tsunami

We present hydrodynamic far-field simulations of the Aleutian tsunami of 1946 April 1, using both a dislocation source representing a slow earthquake and a dipolar one modelling a large landslide. The earthquake source is derived from the recent seismological study by López and Okal, while the landslide source was previously used to explain the exceptional run-up at Scotch Cap in the near field. The simulations are compared to a field data set previously compiled from testimonies of elderly witnesses at 27 far-field locations principally in the Austral and Marquesas Islands, with additional sites at Pitcairn, Easter and Juan Fernández. We find that the data set is modelled satisfactorily by the dislocation source, while the landslide fails to match the measured amplitudes, and to give a proper rendition of the physical interaction of the wavefield with the shore, in particular at Nuku Hiva, Marquesas. The emerging picture is that the event involved both a very slow earthquake, responsible for the far-field tsunami, and a major landslide explaining the near-field run-up, but with a negligible contribution in the far field.     

Gaussian beams in inhomogeneous media: A review

The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method, which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are twins, though they differ by names. The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered medium.     

Numerical Simulation of Andhra Severe Cyclone (2003): Model Sensitivity to the Boundary Layer and Convection Parameterization

The Andhra severe cyclonic storm (2003) is simulated to study its evolution, structure, intensity and movement using the Penn State/NCAR non-hydrostatic mesoscale atmospheric model MM5. The model is used with three interactive nested domains at 81, 27 and 9 km resolutions covering the Bay of Bengal and adjoining Indian Peninsula. The performance of the Planetary Boundary Layer (PBL) and convective parameterization on the simulated features of the cyclone is studied by conducting sensitivity experiments. Results indicate that while the boundary layer processes play a significant role in determining both the intensity and movement, the convective processes especially control the movement of the model storm. The Mellor-Yamada scheme is found to yield the most intensive cyclone. While the combination of Mellor-Yamada (MY) PBL and Kain-Fritsch 2 (KF2) convection schemes gives the most intensive storm, the MRF PBL with KF2 convection scheme produces the best simulation in terms of intensity and track. Results of the simulation with the combination of MRF scheme for PBL and KF2 for convection show the evolution and major features of a mature tropical storm. The model has very nearly simulated the intensity of the storm though slightly overpredicted. Simulated core vertical temperature structure, winds at different heights, vertical winds in and around the core, vorticity and divergence fields at the lower and upper levels—all support the characteristics of a mature storm. The model storm has moved towards the west of the observed track during the development phase although the location of the storm in the initial and final phases agreed with the observations. The simulated rainfall distribution associated with the storm agreed reasonably with observations.     

Study on snowmelt runoff simulation in the Kaidu River basin

Alpine snowmelt is an important generation mode for runoff in the source region of the Tarim River basin, which covers four subbasins characterized by large area, sparse gauge stations, mixed runoff supplied by snowmelt and rainfall, and remarkably spatially heterogeneous precipitation. Taking the Kaidu River basin as a research area, this study analyzes the influence of these characteristics on the variables and parameters of the Snow Runoff Model and discusses the corresponding determination strategy to improve the accuracy of snowmelt simulation and forecast. The results show that: (i) The temperature controls the overall tendency of simulated runoff and is dominant to simulation accuracy, as the measured daily mean temperature cannot represent the average level of the same elevation in the basin and that directly inputting it to model leads to inaccurate simulations. Based on the analysis of remote sensing snow maps and simulation results, it is reasonable to approximate the mean temperature with 0.5 time daily maximum temperature. (ii) For the conflict between the limited gauge sta-tion and remarkably spatial heterogeneity of rainfall, it is not realistic to compute rainfall for each elevation zone. After the measured rainfall is multiplied by a proper coefficient and adjusted with runoff coefficient for rainfall, the measured rainfall data can satisfy the model demands. (iii) Adjusting time lag according to the variation of snowmelt and rainfall position can improve the simulation precision of the flood peak process. (iv) Along with temperature, the rainfall increases but cannot be completely monitored by limited gauge stations, which results in precision deterioration.     

Nonlinear numerical method for earthquake site response analysis II — case studies

This paper presents two numerical case studies of medium and strong motion events, namely Loma-Prieta 1989 and Hyogoken-Nambu (Kobe) 1995. These simulations were performed using CyberQuake model. The cyclic elastoplastic constitutive model is fully detailed in the companion paper. Through these case studies, we demonstrate the importance of using appropriate constitutive modelling when the part played by nonlinear phenomena is preponderant. The need to account for 3D kinematics (i.e. the three components of the input motion), is also demonstrated, even though a 1D geometry is considered, as the plastic coupling existing between components of motion during the earthquake, strongly affects the seismic soil response.