利用尾波包络线反演方法研究伽师强震群区地壳的非均匀结构

以单次散射模型、球形辐射源、一维线性速度模型为基础计算理论包络线,通过尾波归一化方法把理论与观测包络线结合起来形成非线性方程组,用迭代最小二乘方法反演伽师强震群区的散射系数空间分布.结果表明,伽师强震群区下方中上地壳中存在明显的非均匀结构.在10~25 km各深度图像上,研究区中部都存在一低散射系数区,这个低散射系数区在10 km、15 km深度被高散射系数区包围,在20 km深度仅其东侧存在高散射系数区.1997年的7次M_s≥6.0级地震发生在研究区中部的低散射系数区内或低散射系数区向高散射系数区过渡带的低散射系数一侧.     

横向不均匀介质中的面波及地球内部速度结构研究

近十年中地震面波的理论和应用研究有很大进展。由于资料增加、反演方法的改进及利用了振幅和波形等更多的面波信息，对于全球及区域上地幔速度结构及各向异性有了进一步的认识。同时对于显著不均匀三维介质中面波的传播特征，在理论探讨和实际观测方面取得了较多成果。本文对于面波研究的最新成果作了介绍。     

On terrain corrections in terrestrial heat flow

The theory of heat-flow variation over a Lees-type hill is well known for evaluating terrain effects on terrestrial heat flow. However, a hill cannot be converted into a valley by simply changing the sign of one of the terrain-defining parameters, nor can several Lees' hills be superposed to simulate a series of hills and valleys for correcting a terrain-induced disturbance of heat flow. The proper derivations for a Lees-type valley are presented and the superposition is compared with analytic solutions for two parallel semi-cylindrical valleys and two semi-spherical basins embedded in an otherwise planar ground. Generally, heat flux climaxes over central valleys or basins, and dips toward their margins where rapid change in topography occurs. Variation of heat flux induced by a three-dimensional terrain is relatively large, as compared to two-dimensional features, but its areal extent is relatively limited. The two-dimensional effects also extend relatively deeper. In applying two-dimensional algorithms to a three-dimensional terrain, the correction may be over- or underestimated depending on the distance from a borehole to prominent terrain features in the surrounding area.     

An optimized use of limited ground based topographic data for river applications

A prime requirement for hydrological applications,such as sediment budgeting or numerical modelling,is that produced Digital Terrain Models(DTMs)accurately represent the shape of landforms,especially for river reaches where data are not homogeneous.DTM error is a function of data point measurement accuracy and density and also of the field survey strategy when limited amounts of data will be acquired.This paper aims to advance the importance of the field survey strategy for the specific,but common cases,where only limited topographic data will be available.This methodology is based on the idea that any feature can be properly described by a set of cross sections and breaklines describing both main and secondary directions of the flow.Then,a longitudinal linear interpolation can be applied to the defined homogeneous zones.This morphologically oriented(MO)method that includes data acquisition strategy and interpolation,was validated using a reference DTM derived from LiDAR measurements.An estimation of the uncertainties also is suggested based on the distance of the nearest point and the local slope using a geographically weighted regression.The proposed MO method is typically applicable to Alpine river reaches characterized by multiple channels that may always be underwater and not navigable such as an alternate bar system with secondary and transverse channels.     

Three—dimensional P velocity structure in Beijing area

A detail three-dimensional P wave velocity structure of Beijing,Tianjin and Tangshan area(BTT area)was deter-mined by inverting local earthquake data.In total 16 048 Pwave first arrival times from 16048 shallow and mid-depth crustal earthquakes,which occurred in and around the BTT area from 1992to 1999were used.The first arrival times are recorded by Northern China Unived Telemetry Seismic Network and Yanqing-Huailai Digital Seismic Network.Hypocentral parameters of 1 132 earthquakes with magnitude ML=1.7-6.2 and the three-dimensional P wave velocity structure were obtained simultaneously.The inversion result reveals the com-plicated lateral heterogeneity of P wave velocity structure around BTT area.The tomographic images obtained are also found to explain other seismological observations well.     

基于Markov链模型的储层岩相随机模拟

在油气储层随机建模研究中，基于Markov链模型的方法是一类较受欢迎的技术，同时也是一类不成熟的技术，问题的症结之一在于侧向的转移概率矩阵很难求取，针对这种情况，作者在深入理解Walther相律的基础上，借鉴模拟退火算法的相应思路，提出了一种岩相模拟的新方法，该方法依据不同岩相的百分比进行随机模拟得到一幅初始图象，而后以按岩相组织剖面得到的垂向和侧向的岩相转移概率矩阵的相似性作为判别标准对图象进行扰动，直至得到满意的图象，二维模型试算结果表明了这种岩相随机模拟方法的可行性。     

Effect of anelastic and scattering structures of the lithosphere on the shape of local earthquake coda

A simple model of single acoustic scattering is used to study the dependence of the shape of local earthquake coda on the anelastic and scattering structures of the lithosphere. The model is applied to the coda of earthquakes located near Stone Canyon, central California, and provides an explanation for the features observed in the data, which include an interesting temporal variation in the coda shape. A surficial layer with aQ of 50 and thickness of 10 or 25 km underlain by a zone with aQ of 1000 extending to the bottom of the lithosphere, together with a scattering scale length,a, that varies with depthz according to the relationa=0.3 exp[-(z/45)²] are found to constitute the simplest structure of the medium compatible with the coda data and with body and surface wave attenuation data. The profile of heterogeneity sizes implies that the scattering strength increases strongly with depth, a constraint required by the necessity to boost the energy of the later coda without forcing the intrinsicQ to be excessively high in the uppermost mantle. This constraint is viewed as an artifact of the single scattering model which overstimates the scattering coefficient due to the neglect of multiple scattering. The observed temporal variation of the signal is difficult to explain by a simple change of the intrinsicQ at some depth. Rather, it is suggested that the scattering properties at depth changed with time through a variation of the fractional rms velocity fluctuation on the order of one percent.     

岩石破裂过程TMD耦合数值模型研究

从岩石的细观结构层次出发,应用损伤力学和热弹性理论,对热力耦合作用下岩石破裂过程中热-应力-损伤相互作用关系进行了分析。初步建立了细观热-应力-损伤（TMD model,thermal-mechanical-damage）耦合数值模型,并在岩石破裂过程分析系统RFPA2D中加实现。运用该数值模型计算模拟均匀与非均匀材料试样在热力耦合作用下的应力分布及破坏形 态,通过与理论及试验结果对比,证明了该数值模型的合理性和有效性。     

Influence of Heterogeneity of Mechanical Properties on Hydraulic Fracturing in Permeable Rocks

Summary Numerical simulations of circular holes under internal hydraulic pressure are carried out to investigate the hydraulic fracture initiation, propagation and breakdown behavior in rocks. The hydraulic pressure increases at a constant rate. The heterogeneity of the rocks is taken into account in the study by varying the homogeneity index. In addition, the permeability is varied with the states of stress and fracture. The simulations are conducted by using a finite element code, F-RFPA^2D, which couples the flow, stress and damage analyses. The simulation results suggest that the fracture initiation and propagation, the roughness of the fracture path and the breakdown pressure are influenced considerably by the heterogeneity of rocks. The hole diameter elongation and the stress field evolution around the fracture tip during the fracture propagation can also provide useful information for the interpretation of the hydraulic fracturing behaviour.