时间域航空电磁数据的反演

本文用广义逆矩阵理论讨论了时间域航空电磁数据的反演方法.以水平二层大地和球体模型为例进行了计算,结果表明这种方法对时间域电磁数据的反演是行之有效的.对于理论数据,一般仅需迭代数次即可稳定地收敛到真值.通过奇异值分析,可以确定一个模型中的相对重要参数和无关紧要参数;对于求解某一特定参数,可确定哪些数据起决定性的作用.由奇异值分解而得出的信息密度矩阵、模型分辨矩阵,分别给出了模型响应拟合数据和模型参数分辨率的度量.这些信息对于时间域航空电磁系统的设计、野外测量和资料解释都具有重要的参考价值.最后给出了实例.     

唐山余震震源参数和衰减的联合反演

利用唐山余震的数字化地震记录,研究了唐山余震的震源参数.发展了一种 P 波谱的反演方法,以消除用波谱分析求解震源参数过程中的不确定因素的影响.经过数值实验和对实际资料的应用,初步证明了这种方法的实用性.利用5个台站接收的35次地震记录的 P 波频谱资料,反演了唐山滦县地区的震源动力学参数,得到了5个台站的 P 波 Q 值和震源频谱的高频衰减率.Q 值分布在408——847之间,均值为520;值分布在1.54——3.22之间,均值为2.41.在该区域的小震频谱的定标律研究中,发现应力降随地震矩的增加而增加.      

COMPUTER PROCESSING OF AIRBORNE ELECTROMAGNETIC DATA *

With the increasing number of channels in AEM systems, computer data handling is becoming a necessity. The experience gained in processing of seismic and aeromagnetic data cannot be applied directly to low-frequency (100-5000 Hz) AEM methods. A novel scheme has been designed for AEM data processing and tested on 2900 km (1800 miles) of Input surveys. In the first step, the digital flight tapes are merged with digitized flight path recovery to form the primary data set. The validity of the raw data is controlled by the computer, but the interpreter has an option of checking them in perspective plots of channel amplitudes. The primary data set is reduced by processors which determine the location and type of anomalies and discard noise. Unlike the widely used deconvolution, the sequential processor determines first the anomaly location and then estimates parameters, such as peak amplitude, base width, and excess area, which are used as acceptance criteria. Interpretation parameters, such as σt, conductor depth, and dip are estimated by comparison with quantitative models. The recorded channel amplitudes are plotted together with the selected interpretation parameters in a profile form. The secondary data set which includes only the interpretation parameters for selected anomalies is graphically displayed as a schematic map of apparent σt. Elongated features are traced by a fan strike recognition routine and a trend map is automatically compiled and plotted. Disk storage is essential for second pass processing during which parts of the primary data set are searched for undiscovered anomalies matching the analyzed trend. The suggested procedure for AEM data processing is demonstrated on an Input MK V survey, Southern Indian Lake, Manitoba.     

ESTIMATION OF SOIL PARAMETERS FOR ASSESSING POTENTIAL WETNESS: COMPARISON OF MODEL RESPONSES THROUGH GIS

A geographic information system (GIS) is utilized to model wetness potential for a portion of Uwharrie National Forest, North Carolina. The wetness index is derived from TOPMODEL, a hillslope-scale runoff simulation model. The wetness index is a distributed-parameter model, with the input parameters obtained from a digital elevation model (DEM) and Soil Conservation Service (SCS) soils data. The primary objectives of the research are to: (1) compare methods of estimating soil parameters for input into the wetness potential model; and (2) determine how the model outputs vary spatially as a consequence of different methods of estimating soil parameters. Three methods of estimating soil parameters are used: (a) assuming uniform soil properties; (b) using SCS data presented as ranges; and (c) using alternative literature-based estimates of soil parameters. Results indicate that the wetness model responds similarly regardless of how the soil parameters are estimated, but differences in the spatial variability of the wetness potentials occur as a result of estimating soil parameters through alternative approaches. Correlation, pair-wise regression and analysis of regression residuals are used to compare model responses within a GIS environment.     

