浅谈地震保险的必要性及其作用

以防震减灾为视角，针对地震灾害对社会发展及人民生活造成的重大损失，从强化地震保险意识的角度，论述地震保险的必要性和重大意义。同时对地震保险在心理、补偿、救灾、缓冲、防灾及附带等方面的作用加以阐述，从而证明要减少地震灾害的损失，充分认识地震保险的必要性及其作用，建立和完善有效的保险机制，对于减轻地震灾害，促进社会发展和人民安居乐业是一项必要的基础性工作。     

2001年四川省雅江-康定6.0级地震烈度考察及发震构造背景雏议

2001年2月23日四川省发生的6.0级地震，其宏观震中位于雅江县与康定县之间的高山峡谷地带．极震区烈度可达Ⅷ 度，Ⅷ 度区、Ⅶ 度区和Ⅵ 度区面积分别为180km2，1472 km2和3998 km2，等震线分布总体形态呈椭圆形，长轴近南北向．震区建筑物普遍受损，极震区山地灾害较为严重．初步分析认为，此次地震的发生与理塘断裂和玉农希断裂的长期活动有关，主干断裂活动造成的地壳应力调整和集中，最终在块体内部的次级断裂上释放而产生地震．      

唐山地震震源区构造应力场强度的初步分析

利用唐山地区丰富的震源机制解资料反演构造应力场，并结合构造物理研究手段，对唐山地震震源区现代构造应力场的3个主应力绝对量值进行了估算，给出了地壳介质参数取理论平均值的情况下，唐山地震震源区构造应力场３个主应力量值随深度的变化趋势．研究分析了摩擦系数c、孔隙压力P0和应力形因子 等参数与应力值的关系，发现最大主应力1的垂向增长率与摩擦系数c成正比，与孔隙压力P0和应力形因子成反比；而最小应力值3的垂向增长率与孔隙压力P0成正比，与摩擦系数c和应力形因子成反比．本文对应用震源机制解资料反演完整的应力张量进行了一次有意义的尝试．      

由GPS观测结果推导中国大陆现今水平应变场

以中国大陆及周边近400个GPS测站的水平运动速率为基础，给出了现今地壳水平应变场结果表明:①中国大陆水平应变为西强东弱，剪应变数值大于正应变数值(绝对值)，应变量级一般为10-8/a，局部区域达到10-7/a，但应变分布不均匀；②南北向应变最突出的部位为中国西南部西段的喜马拉雅条带、西北部的36N~42N段及柴达木断块的北缘；③东西向应变西边缘变化最大．此外，由西向东还具有正负交替的变化特征；④REN(东-北向剪切应变)与Rmax(最大剪切应变)数值较大的区域分别是喜马拉雅条带、西北部的36N～42N段、柴达木断块的西部、川滇菱形块体，以及阿拉善、祁连及塔里木断块的交界区；⑤青藏块体周边以面收缩为主，内部则以面膨胀为主．其以北的地区以面收缩为主．西界数值最大，东部数值最小(除燕山构造带外)；⑥西部西区主压应变为南北向，主张应变为东西向．西部东缘区主压应变为近东西向，主张应变为近南北向．川滇菱形块体主应变的方向发生了很大的变化，北部地区为东西压南北张，南部地区则恰好相反；⑦中国大陆的应变模式可能是断块模式与连续模式的组合．此外，小尺度优势应变可能是剪切应变．造成上述结果除与印度板块的碰撞及边界耦合有关外，还与深部物质的活动及地壳介质的物性有密切的关系．必须指出，由于GPS测站在空间上分布的不均匀性，那么，由此而来的应变场，其应变尺度也不一样．      

1996年包头6.4级地震的地壳应变特征

根据GPS观测资料求出的水平地应变和由跨断层垂直形变计算出的速率强度累积率，研究了包头-大同地区1992～1995年、1995～1996年和1996～1999年的各时期的应变特征，并对包头6.4级(1996年5月3)地震前后的地应变进行对比，认为以压应变为主导的高值区可能是未来强震孕育的地区．面应变、主压应变、剪应变和趋势累积率同时较高的地区，强震危险性较大．一般低应变区和张应变为主导的地区，孕育强震的可能性小，属于比较稳定的地区．1992～1999年包头-大同地区的GPS水平应变的演变，反映了1996～1998年地震幕的孕育发展及结束的全过程．以压应变为主的高应变区和应变梯度带可作为未来强震危险区的判定标志之一．      

