Response spectra for differential motion of columns paper II: Out-of-plane response

It is shown that the common response spectrum method for synchronous ground motion can be extended to make it applicable for earthquake response analyses of extended structures experiencing differential out-of-plane ground motion. A relative displacement spectrum for design of first-story columns SDC (T, T_T, ζ, ζ_T, τ, δ) is defined. In addition to the natural period of the out-of-plane response, T, and the corresponding fraction of critical damping, ζ, this spectrum also depends on the fundamental period of torsional vibrations, T_T, and the corresponding fraction of critical damping, ζ_T, on the “travel time,” τ (of the waves in the soil over a distance of about one-half the length of the structure), and on a dimensionless factor δ, describing the relative response of the first floor. The new spectrum, SDC, can be estimated by using the empirical scaling equations for relative displacement spectra, SD, and for peak ground velocity, v_max. For recorded strong-motion acceleration, and for symmetric buildings, the new spectrum can be computed from Duhamel's integrals of two uncoupled equations for dynamics equilibrium describing translation and rotation of a two-degree-of-freedom system. This representation is accurate when the energy of the strong-motion is carried by waves in the ground the wavelengths of which are one order of magnitude or more longer than the characteristic length of the structure.     

地震强度指标选择对长周期钢结构余震IDA分析的影响

结构在主震作用下发生损伤,由于刚度和强度退化导致结构周期的延长,再加上高阶周期的影响,在随后的余震对主震损伤结构抗震性能的影响分析中,基于弹性周期的传统谱加速度可能不适用,有必要对地震动强度指标的选择进行进一步研究。本文以9层Benchmark钢框架为研究对象,选取10种地震动强度指标,分别从离散性和相关性两个维度对比分析了其在余震IDA分析中的优劣性。研究结果表明:与最大层间位移角相比,最大残余层间位移角更能够准确地量化主震对结构造成的损伤;对于钢框架结构,基于一段周期的几何平均谱加速度的IM_Boj&Lev指标是最优的地震动强度指标,相对于S_a（T₁）离散性降低了30%左右,相关性提高了6%左右。     

基于核磁共振实验的低渗透砂岩岩电参数分类及应用——以东营凹陷南坡沙四段为例

低渗透砂岩储层孔隙结构复杂, 储层有效性识别及饱和度准确计算难度较大.笔者以东营凹陷南坡沙四段（Es4）低渗透砂岩为研究对象, 根据压汞、物性、薄片及核磁等资料, 将研究区孔隙结构分为三大类、五小类.在岩样孔隙结构分类基础上, 明确了孔隙结构类型与岩电参数之间存在确定的关系, 而核磁共振T2谱定量特征参数在一定程度上能够表征孔隙结构类型及其细节信息, 通过提取T2谱中T2几何平均值（T2g）、T2均值（）、峰度（KG）、可动流体分量（Smf）及区间孔隙分量等孔隙结构参数, 建立了基于核磁T2谱特征参数的孔隙结构识别图版, 显示核磁T2谱孔隙结构参数对不同类型的储层有较好的识别效果, 进而探讨核磁孔隙结构参数和岩电参数之间的关系, 结果表明, T2谱峰度值与孔隙胶结指数（m）值相关性较高, 进一步确定了岩电参数m的核磁计算公式.最终, 将该套方法应用于研究区井筒剖面中, 有效地提高了饱和度计算精度, 也为东营凹陷南坡低渗透砂岩油藏储量估算与高效开发提供了依据.     

Silicate perovskite analogue ScAlO3: temperature dependence of elastic moduli