COMPARISON OF SIGNAL PROCESSING TECHNIQUES FOR ESTIMATING THE EFFECTS OF ANISOTROPY1

Three-component recordings of shear-waves in exploration surveys provide an opportunity to measure crustal anisotropy, which may be important in estimating the geometrical and physical parameters of reservoir rocks. VSPs are particularly important for this purpose as they are less subject to the complex interactions of the shear wavefield with the free surface. The first stage in characterizing the subsurface anisotropy requires that the distinctive phenomenon of shear-wave splitting must be examined for every arrival at each geo-phone. This effect may be defined by two parameters: the polarization of the leading shear-wave and the time-delay between corresponding split shear-waves. A variety of techniques have been designed to estimate these parameters of shear-wave splitting. Here, we classify the published techniques into four main categories and review their properties. Representative procedures from each group are applied to a common synthetic data set contaminated with signal-generated noise. The results allow some general statements to be made about the utility of these methods for processing shear-waves in VSP data.     

综合应用钻孔应变固体潮和应变地震波资料求解地震参数和地层的力学参数

介绍了在井下耦合介质弹性参数等不能直接测定时，应用钻孔应变固体潮和就地震波资料求解地震参数和地层的力学参数的方法，用平衡方程的解来近似代替波动方程解，不但使用钻孔应变仪测到的地震波测地震成为可能，而且比直接应用波动方程求解的方法大大简化。     

APPLICATION AND STATISTICAL ANALYSIS OF RELATIONSHIP BETWEEN SHEAR WAVE VELOCITY AND DEPTH OF SOIL-LAYERS

The shear wave velocity is one of the important parameters in seismic engineering.The common mathematical models of relationship between shear wave velocity and depth of soil-layers are linear function model,quadratic function model,power function model,cubic function model,and quartic function model.It is generally believed that the regression formulae based on aforementioned mathematical models are mainly used for preliminary estimation of the local shear wave velocity.In order to increase the value of test data of wave speed in boreholes,the calculation formulae for the thickness of ground cover layer are derived based on the aforementioned mathematical models and their fitting parameters.The calculation formulae for the mean shear wave velocity of soil-layers are derived by integral mean value theorem.Accordingly,the calculation formulae for the equivalent shear wave velocity of soil-layers are derived.The calculation formulae for the depth of reflective waves in time-depth conversion of the reflection seismic exploration are derived.Through the statistical analysis of test data of shear wave velocity of soil layers in Changyuan County,Henan Province,regression formulae and their fitting parameters of aforementioned mathematical models are obtained.The results show that in the determination of the quality of these regression formulae and their fitting parameters,the adjusted R-square,root mean square error and residual error,the matching on the statistical range between the geometry of function of mathematical models used and the scattergram of the measured data,the application purpose and the simplicity of the regression formulae should be considered.With the aforementioned new formulae,the results show that the calculated values of equivalent shear wave velocity of soil-layers and thickness of ground cover layer meet the engineering needs.The steps for statistics and applications of the relationship between shear wave velocity and depth of soil-layers for a new area are as follows:(1) Analyze the relevant data about the site such as the drilling and wave speed test data,etc.and divide the site into seismic engineering geological units;(2) In a single seismic engineering geological unit,make statistical analysis of the data of borehole wave speed test,comprehensively identify and select mathematical models and their fitting parameters of the relationship between shear wave velocity and depth of soil-layers;(3) Substitute the selected fitting parameters into the formulae,based on their mathematical models for the thickness of ground cover layer,or the equivalent shear wave velocity of soil-layers,or the depth of reflective wave,then the thickness of ground covering layer,equivalent shear wave velocity,and depth of reflective wave are obtained.     

ADAPTIVE PREDICTIVE DECONVOLUTION OF SEISMIC DATA*

