中小震近场地震动估计中地震动衰减关系的适用性分析

核工程场地设计基准地震动确定中涉及弥散地震的地震动估计，这导致中小震近场地震动估计成为被关注的研究课题。本文基于一些中小震近场地震动记录资料，分析和探讨了国内外多种地震动衰减关系对中小震近场地震动估计的适用性问题，并结合我国几个核工程地震问题研究工作的实践对中小震近场地震动衰减关系确定的思路和原则提出了建议。     

Water vapour transport of the pre-monsoon period and the general performance of the Indian summer monsoon

The advective monthly mean transfer of water vapour in the layer below 700 mbar is investigated for India for the years 1962 to 1972 and for the months January to September. The average zonal and meridional components of the transfer of water vapour for India are obtained. They are further averaged for different combinations of the pre-monsoon months from January to May and are correlated with the summer monsoon rainfall. The correlation coefficients for zonal transfer of water vapour are either negative or small positive for different combinations of the months mentioned above. The correlation coefficients for the meridional transfer of water vapour are positive. The maximum value is 0.74 for the March to May combination and is statistically significant at the 1% level.An extensive investigation is, therefore, made for the March to May averages of water vapour transfer for four broad regions of India. The parameters of water vapour transfer for these regions are compared with the threshold values and the prediction category, normal or drought, for the subsequent summer monsoon season is determined for all years. The correlation coefficient between the index of drought, as determined from parameters of water vapour transport and rainfall departure, is statistically significant at the 2% level.     

Determination of GMPE functional form for an active region with limited strong motion data: application to the Himalayan region

Advancement in the seismic networks results in formulation of different functional forms for developing any new ground motion prediction equation (GMPE) for a region. Till date, various guidelines and tools are available for selecting a suitable GMPE for any seismic study area. However, these methods are efficient in quantifying the GMPE but not for determining a proper functional form and capturing the epistemic uncertainty associated with selection of GMPE. In this study, the compatibility of the recent available functional forms for the active region is tested for distance and magnitude scaling. Analysis is carried out by determining the residuals using the recorded and the predicted spectral acceleration values at different periods. Mixed effect regressions are performed on the calculated residuals for determining the intra- and interevent residuals. Additionally, spatial correlation is used in mixed effect regression by changing its likelihood function. Distance scaling and magnitude scaling are respectively examined by studying the trends of intraevent residuals with distance and the trend of the event term with magnitude. Further, these trends are statistically studied for a respective functional form of a ground motion. Additionally, genetic algorithm and Monte Carlo method are used respectively for calculating the hinge point and standard error for magnitude and distance scaling for a newly determined functional form. The whole procedure is applied and tested for the available strong motion data for the Himalayan region. The functional form used for testing are five Himalayan GMPEs, five GMPEs developed under NGA-West 2 project, two from Pan-European, and one from Japan region. It is observed that bilinear functional form with magnitude and distance hinged at 6.5 M _w and 300 km respectively is suitable for the Himalayan region. Finally, a new regression coefficient for peak ground acceleration for a suitable functional form that governs the attenuation characteristic of the Himalayan region is derived.     

On the estimation of accuracy for global models of gravitational field of the earth

The methods and techniques for estimating the accuracy of global models of the Earth’s gravity field in the form of spherical harmonic expansion of the geopotential are analyzed. Various methods for obtaining the a priori and a posteriori estimates for the accuracy are considered and classified. The application of different approaches is illustrated by numerical examples for nine models, including those recently developed using the modern methods of space geodesy. The basic requirements for the database and software for estimating the accuracy are formulated.     

从抗震设防类别探讨室内地震应急避难场所的选择

区分建筑物的抗震设防类别是进行抗震设计的前提和依据。本文针对当前我国室内地震应急避难场所认定面临的现有建筑的抗震安全性问题,根据不同版本的《建筑工程抗震设防分类标准》和《建筑抗震设计规范》,对现有体育、会展、教育建筑等设计、施工时依据的抗震设防类别、标准、规范进行分析研究,以为室内地震应急避难场所的选择提供参考,从而推动室内地震应急避难场所的认定工作。     

Development of marine water quality criteria for the USA