The elastic moduli of ScAlO₃ perovskite, a very close structural analogue for MgSiO₃ perovskite, have been measured between 300 and 600 K using high precision ultrasonic interferometry in an internally heated gas-charged pressure vessel. This new capability for high temperature measurement of elastic wave speeds has been demonstrated on polycrystalline alumina. The temperature derivatives of elastic moduli of Al₂O₃ measured in this study agree within 15% with expectations based on published single-crystal data. For ScAlO₃ perovskite, the value of (∂K_S/∂T)_P is −0.033 GPa K⁻¹ and (∂G/∂T)_P is −0.015 GPa K⁻¹. The relative magnitudes of these derivatives agree with the observation in Duffy and Anderson [Duffy, T.S., Anderson, D.L., 1989. Seismic velocities in mantle minerals and the mineralogy of the upper mantle. J. Geophys. Res. 94, 1895–1912.] that |(∂K_S/∂T)_P| is typically about twice |(∂G/∂T)_P|. The value of (∂K_S/∂T)_P for ScAlO₃ is intermediate between those inferred less directly from V(P,T) studies of Fe-free and Fe- and Al-bearing MgSiO₃ perovskites [Wang, Y., Weidner, D.J., Liebermann, R.C., Zhao, Y., 1994. P–V–T equation of state of (Mg,Fe)SiO₃ perovskite: constraints on composition of the lower mantle. Phys. Earth Planet. Inter. 83, 13–40; Mao, H.K., Hemley, R.J., Shu, J., Chen, L., Jephcoat, A.P., Wu, Y., Bassett, W.A., 1991. Effect of pressure, temperature and composition on the lattice parameters and density of (Mg,Fe) SiO₃ perovskite to 30 GPa. J. Geophys. Res. 91, 8069–8079; Zhang, Weidner, D., 1999. Thermal equation of state of aluminum-enriched silicate perovskite. Science 284, 782–784]. The value of |(∂G/∂T)|_P for ScAlO₃ is similar to those of most other mantle silicate phases but lower than the recent determination for MgSiO₃ perovskite [Sinelnikov, Y., Chen, G., Neuville, D.R., Vaughan, M.T., Liebermann, R.C., 1998. Ultrasonic shear wave velocities of MgSiO₃ perovskite at 8 GPa and 800K and lower mantle composition. Science 281, 677–679].

Combining the results from the previous studies and current measurements on ScAlO₃ perovskite, we extracted the parameters (q and γ₀) needed to fully specify its Mie–Grüneisen–Debye equation-of-state. In this study, we have demonstrated that acoustic measurements of K_S(T), unlike V(P,T) data, tightly constrain the value of q. It is concluded that ScAlO₃ has ‘normal’ γ₀ (1.3) and high q (3.6). The high value of q indicates that ScAlO₃ has very strong intrinsic temperature dependence of the bulk modulus; similar behaviour has been observed in measurements on Fe- and Al-bearing silicate perovskites (Mao et al., 1991; Zhang and Weidner, 1999).     

Liquefaction in silty soils—screening and remediation issues

Current techniques for liquefaction screening, ground modification for liquefaction mitigation, and post-improvement verification rely on knowledge gained from extensive research on clean sands, field observations of liquefied ground, and judicial correlation of normalized penetration resistance [(N₁)₆₀,q_c1N] or shear wave velocity (v_s1) data with field liquefaction observations. Uncertainties prevail on the direct extrapolation of such techniques for silty soil sites. This paper examines laboratory data on liquefaction resistance, strength, and v_s1 of sands and silty soils using grain contact density as the basis. Effect of silt content on cyclic resistance, strength, m_v, and c_v is examined in this light. Rational insights on effects of silt content on the current screening techniques based on (N₁)₆₀, q_c1N, and v_s1 to silty soils are offered. Recent advances and modifications to the traditional densification, drainage, and permeation grouting techniques to make them viable for silty soils are discussed.     

Analysis of the time-scaling behaviour in the sequence of the aftershocks of the Bovec (Slovenia) April 12, 1998 earthquake

Analysis of the temporal fluctuations, described by a power spectral density S(f) decaying as 1/f with exponent >0, in the sequence of aftershocks of the Bovec (Slovenia) earthquake of April 12, 1998 has been performed, using the Allan Factor (AF) method. The sequence of the occurrence times of the events with threshold magnitude M_th≥2.0 is characterised by a scale-invariant behaviour from the time scale τ2·10⁴ s with a scaling coefficient 0.9, evaluated by a least-square method. By gradually increasing the threshold magnitude up to M_th=2.9, the AF curves, associated, respectively with the processes of selected events with magnitude M≥M_th, indicate a monotonic decrease of the value of the scaling exponent , in the range of time scales stretching from 10² to 5·10⁵ s. For all the thresholds considered, the AF curves increase monotonically with the counting time and are well fitted by an increasing power-law function for counting times greater than approximately 2·10⁴ s. This monotonic power-law increase indicates the presence of fluctuations on many time scales and therefore of fractal clustering.     

