Identification of soil degradation during earthquake excitations by Bayesian inference

A Bayesian inference approach is introduced to identify soil degradation behaviours at four downhole array sites. The approach of inference is based on a parametric time‐varying infinite impulse response filter model. The approach is shown to be adaptive to the changes of filter parameters and noise amplitudes. Four sites, including the Lotung (Taiwan), Chiba (Japan), Garner Valley (California), and Treasure Island (California) sites with downhole seismic arrays are analysed. Our results show two major types of soil degradation behaviour: the well‐known strain‐dependent softening, and reduction in stiffness that is not instantaneously recoverable. It is also found that both types of soil degradation are more pronounced in sandy soils than in clayey soils. The mechanism for the second type of soil degradation is not yet clear to the authors and suggested to be further studied. Copyright © 2003 John Wiley & Sons, Ltd.     

A quantitative study of the energy release in the aftershocks of the Bhuj earthquake, 2001, India, using Lg phase

The devastating earthquake on 26 January 2001 at Bhuj, India, resulted in large-scale death and destruction of properties of several million US dollars. The moment magnitude of the earthquake was 7.7 and its maximum focal intensity exceeded X in MM scale. The rate of aftershocks of this earthquake, recorded at Gauribidanur seismic array station (GBA), shows a monotonic decay with time superposed with oscillations. For the Indian continent the Lg phase is a prominent arrival at regional distances. The estimate of Lg amplitude is obtained by optimally fitting the Lg wave train to a exponential decay curve. The logarithm of these amplitudes and logarithm of root mean square (rms) value of actual amplitudes of the Lg are calibrated with USGS m_b to create a local m_bLg magnitude scale. The energy released from these aftershocks is calculated from the rms value of Lg phase. The plot of cumulative energy release with time follows the power law of the form t^p, superposed with oscillations. The exponent of the power law, p, is estimated both by a time-window scanning method and by an interpolation method. The value of p is 0.434 for time-window scanning method and 0.432 for the interpolation method. The predominant periods found in the oscillatory part of the cumulative energy, obtained by differencing the observed from the power law fit, are 10.6, 7.9, 5.4, 4.6 and 3.5 h for time-window scanning method. The corresponding periods for interpolation method are 13.4, 11.5, 7.4, 4.2, 3.5, 2.6 and 2.4 h.     

基于遗传算法的最优证券投资组合模型

为了研究能使风险损失率最小、收益最大的最优证券投资组合问题，首先提出了多目标规划模型，并采用模糊优选法将多目标模型单目标化。随后尝试采用遗传算法来对其求解，给出了模型的遗传算法编码规则和算法步骤。最后通过实例论证了模型的合理性及遗传算法的有效性。     

声波测试中振子衰减的线列阵模型

从声波在岩体中传播产生衰减的基本原理出发，建立了声波衰减振子的线列阵模型。该模型是将传播声波路径的岩体分为若干个振子，从波动概念出发，以振动原理作依据建立波动方程．为探讨岩体声波衰减机制提供了一种新的研究方法，把振动和波动直观地联系起来．单一方向上有限点的声波测试技术有可能被多方向的声测技术所代替．该模型的建立，使模态分析和参数识别方法应用到岩体声波测试技术中成为可能。     

影响上海市地价空间分布的区位因子分析

本文采用ＧＩＳ技术和多元线性回归分析方法，以土地批租地块为样本，分析上海市地价与城市区位因子的线性关系。回归分析结果表明；南京路，市中心和淮海路是影响上海市寺价空间差异的主要区位因子，外滩，次级商业中心和交通结节点对地价的影响很小，表明上海市地价的空间分布具有很高的向心性。     

Phase-array decomposition of a seismic section

A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)     

Noise reduction and detection of weak, coherent signals through phase-weighted stacks