The Wiener prediction filter has been an effective tool for accomplishing dereverberation when the input data are stationary. For non-stationary data, however, the performance of the Wiener filter is often unsatisfactory. This is not surprising since it is derived under the stationarity assumption. Dereverberation of nonstationary seismic data is here accomplished with a difference equation model having time-varying coefficients. These time-varying coefficients are in turn expanded in terms of orthogonal functions. The kernels of these orthogonal functions are then determined according to the adaptive algorithm of Nagumo and Noda. It is demonstrated that the present adaptive predictive deconvolution method, which combines the time-varying difference equation model with the adaptive method of Nagumo and Noda, is a powerful tool for removing both the long- and short-period reverberations. Several examples using both synthetic and field data illustrate the application of adaptive predictive deconvolution. The results of applying the Wiener prediction filter and the adaptive predictive deconvolution on nonstationary data indicate that the adaptive method is much more effective in removing multiples. Furthermore, the criteria for selecting various input parameters are discussed. It has been found that the output trace from the adaptive predictive deconvolution is rather sensitive to some input parameters, and that the prediction distance is by far the most influential parameter.     

MAXIMUM-LIKELIHOOD ESTIMATION OF SEISMIC IMPULSE RESPONSES*

A seismic trace is assumed to consist of a known signal pulse convolved with a reflection coefficient series plus a moving average noise process (colored noise). Multiple reflections and reverberations are assumed to be removed from the trace by conventional means. The method of maximum likelihood (ML) is used to estimate the reflection coefficients and the unknown noise parameters. If the reflection coefficients are known from well logs, the seismic pulse and the noise parameters can be estimated. The maximum likelihood estimation problem is reduced to a nonlinear least-squares problem. When the further assumption is made that the noise is white, the method of maximum likelihood is equivalent to the method of least squares (LS). In that case the sampling rate should be chosen approximately equal to the Nyquist rate of the trace. Statistical and numerical properties of the ML- and the LS-estimates are discussed briefly. Synthetic data examples demonstrate that the ML-method gives better resolution and improved numerical stability compared to the LS-method. A real data example shows the ML- and LS-method applied to stacked seismic data. The results are compared with reflection coefficients obtained from well log data.     

陇西地区黄土的强度特性

本文介绍了对陇西地区黄土的强度特性进行的试验研究结果。陇西地区Q_3和Q_4黄土在偏应力作用下的应力—应变关系基本上可用康德纳双曲线模式来表示。对大量的试验数据采用了最小二乘法计算了黄土的强度参数,结果表明其强度参数是比较低的。本文还对几十组黄土的强度参数与其物理性质指标之间的关系进行了回归分析,建立了它们之间的经验公式,该公式可以直接用于陇西地区黄土的强度和土坡稳定性估算。     

APPLICATION OF TWO SIMPLIFIED NEWMARK MODELS TO THE ASSESSMENT OF LANDSLIDES TRIGGERED BY THE 2008 WENCHUAN EARTHQUAKE

At present, with the wide application of the Newmark method, various Newmark empirical formulas with different ground motion parameters have been fitted by many researchers based on global strong-motion records. However, the existing study about the Wenchuan earthquake does not quantitatively evaluate the applicability of different Newmark models based on the actual landslides distribution. The aim of this paper is to present a comparison between observed landslides from the 2008 Wenchuan earthquake and predicted landslides using Newmark displacement method based on different ground motion parameters. The factor-of-safety map and critical acceleration(a_c)map in the study area are obtained by using the terrain data and geological data. The distribution of Arias intensity(I_a)and PGA in the study area is obtained by using the attenuation formulas of Arias intensity(I_a)and PGA, which is regressed by Wenchuan ground motion records. Based on the distribution of Arias intensity(I_a)and PGA parameters, we obtained the predicted locations of landslide using Newmark regression equations which are generated using global strong-motion records. The results shows that the assessment results can better reflect the macroscopic distribution characteristics of co-seismic landslides, most predicted landslide cells are distributed on the two sides of the Beichuan-Yingxiu Fault, especially the Pengguan complex rock mass in the hanging wall. The abilities to predict landslide occurrence of the two Newmark simplified models are different. On the whole, the evaluated result of simplified model based on parameter I_a is better than that based on PGA parameter. The GFC values obtained by the Newmark model of I_a and PGA parameters are 65.7% and 34.9%respectively. The evaluated result based on I_a can better reflect the macro distribution of coseismic landslides. The Ls_Pred value based on the Newmark model of parameter I_a is 26.5%, and the Ls_Pred value based on the Newmark model of PGA parameter is 10.3%. However the total area of predicted landslides accounts for 2.4% of the study area, which indicates that the predicted landslide cells are greater than the observed landslide cells. This reminds us that depending on the current input of shear strength and ground-motion parameters, we can only conduct landslide hazard assessment in macro areas, the ability to predict landslide can be improved using more accurate topographic data and input parameters.     

