考虑震害因子耦合作用的砌体建筑群破坏概率模型研究

砌体建筑群在地震中往往破坏严重损失巨大,合理评估地震作用对不同种类砌体结构造成破坏的风险变得至关重要。传统基于后验概率的地震危险性分析方法忽略了砌体建筑个体差异性的影响,未深入考虑多种震害因子的耦合作用。本文以华南地区砌体建筑群为例,开发了一种集成概率方法来对城市砌体结构的破坏风险进行建模,考虑建筑年代、层数、使用用途和墙厚四类震害因子的耦合影响,采用(Kolmogorov-Smirnov)K-S检验,在设定地震动参数下选取Gaussian分布、Log-Normal分布、Gumbel分布和Beta分布四种概率分布对该地区砌体建筑物的破坏状态概率分布参数进行拟合。通过均方根误差(Root Mean Square Error)RMSE进行拟合优度评价,最终建立基于Gaussian分布和Log-Normal分布的砌体建筑物破坏联合概率模型。最后,以华南地区三个城市典型砌体建筑物为例进行实例对比验证,将基于本文建立的建筑破坏概率模型推算出的砌体建筑群震害矩阵与基于单体结构分析得到的震害矩阵进行对比,与理论值最大偏差为0.033 3。研究表明：本文构建的集成概率方法能够获得更加合理的城市砌体建筑...     

湖泊蒸发观测与计算方法研究进展

湖泊蒸发是连接湖泊水分循环与能量平衡的关键纽带之一,认识湖泊蒸发过程对于理解湖-气间相互作用机制十分重要。然而,由于湖泊的形态结构、地理位置和气候背景不同,各湖泊在不同时间尺度上的蒸发特征不同。湖泊蒸发存在复杂的物理驱动过程、时间尺度相关的反馈以及空间异质性。准确地捕捉并量化湖泊蒸发过程,仍是湖沼学、水文学和气象学等学科的重要研究内容。本文首先介绍了湖泊蒸发的主要观测手段,然后概述了湖泊蒸发在日内、季节、年际、年代际四种不同时间尺度的变化特征;梳理了不同时间尺度各要素对湖泊蒸发的影响,以及湖泊蒸发与湖泊面积和深度等形态结构特性及纬度和海拔等地理因素之间的关系;对自涡度相关技术等高精度仪器应用以来的湖泊蒸发研究进展做出了总结。     

Matrix expression of thermal radiative characteristics for an open complex

The directionality of thermal radiance of a homogeneous isothermal non-black plane surface is totally decided by its directional emissivity, which depends on the complex dielectric constant and roughness of surface. It can be expressed by This paper proves that it is necessary to express emissivity by a matrix when a target becomes an inhomogeneous non-isothermal open complex with complicated inner geometric structure. The matrix describes the inner radiative interaction among components accurately and also expresses its thermal radiative directionality and structural characteristics completely. Advantages of matrix expression are as follows: first, the physical mechanics of effective emissivity of an open complex is described in a simple and perfect way; second, it becomes easy to understand the principle and method to retrieve components temperature from multi-angle thermal remotely sensed data; and third, the differences of directionalities between an open complex and a homogeneous isothermal non-black plane body are expressed by just using an effective emissivity matrix instead of an emissivity vector. Formula in classic physics is only the special case of matrix expression; therefore, the matrix is a universal unconditional expression to describe the directionality of thermal radiance.     

Monitoring volcanic hazard using eddy covariance at Solfatara volcano, Naples, Italy

An eddy covariance (EC) station was deployed at Solfatara crater, Italy, June 8–25, 2001 to assess if EC could reliably monitor CO₂ fluxes continuously at this site. Deployment at six different locations within the crater allowed areas of focused gas venting to be variably included in the measured flux. Turbulent (EC) fluxes calculated in 30-min averages varied between 950 and 4460 g CO₂ m⁻² d⁻¹; the highest measurements were made downwind of degassing pools. Comparing turbulent fluxes with chamber measurements of surface fluxes using footprint models in diffuse degassing regions yielded an average difference of 0% (±4%), indicating that EC measurements are representative of surface fluxes at this volcanic site. Similar comparisons made downwind of degassing pools yielded emission rates from 12 to 27 t CO₂ d⁻¹ for these features. Reliable EC measurements (i.e. measurements with sufficient and stationary turbulence) were obtained primarily during daytime hours (08:00 and 20:00 local time) when the wind speed exceeded 2 m s⁻¹. Daily average EC fluxes varied by ±50% and variations were likely correlated to changes in atmospheric pressure. Variations in CO₂ emissions due to volcanic processes at depth would have to be on the same order of magnitude as the measured diurnal variability in order to be useful in predicting volcanic hazard. First-order models of magma emplacement suggest that emissions could exceed this rate for reasonable assumptions of magma movement. EC therefore provides a useful method of monitoring volcanic hazard at Solfatara. Further, EC can monitor significantly larger areas than can be monitored by previous methods.     

Weak-contrast approximation of the elastic scattering matrix in anisotropic media