易县地震台分量钻孔应变观测数据分析

选取易县地震台YRY-4型分量式钻孔应变仪2012—2020年观测数据,进行Venedikov调和分析,并对2015—2020年四分量数据S₁、S₂、S₃、S₄进行自洽检验,发现应变观测数据连续性、可靠性较高。选取该台2018—2019年钻孔应变原始数据S₁、S₂、S₃、S₄,利用面应变S₁ + S₃、S₂ + S₄相关系数,分析河北省3次M_L≥4.0地震前数据异常变化,并利用剪应变S₁ - S₃、S₂ - S₄的大小,分析受抽水干扰的NE分量主应力方位角,判断干扰源方位,结果发现：①面应变与剪应变在地震发生前均会出现异常变化,且剪应变异常变化幅度较大,应与震中相对于观测井的位置有关;②对于受抽水干扰的钻孔应变NE分量,利用剪应变固体潮数据推算主应力方位角,可得干扰源在钻孔应变观测室NW80.96°方位,与实际位置存在一定偏差。     

Dynamic properties and liquefaction resistance of two soil materials in an earthfill dam—Laboratory test results

Results are presented of laboratory resonant column and cyclic triaxial tests on specimens of two compacted soils (a sandy–silty clay and a sand–gravel mixture), planned to be used in the core and the shells, respectively, of a proposed earthfill dam. The values of low-amplitude shear modulus of the clayey material were found to increase with increasing confining pressure and decreasing water content, with deviations of ±20% from the predictions of the “Hardin equation”. On the other hand, the low-amplitude damping ratio was found to be at least four times higher than the values corresponding to natural undisturbed cohesive soils. The proposed G/G₀−γ_c curve for the compacted cohesive soil was found to be independent of confining pressure and small variations of the water content on either side of the optimum value and showed a remarkable agreement with recently published similar curves for natural cohesive soils. For the case of the sand–gravel mixture, normalized G/G₀–γ_c and D–γ_c curves are proposed, based on recently published results for gravelly soils and the limited data of the present study. The liquefaction resistance of the saturated sand–gravel mixture was found to be strongly dependent on its relative density, especially for high values of cyclic stress ratio. A relative density of at least 55% was found to be necessary to assure safety against earthquake-induced liquefaction of the material. The results presented herein may be used (directly or as guide) in the seismic analysis of (new or existing) earth dams constructed from similar soil materials and in addition they provide insight into the dynamic behavior of compacted soils.     

中国东北及邻区竖直向地震动强度包络函数参数的确定及其衰减关系

选用2008—2019年中国东北及邻区27次地震竖直向105条强震动记录,以三段式强度包络函数模型为目标,利用加权最小二乘回归方法,尝试建立了东北及邻区竖直向地震动时程强度包络函数参数的衰减公式,以使我国的地震动区域划分更具有普遍性。同以往的研究类似,上升段持时t₁和平稳段持时t_s的确定方法采用了70%能量持时法。所选地震的震级为3.3—5.8级,最大震中距为371 km。基于土层场地的竖直向记录,通过与我国目前常用的时程强度包络函数参数的衰减关系进行对比分析发现,t₁与震级之间有明显的相关性,实际应用中应考虑震级对t₁的影响;地震整体持时较短造成t_s上升趋势随距离的增加而变缓;下降段由衰减系数c的变化可见,记录整体衰减较快。     

Understanding the effect of a static driving shear stress on the liquefaction resistance of medium dense granular soils

The influence of a driving static shear stress on the liquefaction resistance of medium dense granular soils is investigated. A laboratory study of the behavior of five sands (percentage of silt varying between 2 and 42%, and D₅₀ varying between 0.55 and 0.09 mm) was undertaken. These five sands were sampled from the foundations of two dams in Quebec (Canada). Sixty cyclic direct simple shear (DSS) tests were performed with and without a driving static shear stress. Relative densities D_r varied between 60 and 70%, and static shear stress ratios =τ_st/σ′_vc varied between 0.2 and 0.6. The presence of static shear stress increased the cyclic strength of these soils. A relation between and K() factor giving the effect of preshearing on cyclic shear resistance is obtained. The mechanical competence is examined as a function of D₅₀ and percentage of fines. The observed dilatancy yielded rather cyclic mobility than liquefaction.     