用应变积累释放模型研究中国大陆地块分区地震活动

依据活动地块假说及活动边界的划分，对中国大陆地震活动进行分区．在此基础上，应用应变积累释放模型，通过地震应变积累释放的计算，研究了中国大陆各活动地块分区的地震活动性，初步探讨了各地块的地震活动阶段划分．此外，还对应变积累释放模型的应用前提条件进行了必要的讨论，并对模型结果中可能出现的问题给予一定分析．      

Efficient methods for solving water flow in variably saturated soils under prescribed flux infiltration

The non-linear solvers in numerical solutions of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties and in this paper we address the gravity term and the prescribed-flux boundary in the Picard iteration. The problem of the gravity term in the Picard iteration is iteration-to-iteration oscillation as the gravity term is treated, by analogy with the time-step advance technique, ‘explicitly’ in the iteration. The proposed method for the gravity term is an improvement of the ‘implicit’ approach of Zhang and Ewen [Water Resour. Res. 36 (2000) 2777] by extending it to heterogeneous soil and approximating the inter-nodal hydraulic conductivity in the diffusive term and the gravity term with the same scheme. The prescribed-flux boundary in traditional methods also gives rise to iteration-to-iteration oscillation because there is no feedback to the flux in the solution at the new iteration. To reduce such oscillation, a new method is proposed to provide such a feedback to the flux. Comparison with traditional Picard and Newton iteration methods for a wide range of problems show that a combination of these two proposed methods greatly improves the stability and consequently the computational efficiency, making the use of small time step and/or under-relaxation solely for convergence unnecessary.     

Screening model for volatile pollutants in dual porosity soils

This paper develops mass fraction models for transport and fate of agricultural pollutants in structured two-region soils. Mass fraction index models, based on a semi-infinite domain solution, are derived that describe leaching at depth, vapor losses through soil surface, absorption, and degradation in the dynamic- and stagnant-water soil regions. The models predict that leaching is the result of the combined effect of the upward vapor-phase transport relative to downward advection, residence time relative to half-life, dispersion, and lateral diffusive mass transfer. Simulations show that leached fraction of volatile compounds does not always decrease monotonically with increased residence time relative to the pollutant half-life, as a result of complex interactions among the different physical and biochemical processes. The results show that leaching, volatilization, and degradation losses can be affected significantly by lateral diffusive mass transfer into immobile-water regions and advection relative to dispersion (i.e. Peclet number) in the mobile-water regions. It is shown that solute diffusion into the immobile phase and subsequent biochemical decay reduces leaching and vapor losses through soil surface. Potential use of the modified leaching index for the screening of selected pesticides is illustrated for different soil textures and infiltration rates. The analysis may be useful to the management of pesticides and the design of landfills.     

Common-Reflection-Surface Stack for Converted Waves

The finite-offset (FO) common-reflection-surface (CRS) stack has been shown to be able to handle not only P-P or S-S but also arbitrarily converted reflections. It can provide different stack sections such as common-offset (CO), common-midpoint (CMP) and common-shot (CS) sections with significantly increased signal-to-noise ratio from the multi-coverage pre-stack seismic data in a data-driven way. It is our purpose in this paper to demonstrate the performance of the FO CRS stack on data involving converted waves in inhomogeneous layered media. In order to do this we apply the FO CRS stack for common-offset to a synthetic seismic data set involving P-P as well as P-S converted primary reflections. We show that the FO CRS stack yields convincing improvement of the image quality in the presence of noisy data and successfully extracts kinematic wavefield attributes useful for further analyses. The extracted emergence angle information is used to achieve a complete separation of the wavefield into its P-P and P-S wave components, given the FO CRS stacked horizontal and vertical component sections.     

Optimal placement of dampers and actuators based on stochastic approach

A general method is developed for optimal application of dampers and actuators by installing them at optimal location on seismic-resistant structures. The study includes development of a statistical criterion, formulation of a general optimization problem and establishment of a solution procedure. Numerical analysis of the seismic response in time-history of controlled structures is used to verify the proposed method for optimal device application and to demonstrate the effectiveness of seismic response control with optimal device location. This study shows that the proposed method for the optimal device application is simple and general, and that the optimally applied dampers and actuators are very efficient for seismic response reduction.