The US Environmental Protection Agency has developed guidelines for deriving numerical national water quality criteria for the protection of aquatic organisms and their uses. These guidelines provide the method for deriving water quality criteria, including minimum data base requirements, data evaluation procedures, and calculations. The guidelines have been in place for a long time, and states have used them to derive water quality standards for their water bodies. More recent efforts have been directed towards the development of technical guidance based on the concept that bioassessment and biocriteria programs for estuaries and near coastal waters are interrelated and are critical components of comprehensive water resource protection and management. This is a holistic approach to protection and management, integrating biological assessments into traditional chemical and physical evaluations and augmenting the established water quality criteria. The method for deriving water quality criteria and the approach for biocriteria development for marine systems are described.     

Changes in the scope for growth of blood cockles (Anadara granosa) exposed to industrial discharge

Blood cockles (Anadara granosa) were placed in plastic cages and deployed at eight locations along a pollution gradient. A control site was chosen a few kilometres away. At 2-week intervals for 8 weeks, the scope for growth of the cockles was determined. A total of 11 environmental parameters were measured during each sampling exercise for correlation with the growth index. Values of scope for growth for all the eight stations, except for stations 7 and 8, were significantly lower than the value for the control site; station 1 (nearest the discharge point) had the lowest value. Strong correlations were found between the scope for growth of the cockles and seven of the environmental parameters. The pattern of the impact of the discharge on the scope for growth of the cockles suggests that longshore currents may play a dominant role in the distribution of the industrial wastes.     

Engineering-geological substantiation of water resource protection in the zones of oil and gas fields

The hydrochemical background and the specific features of pollution of geological environmental components by hydrocarbons and chemical components in the industrial sites of oil and gas fields is characterized for two types of natural-anthropogenic environments with different hydrogeological conditions. A methodology is proposed for comprehensive environmental-hydrogeological studies for the development of information support for the assessment of the hazard and risk for the environmental conditions, organization of complex environmental monitoring, and projecting of measures for protection of water resources.     

Unified fractional differential approach for transient interporosity flow in naturally fractured media

A unified approach to modeling flows of slightly compressible fluids through naturally fractured media is presented. The unified fractional differential model is derived by combining the flow at micro scale for matrix blocks and macro scale for fractures, using the transient interporosity flow behavior at the interface between matrix blocks and fractures. The derived model is able to unify existing transient interporosity flow models formulated for different shapes of matrix blocks in any medium dimensions. The model is formulated in the form of a fractional order partial differential equation that involves Caputo derivative of order 1/2 with respect to time. Explicit solutions for the unified model are derived for different axisymmetrical spatial domains using Hankel or Hankel–Weber finite or infinite transforms. Comparisons between the predictions of the unified model and those obtained from existing transient interporosity flow models for matrix blocks in the form of slabs, spheres and cylinders are presented. It is shown that the unified fractional derivative model leads to solutions that are very close to those of transient interporosity flow models for fracture-dominant and transitional fracture-to-matrix dominant flow regimes. An analysis of the results of the unified model reveals that the pressure varies linearly with the logarithm of time for different flow regimes, with half slope for the transitional fracture-to-matrix dominant flow regime vs. the fracture and matrix dominant flow regimes. In addition, a new re-scaling that involves the characteristic length in the form of matrix block volume to surface area ratio is derived for the transient interporosity flow models for matrix blocks of different shapes. It is shown that the re-scaled transient interporosity flow models are governed by two dimensionless parameters Θ and Λ compared to only one dimensionless parameter Θ for the unified model. It is shown that the solutions of the transient interporosity flow models for different shapes of matrix blocks are almost identical for the re-scaled variables. Furthermore, the driving parameters for solution behavior are identified based on asymptotic approximations for different flow regimes. It is found that the matrix diffusion and the matrix area-to-volume ratio affect the solution behavior only for the transitional fracture-to-matrix dominant flow regime, that the capacitance ratio affects the solution behavior only for transitional and matrix dominant flow regimes and that the fracture diffusion is involved in all three flow regimes. Similar identification of the driving parameters is also presented in the re-scaled case.     

Magnitude assessment of northwestern european earthquakes

A homogeneous and internally consistent body of regional magnitudes is needed not only for the statistical study of seismicity but also for the modelling of seismic sources and for the hazard assessment of engineering sites. Also the assessment of the design ground motions for engineering projects requires homogeneous magnitude estimates. The nature of the available data for U.K. and Northwestern European earthquakes, and of the events themselves, is such that we are concerned with the assessment of magnitudes at the lower end of the scale, and thus great importance focuses on the uniform calculation of magnitudes within a relatively narrow range of M values. The purpose of this study is to present the uniform re-assessment of magnitudes for British and Northwestern European earthquakes for the period 1900 to 1984. It is shown that, for the larger events, surface-wave magnitudes can be estimated uniformly and that neither a period constraint at 20 seconds nor distance effects for the regional conditions considered appear to play an important role. For smaller events, crustal phases at short distances have been used to derive calibration curves for a magnitude which, at larger distances and longer periods would fit the magnitude scale of surface waves.     

