基于全波形激光雷达数据分解方法的研究

波形分解是全波形激光雷达数据处理的关键技术之一,通过雅克比矩阵更新与否以及阻尼系数的调整,对优化LM(Levenberg-Marquard)算法进行了研究。实验表明,该方法可最优化求解足印内各波形分量的振幅、中心位置和半宽,而不必每次迭代时都需计算雅克比矩阵,简化了计算过程,提高了计算效率,最终得到的模拟波形能够非常好地拟合原始波形。     

Bivariate Gaussian mixture–based equivalent linearization method for stochastic seismic analysis of nonlinear structures

To address challenges in stochastic seismic analysis of nonlinear structures, this paper further develops a recently proposed Gaussian mixture–based equivalent linearization method (GM‐ELM). The GM‐ELM uses a Gaussian mixture distribution model to approximate the probabilistic distribution of a nonlinear system response. Using properties of the Gaussian mixture model, GM‐ELM can decompose the non‐Gaussian response of a nonlinear system into multiple Gaussian responses of linear single–degree of freedom oscillators. With the set of the equivalent linear systems identified by GM‐ELM, response statistics as crossing rate and first‐passage probability can be computed conveniently using theories of linear random vibration analysis. However, the original version of GM‐ELM may lead to an inaccurate estimate because of the heuristic parameters of the linear system introduced to supplement insufficient information. To overcome this limitation and define unique equivalent linear systems, this paper proposes a further developed version of GM‐ELM, which uses a mixture of bivariate Gaussian densities instead of univariate models. Moreover, to facilitate the use of elastic response spectra for estimating the mean peak responses of a nonlinear structure, a new response spectrum combination rule is proposed for GM‐ELM. Two numerical examples of hysteretic structural systems are presented in this paper to illustrate the application of the bivariate GM‐ELM to nonlinear stochastic seismic analysis. The analysis results obtained by the bivariate GM‐ELM are compared with those obtained by the univariate GM‐ELM, the conventional equivalent linearization method, the tail equivalent linearization method, and Monte Carlo simulation. The supporting source code and data are available for download at https://github.com/yisangri/GitHub‐bGM‐ELM‐code.git     

基于误差分布的震源区波速比反演及其应用:乳山震群源区介质性质变化研究

本文利用误差分布和概率统计分析改进了波速比计算方法.对于震源位置相对集中的震群活动,对台站震相到时进行两次差分,通过对差分后的震相数据对进行二维高斯分布拟合,可以稳健地估计震群活动震源区波速比.该方法充分利用了不同台站Pg、Sg到时差的所有信息,其优势是不需要地震事件的震源位置,并且不依赖震源区以外的速度变化,有效消除了震源区到台站的传播路径效应的影响;相对于传统的平均波速比,本文方法得到的震源区波速比,更能真实地反映震源区介质的性质.我们将该方法应用到2013—2016年的乳山震群,结果显示:震源区波速比的变化与震群活动过程密切相关,波速比的变化反映了序列活动的阶段性特征.     

Modeling of an actively braced full‐scale building considering control–structure interaction

Recently, the application of active control to seismic‐excited buildings has attracted international attention. To demonstrate the practical applicability of active control, we have conducted experimental tests using a full‐scale three‐storey building equipped with active bracing systems on the shake table at the National Center for Research on Earthquake Engineering (NCREE), Taiwan. Experimental results indicate that the control–structure interaction (CSI) effect is significant. A state‐space analytical model of this actively controlled building taking into account the CSI effect is established in this paper using a system identification technique based on curve‐fitting of transfer functions. To verify the accuracy of the analytical model for simulating the controlled response, four sets of linear quadratic Gaussian (LQG) controllers using acceleration feedback are designed and further experimental tests are conducted for comparison. It is demonstrated that the correlations between the simulation and experimental results are remarkable. The construction of an accurate analytical model is important for active control, and such an analytical model can be used for future benchmark studies of different control algorithms based on numerical simulations. Copyright © 2000 John Wiley & Sons, Ltd.     

基于高斯距离加权格网法研究华北地区地壳形变特征