基于等效源的总强度磁异常非线性处理方法

总强度磁异常（ΔT₁）常规处理方法通常将其近似当作磁异常矢量在地磁正常场方向上的投影（ΔT₂）.然而对于高磁环境如磁铁矿处的磁异常场幅值可达10⁴ nT甚至更大的情况,上述近似假设则不再成立,如采用常规处理方法可能会带来明显的误差.本文针对该问题,提出一种基于等效源的总强度磁异常非线性处理方法;该方法根据总强度磁异常获取流程,直接反演实测地磁场总强度幅值与正常场强度幅值之差来求取等效场源.本文首先通过模型试验,分析ΔT₁与ΔT₂的差异;然后将ΔT₁当作ΔT₂采用常规处理方法以及ΔT₁采用新方法得到的结果作对比分析,结果表明新方法处理结果与理论值的差异为常规方法处理结果的1/5甚至更小,充分说明了新方法的有效性;最后将该方法应用于铁矿实测总强度磁异常处理实例中来转换计算磁异常总模量,其实际应用效果进一步体现了高幅值总强度磁异常数据处理过程中采用新方法的必要性.

     

川滇地块西部差异性旋转的构造意义:青藏高原东南缘白垩纪红层古地磁学新证据

印度-欧亚板块碰撞以来青藏高原内部及其周缘地区经历了复杂的构造演化,复杂构造变形区的复合构造使得古地磁的数据解释究竟代表区域的构造旋转还是只能反映局部的构造变形一直是备受关注的问题.本文通过采集川滇地块西缘渔泡江断裂东侧三岔河地区白垩纪红层古地磁样品,揭示采样区差异性旋转并探讨川滇地块西部自中新世以来的构造演化规律.前人的地质调查表明川滇地块渔泡江断裂东侧上白垩统赵家店组地层发育倾伏褶皱.三岔河剖面以三岔河镇为界分为南北两段,三岔河南段剖面高温剩磁分量平均方向在倾斜校正后D_s=29.3°,I_s=45.7°,k_s=54.3,α₉₅=6.6°,倾伏地层产状校正后D_s=30.6°,I_s=46.6°,k_s=69.3,α₉₅=5.8°;而三岔河北侧剖面高温剩磁分量平均方向在倾斜校正后D_s=350.4°,I_s=42.1°,k_s=69.4,α₉₅=9.2°,倾伏地层产状校正后D_s=347.4°,I_s=41.9°,k_s=96.6,α₉₅=7.8°;两组高温剩磁分量均通过了褶皱检验,表明其获得于褶皱形成之前.相对于东亚稳定区80 Ma古地磁极,三岔河南侧剖面发生了20.5°±4.8°的顺时针构造旋转量,与楚雄盆地核部之间不存在差异性旋转;但三岔河镇以北剖面却发生了22.7°±6.6°的逆时针旋转.综合分析川滇地块内部的古地磁数据表明自中新世以来川滇地块南部楚雄盆地经历了约20°的顺时针构造旋转,而三岔河镇北侧经历了约20°逆时针旋转.进一步分析表明三岔河北侧剖面相对于南侧剖面经历了约40°的逆时针旋转,可能由于研究区的滑脱构造导致岩石薄弱层拆离滑脱所引起.

     

基于IDA方法的交河故城金刚宝座式塔土建筑遗址地震易损性分析

增量动力分析（IDA）作为动力形式的静力推覆法,被广泛应用于建筑结构的抗震能力分析和性能评估中,但其在土建筑遗址地震易损性分析中的应用甚少。基于IDA方法尝试对交河故城塔林区金刚宝座式塔土建筑遗址进行在地震作用下的易损性和抗震能力分析,以谱加速度S_a（T₁,5%）作为地震动强度指标,最大层间位移角θ_max作为结构损伤指标进行IDA分析,得到土建筑遗址在一系列调幅地震动作用下的IDA曲线、地震需求概率模型和易损性曲线,计算出模型的倒塌率和CMR值。通过以上分析得出金刚宝座式塔的易损性和抗震能力,为今后土建筑遗址抗震保护提供理论基础和分析思路。     

