地表一致性静校正量误差分析

目前,在生产中广泛使用的静校正方法几乎都是基于地表一致性假设的。在地表复杂地区,地形起伏较大,低速带厚度和横向速度变化较大,地表一致性假设往往与实际地质情况不符,在这种情况下再运用地表一致性静校正,势必产生很大的误差。笔者利用模型正演方法着重讨论复杂地区地表一致性静校正误差存在的原因及这种原因跟基准面、炮检距、低速带和基岩速度、低速带厚度、反射面深度等之间的关系,分析其变化规律,提出解决复杂地区静校正问题的建议,以期改进现有的静校正方法,达到减小误差的目的。     

北大别造山后伸展期中-下地壳流动形式数值模拟研究

造山后伸展过程中,深部地壳的变形形式反映深部岩石圈的构造活动方式,是造山后伸展深部动力学研究的主要内容。本次工作以大别造山带造山后伸展变形构造为研究对象,在详细的构造观察和组构分析基础上,以实际变形组构为约束开展数值模拟,系统地研究伸展变形构造在垂向剖面上岩石变形机制的空间分布,进而揭示大别造山带造山后伸展过程中深部地壳的变形形式。野外观察表明,北大别穹隆边缘区域拉伸线理优选方位显著,而到穹隆中部拉伸线理定向性减弱,再到穹隆核部线理定向性有所增强。运动学上,北大别穹隆内呈现了一致的上盘向NW的剪切指向。以此为约束的数值模拟研究表明,大别造山带造山后伸展期中-下地壳流动方式总体上是介于透入性简单剪切形式与纯剪切形式之间的组合形式;垂直剖面上呈现为由上、下边界向中心简单剪切分量显著降低,中-下地壳流变带中部以纯剪切变形为主。我们用造山带加厚地壳的韧性垮塌和深俯冲的太平洋板块"后撤"来解释其动力学来源。     

复杂断层网络环境下地下洞室群三维地质分析与应用

断层构造是制约地下工程设计与施工的关键因素。综合考虑复杂地质构造与地下洞室群结构的耦合关系,采用三维地质建模方法建立实现了地层体、断层体与地下洞室群结构的统一模型。基于三维模型提出了断层复杂度的概念,为水利水电工程地下厂房洞室群的选址方案分析提供了一种新的手段。针对地下围岩稳定性问题,提出基于断层与断层、断层与开挖面相互切割的曲面块体识别方法,搜索地下厂房洞室群内可能的失稳区域。上述方法已成功应用在某水电工程地下厂房洞室群分析中,在设计方案优化分析中取得了良好的效果。     

谱模拟方法在高分辨率地震资料处理中的应用

为了克服地震子波的时变影响,避免复杂的时窗大小选取问题,使得谱模拟方法能更好地适应地震信号是平稳信号的假设条件,将改进S变换与谱模拟方法相结合,形成时频域谱模拟方法。在二维时频谱中进行谱模拟处理,该方法能够克服地震子波随地层深度变化的限制。理论模型及实际资料处理效果显示,时频域谱模拟方法能够有效提高地震分辨率,实现地震资料的高分辨处理。     

Fault type analysis along the San Andreas Fault zone: A numerical approach

Finite Element (FE) modeling under plane stress condition is used to analyze the fault type variation with depth along and around the San Andreas Fault (SAF) zone. In this simulation elastic rheology was used and was thought justifiable as the variation in depth from 0.5 km to 20 km was considered. Series of calculations were performed with the variation in domain properties. Three types of models were created based on simple geological map of California, namely, 1) single domain model considering whole California as one homogeneous domain, 2) three domains model including the North American plate, Pacific plate, and SAF zone as separate domains, and 3) Four domains model including the three above plus the Garlock Fault zone. Mohr-Coulomb failure criterion and Byerlee's law were used for the calculation of failure state. All the models were driven by displacement boundary condition imposing the fixed North American plate and Pacific plate motion along N34°W vector up to the northern terminus of SAF and N50°E vector motion for the subducting the Gorda and Juan de Fuca plates. Our simulated results revealed that as the depth increased, the fault types were generally normal, and at shallow depth greater strike slip and some thrust faults were formed. It is concluded that SAF may be terminated as normal fault at depth although the surface expression is clearly strike slip.     

Metabasites with eclogite facies relics from Variscides in Sardinia,Italy: a review

