永登5.8级地震极震区黄土微结构孔隙性的定量研究

通过电镜扫描手段和计算机图像处理,对1995年永登5．8级地震极震区遭受地震影响和未遭受地震影响的黄土微结构变化的各项指标进行了对比试验和分析,并将分析结果绘成受震与未受震两种孔隙分布曲线。从这些曲线可以看出,不同地点的黄土孔隙分布特征不同;同一深度不同地点的孔隙分布亦有明显差异。对未受震(粒状架空孔隙结构)与受震(粒状镶嵌孔隙结构)黄土的微观结构进行了定量处理,其分析结果较好地解释了土体的工程性质。同时介绍了黄土微结构研究的方法和要求,通过对黄土微观结构的定量处理,揭示了黄土震陷性与黄土孔隙含量所占比例大小的内在关系,达到了确定黄土场地震陷量的目的,验证了利用微结构研究黄土结构变化的可行性。     

强震作用下抗震陷黄土改良地基的微观特征分析

基于动三轴试验和SEM细观结构测试试验,结合图像分析处理软件,对强震荷载前后宝兰客运专线沿线典型震陷性黄土及其经物理、化学和复合改良方法处理后的试样进行微观尺度的结构变形演化规律研究,探讨不同改良方法对土体微结构的影响及其与震陷系数之间的内在联系。结果表明:(1)物理改良方法对于消除大孔隙和架空孔隙结构的效果最为明显,同时对颗粒级配和结构也有调整;(2)化学改性方法则从颗粒接触方式、粒间胶结程度等方面影响土体强度,且不同化学反应的参与可生成独特的玻璃微珠或絮凝状细结构,从而大大提升土体某项参数指标,进而在强震中分别起到填充、胶结或缓冲作用;(3)强震作用前后不同改良方法对各微观要素的改变与相应的残余应变现象吻合较好,说明微结构能有效反映改良黄土残余变形的强弱。     

我国黄土的微结构类型与震陷区域划分

对比我国中西部成土年代相近的浅层风成黄土（不考虑地形地貌）,针对不同成土环境下微观结构及其震陷性大小,初步将中西部黄土的微结构类型划分为五类：Ⅰ,近砂源快速降砂微结构;Ⅱ,强降尘弱成壤微结构;Ⅲ,冷干慢速降尘弱成壤微结构;Ⅳ,中湿成壤微结构;Ⅵ,温湿成壤微结构.微结构类型对于震陷性等工程性质分析有借鉴作用,具体体现在颗粒大小、粒径分布、排列方式,颗粒之间的胶结程度、接触方式等.通过分析黄土地区气候变化引起的地域性微观结构类型差异,区域性微观结构的气候形成机制来研究黄土的致灾特征,进一步获得黄土地区土层微结构的地域性差别导致的震陷变形强弱信息.这对于黄土建筑工程场地的设计以及提出针对性的防护措施有重要意义.     

黄土的震陷性与其微结构特征的关系研究

借助电镜扫描分析手段，对Q3震陷性黄土的微结构特征、颗粒形态等进行了分析。结果表明黄土的显微结构对震陷系数有较大影响，并与震陷的临界动应力关系密切。通过孔隙含量的定量测试分析，探讨了中大孔隙对黄土震陷影响，并对土样的微结构某些指标进行了比较和分析。     

黄土场地震陷的灾害特征及成因研究综述

基于文献史料记载和现场科考成果,考证了黄土地区几次中强以上地震引发的黄土震陷灾害,分析了黄土震陷灾害的破坏模式、发育特征及其形成条件。依据大量不同地区黄土的动三轴试验结果,比较分析了黄土震陷形成的物性指标的区域性变化规律,并进一步研究了震陷临界动应力的影响因素。研究显示,黄土地区的几次中强以上地震的极震区均有黄土震陷事例;黄土震陷的灾害模式可区分为振密型震陷、震陷型滑坡和液化型震陷,并由土性条件、地形地貌和地震动作用形式决定;孔隙比大于0.8时,黄土高原地区的黄土需要考虑不同含水量条件下的震陷性,而黄土高原地区东南部的河南、河北和晋东南地区的黄土则可根据缩限含水量判定其震陷性;临界动应力的大小由黄土的结构强度所主导,临界动应力比与含水量之间存在很好的幂函数关系。     

临洮黄土土层微观差异及动力学指标分析

取甘肃临洮黄土不同深度土层4~16 m样品,分别研究微观结构类型差异及震陷性、剪切波速和元素比值等指标的关系。由于不同深度土层样品形成年代不同,气候条件变化很大,同一场地不同深度剖面所形成的黄土微观结构也存在较大差异,各深度(4 m,8 m,12 m,16 m)土层样品的微结构类型依次为弱剪切结构,易压密结构,密实结构,中强胶结结构;同时震陷系数依次降低,剪切波速依次增大,其中的弱剪切结构的剪切波速仅为142 m/s。研究表明通过土层样品的微观结构类型特征及指标参数分析,可较为全面地判断土层场地,找到震害分析和土层处理的有效指标和参数,并据此进行合理的工程设计,可为黄土地区震害的成因和预防提供理论依据。     