The equation of state of CaSiO3 perovskite to 108 GPa at 300 K

Pressure–volume measurements have been performed for CaSiO₃ perovskite to 108 GPa at 300 K using NaCl and argon pressure-transmitting media, and energy dispersive X-ray diffraction (EDXD) in a diamond-anvil cell (DAC). By determining a parameter that is the product of the elastic anisotropy, S, and the uniaxial stress component, t, for each data point, we define the stress condition of the sample. For different points at the same pressure in a temperature-quenched sample, the St value can differ by as much as a factor of 5, indicating heterogeneity in the stress condition. This may be responsible for the large scatter of earlier P–V measurements in the DAC which in general used a large diameter X-ray beam. Also, the St value provides insight into the elastic anisotropy, S, of CaSiO₃ perovskite and platinum. The sign of S (positive) for CaSiO₃ perovskite agrees with first principles calculations but the magnitude may be inconsistent. A new compression curve at 300 K was obtained for CaSiO₃ perovskite by using those data points which represent the most nearly hydrostatic conditions. It is observed that the data points with high St values yield larger volumes than the points with small St values at a given pressure. By selecting the data points having low St values (St≤0.005), combining with lower pressure large volume press (LVP) measurements and fitting to third order Birch–Murnaghan equation of state (EOS), we find that CaSiO₃ perovskite is more compressible (V₀=45.58±0.05 ų, K_T0=236±4 GPa, and K_T0′=3.9±0.2 GPa) than suggested by previous studies. The density and bulk modulus of CaSiO₃ perovskite at lower mantle pressures and 300 K are 1–3% greater and 5–15% smaller, respectively, than found in previous studies. This study demonstrates that defining the stress state of the sample is crucial to obtain an accurate 300 K compression curve for unquenchable high-pressure phases.     

Apparent stress scaling and statistical trends

Statistical relations between different earthquake parameters, such as M₀ (seismic moment), E_S (seismic energy), τ_a (apparent stress), A (rupture area), g (average slip acceleration), are investigated. For this purpose, a kinematic earthquake model representing averaged earthquake rupture process is formulated. The model implies a scaling relationship for τ_a as a function of three other parameters, related to kinematic (M₀), geometric (A) and material (g) source characteristics, which, according to the model, can change independently. This scaling relation is used to explain statistical trends that characterize different earthquake data sets (including micro-, small, moderate and large events) plotted in the log τ_a − log M₀ space, and to determine the area in this space, where typical earthquakes occur. The scaling relationship is interpreted in terms of the apparent stress minimum (i.e., the most uniform among the possible earthquake rupture patterns). It is concluded that, although the apparent stress increases on an average with increasing seismic moment, small and large earthquakes are essentially similar.     

岩石润湿性的核磁共振表征方法与初步实验结果

随着石油勘探开发领域的不断扩大,某些复杂油气藏尤其是致密油气的岩石物理响应,受岩石润湿性变化的影响越来越明显,影响油气藏的评价以及后续勘探开发,识别和评价岩石润湿性成为重要的研究课题.核磁共振弛豫可以有效表征多孔介质孔隙及其流体分布,特别是通过T₁-T₂弛豫图谱能够获得流体分子与孔隙表面相互作用的信息,进而为孔隙表面润湿性质的判断提供了可能性.本文首先讨论了T₁-T₂弛豫谱表征孔隙介质润湿强度的理论基础,然后,采用一阶近似方法对孔隙表面流体分子的核磁共振弛豫分布进行正演模拟,建立T₁/T₂比值法判别孔隙介质润湿性的基本方法.最后,选取玻璃珠和天然砂岩样品进行初步实验,验证方法的合理性.结果表明,当岩石饱和双相流体时,润湿相流体T₁/T₂值高于非润湿相流体,通过观察不同饱和度状态T₁-T₂二维图谱的变化,能够判断岩石的润湿性.本文的理论分析与初步实验结果为岩石润湿性研究提供了新的手段,对核磁共振测井探测储层润湿性也具有一定的应用价值.

     