The dynamics of a rigid foundation on the surface of an elastic half-space

The response of a rigid rectangular foundation block resting on an elastic half-space has been determined by considering first the displacement functions for any position on the surface of an unloaded half-space due to a harmonic point force. The influence of the foundation has been taken into account by assuming a relaxed condition at the interface, i.e. a uniform displacement under the foundation and that the sum of the point forces must equal the total applied force. The three motions of vertical, horizontal and rocking have been considered and numerical values for the in-phase and the quadrature components of the displacement functions are presented for a Poisson's ratio of 0.25. The effect of the mass and inertia of the foundation can be allowed for by an impedance matching technique. Response curves and non-dimensional resonant frequency curves are given for a square and a rectangular foundation for different mass and inertia ratios and for several values of Poisson's ratio. These curves are for design purposes and are an addition to similar published curves for circular and infinitely long rectangular foundations. Some of the calculated results have been verified by a laboratory experiment.     

首都减灾圈的思考与建议

首都减灾圈,系首都减轻自然灾害预测防治圈。1991年12月20—21日在北京召开了《首都圈自然灾害及其减灾对策研讨会》。本文根据这次会议所提供的材料,在从整体上实现减灾的思想指导下,就首都减灾圈的“成灾背景”、“首都减灾圈的组成”、“首都减灾圈的灾害预测与防治状况”、“首都减灾圈的灾害关联性分析”、“首都减灾圈的减灾实效预估”和“首都减灾圈的减灾对策与实施”等六个带有共同性的问题,进行了讨论。供制定首都圈减灾方案时参考。     

Attenuation relations for horizontal peak ground acceleration and response spectrum in northeastern Tibetan Plateau region

Introduction Developing local attenuation relations of ground motion is one of the key steps in seismic hazard assessment. Because of inadequate strong ground motion records in China, the attenuation relations used in China are usually developed by using the transforming method (Hu, Zhang, 1983; HU, ZHANG, 1984). To use this method, we need to have both the attenuation relation of seismic intensity for the studied region and the attenuation relations of seismic intensity and ground mo-tion…     

A3M联测技术在岩芯分析中的研究及应用

总结国内外岩芯分析技术，综合应用现代控制论、信息论、系统工程论的基本原理，提出A3M理论性认识和联测技术的新概念，研制出一套具有多学科、多目的（任务）、多项目（参数）的岩芯综合分析技术，即岩性、物性、岩石电性、声波等联测的技术．完成了一套与联测技术流程相配套的实用软件包，这套技术提供了系统的、有代表性的、可对比的岩芯分析资料，可直接用于油田现场与实验室岩芯分析生产、科研和管理及油田的开发中，有利于认识非均质储层岩石地质、物理特征以及各种信息之间的内在规律．     

Logic-tree Approach for Probabilistic Tsunami Hazard Analysis and its Applications to the Japanese Coasts

For Probabilistic Tsunami Hazard Analysis (PTHA), we propose a logic-tree approach to construct tsunami hazard curves (relationship between tsunami height and probability of exceedance) and present some examples for Japan for the purpose of quantitative assessments of tsunami risk for important coastal facilities. A hazard curve is obtained by integration over the aleatory uncertainties, and numerous hazard curves are obtained for different branches of logic-tree representing epistemic uncertainty. A PTHA consists of a tsunami source model and coastal tsunami height estimation. We developed the logic-tree models for local tsunami sources around Japan and for distant tsunami sources along the South American subduction zones. Logic-trees were made for tsunami source zones, size and frequency of tsunamigenic earthquakes, fault models, and standard error of estimated tsunami heights. Numerical simulation rather than empirical relation was used for estimating the median tsunami heights. Weights of discrete branches that represent alternative hypotheses and interpretations were determined by the questionnaire survey for tsunami and earthquake experts, whereas those representing the error of estimated value were determined on the basis of historical data. Examples of tsunami hazard curves were illustrated for the coastal sites, and uncertainty in the tsunami hazard was displayed by 5-, 16-, 50-, 84- and 95-percentile and mean hazard curves.     

Streamflow simulation and skewness preservation based on the bootstrapped stochastic models

