高温冻土的蠕变特性试验及蠕变模型研究

为了研究高温冻土蠕变变形特征以及各影响因素对蠕变的作用, 分别在含水量15%、 25%及35%, 荷载100 kPa、 200 kPa及300 kPa, 温度-1.5 ℃、 -0.7 ℃及-0.3 ℃的条件下开展了室内单轴蠕变试验, 分析在无侧限条件下高温冻土在不同温度、 荷载及含冰量条件下的蠕变变形特性。结果表明： 在当前试验条件下, 冻土蠕变变形非常可观, 且蠕变曲线都没有出现渐进流阶段; 温度是影响冻土蠕变的最重要的外在因素, 而含冰量是影响冻土蠕变的关键内在因素; 在高含冰量条件下温度及荷载的改变对蠕变速率的影响非常显著, 甚至引起量级上的差别。在现有试验条件下, 高温冻土蠕变过程可利用Burgers黏弹性模型来较好地描述。     

Temporal and spatial variations of post-seismic deformation following the 1999 Chi-Chi, Taiwan earthquake

We use GPS displacements collected in the 15 months after the 1999 Chi-Chi, Taiwan earthquake  ( M _w 7.6)  to evaluate whether post-seismic deformation is better explained by afterslip or viscoelastic relaxation of the lower crust and upper mantle. We find that all viscoelastic models tested fail to fit the general features in the post-seismic GPS displacements, in contrast to the satisfactory fit obtained with afterslip models. We conclude that afterslip is the dominant mechanism in the 15-month period, and invert for the space–time distribution of afterslip, using the Extended Network Inversion Filter. Our results show high slip rates surrounding the region of greatest coseismic slip. The slip-rate distribution remains roughly stationary over the 15-month period. In contrast to the limited coseismic slip on the décollement, afterslip is prominent there. Maximum afterslip of 0.57 m occurs downdip and to the east of the hypocentral region. Afterslip at hypocentral depths is limited to the southern part of the main shock rupture, with little or no slip on the northern section where coseismic slip was greatest. Whether this results from along strike variations in frictional properties or dynamic conditions that locally favour stable sliding is not clear. In general, afterslip surrounds the area of greatest coseismic slip, consistent with post-seismic slip driven by the main shock stress change. The total accumulated geodetic afterslip moment is  3.8 × 10¹⁹ N m  , significantly more than the seismic moment released by aftershocks,  6.6 × 10¹⁸ N m  . Afterslip and aftershocks appear to have different temporal evolutions and some spatial correlations, suggesting that aftershock rates may not be completely controlled by the rate of afterslip.     

Relaxation time distribution from time domain induced polarization measurements

We propose an algorithm for inverting time domain induced polarization data to a relaxation time distribution. The algorithm is based on the (Tikhonov) regularized solution of the 1^st kind Fredholm integral equation. We test the algorithm on synthetic data, and show its robustness for a noise level, typical of laboratory time domain measurements. We also show that, for the inversion purpose, the time domain data must be obtained with the different current wavelengths. We then test the algorithm on the experimental data recently obtained on brine-saturated medium-grained quartz sand (average grain diameter of 4 × 10⁻⁴ m), and on sand mixtures. For the medium-grained sand, relaxation time distribution contains a main peak at 25 s. Different amounts (3%, 8% and 12%) of fine-grained quartz sand (average grain diameter of 1.12 × 10⁻⁴ m) were added to the medium-grained quartz sand. For the sand mixture, an additional peak is observed in the relaxation time distributions, in the time range from 1.0 to 2.5 s. The magnitude of the second peak increases with the increase of the fine-grained sand content. Therefore, the experimental data show that peaks in the relaxation time distributions are related to the grain size. On the basis of the known peak location, and of the known grain size value, we obtained the values of the diffusion coefficient, which were found to be of the same order of magnitude as those in the bulk solution.     

Construction of constant-Q viscoelastic model with three parameters

