井下地电阻率观测中地表电流干扰影响计算

本文把地电阻率台站地下介质简化为水平三层均匀介质模型,以点电流模拟地表干扰电流,针对对称四极观测装置,计算了在不同电性结构中的不同深度观测时,地表电流干扰源对对称四极装置地电阻率观测的影响,计算得到:地表干扰源对电阻率观测的影响取决于电性结构的类型和层参数、供电电极和测量电极的埋深以及避开干扰源的距离.本文研究结果对实施井下电阻率观测中台址电性结构的选择、电极埋深、干扰源避让距离等有参考意义.     

使用地面单台三分向地震记录资料评估场地效应方法的研究

利用唐山强震观测台阵提供的井下和地面数字测震资料，研究了用单台地震记录评估场地效应的方法。利用井下基岩地震记录作参考台，用地面和井下的记录求取经验传递函数。同时应用Nakamura单台评估方法求解Nakamura谱比率。通过对比，证实了对场地效应的评估可以使用地面单台测震资料进行，不必参考基岩台的资料。因此，该方法是评估重要工程场地效应的一种简便可行的方法。使用该方法可以降低地震小区划的成本。     

地表点电流对地电阻率观测的影响

将地电阻率台站地下介质简化为水平层状均匀介质模型,以点电流模拟地表干扰电流源,针对对称四极观测装置,计算干扰源电流增减及位置变化对地电阻率观测的影响。结果表明,地电阻率的受干扰程度取决于地表干扰点电流的大小、在测区的方位及距中心点的距离。     

凹凸体位错非均匀性对断层引起地表位移的影响

在已有的凹凸体震源模型基础上考虑凹凸体位错非均匀性,提出改进的凹凸体位错模式.新的位错模式随机设置凹凸体上各子断层位错量,同时在凹凸体与背景区交界处设置位错平缓变化的区域.以1989年美国 LomaPrieta地震为例,计算不同断层位错模式下引起的地表位移,并与反演得到真实情况下地表的位移作对比,验证改进凹凸体位错模式的可靠性.结果显示,使用改进的渐变凹凸体非均匀位错模式计算时,与真实情况下竖向地表位移差大于8cm 的区域面积为 43km2,相较凹凸体均匀位错模式缩减25%;水平地表位移差大于8cm 区域面积为117km2,相较凹凸体均匀位错模式缩减达31%.利用改进的渐变凹凸体非均匀位错模式,计算1679年三河—平谷大地震在北京地区形成的地表位移场.北京市外,最大竖向地表位移大于4．8m,最大水平地表位移大于2．6m,均出现在三河附近;北京市内,通州区至平谷区一带地表位移最大,最大水平地表位移大于 2．6m,最大竖向地表位移大于1m.研究结果可为今后北京地区的抗震设防提供参考.     

水平螺旋度与沙尘暴的动力学关系研究

应用2002～2006年高空流场和地面观测资料,计算近地面至500 hPa水平螺旋度的大小.结果表明：水平螺旋度负值中心值越大,500 hPa到近地面风速越大,西风增强,风速垂直切变越大,辐合上升运动越强,形成沙尘暴的强度就越强.水平螺旋度负值中心常常在河西走廊附近最强,导致其下游东南方发生沙尘暴.水平螺旋度负值中心与其下游沙尘暴发生强度有一致的对应关系：当水平螺旋度负值中心≤-200 m²/s²时,未来24小时内该区下游将有沙尘天气出现,当≤-600 m²/s²时,6小时内该区下游将有能见度小于500 m的强沙尘暴天气出现,当≤-1000 m²/s²时,6小时内该区下游将有能见度小于50 m的特强沙尘暴天气出现.     

Estimates of the energy of surface waves from atmospheric explosions and the source parameters of the Tunguska event

The results of the numerical modeling of experimental concentrated atmospheric explosions, both nuclear and chemical, are presented. The surface pressure spectra are calculated and constraints are gained on their relation to the energy of surface seismic waves and to earthquake magnitudes. For a more accurate determination of the source parameters of the Tunguska event (June 30, 1908), atmospheric explosions at various altitudes and vertical falls of meteoroids are modeled numerically and their seismic effect is estimated.     

Near-surface layer traveltime inversion: a synthetic example1