Non-Darcian flow in a single confined vertical fracture toward a well

Non-Darcian flow towards a well which fully penetrates a confined single vertical fracture is presented in this paper on the basis of the Izbash equation. We have obtained semi-analytical solutions for non-Darcian flow by using the Boltzmann transform and developed the non-Darcian flow well functions for cases with and without the effect of wellbore storage. The results show that the non-Darcian flow type curves are more or less deviated from the Darcian flow type curve. The non-linear effect is mainly attributable to the turbulent factor, v, a dimensionless parameter related to the pumping rate, the fracture aperture, the fracture thickness, and two constants k′ and n used in the Izbash power-law. The non-linear effect appears to be less sensitive to the power-law index, n. When excluding wellbore storage, the well function at early times is proportional to v^−1/(n−1)u^−n/(n−1), where u is a dimensionless term inversely proportional to time; whereas the well function at late times is approximated as , where A₀(n) is a finite term depending on n. When considering wellbore storage, drawdowns inside the well with different v values approach the same asymptotic value at small times, and the effect of wellbore storage is only found at the early stage of pumping.     

内蒙古东部地区一维地壳速度结构分析

采用内蒙古测震台网2009—2016年记录的内蒙古东部地区131个地震资料,使用速度拟合、分区扫面、折合走时方法,反演得到该区域速度模型：v₁=6.10 km/s、v_Pb=6.72 km/s、v_n=8.05 km/s、H₁=23 km和H₂=16 km。东部模型检验结果显示,定位残差均值较华南模型和2015内蒙最优模型有明显的降低,且更加均匀稳定;东部模型与编目定位震中差较华南与编目、2015内蒙最优模型与编目和华南与编目震中差均值降低1 km左右。可见,东部模型更适合内蒙古东部地区。     

四川盆地早白垩世红层磁倾角偏低研究

对四川盆地东北部巴中地区和西部雅安地区早白垩世红层分别采集9块手标本,进行了详细的古地磁学研究.系统退磁获得两地的特征剩磁,其中巴中地区平均方向为,倾斜校正之前：D_g/I_g=25.3°/19.0°,k=18.6,α₉₅=8.8°;倾斜校正之后：D_s/I_s=25.8°/18.9°,k=24.3,α₉₅=7.6°.雅安地区平均方向为,倾斜校正之前：D_g/I_g=24.5°/45.0°,k=15.7,α₉₅=9.0°;倾斜校正之后：D_s/I_s=356.7°/35.6°,k=28.5,α₉₅=6.6°.磁化率各向异性实验结果显示两地均未受到显著的构造应力影响.等温剩磁各向异性（AIR）实验结果,巴中地区IRMz/IRMx平均值为0.8194,表明18%的压实率;雅安地区IRMz/IRMx平均值为0.8909,表明11%的压实率,计算得到巴中和雅安地区校正后的磁倾角分别为22.7°和38.8°.根据等温剩磁各向异性（AIR）实验和Tauxe and Kent（2004）提出的EI校正法得到的结果表明,四川盆地早白垩世陆相碎屑岩层中存在沉积压实作用造成的磁倾角偏低现象,而且川东北巴中地区的偏低程度强于川西雅安地区.

     

Physical basis of earthquake magnitudes: an extreme value of seismic amplitudes from incoherent fracture of random fault patches

A relation between surface-wave magnitude M_s and fault area S for great earthquakes has been proposed

log S 2M_{s − 11.5 (}M_{s > 7.5)}

A similar formula has been also derived for body-wave magnitude m_b^* redetermined from maximum amplitudes of short-period P-wavetrains

log S ∝ 1.7 m_b^*

These are quite different from a theoretical relation expected on the basis of long-wave approximation. Because wavelengths of seismic waves used for the above magnitudes are very short compared to the size of earthquake sources, these relations represent the short-period nature of the earthquake process. The statistical theory of extreme values has been applied to understand the relations considering that the component waves which constitute the wavetrains for m_b^* and M_s determinations originate from the random fracture of fault heterogeneities.