Metabasites with eclogite facies relics occur in northern Sardinia as massive to strongly foliated lenses or boudins embedded within low- to medium-grade rocks (Anglona) and migmatites (NE Sardinia). U–Pb zircon dating yielded 453 ± 14, 457 ± 2 and 460 ± 5 Ma as the protolith ages; 400 ± 10 and 403 ± 4 Ma have been interpreted as the ages of the HP event and 352 ± 3 and 327 ± 7 Ma as the ages of the main Variscan retrograde events. A pre-eclogite stage is documented by the occurrence of tschermakite, zoisite relics within garnet porphyroblasts (Punta de li Tulchi) and an edenite–andesine inclusion within a relict kyanite porphyroblast (Golfo Aranci). Four main metamorphic stages have been distinguished in the eclogite evolution: (1) eclogite stage, revealed by the occurrence of armoured omphacite relics within garnet porphyroblasts. The Golfo Aranci eclogites also include kyanite, Mg-rich garnet and pargasite; (2) granulite stage, producing orthopyroxene and clinopyroxene–plagioclase symplectites replacing omphacite. At Golfo Aranci, the symplectitic rims around relict kyanite consist of sapphirine, anorthite, corundum and spinel; (3) amphibolite stage, leading to the formation of amphibole–plagioclase kelyphites between garnet porphyroblasts and pyroxene–plagioclase symplectites and to the growth of cummingtonite on orthopyroxene. Tschermakite to Mg-hornblende, plagioclase, cummingtonite, ilmenite, titanite and biotite are coexisting phases; (4) greenschist to sub-greenschist stage, defined by the appearance of actinolite, chlorite, epidote ss, titanite, sericite and prehnite. The following P–T ranges have been estimated for the different stages. Eclogite stage 550–700°C; 1.3–1.7 GPa; granulite stage 650–900°C; 0.8–1.2 GPa, clustering in the range 1.0–1.2 GPa; amphibolite stage 550–740°C; 0.3–0.7 GPa; greenschist stage 300–400°C; 0.2–0.3 GPa. Comparable ranges characterise the other Variscan massifs in Europe; eclogite stage: T = 530–800°C; P from 0.7–1.1 to 1.7 ± 0.3 GPa; granulite stage T = 760–870°C and P from 1.1–1.4 to 7.2–9.9 GPa, clustering around 1.0–1.2 GPa. Whole-rock chemistry: Sardinian eclogites are N- to T-MORB; European ones N- to E-MORB or calc-alkaline.     

有限差分的时间和空间频散及其对逆时偏移和全波形反演的影响

有限差分是最常用的地震波方程数值模拟方法,但时间和空间离散会产生数值频散.正演是逆时偏移和全波形反演的基本单元,成像和反演的精度很大程度依赖于所采用的数值模拟算法.本文研究了有限差分的时间和空间频散特性及其对逆时偏移和全波形反演的影响.通过时间有限差分+伪谱法、时间频散校正+空间有限差分、时间频散校正+伪谱法获取时间频散、空间频散和无频散数据;发展了抗时间频散、抗空间频散、抗时间+空间频散的逆时偏移和全波形反演方法;采用理论模型和实际资料对提出的方法进行了测试.数值结果表明：逆时偏移同时受时间和空间频散影响,时间频散导致同相轴不聚焦、成像位置偏离,空间频散会产生高频噪声和虚假反射界面;全波形反演在低频大尺度反演中几乎不受时间和空间频散影响,高频精细反演中时间频散引起波形相移、降低反演精度,空间频散增加多解性、导致反演不收敛;抗频散方法可以有效缓解时间和空间频散影响,获得高质量的偏移剖面和反演结果.

     

Characterizing landscape spatial heterogeneity in multisensor images with variogram models

Most evaluation of the consistency of multisensor images have focused on Normalized Difference Vegetation Index （NDVI） products for natural landscapes, often neglecting less vegetated urban landscapes. This gap has been filled through quantifying and evaluating spatial heterogeneity of urban and natural landscapes from QuickBird, Satellite pour l＇observation de la Terre （SPOT）, Ad- vanced Spacebome Thermal Emission and Reflection Radiometer （ASTER） and Landsat Thematic Mapper （TM） images with variogram analysis. Instead of a logarithmic relationship with pixel size observed in the corresponding aggregated images, the spatial variability decayed and the spatial structures decomposed more slowly and complexly with spatial resolution for real multisensor im- ages. As the spatial resolution increased, the proportion of spatial variability of the smaller spatial structure decreased quickly and only a larger spatial structure was observed at very coarse scales. Compared with visible band, greater spatial variability was observed in near infrared band for both densely and less densely vegetated landscapes. The influence of image size on spatial heterogeneity was highly dependent on whether the empirical sernivariogram reached its sill within the original image size. When the empirical semivariogram did not reach its sill at the original observation scale, spatial variability and mean characteristic length scale would increase with image size; otherwise they might decrease. This study could provide new insights into the knowledge of spatial heterogeneity in real multisen- sor images with consideration of their nominal spatial resolution, image size and spectral bands.     

Multi-physics modeling of flow and transport in porous media using a downscaling approach

Multi-phase flow and transport processes generally occur on different spatial and temporal scales. Very often also, within a physical system, they vary in space meaning that different kinds of processes might take place in different parts of the system. In order to account for the variety of processes and to take their scale-dependence into account, the development of multi-scale multi-physics techniques can be envisaged.     

Parametric study for optimal design of large piled raft foundations on sand

In designing piled raft foundations, controlling the total and differential settlements as well as the induced bending moments of the raft is crucial. The majority of piled raft foundations have been designed by placing piles uniformly. In such a design method, the settlements of the piled rafts are likely to be large, which leads to an increase of the pile length and/or number of piles required to reduce the settlements. However, this increase does not satisfy the requirement for economical design. On the basis of a parametric study, this paper contributes a framework for considering an economical design methodology in which piles are placed more densely beneath the column positions when the piled raft is subjected to column loads. The analysis uses PLAXIS 3D software, and the validity of the parametric study is examined through the results of centrifuge model tests conducted by the authors. The study shows that the concentrated pile arrangement method can help to considerably reduce the total and differential settlements as well as the induced bending moments of the raft. Moreover, the effects of parameters, such as pile length, pile number, raft thickness and load types, on the piled raft behavior are investigated. This study can help practicing engineers choose pile and raft parameters in combination with the concentrated pile arrangement method to produce an economical design.