The near-surface layer is modelled as a constant-velocity layer with varying thickness. The base of the layer is described by a B-spline curve. The optimum model is calculated by minimizing, with respect to the model parameters, the difference between traveltimes predicted by the model and those observed in the data. Once a model has been produced, corrections that are dependent on the raypath geometry through the near-surface layer can be calculated. The effect of the near-surface layer is normally considered to be consistent at each shot or geophone station for all traveltimes arriving at that location (the surface-consistent approximation). This assumption linearizes the problem, allowing timeshifts to be calculated and the traveltimes corrected to a chosen datum, representing static corrections. The single correction at each point is an averaged correction, based on an assumption that is particularly inaccurate in the presence of lateral variations of velocity or thickness of the surface layer, in the presence of large surface layer velocities or in the presence of a thick surface layer. The method presented considers the non-linear relationship between data and model explicitly, hence the correction that is dependent on the raypath. Linearization removes this dependence and reduces the problem to a surface-consistent approximation. The method is applied to synthetic data calculated from a model with surface layer variations. Comparisons are made between the corrected data resulting from the method described here and the conventional surface-consistent approach. From these results it becomes apparent that the near-surface layer inversion method presented here can reproduce accurate models and correct for near-surface layer effects in cases where conventional methods encounter difficulties. Additionally the method can be readily extended to 3D.     

INVERSION OF GRAVITY AND MAGNETIC DATA FOR THE LOWER SURFACE OF A 2.5 DIMENSIONAL SEDIMENTARY BASIN1

Gravity and magnetic data have been inverted to obtain the continuous lower surface of a 2.5 dimensional sedimentary basin. The non-linear problem is linearized and a solution is calculated through a recursive process until the predicted data matches the observed data. An average model is then calculated and a resolution analysis shows which features are uniquely determined. The results of individual inversion indicate that a final solution is initial model dependent but the average models are independent of the initial model except at the margins. The average model for the magnetic solutions have uniformly smaller spreads than the gravity solutions. The algorithms were applied to data from the Sanford Basin in North Carolina. The results indicate that the basin is asymmetrical in shape with a maximum depth of 3.2 km. Comparing these results with those obtained from a generalized linear inverse (GLI) algorithm indicate that the higher-frequency features determined from the GLI algorithm are not resolved.     

基于数字图像测量系统的三轴试样变形研究

基于亚像素角点检测的图像测量技术能够记录试样表面方块角点每一时刻的位移,从而获得试样的表面应变场。将基于试样表面6行方块分别计算的应变与基于全部6行方块的整体平均计算得到的结果进行比较,显示中间4行方块区域的变形基本不受端部约束影响。还对采用试样上、中、下三个部分与整个试样的测量结果获得的应力与应变特性进行了对比分析。该方法能够测量任意时刻试样的轴向应变场,通过轴向应变场的变化可以准确地描述试样剪切带的产生、发展以及完全形成的全过程。     

具有低频超宽衰减域的部分埋入式表面波屏障

针对低频Rayleigh表面波,设计了部分埋入式工字形截面周期波屏障。利用有限元方法计算了结构的频散曲线,分析了带隙的形成机理,讨论了屏障埋入土体深度和截面参数对带隙的影响,在此基础上设计了具有低频超宽衰减域的梯度及分段梯度波屏障并计算了其传输谱。结果表明:周期波屏障存在较宽带隙,板埋入深度和端部尺寸是影响带隙的关键参数,通过参数调节可实现不同频段Rayleigh波的调控。工字形变截面波屏障比等截面具有更优越的隔震性能且节省材料。梯度及分段梯度波屏障显著拓宽了衰减域的频率范围,对1.5~20 Hz范围内的Rayleigh表面波实现了全覆盖,用小尺寸控制了大波长。     

南北地震带北段地温场的初步探讨

本文利用39个地温梯度值及估算的热流值讨论了本区地温场的分布特征。并求取了居里等温面。结果表明:(1)全区平均地温梯度为3.2℃/100m,平均热流值为1.5HFU;(2)居里面的深度为14—34公里,强震主要发生在居里等温线的梯度带上。最后就地温场与地震的关系做了初步分析。     

Vibrations of the elastic half-space due to vertical surface loads

The displacements of the surface of the homogeneous elastic half-space subjected to a vertical concentrated load acting at the surface are calculated. The solution is shown to approach Lamb's simple formula for the far-field. The solution for the concentrated load can be used to calculate displacements due to distributed loads with high accuracy and economically. The results are applied to define the regions ‘near-field’ and ‘far-field’. Some near-field phenomena are described.     

Dynamic shear modulus of undisturbed soil under different consolidation ratios and its effects on surface ground motion

The dynamic shear modulus for three types of undisturbed soil under different consolidation ratios is presented by using the resonant column test method. Its effects on surface ground motion is illustrated by calculation. The test results indicate that the power function is a suitable form for describing the relationship between the ratio of the maximum dynamic shear modulus due to anisotropic and isotropic consolidations and the increment of the consolidation ratio. When compared to sand, the increment of the maximum dynamic shear modulus for undisturbed soil due to anisotropic consolidation is much larger. Using a one-dimensional equivalent linearization method, the earthquake influence factor and the characteristic period of the surface acceleration are calculated for two soil layers subjected to several typical earthquake waves. The calculated results show that the difference in nonlinear properties due to different consolidation ratios is generally not very notable, but the degree of its influence on the surface acceleration spectrum is remarkable for the occurrence of strong earthquakes. When compared to isotropic consolidation, the consideration of actual anisotropic consolidation causes the characteristic period to decrease and the earthquake influence factor to increase.     