含水量对黄土震陷性定量影响研究

含水量是影响黄土震陷性非常重要的参数。本文以洮河三级阶地上典型的Q3黄土为研究对象,对其震陷性特征进行了全面、系统的研究。通过在20 m深探井中依次取得4 m、8 m、12 m、16 m、20 m处黄土原状试样进行室内试验,得到以下结果:(1)天然黄土存在震陷性随深度增加而逐渐减弱的规律;(2)通过增湿、减湿获得5%,10%,15%和20%四种不同含水量黄土试样,同一深度黄土随着含水量增加其震陷性增长的规律明显;(3)通过实验数据分析得到了含水量对黄土震陷性定量影响规律,并根据其规律给出了不同动应力作用下的黄土震陷临界含水量。     

不同加载方式下黄土累积变形特征研究

土体动力特性主要与其受力状态、土体加载类型、试验控制方式以及荷载等参数有关。以往非饱和黄土累积变形试验研究通过单向动应力加载方式开展相关研究,并不能真实反映土体受力状态。现运用WF-12440型动三轴-空心圆柱扭剪试验系统,通过单向荷载和双向荷载两种不同的动应力加载方式,进行黄土的长期重复载荷试验,对比不同动应力加载方式对非饱和黄土累积变形发展特征的影响。试验结果发现：黄土在单向荷载或双向荷载作用下,其软化指数均随着循环次数的增加,呈非线性减小的趋势,在相同循环振次时,土体软化指数随着动应力幅值的增大而减小;黄土在不同加载方式作用下,软化指数减小速率在循环振次100次以内较快,后期随振次增加,土体软化速率逐渐趋于平稳;黄土在双向荷载作用下,当轴向动荷载较小时,径向荷载的施加,将加速黄土的软化程度,随着轴向动荷载的增大,径向荷载对黄土软化的影响逐渐减弱,其对于黄土变形的影响存在一临界轴向循环动应力。本研究的开展,可更准确地进行不同荷载耦合变化下,黄土场地震陷评估或路基土体变形计算,实现有效控制路基的整体稳定和工后沉降量,以期根据路基土参数预测后期沉降,极具重要的科研和工程意义。     

黄土场地震陷量的试验预测

通过不同地震荷载作用下黄土震陷的试验研究对兰州地区几个黄土场地的震陷量进行了预测．发现黄土在随机地震荷载和等幅正弦荷载作用下的震陷量预测结果存在一定的差别     

黄土的湿陷和震陷

黄土的湿陷性和震陷性分别是黄土在遇水和地震作用下产生的岩土灾害。黄土湿陷性评价及防治已经相对成熟。有些专家试图联系二者,以便解决黄土震陷性的工程评价和处理。本研究中,作者分析了黄土湿陷性与震陷性之间的联系与差异。研究认为,二者的产生都基于黄土的亚稳态结构,都是黄土架空结构破坏所引起,但其触发机制却是不同的。湿陷性产生的主要原因是吸力丧失,而震陷性主要由黄土在动荷载作用下的土体剪切破坏所导致。因此,关于二者之间的联系还不存在坚实的物理基础。     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Solute concentrations in the lower niger river and the source rock contribution

Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO₃ > SO₄ > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1^?1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Late Permian post‐ophiolitic trondhjemites from Central Iran: a mark of subduction role in growth of Paleozoic continental crust

Late Permian trondhjemites in the Anarak area occur as stocks and dykes, which cross cut the Anarak ophiolite and its overlying metasedimentary rocks, and are exposed along the northern Anarak east–west main faults. These leucocratic intrusive bodies have enclaves of all ophiolitic units and metamorphic rocks. They are composed of amphibole, plagioclase (oligoclase), quartz, zircon and muscovite. Secondary minerals are chlorite (pycnochlorite), epidote, albite, magnetite and calcite. Whole‐rock major‐ and trace‐element analyses reveal that they are characterized by high SiO ₂ (67.8–71.0 wt%), Al ₂ O ₃ (14.9–17.1 wt%) and Na ₂ O (5.3–8.6 wt%), low K ₂ O (0.1–1.5 wt%; average: 0.8 wt%), low Rb/Sr ratio (0.01–0.40; average: 0.09), low Y (3–6 ppm), negative Ti, Nb and Ta anomalies, slightly negative or positive Eu anomaly, LREE enrichment and fractionated HREE. These rocks present 2 to 40 times enrichment in inclined chondrite‐normalized REE patterns. Geochemical characteristics of the Anarak trondhjemites all reflect melting of a mafic protolith at more than 10 kbar. The field evidence and whole‐rock chemistry reveal that these rocks have been crystallized from magmas derived from melting of subducted Anarak oceanic crust. This study reveals that melting of garnet amphibolite was an important element of continent formation in the study area.