Discussion on the feature of strong earthquake: Orderly distribution in time, space and intensity before the Western Kunlun Mountain Pass M=8.1 earthquake

Introduction The Western Kunlun Mountain Pass M=8.1 earthquake occurred on November 14, 2001 is the other M=8 earthquake occurred 50 years after Dangxiong, Tibet M=8.0 earthquake in Chinese mainland. The earthquake has caused the attention of the seismologists in the following aspects: 1) The fracture length is more than 400 km, which is far away from the estimated length by the statistic empirical function between the magnitude and the fracture length (WANG, et al, 2002); 2) The aftersh…     

Petrography and uplift history of the Quaternary Takidani Granodiorite: could it have hosted a supercritical (HDR) geothermal reservoir?

The Quaternary Takidani Granodiorite (Japan Alps) is analogous to the type of deep-seated (3–5 km deep) intrusive-hosted fracture network system that might support (supercritical) hot dry/wet rock (HDR/HWR) energy extraction. The I-type Takidani Granodiorite comprises: porphyritic granodiorite, porphyritic granite, biotite-hornblende granodiorite, hornblende-biotite granodiorite, biotite-hornblende granite and biotite granite facies; the intrusion has a reverse chemical zonation, characterized by >70 wt% SiO₂ at its inferred margin and <67 wt% SiO₂ at the core. Fluid inclusion evidence indicates that fractured Takidani Granodiorite at one time hosted a liquid-dominated, convective hydrothermal system, with <380°C, low-salinity reservoir fluids at hydrostatic (mesothermal) pressure conditions. ‘Healed’ microfractures also trapped >600°C, hypersaline (35 wt% NaCl_eq) fluids of magmatic origin, with inferred minimum pressures of formation being 600–750 bar, which corresponds to fluid entrapment at 2.4–3.0 km depth. Al-in-hornblende geobarometry indicates that hornblende crystallization occurred at about 1.45 Ma (7.7–9.4 km depth) in the (marginal) eastern Takidani Granodiorite, but later (at 1.25 Ma) and shallower (6.5–7.0 km) near the core of the intrusion. The average rate of uplift across the Takidani Granodiorite from the time of hornblende crystallization has been 5.1–5.9 mm/yr (although uplift was about 7.5 mm/yr prior to 1.2 Ma), which is faster than average uplift rates in the Japan Alps (3 mm/yr during the last 2 million years). A temperature–depth–time window, when the Takidani Granodiorite had potential to host an HDR system, would have been when the internal temperature of the intrusive was cooling from 500°C to 400°C. Taking into account the initial (7.5 mm/yr) rate of uplift and effects of erosion, an optimal temperature–time–depth window is proposed: for 500°C at 1.54–1.57 Ma and 5.2±0.9 km (drilling) depth; and 400°C at 1.36–1.38 Ma and 3.3±0.8 km (drilling) depth, which is within the capabilities of modern drilling technologies, and similar to measured temperature–depth profiles in other active hydrothermal systems (e.g. at Kakkonda, Japan).     

Tensor time domain electromagnetic resistivity measurements at Ngatamariki geothermal field, New Zealand

Experimental measurements in the Ngatamariki geothermal field, North Island, New Zealand were made to test the applicability of the time domain electromagnetic method for detailed investigation of the resistivity structure within a geothermal field. Low-frequency square wave signals were transmitted through three grounded bipole current sources sited about 8 km from the measurement lines. Despite high levels of electrical noise, transient electric field vectors could be determined reliably for times between 0.02 and 3.3 s after each step in the source current. Instantaneous apparent resistivity tensors were then calculated. Apparent resistivity pseudosections along the two measurement lines show smooth variations of resistivity from site to site. Over most of the field the images consistently show a three-layer resistivity structure with a conductive middle layer (3–10 Ωm) representing the conductive upper part of the thermal reservoir. A deep-seated region of low resistivity in the northwest of the field may indicate a conductive structure at about 1 km associated with a deeper diorite intrusion. Measurements sited closer than about 100 m to drillholes appear to have been disturbed by metallic casing in the holes. A change in resistivity structure in the east of the field may indicate a major geological or hydrothermal boundary.     

Optimum multiple tuned mass dampers for structures under the ground acceleration based on the uniform distribution of system parameters

The five MTMD models, with natural frequencies being uniformly distributed around their mean frequency, have been recently presented by the first author. They are shown to have the near‐zero optimum average damping ratio (more precisely, for a given mass ratio there is an upper limit on the total number, beyond which the near‐zero optimum average damping ratio occurs). In this paper, the eight new MTMD models (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1～US‐MTMD3, UD‐MTMD1 and UD‐MTMD2), with the system parameters (mass, stiffness and damping coefficient) being, respectively, uniformly distributed around their average values, have been, for the first time here, proposed to seek for the MTMD models without the near‐zero optimum average damping ratio. The structure is represented by the mode‐generalized system corresponding to the specific vibration mode that needs to be controlled. Through minimization of the minimum values of the maximum dynamic magnification factors (DMF) of the structure with the eight MTMD models (i.e. through the implementation of Min.Min.Max.DMF), the optimum parameters and values of Min.Min.Max.DMF for these eight MTMD models are investigated to evaluate and compare their control performance. The optimum parameters include the optimum mass spacing, stiffness spacing, damping coefficient spacing, frequency spacing, average damping ratio and tuning frequency ratio. The six MTMD models without the near‐zero optimum average damping ratio (i.e. the UM‐MTMD1～UM‐MTMD3, US‐MTMD1, US‐MTMD2 and UD‐MTMD2) are found through extensive numerical analyses. Likewise, the optimum UM‐MTMD3 offers the higher effectiveness and robustness and requires the smaller damping with respect to the rest of the MTMD models in reducing the responses of structures subjected to earthquakes. Additionally, it is interesting to note, by comparing the optimum UM‐MTMD3 with the optimum MTMD‐1 recently investigated by the first author, that the effectiveness and robustness for the optimum UM‐MTMD3 is almost identical to that for the optimum MTMD‐1 (without inclusion of the optimum MTMD‐1 with the near‐zero optimum average damping ratio). Recognizing these performance benefits, it is preferable to employ the optimum UM‐MTMD3 or the optimum MTMD‐1 without the near‐zero optimum average damping ratio, when installing the MTMD for the suppression of undesirable oscillations of structures under earthquakes. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Predictive instantaneous optimal control of elastic structures during earthquakes