山东地区地震加速度反应谱特征周期的统计研究

本文收集整理了2010年至2013年期间山东地区830个地震安全性评价工程的场地类型、覆盖层厚度、钻孔剪切波速、自由基岩地震动峰值加速度和地震加速度反应谱特征周期等资料,在计算出场地指数的基础上,采用回归分析方法,对特征周期与基岩峰值加速度、场地指数之间的关系进行了统计研究,并得出了三者之间的定量表达式。结果表明,特征周期与基岩峰值加速度的自然对数呈线性相关,与场地指数呈非线性相关。     

Numerical calculations of the perturbation of normal heat flow by topographic and local conductivity structure

A numerical method is used to investigate the effect of topographic and local thermal conductivity anomalies on near-surface heat flow for two-dimensional models. Heat flow associated with a sloping topographic structure is calculated. Also, the effects of a fault structure associated with the sloping topography are considered. Vertical and horizontal heat flow components are calculated alone; the surface of the earth as well as throughout the whole region of interest. The results indicate that surface heat flow is substantially affected by topographic relief and the horizontal heat flow component associated with topographic features can be large. Also, regional heat flow is greatly perturbed by local thermal conductivity anomalies and the effect of topographic features may be considerably modified by the subsurface structure.     

抚仙湖集水域地表径流入湖水量模拟

采用适用于无资料流域、参数较少的SCS模型计算抚仙湖集水域地表径流量.模型考虑了集水域下垫面条件的空间差异,利用Maplnfo/Arc view软件按照土地利用方式与土壤类型的不同,把集水域划分为若干个水文响应单元,分别计算产流量,较准确地模拟了入湖径流量.通过对梁王河流域和大鲫鱼沟流域实测降雨径流资料的分析与反演,提出了适合该区域的产流计算CN值.在对CN值作坡度修正后再应用到其它无观测数据区域.通过模型计算得到的抚仙湖集水域2005年3月1日-2006年2月28日地表径流量为1.74×108 m3,陆面径流系数为0.395.模型为指导抚仙湖集水域径流观测及入湖污染物负荷的计算提供依据.     

Determination of the near field distance for point and line sources acting on the surface of an homogeneous and viscoelastic half-space

A methodology to determine the near field region, defined as the region where volumetric waves significantly influence the surface vibration levels, is presented. This interference region between Rayleigh and volumetric waves is calculated for the case of an homogeneous and viscoelastic half-space excited by a surface point source and a surface infinite line source. This calculation is performed following an analytical approach. The effects of the mechanical properties of the ground on the dimensions of this interference region, characterized here by the near field distance, are also investigated. The results show near field distances for both infinite line and point sources and for five different real grounds, typical of the Barcelona (Spain) metropolitan area.     

核－幔边界的动力学背景

根据传播矩阵方法，并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷，以求取核－幔边界的动力学背景．计算结果显示：１．所求的核－幔边界起伏图像与Ｈａｇｅｒ等根据格林函数方法所求得的核－幔边界的起伏在全球范围内基本相符．２．核－幔边界处的环型场仅在数量上降为地表处的环型场的１／８左右，而极型场较地表处的极型场的流动图像有显著变化，数值也增为地表处极型场的３倍左右．     

核－幔边界的动力学背景

根据传播矩阵方法，并把由联合反演得到的同时满足长波地形起伏、板块运动速度、重力位异常资料以及地震层析先验知识的全地幔三维异常密度作为载荷，以求取核－幔边界的动力学背景．计算结果显示：１．所求的核－幔边界起伏图像与Ｈａｇｅｒ等根据格林函数方法所求得的核－幔边界的起伏在全球范围内基本相符．２．核－幔边界处的环型场仅在数量上降为地表处的环型场的１／８左右，而极型场较地表处的极型场的流动图像有显著变化，数值也增为地表处极型场的３倍左右．     

球层岩石层模型对表面荷载的响应

本文给出了在球层模型下,地球岩石层对表面荷载响应的解析解。计算了地球岩石层对圆盘形荷载和正弦型荷载的响应,及由此产生的重力异常。结果表明:1.岩石层弯曲不仅与其厚度有密切关系,而且与荷载尺度有关。当岩石层厚度与荷载尺度相比很小时,其弯曲与厚度关系不大,并呈现薄板效应。2.弹性模型岩石层在外加荷载下能承受数千巴的弯曲应力。3.弯曲应力在岩石层中的分布较为规则。由于大尺度荷载下岩石层出现薄板效应,传统的由弯曲刚度估算等效岩石层厚度的部分结果需要重新估价。