We present a new tool for efficient incoherent noise reduction for array data employing complex trace analysis. An amplitude-unbiased coherency measure is designed based on the instantaneous phase, which is used to weight the samples of an ordinary, linear stack. The result is called the phase-weighted stack (PWS) and is cleaned from incoherent noise. PWS thus permits detection of weak but coherent arrivals. The method presented can easily be extended to phase-weighted cross-correlations or be applied in the τ p domain. We illustrate and discuss the advantages and disadvantages of PWS in comparison with other coherency measures and present examples. We further show that our non-linear stacking technique enables us to detect a weak lower-mantle P -to- S conversion from a depth of approximately 840 km on array data. Hints of an 840 km discontinuity have been reported; however, such a discontinuity is not yet established due to the lack of further evidence.     

运用X射线衍射和多重峰分离程序解析高岭石的结构缺陷

本文尝试运用Ｘ射线衍射与多重峰分离方法研究高岭石的结构缺陷，重新评价了Ｈｉｎｃｋｌｅｙ（１９６３）结晶指数（ＨＩ），运用０２，１１区域对（１１１）和（１１１）峰的分峰结果，定义了新的结晶度指数：ＣＩ＝（Ｉ１１１＋Ｉ１１１）／Ｉ１１０，它完全以衍射峰的强度进行计算，较灵敏地反映出高岭石真实的结构信息，探讨了高岭石结晶度与高峰石对称性，有序度以及与高岭石Ｐｂ／３滑移缺陷的关系，结果证明高岭石结晶度越高     

Site effects at a vertical accelerometer array near Paducah, Kentucky

The downhole vertical accelerometer array VSAP near Paducah, KY, consists of three-component accelerometers at the surface, the top of the McNairy Formation (−41 m), and the top of the Paleozoic bedrock (—99 m). The array is at the northern end of the Mississippi Embayment, and it was installed to verify the ground-motion modeling for the site, assuming a significant earthquake in the New Madrid Seismic Zone. Accelerograms from 4.2 and 2.0 m_b earthquakes were used to check aspects of the modeling pertaining to linear behavior of the soil column, and to review the soil column models derived by drilling and geotechnical methods and through the use of high-resolution P- and SH-wave seismic refraction in reflection techniques. Results of the study indicate that for the linear case the soil column models derived by the two techniques are equivalent, and that the most important boundary in the soil column, with respect to amplification of the ground motions, is the interface between the limestone bedrock and soil.     

Multiple episodes of tectonothermal processes in the Eastern Ghats granulite belt

A suite of rocks from Borra Carbonate Granulite Complex (BCGC) in the Eastern Ghats granulite belt displays superposed structures and overprinted mineral assemblages that reveal multiple episodes of tectonothermal reworking of the complex under granulite facies condition. Five distinct episodes of deformation (D₁, D₂, D₃, D₄ and D₅) and four phases of metamorphism (M₁, M₂, M₃ and M₄) are recorded. The signature of the earliest tectonothermal event, D₁ is a gneissic foliation (S₁) denned by segregation of peak granulite facies mineral assemblages corresponding to prograde M₁ metamorphism. M₂ metamorphic overprint represents an episode of near-isobaric cooling of the complex under a static condition. D₂ represents an episode of ductile deformation manifested by isoclinal folding (F₂) and associated extensional structures, within a broad framework of coaxial bulk deformation. The present study reveals that D₂ took place subsequent to M₂ - Subsequent deformation, D₃, produced F₃ folds and also deformations of boudins formed during D₂. M₃, which is synchronous with F₃, represents a near isothermal decompression of the BCGC. This was followed by a weak structural readjustment (D₄), producing E-W cross folds. The latter was not, however, associated with any recognizable petrological reworking. In the terminal events, deformation (D₅) and mineral reactions (M₄) were localized along narrow intersecting shear zones. The latter acted as channelways for carbonic and still later hydrous fluid infiltration. The available thermobarometric data from BCGC and other areas of the Eastern Ghats belt reveal that reworking during M₂ and M₃ ensued in a thermally perturbed regime. The high thermal regime might also have persisted during carbonic fluid infiltration related to terminal reworking (M₄).