AN APPROACH TO THE MEASUREMENT OF THE POTENTIAL STRUCTURAL DAMAGE OF EARTHQUAKE GROUND MOTIONS

The quantification and prediction of damage due to different seismic actions to structure types of different strength is an important problem not yet solved in the Earthquake Engineering field. In addition, owing to the fact that macroseismic information cannot be used directly in dynamic calculations, a new problem appears when these are the only kind of data available. Thus, there is a need to estimate a parameter to relate the energy of the ground motion and the damage occurrence, and eventually achieve a better seismic risk assessment. After the study and review of some representative potential damage parameters, attention has been paid to the Arias intensity (unfiltered and filtered in certain frequency ranges) and the Cumulative Absolute velocity (CAV) as the parameters to evaluate the energy of movement, and to relate them with the observed damage. The data used to infer these correlations have been provided by the ENEA-ENEL (Italy). The information consists of strong motion records from the Campano Lucano (1980), Umbria (1984) and Lazio-Abruzzo (1984) earthquakes, and data of damage to buildings in the vicinity of recording instruments (within a maximum radius of 300 m, where the soil conditions remain constant). In this paper, some relations have been obtained to quantify the damage level for different seismic inputs. The results suggest that unfiltered Arias intensity and CAV (for calculation threshold 20 cm/s²) correlate well with the macroseismic information used. Best fits are obtained between the quoted parameters and the observed damage in type A structures. © 1997 by John Wiley & Sons, Ltd.     

APPLICATION OF OFF-LINE AND ON-LINE IDENTIFICATION TECHNIQUES TO BUILDING SEISMIC RESPONSE DATA

The objectives of this paper are to present a comparison of the dynamic characteristics of a seven-storey reinforced concrete building (Van Nuys–Holiday Inn) identified from four recorded strong-motion response data (Whittier earthquake, Landers earthquake, Big Bear earthquake and Northridge earthquake). In the analysis, time-domain methods for estimating the system parameters and the modal properties of the building are studied. Both off-line and on-line identification algorithms are applied to these seismic response data. Under the assumption of a linear time-invariant system the ARX model and ARMAX model are used. Comparison of the identification results using different models are made. In addition, recursive procedures are adapted as on-line identification and the time-varying modal parameters are estimated. For structural systems under strong earthquake excitation, a recursive identification method, adaptive forgetting through multiple models (AFMM), is introduced to identify systems with rapidly changing parameters. Through the analysis of the seismic response data of the building subjected to four earthquakes the identification algorithm and the identification results are discussed.     

用地震面波研究上地幔各向异性

本文论述了上地幔各向异性的参数描述和主要证据,以及与地球动力学的关系。并就近几年来国外利用地震面波研究上地幔各向异性的重要成果和主要方法作了综合评述。     

从雷电脉冲谱推算出低电离层等效结构

从低电离层-地面波导中电波传播的波导模理论出发,分析了雷电产生的电磁脉冲（LEP）谱随距离变化与低电离层结构之间关系。利用LEP谱中干扰最小值时的频率与指数电离层模式参量的关系得出低电离层等效模拟参量。从我们观测出的LEP数据得出指数电离层等效参量是:在白天β=0.3km^-1,h_i=70km;在夜间β=0.5km^-1,h_i=88km。并指出使用这种方法可以观测出低电离层在时间上和地理位置上的变化。     