分析了高斯距离加权格网法计算应变率的原理和特征,研制了相应的计算程序,并用于华北地区现今地壳应变特征研究。基于华北地区1999～2009年的GPS观测结果,研究了确定高斯平滑半径的原则和方法,并给出华北地区地壳形变多个应变率物理量的分布。结果表明,基于150 km高斯平滑半径计算水平速度梯度能够较好地分辨华北地区的应变率空间变化特征,首都圈地区地壳以张渤带上的左旋剪切为主,其中唐山-秦皇岛地区同时存在南北向的拉张变形;山西带整体表现为沿断裂的右旋剪切变形,同时兼具“南挤北张”特性;郯庐断裂带整体形变特征不显著,在南北两端局部区域存在较小的剪切形变;鄂尔多斯块体内部比较稳定,北缘以拉张和左旋剪切变形为主,西缘表现为左旋剪切以及EW向的挤压特征,南缘的秦岭块体南北两侧分别具有左旋剪切和右旋剪切的特性,东西两侧分别表现为正断拉张与逆断挤压特征。     

基于有效邻域波场近似的起伏地表保幅高斯束偏移

随着我国陆上地震勘探向复杂地表探区的转移,高精度、适应性强的地震成像方法在地震资料的处理、解释及后续属性分析、储层预测中具有重要意义.本文基于有效邻域波场近似理论发展了一种成像精度更高且适用于复杂起伏地表条件的叠前保幅高斯束偏移方法.在传统水平地表高斯束偏移的基础上,本文根据中心射线附近有效邻域内高斯束表征的近似波场,导出了起伏地表条件下具有相对振幅保持的高斯束偏移公式,并给出了一种精度更高的旁轴射线传播角度计算方法.同现有的高斯束偏移方法相比,本文方法不仅考虑了起伏地表对高斯束走时的线性影响,而且首次引入了由地表高程差异和近地表速度变化引起的二次时差校正项和振幅校正项,使得成像结果更加准确可靠.两个典型模型算例验证了本文方法的正确性和有效性.     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

利用GRACE卫星数据提取苏门答腊地震同震和震后重力变化

基于GRACE重力卫星数据,采用P3M6去相关滤波技术、350 km半径的高斯平滑技术与差分方法,提取2004年苏门答腊Mw9.3地震的同震和震后重力场时空演化过程。该地震断层下盘同震重力变化达到4.6 μGal,上盘同震重力变化最大值为-8.3 μGal,该结果与球体位错理论计算结果在量级和分布形态上具有较好的一致性。震后7 a断层下盘重力变化累计达到5.6 μGal,超过同震变化信号|而断层上盘震后重力变化只有0.6 μGal,仅为同震变化信号的11%。该现象表明,断层下盘所在海洋地层的粘滞性较低,地层相对柔软,其地层在震后由于较大的粘滞性松弛效应而持续变形,并伴随较大的重力变化|而断层上盘所在的地层接近陆地地壳特性,地层相对脆硬,震后受地层的粘滞性响应而导致的重力变化相对较小。     

Enrichment of seismic ground motion databases using Karhunen–Loève expansion

Best estimate seismic analysis are generally based on time‐domain simulations of structural responses. The seismic load is then modeled by a stochastic process representing ground motion. For this purpose, the analyst can use recorded accelerograms or work with synthetically generated ones. The number of ground motion time‐histories available for a given scenario and site condition is limited and generally not sufficient for carrying out more advanced probabilistic structural response analysis. It is then necessary to have at our disposal methods that allow for generating synthetic accelerograms that realistically characterize earthquake ground motions. However, most of the methods proposed in literature for generating synthetic accelerograms do not accurately reproduce the natural variability of ground motion parameters (such as PGA, cumulative absolute velocity, and Arias intensity) observed for recorded time histories. In this paper, we introduce a new method for generating synthetic ground motion, based on Karhunen‐Loève decomposition and a non‐Gaussian stochastic model. The proposed method enables the structural analyst to simulate ground motion time histories featuring the properties mentioned above. To demonstrate its capability, we study the influence of the simulation method on different ground motion parameters and on soil response spectra. We finally compute fragility curves to illustrate the practical application of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Improved estimators of correlation and R2 for skewed hydrologic data

ABSTRACT

The coefficient of determination R² and Pearson correlation coefficient ρ = R are standard metrics in hydrology for the evaluation of the goodness of fit between model simulations and observations, and as measures of the degree of dependence of one variable upon another. We show that the standard product moment estimator of ρ, termed r, while well-behaved for bivariate normal data, is upward biased and highly variable for bivariate non-normal data. We introduce three alternative estimators of ρ which are nearly unbiased and exhibit much less variability than r for non-normal data. We also document remarkable upward bias and tremendous increases in variability associated with r using both synthetic data and daily streamflow simulations from 905 calibrated rainfall–runoff models. We show that estimators of ρ = R accounting for skewness are needed for daily streamflow series because they exhibit high variability and skewness compared to, for example, monthly/annual series, where r should perform well.