软岩巷道弹塑性变形的理论分析

本文采用莫尔－库仑屈服准则和非关联塑性流动法则分析了轴对称软岩巷道的变形规律。从轴对称平面应变问题的基本方程出发，导出了软岩巷道弹塑性变形的位移理论解答，并给出了现有应力与位移理论解的适用条件。通过实例计算结果表明，本文得出的位移理论解与有限元数值解相当吻合。     

Nonequilibrium thermodynamics of pressure solution

This paper is concerned with the thermodynamic theory of solution and precipitation processes in wet crustal rocks and with the mechanism of steady pressure-solution slip in contact zones, such as grain-to-grain contacts, fracture surfaces, and permeable gouge layers, that are infiltrated by a mobile aqueous solution phase. A local dissipation jump condition at the phase boundary is fundamental to identifying the thermodynamic force driving the solution and precipitation process and is used here in setting up linear phenomenological relations to model near-equilibrium phase transformation kinetics. The local thermodynamic equilibrium of a stressed pure solid in contact with its melt or solution phase is governed by Gibbs's relation, which is rederived here, in a manner emphasizing its independence of constitutive assumptions for the solid while neglecting surface tension and diffusion in the solid. Fluid-infiltrated contact zones, such as those formed by rough surfaces, cannot generally be in thermodynamic equilibrium, especially during an ongoing process of pressure-solution slip, and the existing equilibrium formulations are incorrect in overlooking dissipative processes tending to eliminate fluctuations in superficial free energies due to stress concentrations near asperities, defects, or impurities. Steady pressure-solution slip is likely to exhibit a nonlinear dependence of slip rate on shear stress and effective normal stress, due to a dependence of the contact-zone state on the latter. Given that this dependence is negligible within some range, linear relations for pressure-solution slip can be derived for the limiting cases of diffusion-controlled and interface-reaction-controlled rates. A criterion for rate control by one of these mechanisms is set by the magnitude of the dimensionless quantityk/2C _pD, wherek is the interfacial transfer coefficient, is the mean diffusion path length,C _p is the solubility at pressurep, andD is the mass diffusivity.     

Diffusion models of slope development

This paper represents a systematic investigation of slope evolution diffusion models and has the following sections: (1) The model of slope development with linear coefficient k = k₀x; (2) The model of slope development with quadratic coefficient in x; (3) Slope development model with vertical lowering of base level (downcutting); (4) Slope development model with the base level a horizontal variable (undercutting); (5) Steady-state regime of undercut slopes; (6) Model of a pediment and scree slope formation. The comparison is made between mathematical and classical methods of slope evolution analysis.     

Simulation of the hydrocarbon structures using the P-SV wave solution

Numerical solution of the scalar and elastic wave equations has greatly aided geophysicists in modeling seismic wave fields in the complicated geologic structures containing hydrocarbons and hence increases the geologic interpretation. Finite-difference method offers a versatile approach to compute synthetic seismograms numerically for given subsurface complex geological structures. To avoid the spatial derivative of the elastic parameters and density, elastodynamic equation (first-order hyperbolic equation) has been solved using the Lax-Wendroff scheme. A numerical finite-difference modeling program has been developed for the P-SV wave using the above solution. A line source with a time delay of 0.015s and dominant frequency of 120 Hz has been utilized in the simulation. In order to avoid the large values of the displacement vector in the source region,Alterman andKaral's method (1968) has been utilized. Horizontal and vertical component synthetic seismograms have been computed for two different geological models with and without oil and gas bearing zones. It has been concluded from the response that a finite-difference technique not only yields the relative arrival times but also accounts for the variation in amplitude and phase according to the elastic impedance contrast across the interfaces. It should come as no surprise to learn that in spite of the limitation of this numerical method, the scheme has provided a valid response for the thin layer, high acoustic impedance contrast and the pinch out.     

求震源机制P波初动解的格点尝试概率法