A predictive instantaneous optimal control (PIOC) algorithm is proposed for controlling the seismic responses of elastic structures. This algorithm compensates for the time delay that happens in practical control applications by predicting the structural response over a period that equals the time delay, and by substituting the predicted response in the instantaneous optimal control (IOC) algorithm. The unique feature of this proposed PIOC algorithm is that it is simple and at the same time compensates for the time delay very effectively. Numerical examples of single degree of freedom structures are presented to compare the performance of PIOC and IOC systems for various time delay magnitudes. Results show that a time delay always causes degradation of control efficiency, but PIOC can greatly reduce this degradation compared to IOC. The effects of the structure's natural periods and the choice of control gains on the degradation induced by the time delay are also analyzed. Results show that shorter natural periods and larger control gains are both more sensitive and more serious to the degradation of control efficiency. Finally, a practical application of PIOC is performed on a six‐story moment‐resisting steel frame. It is demonstrated that PIOC contributes significantly to maintain stability in multiple degree of freedom structures, and at the same time PIOC has a satisfactory control performance. Copyright © 2003 John Wiley & Sons, Ltd.     

Predictive instantaneous optimal control of inelastic structures during earthquakes

A new inelastic structural control algorithm is proposed by incorporating the force analogy method (FAM) with the predictive instantaneous optimal control (PIOC) algorithm. While PIOC is very effective in compensating for the time delay for elastic structures, the FAM is highly efficient in performing the inelastic analysis. Unlike conventional inelastic analysis methods of changing stiffness, the FAM analyzes structures by varying the structural displacement field, and therefore the state transition matrix needs to be computed only once. This greatly simplifies the computation and makes inelastic analysis readily applicable to the PIOC algorithm. The proposed algorithm compensates for the time delay that happens in practical control systems by predicting the inelastic structural response over a period that equals the magnitude of the time delay. A one‐story frame with both strain‐hardening and strain‐softening inelastic characteristics is analyzed using this algorithm. Results show that the proposed control algorithm is feasibile for any inelastic structures. While the control efficiency deteriorates with the increase in magnitude of the time delay, the PIOC maintains acceptable performance within a wide range of time delay magnitudes. Finally, a computer model of a six‐story moment‐resisting steel frame is analyzed to show that PIOC has good control results for real inelastic structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of predictors of non‐linear seismic demands using ‘fishbone’ models of SMRF buildings

Predictors (or estimates) of seismic structural demands that are less computationally time‐consuming than non‐linear dynamic analysis can be useful for structural performance assessment and for design. In this paper, we evaluate the bias and precision of predictors that make use of, at most, (i) elastic modal vibration properties of the given structure, (ii) the results of a non‐linear static pushover analysis of the structure, and (iii) elastic and inelastic single‐degree‐of‐freedom time‐history analyses for the specified ground motion record. The main predictor of interest is an extension of first‐mode elastic spectral acceleration that additionally takes into account both the second‐mode contribution to (elastic) structural response and the effects of inelasticity. This predictor is evaluated with respect to non‐linear dynamic analysis results for ‘fishbone’ models of steel moment‐resisting frame (SMRF) buildings. The relatively small number of degrees of freedom for each fishbone model allows us to consider several short‐to‐long period buildings and numerous near‐ and far‐field earthquake ground motions of interest in both Japan and the U.S. Before doing so, though, we verify that estimates of the bias and precision of the predictor obtained using fishbone models are effectively equivalent to those based on typical ‘full‐frame’ models of the same buildings. Copyright © 2003 John Wiley & Sons, Ltd.     

时间谱电阻率法中的剩余电磁效应研究

根据频谱激电（SIP）法或复电阻率（CR）法中研究剩余电磁（REM）效应的思想，通过分析时间谱电阻率（TSR）法中场的基本特征，在TSR法中建立了与SIP（或CR）法中参数Φmax/Φ0max相类似的参数－Emax/E0，用来描述了TSR法中的REM效应。通过对一维可极化大地模型的REM效应的理论计算，验证了用Emax/E0表示的REM异常可用来形象地反映地电断面电阻率的高低。     

时频测控技术的发展

目前，在所有的物理量中时间与频率标准的物理实现具有最高的准确度和稳定度，它们的发展不但在计量技术中，而且在几乎整个高科技领域中起着十分重要的作用。这方面的发展一直受到发达国家的高度重视。根据近年来与国外同行的技术交流和了解的情况，对时频测控技术的发展和影响等作了介绍。