The plane-wave reflection and transmission coefficients at a plane interface between two anisotropic media constitute the elements of the elastic scattering matrix. For a 1-D anisotropic medium the eigenvector decomposition of the system matrix of the transformed elasto-dynamic equations is used to derive a general expression for the scattering matrix. Depending on the normalization of the eigenvectors, the expressions give scattering coefficients for amplitudes or for vertical energy flux.Computing the vertical slownesses and the corresponding polarizations, the eigenvector matrix and its inverse can be found. We give a simple formula for the inverse, regardless of the normalization of the eigenvectors. When the eigenvectors are normalized with respect to amplitudes of displacement (or velocity), the calculation of the scattering matrix for amplitudes is simplified.When the relative changes in all parameters are small, a weak-contrast approximation of the scattering matrix, based on the exactly determined polarization vectors in an average medium, is obtained. The same approximation is also derived directly from the transformed elasto-dynamic equations for a smooth vertically inhomogeneous medium, proving the consistency of the approximation.For monoclinic media, with the mirror symmetry plane parallel to the interface, the approximative scattering matrix is given in terms of analytic expressions for the non-normalized eigenvectors and vertical slownesses. For transversely isotropic media with a vertical axis of symmetry (VTI) and isotropic media, explicit solutions for the weak-contrast approximations of the scattering matrices have been obtained. The scattering matrix for amplitudes for isotropic media is well known. The scattering matrix for vertical energy flux may have applications in AVO analysis and inversion due to the reciprocity of the reflection coefficients for converted waves.Numerical examples for monoclinic and VTI media provide good agreement between the approximative and the exact reflection matrices. It is, however, expected that the approximations cannot be used when the symmetry properties of the two media are very different. This is because the approximation relies on a small relative contrast between the eigenvectors in the two media.Presented at the Workshop Meeting on Seismic Waves in Laterally Inhomogeneous Media, Castle of Trest, Czech Republic, May 22–27, 1995.     

动载作用下混凝土细观损伤的CT图像分割和灰度共生矩阵特征值研究

动载作用下混凝土随应力变化的系列CT图像完整地记录了混凝土内部破损过程,在CT图像中能直接观察到细观裂纹时混凝土往往已发生宏观破坏,挖掘混凝土材料损伤的细观信息是混凝土细观力学分析的核心。本文基于混凝土单轴动态压缩CT试验,获得不同应力阶段的横断面CT图像,分别应用图像分割技术和灰度共生矩阵法提取孔隙率和4个特征值,研究孔隙率和4个特征值随应力的变化规律。结果表明孔隙率随应力的增加显示出单调增加的总趋势,个别部位在应力水平较小时孔隙率随应力增加变化不大,反映出混凝土的压密效应。受压密效应影响,对比度在低应力水平时随应力增加有升有降,能量对应力变化不够敏感,相关性和同质性随应力的增加单调增加,对应力变化较为敏感,能反映混凝土整体损伤过程。对混凝土CT图像进行分区后再研究孔隙率、同质性和相关性随应力的变化,可以更好地反映动载作用下混凝土的细观损伤过程。     

图像重建块迭代算法中松弛参数选取的研究

在图像重建中,Landweber迭代算法是图像重建算法中的重要方法。本文将针对Landweber分块迭代算法中松弛参数的选取进行研究。在重建过程中采取对投影矩阵按投影角度分块的方法,选取特定的松弛参数。通过数值实验得出结论:对于按角度分块的块迭代算法,松弛参数选取为λ乘以块矩阵与其共轭转置矩阵乘积的最大特征值分之一,当采集完全投影数据,且λ接近(1/6)~(1/7)时效果最好。另外,本文按角度分块的做法和松弛参数的选取方法对于有限角度图像重建问题也是可行的。     

四川省房屋建筑易损性研究

研究房屋建筑易损性特征是防震减灾工作必不可少的环节之一，本文通过对多年工作的系统总结，得出四川省不同地区的房屋建筑易损性矩阵，全省房屋建筑抗震性能总体水平不高，这是造成地震灾害加重的主要原因之一，随着经济发展实力的提高，适当提高全省房屋建筑抗震设防水平可大大减轻特别是中强地震所造成的经济损失。     

地震人员伤亡的动态评估

引入了初始人员伤亡矩阵的概念,并借助于几个概念分布函数确定了砖砌体系结构的初始人员伤亡矩阵,建立了地震人员伤亡动态量化以在方法,并以唐山地震为例使用该方法进行再评估,得到了一些应用结果。     

公路桥梁地震易损性和震后恢复过程

采用专家调查的方式,统计的我国各种主要公路桥梁在不同烈度的地震发生后可能遭受的结构破坏和功能损失,以及地震发生后所需要的恢复时间。通过分析,给出了大、中、小种规模的斜拉、悬索、板梁式和拱式桥梁的结构破坏概率矩阵、设施损失率、功能损失率分布、损失功能恢复所需要的时间和桥梁在遭受破坏期间所损失的工作日数,并将研究结果与唐山地震之后的实际调查结果进行了对比,得到了一些有益的结果。