印度板块岩石圈地幔向北俯冲到羌塘地体之下的远震P波层析成像证据

利用PASSCAL、INDEPTH Ⅱ、INDEPTH Ⅲ、HIMNT等研究计划,及中国新疆地学断面和国家973项目在青藏高原布设的流动台站记录的到时数据,以及自1990年1月到2004年2月全球地震事件的震相报告,作者对覆盖印度大陆的恒河平原和整个青藏高原的305个地震台站记录的9649个远震事件,共139021条P波初至到时资料进行了层析成像反演．结果表明：印度岩石圈地幔在不同的位置向北俯冲的形态不同,但俯冲前缘都到达羌塘地体之下．沿88°E剖面显示,厚约100 km的印度岩石圈地幔从南部的恒河平原向北一直俯冲到青藏高原之下．在主边界逆冲断裂之下100 km深度处以约22°角度开始向北俯冲,俯冲最前缘到达羌塘地体的中部地区约34°N,之后进入上地幔深处．而沿北东方向的剖面则显示,印度岩石圈地幔以近水平的角度俯冲到青藏高原之下,向北越过班公－怒江缝合带,到达33°N附近,然后以大角度近乎垂直地向下俯冲断离,并引起地幔热物质的上涌,形成羌塘地体之下大规模的低速带．     

1992年美国加州兰德斯地震——地壳结构不均匀性对地震发生的影响

为了理解地震的发生和地壳结构不均匀性的关系,利用南加州地震台网的固定台和临时台所记录的2863个兰德斯余震和区域地震,共计107401个P波和19624个S波高质量的到时数据,采用地震层析成像方法得到了兰德斯地震区P波和S波的精细的三维速度结构和泊松比分布．结果显示,地震的发生和分布与地壳结构的横向不均匀性有密切的关系．总体上看,兰德斯地震区余震成丛分布,并被低速块体截断,其中4级以上地震大多分布于P波高、低速异常过渡区域或偏向高速块体一侧,这可能是因为高速区多属地壳脆性介质,易于造成应力集中,导致地震;反之,低速度区则可能代表破碎程度较高、富含流体或温度较高区域,因而更倾向于产生无震变形．基于兰德斯地震区强震震源位置、地震区P波、S波速度异常与泊松比分布推断,兰德斯地震区可能有流体存在．地壳流体易使地壳岩石弱化,从而引发大地震．     

使用纵波和横波走时层析成像研究菲律宾海板块西边缘带和南海地区的深部结构

使用纵波震相及同台同源纵波和横波震相的走时资料及层析成像反演方法,分别给出了菲律宾海板块(PHP)和南海地区的纵波速度V_pf及同台同源的纵波和横波速度V_p和V_s结构. 结果表明,（1） PHP与欧亚板块(EUAP)的俯冲接触关系随地段而异,在琉球海沟,PHP向EUAP之下俯冲深达400 km;在台湾岛,EUAP先近陡直地俯冲到深度400 km,然后折向PHP之下达到660 km左右;在马尼拉海沟北段,俯冲板片几乎垂直达到660 km附近;在菲律宾海沟,PHP向EUAP之下的俯冲深度超过660 km．（2）南海地区之下是一个深达300～400 km的宽阔低速区,并且大致在莺琼海盆与700 km深处另一低速区曲折相通;在该宽阔低速区内部,有两个小而明显的低速区分别在海口火山和珠江口盆地下方．（3）对V_p和V_s及据其算出的容声速度V_b作分析发现,V_s和V_b的平均扰动量对深度的变化在一些深度范围内是反向的;年龄较大的太平洋板块俯冲体的V_s相对扰动量的平均值大于V_b的,而在较年轻的PHP俯冲体中则正好相反．     

中国大陆中上地壳剪切波速结构

为使已获得的中国大陆中上地壳结构更为可靠,本文搜集了很多对浅部结构分辨率较好的短周期面波资料. 与传统面波层析成像反演方法不同,本文在第二步由面波频散得到剪切波速的过程中不再对每个结点进行一维波速模型分别反演,而是直接将所有结点上的区域化频散转换成三维波速结构的线性化约束,实现了直接的三维反演. 检测板测试结果显示该方法可以得到理想的反演结果. 本文得到的波速模型显示,中国大陆中上地壳的速度分布存在明显的横向变化和分区特征. 较低的波速异常很好地勾勒出我国主要的沉积盆地,波速异常在不同深度上的变化在一定程度上反映了各盆地结晶基底的深度. 以东经95deg;为界,特提斯构造域西部具有明显的低速异常,而东部基本没有低速异常. 基于油气资源多存在于沉积层中,而沉积层表现低速异常,我们推测特提斯构造域西部油气前景比东部好. 另外, 由于特提斯构造域西部低速非常明显,这可能也说明了其地壳温度较高. 兴安造山带的低速异常可能也说明了其地壳温度较高. 20 km深度上鄂尔多斯盆地西侧的弱低速带, 很好地勾勒出中国大陆近似沿105deg;经线的强震带的走势.      