The popularly used viscoelastic models have some shortcomings in describing relationship between quality factor (Q) and frequency, which is not consistent with the observation data. Based on the theory of viscoelasticity, a new approach to construct constant-Q viscoelastic model in given frequency band with three parameters is developed. The designed model describes the frequency-independence feature of quality factor very well, and the effect of viscoelasticity on seismic wave field can be studied relatively accurate in theory with this model. Furthermore, the number of required parameters in this model has been reduced fewer than that of other constant-Q models, this can simplify the solution of the viscoelastic problems to some extent. At last, the accuracy and application range have been analyzed through numerical tests. The effect of viscoelasticity on wave propagation has been briefly illus-trated through the change of frequency spectra and waveform in several different viscoelastic models.     

Reflection/transmission laws for slowness vectors in viscoelastic anisotropic media

The reflection/transmission laws (R/T laws) of plane waves at a plane interface between two homogeneous anisotropic viscoelastic (dissipative) halfspaces are discussed. Algorithms for determining the slowness vectors of reflected/transmitted plane waves from the known slowness vector of the incident wave are proposed. In viscoelastic media, the slowness vectors of plane waves are complex-valued, p = P + iA, where P is the propagation vector, and A the attenuation vector. The proposed algorithms may be applied to bulk plane waves (A = 0), homogeneous plane waves (A ≠ 0, P and A parallel), and inhomogeneous plane waves (A ≠ 0, P and A non-parallel). The manner, in which the slowness vector is specified, plays an important role in the algorithms. For unrestricted anisotropy and viscoelasticity, the algorithms require an algebraic equation of the sixth degree to be solved in each halfspace. The degree of the algebraic equation decreases to four or two for simpler cases (isotropic media, plane waves in symmetry planes of anisotropic media). The physical consequences of the proposed algorithms are discussed in detail. vcerveny@seis.karlov.mff.cuni.cz     

根据因果岸三维地震资料研究地裂缝形成机理

塌隐地裂缝，改变了地层的几何结构，破坏了三维地震反射层的连续性，出现了反射波同相轴错断、反射层产状变化、反射波能量转换、反射波增多或减少等特征。据此利用石塘湾因果岸三维地震资料对其地裂缝地质灾害形成进行了研究。通过对实际资料的分析，指出本区地裂缝类型为张跌型，并以较多的资料翔实地说明了诱发本区地裂地质灾害不可分割的五个因素，从而为治理地裂缝地质灾害提供了丰富的地质资料和治理方向。     

随机地震下阻尼器的优化设置

根据粘弹性阻尼结构的减震原理,利用随机振动理论中的李雅普诺夫方程求解粘弹性阻尼结构反应的方差矩阵,再以层间位移的方差矩阵为控制函数进行阻尼器的优化设置,并给出发实例分析。     

微粒对多晶冰流变行为的影响(Ⅰ):应力下的流变行为

研究了应力下微粒对多晶冰流变行为的影响.结果表明:微粒分布在晶界或同时分布在晶界和晶内都将增加多晶冰的流变速率.当微粒只分布在晶界时,含微粒浓度为1 wt.%的多晶冰有最高的最小流变速率;当微粒同时分布在晶内与晶界时,多晶冰的最小流变速率基本与微粒浓度无关.微粒提高流变速率的原因来源于提高了流变过程中的位错密度,分布在晶界和晶内的微粒可以分别通过阻碍晶界滑移,发展内应力及Frank-Read位错增值机制提高流变过程中的位错密度.     

土体动非线性黏弹性模型及其ABAQUS软件的实现

基于土的动应力-应变关系--Davidenkov骨架曲线,采用破坏剪应变幅上限值作为分界点,对Davidenkov骨架曲线进行了修正,即当剪应力值大于破坏剪应力值时,土体产生破坏,土体的动剪切模量采用破坏后的动剪切模量;参考Mashing法则,构造了修正后的Davidenkov骨架曲线土体加卸载对应的应力-应变滞回圈曲线。基于ABAQUS软件的操作平台,开发了土体动弹塑性本构模型的子程序。选择南京某个典型软弱场地为研究对象,首先输入余弦加速度时程曲线,然后选取一条强震加速度记录作为基岩输入地震动,对不同峰值加速度水平下该场地的地震反应进行了二维有限元非线性分析,计算结果验证了所建模型的准确性和可用性。