考虑到Ｐ波初动数据在震源球面上分布的不均匀性对节面解的可能影响，不用矛盾符号所占比例的大小作为选择可取节面解的标准，而是将震源球面分成许多小块，再确定有Ｐ波初动数字的小块为正号区或负号区的概率，最后以全部有数据小块的平均概率的高低作为选择解答的标准。用此改进方法重新测定了徐纪仁等（Ｘｕ ｅｔ．ａｌ，，１９８８）发表的青藏高原地区９９个地震的震源机制解，从中筛选出了２０个可信的解答。     

千岛岛弧大震前哈佛大学矩心矩张量（CMT）解一致性的预测意义s

千岛岛弧地区属于全球地震活动最为活跃的地区之一. 本文利用哈佛大学测定的千岛岛弧地区地震的矩心矩张量（CMT）解， 分析该地区震源机制的一致性特征， 提出利用震源机制和构造应力场的一致性参数a进行地震预测的思路. 研究结果表明，MＷge;7.5地震之前， 都有一致性参数a降低的现象，a的低值起始的时间在发生大震之前的10多天至110多天， a的低值截止的日期距大震在30多天至2天. 相互之间虽然并不完全一致， 但是差别不大. 这种现象的稳定性， 尚需时间的检验， 但是设想在长达数百公里的地区， 连续发生MＷge;5.3的地震的震源机制都与构造应力场一致， 应当不是随机的现象， 而是一种具有预测意义的现象. 当积累的震例足够多时， 则有可能确定统一的评判标准和预测准则.      

Buckling of Piles under Liquefied Soil Conditions

The paper pertains to the analysis of piles embedded in liquefiable soils to predict its’ critical buckling load under partial to full loss of lateral support over a portion of the pile length. The analysis is based on extension of Mindlin solution for a point load acting inside a semi infinite elastic half space. Degenerated solutions obtained by using the developed method compares very well with reported results. Parametric studies showed that the depth of liquefiable soil, degradation of soil strength on liquefaction, slenderness ratio, pile stiffness factor and end conditions have significant influence on the buckling behavior of the piles.     

Focal mechanisms of recent earthquakes in the Southern Korean Peninsula

We evaluate the stress field in and around the southern Korean Peninsula with focal mechanism solutions, using the data collected from 71 earthquakes ( M_L = 1.9–5.2) between 1999 and 2004. For this, the hypocentres were relocated and well-constrained fault plane solutions were obtained from the data set of 1270 clear P -wave polarities and 46 SH / P amplitude ratios. The focal mechanism solutions indicate that the prevailing faulting types in South Korea are strike-slip-dominant-oblique-slip faultings with minor reverse-slip component. The maximum principal stresses (σ₁) estimated from fault-slip inversion analysis of the focal mechanism solutions show a similar orientation with E–W trend (269°–275°) and low-angle plunge (10°–25°) for all tectonic provinces in South Korea, consistent with the E–W trending maximum horizontal stress (σ_Hmax) of the Amurian microplate reported from in situ stress measurements and earthquake focal mechanisms. The directions of the intermediate (σ₂) and minimum (σ₃) principal stresses of the Gyeongsang Basin are, however, about 90 deg off from those of the other tectonic provinces on a common σ₂–σ₃ plane, suggesting a permutation of σ₂ and σ₃. Our results incorporated with those from the kinematic studies of the Quaternary faults imply that NNW- to NE-striking faults (dextral strike-slip or oblique-slip with a reverse-slip component) are highly likely to generate earthquakes in South Korea.     

Two-dimensional power series solution for non-axisymmetrical transport in a radially convergent tracer test with scale-dependent dispersion

It has been known for many years that dispersivity increases with solute travel distance in a subsurface environment. The increase of dispersivity with solute travel distance results from the significant variation of hydraulic properties of heterogeneous media and was identified in the literature as scale-dependent dispersion. This study presents an analytical solution for describing two-dimensional non-axisymmetrical solute transport in a radially convergent flow tracer test with scale-dependent dispersion. The power series technique coupling with the Laplace and finite Fourier cosine transform has been applied to yield the analytical solution to the two-dimensional, scale-dependent advection–dispersion equation in cylindrical coordinates with variable-dependent coefficients. Comparison between the breakthrough curves of the power series solution and the numerical solutions shows excellent agreement at different observation points and for various ranges of scale-related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at any observation point.