封闭系统中多孔介质甲烷水合物的CT实验研究

阐述了在封闭系统中的查理模型,建立了利用查理定律计算分析甲烷水合物生成和分解的方法。通过实验,对封闭系统中不同含水量多孔介质甲烷水合物生成、发育和分解的实验方法、相平衡条件、生成量关系等进行了探讨。在一定温度差的介质里,水分克服重力向冷端迁移并生成水合物。实验中利用查理数的改变计算水合物的生成量变化,连通性好的含水多孔介质有很好的水—气结合条件,达到平衡条件后迅速生成水合物,在240分钟内完成,含水量线性地影响甲烷水合物的生成量。同时,利用CT扫描方法直接观测含水粗砂在低温—高压环境下的甲烷水合物生成和分解过程,并可以利用CT图像数据计算出多孔介质中水分迁移、区域密度改变和砂体的移动、及分解后介质密度分布变化特征。     

鄂尔多斯地块周缘三维地壳速度结构研究

鄂尔多斯地块及周缘位于青藏高原与华北地块之间,其内部无大断裂分布,而周缘却由活跃的地震带——活动断陷带和弧形断裂带包围.因此,鄂尔多斯地块内部无6级以上大地震分布,而周缘地震活动则频度高、强度大(图1),历史上发生多次破坏性大地震,如:1654 年天水南8.0 级、1879 年武都南8.0级、1920 年宁夏海原81/...     

Anticipating the effects of groundwater withdrawal on seawater intrusion and soil settlement in urban coastal areas

Intensive pumping in urban coastal areas is a common threat to water resource quality due to seawater intrusion. In those areas where subsurface water resources are not usually used for human consumption or irrigation, intensive pumping is associated with other activities like the lowering of the water table necessary to support underground structures and building foundations. This activity also increases the likelihood of soil settlement that affects building stability and the corrosion of concrete structures due to groundwater salinity. Under these circumstances, the awareness of a certain municipality (Calonge, NE Spain) of the potential effects of groundwater withdrawal upon foundations has led to an integrated approach to anticipate seawater intrusion related to urban development. Geological mapping and correlation of borehole logs, electrical resistivity tomography, and hydrochemical data provide comprehensive knowledge of the geology and hydrogeology of the area and act as screening tools necessary to discern the influence of hydrological processes in coastal areas. Developing Strack's analytical solution, new comprehensive, dimensionless expressions are herein derived to determine the critical pumping rate necessary to prevent seawater intrusion, as well as to reproduce the evolution of the wedge toe and the water table stagnation point under different withdrawal rates. Furthermore, the Dupuit–Forchheimer well discharge formula allows the estimation of the effects of the water table lowering due to such critical pumping in the surrounding building foundations. Field data from the Calonge coastal plain illustrate this approach and provide assessment criteria for future urban development and planning. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

Integration of remote sensing data into hydrological models for reservoir management

Abstract

The purpose of this paper is to present the methodology set up to derive catchment soil moisture from Earth Observation (EO) data using microwave spaceborne Synthetic Aperture Radar (SAR) images from ERS satellites and to study the improvements brought about by an assimilation of this information into hydrological models. The methodology used to derive EO data is based on the appropriate selection of land cover types for which the radar signal is mainly sensitive to soil moisture variations. Then a hydrological model is chosen, which can take advantage of the new information brought by remote sensing. The assimilation of soil moisture deduced from EO data into hydrological models is based principally on model parameter updating. The main assumption of this method is that the better the model simulates the current hydrological system, the better the following forecast will be. Another methodology used is a sequential one based on Kalman filtering. These methods have been put forward for use in the European AIMWATER project on the Seine catchment upstream of Paris (France) where dams are operated to alleviate floods in the Paris area.     

Inversion of multiple thermal indicators: Quantitative methods of determining paleoheat flux and geological parameters. II. Theoretical development for chemical, physical, and geological parameters

Using the logic for quantitative inversion of present-day downhole thermal indicators, the inversion procedure can be used to determine unknown, or poorly known, chemical and physical parameters as well as other geological quantities of interest which impact on burial history and thermal history of an evolving sedimentary basin. Some such quantities are: amount of erosion and timing of unconformities, paleo-overpressuring, stratigraphic age, timing of igneous intrusion and insertion temperature, overthrust timing and frictional heating, fault and slump timing, effects due to emplacement of a radiocative layer, and salt emplacement and dissolution timing. Combining a priori unknown values of these chemical/physical and geological parameters with unknown (a priori) paleoheat flux variations, a theoretical scheme, called thermal indicator tomography, is developed for the systematic determination of all parameters at the same time and on the same footing. Case histories will be discussed in subsequent papers in this series.