The stochastic model has been widely used for the simulation study. However, there was a difficulty in the reproduction of the skewness of observed series and so the stochastic model for the skewness preservation was appeared. While the skewness in the residuals of the stochastic model has been considered for the skewness preservation this study uses a random resampling technique of residuals from the stochastic models for the simulation study and for the investigation of the skewness coefficient. The main advantage of this resampling scheme, called the bootstrap method is that it does not rely on the assumption of population distribution and this study uses the combined model of the stochastic and bootstrapped models. The stochastic and bootstrapped stochastic (or combined) models are used for the investigations of skewness preservation and of the reproduction of probability density function between the simulated series. The models are applied to the annual and monthly streamflows of Yongdam site in Korea and Yakima river, Washington, USA for the streamflow simulation study then the statistics and probability density functions for the observed and simulated streamflows are compared. As the results the bootstrapped stochastic model reproduces the skewness and probability density function much better than the stochastic model. This evidences suggest that the bootstrapped stochastic model might be more appropriate than the stochastic model for the preservation of skewness and for simulation purposes of the series.     

Parallel distributed computing using Python

This work presents two software components aimed to relieve the costs of accessing high-performance parallel computing resources within a Python programming environment: MPI for Python and PETSc for Python.MPI for Python is a general-purpose Python package that provides bindings for the Message Passing Interface (MPI) standard using any back-end MPI implementation. Its facilities allow parallel Python programs to easily exploit multiple processors using the message passing paradigm. PETSc for Python provides access to the Portable, Extensible Toolkit for Scientific Computation (PETSc) libraries. Its facilities allow sequential and parallel Python applications to exploit state of the art algorithms and data structures readily available in PETSc for the solution of large-scale problems in science and engineering.MPI for Python and PETSc for Python are fully integrated to PETSc-FEM, an MPI and PETSc based parallel, multiphysics, finite elements code developed at CIMEC laboratory. This software infrastructure supports research activities related to simulation of fluid flows with applications ranging from the design of microfluidic devices for biochemical analysis to modeling of large-scale stream/aquifer interactions.     

围压对低渗砂岩渗透率有效应力系数的影响

多孔岩石渗透率有效应力系数的取值至今仍没有一致的认识,影响取值的因素较为复杂.为了探索围压对低渗砂岩渗透率有效应力系数的影响,实验上采用较大的围压范围(10MPa～50MPa),选取两块低渗致密砂岩岩样,进行渗透率测试,并利用二元线性回归方法,计算了不同围压组合下的渗透率有效应力系数值,获取了该系数随围压的变化趋势.实验中,两块岩样的渗透率有效应力系数随围压的增加而减小,但各自的减小趋势略有不同.利用Bernabé和Berryman提出的模型,理论上解释了围压对渗透率有效应力系数的影响,并得到渗透率有效应力系数的下限值为岩样孔隙度.尽管渗透率有效应力系数不是一个常数,但可以表示为压差的函数,进而起到简化分析的作用.     

太湖流域水文数学模型

本文针对太湖流域研制了全流域水文数学模型,该模型包括河流、湖泊模拟、边界条件模拟、降雨径流模拟、工程情况及控制运行方式模拟以及骨干河网中水流运动模拟五个方面。模型经1984、1985两年资料进行了上述五个方面的全面率定和检验,模拟结果与实测基本吻合。本模型全面通用,因此可以用来研究太湖流域洪水、枯水及调度等诸方面课题,本文简要地介绍应用该模型研究围垦的影响、规划工程的防洪效果及设计典型年选择三方面问题。     

CALCULATING THE DYNAMIC STIFFNESS MATRIX OF 2-D FOUNDATIONS BY DISCRETE WAVE NUMBER INDIRECT BOUNDARY ELEMENT METHODS

Apart from some special cases, calculating the dynamic stiffness matrix of foundations on a layered half-space, especially for embedded foundations, is computationally expensive. An efficient method for two-dimensional foundations in a horizontally layered soil media is presented in this paper. This method is based on indirect boundary element methods and uses discrete wave number solution methods for calculating Green's functions for displacements and analytical methods for the integrations over the boundary. For surface foundations, the present method applies at all frequencies. For embedded foundations or for constructing energy transmitting boundaries, because the free-field part is modelled by boundary elements and the excavated part is modelled by finite elements, the present method applies only at low frequencies for the spring coefficients (the real parts of the dynamic stiffness matrix) but applies at all frequencies for the damping coefficients (the imaginary part of the dynamic stiffness matrix) for undamped sites. The novelty of the method can be used for three-dimensional foundations. © 1997 by John Wiley & Sons, Ltd.