Joint multivariate statistical model and its applications to synthetic earthquake predic-tion

IntroductionSeismologistsinChinahaveestablishedmanyseismologicalmethodsandearthquakepredic-tionfactorsintheearthquakeforecastpracticesinthelastthirtyyears.Atpresent,morethanonehundredseismometricfactorsareusedinthemedium-termandmedium-short-termearthquakepredictions,suchastheb-valuetomirrorthestressstateandthedistributiondegreeofthemedium,themf-valuetoexpressiftheseismicactivityisincreasedornot(WANG,etal,1994),thequanti-tativeparameterA(b)todescribetheseismicactivitiesinadistrict(WU,CAO,19…     

多元统计组合模型在地震综合预报中的应用

针对目前地震综合预报中的一些问题，利用近30年来迅速发展的多元统计分析中主成分分析、判别分析组成多元统计组合模型，在众多的地震预报指标（预报因子）中采用信息最大化方法，选择对中期预测信息累积贡献率大于90%地震预报指标，分别进行相关分析、预测、检验，最终应用马氏距离判别作外推综合预报；并以华北地区（30~42N，108~125E）为例进行模型的应用检验，初步研究已取得了较好的效果.      

Application of support vector machine to synthetic earthquake prediction

This paper introduces the method of support vector machine (SVM) into the field of synthetic earthquake pre-diction, which is a non-linear and complex seismogenic system. As an example, we apply this method to predict the largest annual magnitude for the North China area (30°E-42°E, 108°N-125°N) and the capital region (38°E-41.5°E, 114°N-120°N) on the basis of seismicity parameters and observed precursory data. The corresponding prediction rates for the North China area and the capital region are 64.1% and ...     

青海地区强震活动的统计分析及预报探讨

从统计的角度出发，研究了青海地区１９９３年Ｍ≥６强震的发生概率，在此基础上应用非线性最小二乘法拟合的Ｇｏｍｐｅｒｔｚ与Ｌｏｇｉｓｔｉｃ曲线模型以及可公度模型对该区Ｍ≥６．５强震进行了预测探索。结果表明，１９９３年该地区Ｍ≥６强震的发生概率将达到０．７６－０．８７，１９９３－１９９４年有可能发生Ｍ≥６．５强震。     

川滇地区地震动预测模型及其对2021年漾濞6.4级地震的适用性

随着川滇地区强震记录的不断增加,为了建立更符合该区域地震动特征的预测模型,文中基于该区域现有的地震动数据,通过随机效应回归模型建立适用于川滇地区的地震动预测模型;2021年5月21日,云南省大理州漾濞县发生6.4级地震,为了分析文中预测模型对漾濞地震的适用性,首先根据预测模型的适用范围选取合适的漾濞地震数据,计算真实记...     

A multivariate non‐parametric model for synthetic generation of daily streamflow

A p‐order multivariate kernel density model based on kernel density theory has been developed for synthetic generation of multivariate variables. It belongs to a kind of data‐driven approach and is able to avoid prior assumptions as to the form of probability distribution (normal or Pearson III) and the form of dependence (linear or non‐linear). The p‐order multivariate kernel density model is a non‐parametric method for synthesis of streamflow. The model is more flexible than conventional parametric models used in stochastic hydrology. The effectiveness and satisfactoriness of this model are illustrated through its application to the simultaneous synthetic generation of daily streamflow from Pingshan station and Yibin‐Pingshan region (Yi‐Ping region) of the Jinsha River in China. Copyright © 2007 John Wiley & Sons, Ltd.     

基于PCA-PSO-SVM的地震死亡人数预测模型研究

为准确预测地震死亡人数,提出了基于主成分分析法(PCA)和粒子群算法(PSO)优化的支持向量机(SVM)模型。首先利用主成分分析法对地震死亡人数7个影响因子中的6个进行数据降维,同时对第7个发震时刻因子单独进行区间分类,然后对提取出的主成分进行归一化处理,将归一化的主成分数据作为支持向量机的输入向量,通过粒子群算法寻优获得最优支持向量机模型参数,最终建立基于PCA-PSO-SVM的地震死亡人数预测模型,并对5组样本进行死亡人数预测,同时对比分析包含和不包含发震时刻因子的2种情况下的模型预测效果。结果表明:在不考虑发震时刻因子的情况下,使用PCA-PSO-SVM模型的最小误差、最大误差和平均误差分别为0.85%、20%、10%,其平均误差相比PSO-SVM、SVM模型分别降低2.08%、2.28%;输入向量加入发震时刻因子分类数据后,PCA-PSO-SVM模型的最小误差、最大误差和平均误差分别为0.25%、20%、7.18%,其平均误差相比PSO-SVM、SVM模型分别降低3.34%、3.50%。因此,加入发震时刻因子后3种模型的平均误差明显降低,同时由于PCA-PSO-SVM模型进行主成分降维处理,能够明显提高运行效率和预测精度,故降低了模型复杂度。     

The method of AF synthetic detection on crust deformation anomaly and its application in earthquake prediction

In this paper, the unstable model of crust deformation and the anomaly F-detection are discussed. On this basis, the quantitative parameters of deformation anomaly and the detection criteria for identifying anomaly have also been given, and we detected the anomaly of simulated numerical sequences. Using the observed multi-dimension data, we computed the anomaly intensity and analyzed the possible relationship between the AF anomaly and the accumulated energy released by earthquakes. Results demonstrated that the AF index can be used to illustrate an anomaly intensity for some complicated processes and the anomaly detection has reliability and stability in some degree. Sponsored by The Chinese Joint Seismological Science Foundation (197066) and Project Nos. 95-04-05-03-01, 96-913-01-01-05.     

地震统计力学—几个非线性地震模型的相互关系及其物理含义

近年来地震统计力学研究的进展已经使我们有可能以一种简洁的方式将一些非线性地震模型物理地联系起来,并从地震学角度阐明这些模型的意义和限度。本文结合逾渗模型、热力学弹性回跳模型及弹簧滑块模型,系统地评述了这方面研究的进展。     

Application of mechanical and statistical models to the study of seismicity of synthetic earthquakes and the prediction of natural ones

Based on a cellular automata model, which obeys rules produced by a spring-damper-slider model of fault network, a synthetic earthquake catalog of deterministic chaos can be produced. It is a challenge to inverse the model parameters and predict the synthetic earthquakes from the synthetic catalog. We apply an improved statistical model of coupled stress release to the synthetic catalog, it does successfully inverse the stress evolution of the system, but effectively prediction of the earthquakes from the stresses remains a problem. This research is partially supported by the Chinese Joint Seismological Science Foundation and Asian 2000 Foundation of New Zealand.     

电阻率各向异性及在地震研究中的应用

地电阻率各向异性广泛存在地壳中,尤其是浅部.因地震、矿山开采、火山喷发以及大型地质构造运动会明显改变地应力场,从而导致裂隙优势排列及裂隙中充填的流体重新分布,最终导致明显的电阻率各向异性.通过监测其变化,不仅可研究和了解地应力场的动态演化过程,还可以作为预测地震,矿震、火山活动和灾害性地质构造运动的依据,并研究它们形成的物理过程和物理机制.本文从室内、野外实验和数值模拟方面对电阻率各向异性的研究进行了概述,阐述了其在地震监测预报研究中的应用.探讨了室内电阻率各向异性岩石实验研究的发展特点,提出了新手段--电成像技术实验研究岩石电阻率各向异性的思路,期望建立研究电阻率各向异性的阵列式观测系统.     

Seismicity acceleration model and its application to several earthquake regions in China

With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates: North China Subplate, Chuan-Dian Block and Xinjiang Subplate, and divide the three subplates into seven researched regions by the difference of seismicity and tectonic conditions. With the modified equation given by Sornette and Sammis (1995), we analysis the seismicity of each region. To those strong earthquakes already occurred in these region, the model can give close fitting of magnitude and occurrence time, and the result in this article indicates that the seismicity acceleration model can also be used for describing the seismicity of intraplate. In the article, we give the magnitude and occurrence time of possible strong earthquakes in Shanxi, Ordos, Bole-Tuokexun, Ayinke-Wuqia earthquake regions. In the same subplate or block, the earthquake periods for each earthquake region are similar in time interval. The constant αin model can be used to describe the intensity of regional seismicity, and for the Chinese Mainland, α is 0.4 generally. To the seismicity in Taiwan and other regions with complex tectonic conditions, the model does not fit well at present.     

Strong ground motion record selection for the reliable prediction of the mean seismic collapse capacity of a structure group

How to select a limited number of strong ground motion records (SGMRs) is an important challenge for the seismic collapse capacity assessment of structures. The collapse capacity is considered as the ground motion intensity measure corresponding to the drift‐related dynamic instability in the structural system. The goal of this paper is to select, from a general set of SGMRs, a small number of subsets such that each can be used for the reliable prediction of the mean collapse capacity of a particular group of structures, i.e. of single degree‐of‐freedom systems with a typical behaviour range. In order to achieve this goal, multivariate statistical analysis is first applied, to determine what degree of similarity exists between each selected small subset and the general set of SGMRs. Principal Component analysis is applied to identify the best way to group structures, resulting in a minimum number of SGMRs in a proposed subset. The structures were classified into six groups, and for each group a subset of eight SGMRs has been proposed. The methodology has been validated by analysing a first‐mode‐dominated three‐storey‐reinforced concrete structure by means of the proposed subsets, as well as the general set of SGMRs. The results of this analysis show that the mean seismic collapse capacity can be predicted by the proposed subsets with less dispersion than by the recently developed improved approach, which is based on scaling the response spectra of the records to match the conditional mean spectrum. Copyright © 2010 John Wiley & Sons, Ltd.     

Seismicity anomalies before the great earthquake of M_S=8.1 in the Kunlun Pass and its significance to earthquake prediction

Introduction The MS=8.1 earthquake occurred in west of the Kunlun Pass on November 14, 2001. It is the greatest earthquake occurred in China since the last half of the century and is an important event in recent seismic history of China. Some specialists consider that the earthquake occurred in the area where the earthquake monitoring capability is lowest in Chinese mainland; no striking precursory seismicity was found. The study on the precursory seismicity before the earthquake has not b…     

随机地震动模拟的时间序列法及其工程应用

系统研究了AR和ARMA时间序列方法生成给定加速度功率谱人工地震波的新方法，根据相应理论编制了程序，并生成了人工地震加速时程曲线，然后用生成的人工地震波进行结构动力时程分析，与实际地震波分析的结果进行了对比，计算结果表明，本文用于合成人工地震波的方法是可行的，生成的人工波可以用于实际结构的地震反应分析。     

诱发地震主震发生时间的非线性灰色预测及其软件系统

在深入分析非线性灰色模型基本原理的基础上,开发了诱发地震主震发生时间的预测软件,并通过大量实例验证了本模型用于预测主震发生时间的适用性以及本软件的可靠性。     

Simulation of synthetic ground motions for specified earthquake and site characteristics

A method for generating a suite of synthetic ground motion time‐histories for specified earthquake and site characteristics defining a design scenario is presented. The method employs a parameterized stochastic model that is based on a modulated, filtered white‐noise process. The model parameters characterize the evolving intensity, predominant frequency, and bandwidth of the acceleration time‐history, and can be identified by matching the statistics of the model to the statistics of a target‐recorded accelerogram. Sample ‘observations’ of the parameters are obtained by fitting the model to a subset of the NGA database for far‐field strong ground motion records on firm ground. Using this sample, predictive equations are developed for the model parameters in terms of the faulting mechanism, earthquake magnitude, source‐to‐site distance, and the site shear‐wave velocity. For any specified set of these earthquake and site characteristics, sets of the model parameters are generated, which are in turn used in the stochastic model to generate the ensemble of synthetic ground motions. The resulting synthetic acceleration as well as corresponding velocity and displacement time‐histories capture the main features of real earthquake ground motions, including the intensity, duration, spectral content, and peak values. Furthermore, the statistics of their resulting elastic response spectra closely agree with both the median and the variability of response spectra of recorded ground motions, as reflected in the existing prediction equations based on the NGA database. The proposed method can be used in seismic design and analysis in conjunction with or instead of recorded ground motions. Copyright © 2010 John Wiley & Sons, Ltd.     

Using multivariate statistical techniques and geochemical modelling to identify factors controlling the evolution of groundwater chemistry in a typical transitional area between Taihang Mountains and North China Plain

Identifying the key factors controlling groundwater chemical evolution in mountain-plain transitional areas is crucial for the security of groundwater resources in both headwater basins and downstream plains. In this study, multivariate statistical techniques and geochemical modelling were used to analyse the groundwater chemical data from a typical headwater basin of the North China Plain. Groundwater samples were divided into three groups, which evolved from Group A with low mineralized Ca-HCO₃ water, through Group B with moderate mineralized Ca-SO₄-HCO₃ water, to Group C with highly saline Ca-SO₄ and Ca-Cl water. Water-rock interaction and nitrate contamination were mainly responsible for the variation in groundwater chemistry. Groundwater chemical compositions in Group A were mainly influenced by dissolution of carbonates and cation exchange, and suffered less nitrate contamination, closely relating to their locations in woodland and grassland with less pronounced human interference. Chemical evolution of groundwater in Groups B and C was gradually predominated by the dissolution of evaporites, reverse ion exchange, and anthropogenic factors. Additionally, the results of the inverse geochemical model showed that dedolomitization caused by gypsum dissolution, played a key role in the geochemical evolution from Group A to Group B. Heavy nitrate enrichment in most groundwater samples of Groups B and C was closely associated with the land-use patterns of farmland and residential areas. Apart from the high loads of chemical fertilizers in irrigation return flow as the main source for nitrate contamination, the stagnant zones, flood irrigation pattern, mine drainage, and groundwater-exploitation reduction program were also important contributors for such high mineralization and heavy NO₃⁻ contents in Group C. The important findings of this work not only provide the conceptual framework for the headwater basin but also have important implications for sustainable management of groundwater resources in other headwater basins of the North China Plain.     

平稳随机地震地面运动过程模型及其统计特征

地震地面运动过程具有强烈的随机性,应用随机理论对实际工程结构进行地震可靠性分析和抗震设计与加固时都需要建立合理的随机地震地面运动模型,本文选择3种典型的随机地震动模型,即理想白噪声模型、金井清模型和改进的金井消模型,分析了它们的物理概念、频域特征以及适用范围。引入状态向量,建立状态方程．通过复振型叠加法分析了地震地面运动过程的时域统计特性,推导出3种随机地震动模型的相关函数的解析表达式．这些结果可为结构随机地震反应时域分析和抗震可靠性评估提供基础。     

Application of fuzzy clustering method to determining sub-fault planes of earthquake from aftershocks sequence

Earthquake rupture process generally involves several faults activities instead of a single fault. A new method using both fuzzy clustering and principal component analysis makes it possible to reconstruct three dimensional structure of involved faults in earthquake if the aftershocks around the active fault planes distribute uniformly. When seismic events are given, the optimal faults structures can be determined by our new method. Each of sub-fault planes is fully characterized by its central location, length, width, strike and dip. The resolution determines the number of fault segments needed to describe the earthquake catalog. The higher the resolution, the finer the structure of the reconstructed fault segments. The new method successfully reconstructs the fault segments using synthetic earthquake catalogs. By taking the 28 June 1992 Landers earthquake occured in southern California as an example, the reconstructed fault segments are consistent with the faults already known on geological maps or blind faults that appeared quite frequently in longer-term catalogs.