NON-LINEAR AVO INVERSION FOR A STACK OF ANELASTIC LAYERS1

Parameters in a stack of homogeneous anelastic layers are estimated from seismic data, using the amplitude versus offset (AVO) variations and the travel-times. The unknown parameters in each layer are the layer thickness, the P-wave velocity, the S-wave velocity, the density and the quality factor. Dynamic ray tracing is used to solve the forward problem. Multiple reflections are included, but wave-mode conversions are not considered. The S-wave velocities are estimated from the PP reflection and transmission coefficients. The inverse problem is solved using a stabilized least-squares procedure. The Gauss-Newton approximation to the Hessian matrix is used, and the derivatives of the dynamic ray-tracing equation are calculated analytically for each iteration. A conventional velocity analysis, the common mid-point (CMP) stack and a set of CMP gathers are used to identify the number of layers and to establish initial estimates for the P-wave velocities and the layer thicknesses. The inversion is carried out globally for all parameters simultaneously or by a stepwise approach where a smaller number of parameters is considered in each step. We discuss several practical problems related to inversion of real data. The performance of the algorithm is tested on one synthetic and two real data sets. For the real data inversion, we explained up to 90% of the energy in the data. However, the reliability of the parameter estimates must at this stage be considered as uncertain.     

小震震源参数与震源机制的联合测定

本文采用一个非均匀破裂圆盘位错模式,考虑到介质-仪器系统的影响,计算了它在不同幅射方向的P波初动位移与半周期,以及几组用于测定震源参数的理论曲线。在采用1/4周期测定介质的Q值和球面拟合震源机制解的基础上,给出了一个小震震源参数与震源机制的联合测定方法。该方法的实质是,将震源断层面作为被测参数的一部分,并同时考虑它对其它被测参数的影响,从而获得单个小震震源参数与震源机制的相应解。最后,根据新丰江水库区四个台的小震资料作了实测与分析,所得结果为:介质的平均Q值约450,地震矩1019——1020达因厘米,平均应力降约35巴;其震源机制解亦反映了本区构造应力场的主要特征。      

地壳主要岩石流变参数及华北地壳流变性质研究

岩石流变参数和变形机制是根据断层摩擦和岩石幂次流动本构关系建立岩石圈强度剖面的基础。近 30年来 ,高温高压实验取得了很大进展 ,获得了大量地壳矿物和岩石流变资料。本文系统总结了这些流变实验资料 ,并应用流变数据结合地震震源深度分布 ,对华北地壳流变性质进行了研究。结果表明 ,以花岗岩和低级变质岩为代表的上地壳为脆性破裂 ,其强度受断层摩擦约束 ,以长英质片麻岩为主的中地壳和以中性麻粒岩为主的下地壳上层处于塑性流变状态 ,由干的基性麻粒岩组成的下地壳下层处于脆性向塑性流变的过渡状态。华北地壳的这种物质组成和流变为地壳不同层次的解耦和强震孕育提供了力学条件 ,也构成了不同尺度块体的底边界     

2004年首都圈地区中小地震的矩张量反演

通过首都圈地区宽频带数字地震波形反演, 获得了2004年发生在该地区的51次中小地震的矩张量解, 由此确定了这些地震的标量地震矩、 矩震级、 双力偶分量和补偿线性矢量偶极分量的大小以及断层面参数和应力轴参数, 并通过数值试验对反演结果进行了评价.      

安庆4.8级地震序列的地震波谱变化研究

采用安徽区域数字地震台网记录的安庆4.8级地震序列波形资料,开展地震波谱研究,求取了零频位移谱值、应力降和地震矩等震源参数。分析了纵、横波零频谱值比和应力降随时间的变化趋势,并讨论了震源参数之间的关系。结果表明,安庆4.8级主震后,纵、横波零频位移谱值比和应力降呈现一定的变化规律,可为主震—余震序列的最大余震的预测提供一定的参考依据,同时发现震源参数与震级有一定相关。